This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include "bits/stdc++.h"
#define MAXN 300009
#define INF 1000000007
#define mp(x,y) make_pair(x,y)
#define all(v) v.begin(),v.end()
#define pb(x) push_back(x)
#define wr cout<<"----------------"<<endl;
#define ppb() pop_back()
#define tr(ii,c) for(__typeof((c).begin()) ii=(c).begin();ii!=(c).end();ii++)
#define ff first
#define ss second
#define my_little_dodge 46
#define debug(x) cerr<< #x <<" = "<< x<<endl;
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
template<class T>bool umin(T& a,T b){if(a>b){a=b;return 1;}return 0;}
template<class T>bool umax(T& a,T b){if(a<b){a=b;return 1;}return 0;}
const int LOG=18;
int S[MAXN][LOG],c[MAXN],ord[MAXN],lcp[MAXN],rad[MAXN<<1];
PII C[MAXN];
char s[MAXN],tmp[MAXN<<1];
bool cmp(int x,int y){
return (C[x]<C[y]);
}
int ata[MAXN],sz[MAXN];
int tap(int x){
if(ata[x]==x)return x;
return ata[x]=tap(ata[x]);
}
void merge(int x,int y){
if((x=tap(x))==(y=tap(y)))
return;
ata[y]=x;
}
int main(){
//freopen("file.in", "r", stdin);
scanf("%s",s+1);
int n=strlen(s+1);
for(int i=1;i<=n;i++)
S[i][0]=s[i]-'a'+1;
for(int j=1;j<LOG;j++){
for(int i=1;i<=n;i++)
C[i]=mp(S[i][j-1],(i+(1<<(j-1))<=n?S[i+(1<<(j-1))][j-1]:-1)),c[i]=i;
sort(c+1,c+n+1,cmp);
for(int i=1;i<=n;i++)
S[c[i]][j]=S[c[i-1]][j]+(C[c[i]]!=C[c[i-1]]);
}
for(int i=1;i<=n;i++)
ord[S[i][LOG-1]]=i;
vector<pair<int,PII> >v;
vector<PII>g;
for(int i=1;i<n;i++){
int a=ord[i],b=ord[i+1];
for(int j=LOG-1;j>=0;j--)
if(a<=n and b<=n and S[a][j]==S[b][j])
a+=(1<<j),b+=(1<<j),lcp[i]+=(1<<j);
v.pb(mp(lcp[i],mp(ord[i],ord[i+1])));
}
sort(all(v),greater<pair<int,PII>>());
for(int i=0;i<2*n+1;i++)tmp[i]='#';
for(int i=0;i<n;i++)tmp[i*2+1]=s[i+1];
for(int i=0,j=0,k;i<n*2+1;i+=k){
while(i-j>=0&&i+j<n*2+1&&tmp[i-j]==tmp[i+j])j++;
rad[i]=j;k=1;
while(rad[i-k]<rad[i]-k)rad[i+k]=rad[i-k],k++;
j=max(0,j-k);
}
//odd
for(int i=1;i<=n;i++)g.pb(mp(rad[i*2-1]/2,i));
sort(all(g),greater<PII>());
int p=0;ll ans=0;
for(int i=1;i<=n;i++)ata[i]=i,sz[i]=0;
for(int i=0;i<n;i++){
while(p<int(v.size()) and v[p].ff>=g[i].ff)
merge(v[p].ss.ff,v[p].ss.ss),p++;
sz[tap(g[i].ss)]++;
umax(ans,sz[tap(g[i].ss)]*1LL*(g[i].ff*2-1));
}
//even
g.clear();
for(int i=1;i<n;i++)g.pb(mp(rad[i*2]/2,i+1));
sort(all(g),greater<PII>());p=0;
for(int i=1;i<=n;i++)ata[i]=i,sz[i]=0;
for(int i=0;i<n-1;i++){
while(p<int(v.size()) and v[p].ff>=g[i].ff)
merge(v[p].ss.ff,v[p].ss.ss),p++;
sz[tap(g[i].ss)]++;
umax(ans,sz[tap(g[i].ss)]*1LL*(g[i].ff*2));
}
printf("%lld\n",ans);
return 0;
}
Compilation message (stderr)
palindrome.cpp: In function 'int main()':
palindrome.cpp:39:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
39 | scanf("%s",s+1);
| ~~~~~^~~~~~~~~~
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |