BZOJ 3676: [Apio2014]回文串 后缀自动机 Manacher 倍增

2021-04-21 14:31

阅读:389

http://www.lydsy.com/JudgeOnline/problem.php?id=3676

过程很艰难了,第一次提交Manacher忘了更新p数组,超时,第二次是倍增的第0维直接在自动机里完成,但是忽略了增加新点时fa变动的情况,还是肉眼查错最管用。

得到的教训是既然倍增就在倍增的函数里完成,自动机就在自动机里完成,不要随便乱搞赋值。

技术分享图片技术分享图片
 1 #include 2 #include 3 #include
 4 #include 5 #include 6 #include 7 using namespace std;
 8 const int maxn=300010;
 9 char ch[maxn]={};
10 char ch1[maxn*2]={};
11 int siz,siz1;
12 struct sam{
13     int sig[26];
14     int f,len,d[21],num;
15 }t[maxn*2];
16 int tot=1,la=1;
17 int loc[maxn]={},cnt[maxn*2]={},b[maxn*2]={},p[maxn*2]={};
18 long long ans=0;
19 void add(int z){
20     int x=++tot,i=la;
21     t[x].len=t[la].len+1;
22     for(;i&&!t[i].sig[z];i=t[i].f)
23         t[i].sig[z]=x;
24     if(!i)t[x].f=1;
25     else{
26         int p=t[i].sig[z];
27         if(t[p].len==t[i].len+1)t[x].f=p;
28         else{
29             int y=++tot;
30             t[y]=t[p];t[y].len=t[i].len+1;
31             t[p].f=t[x].f=y;
32             for(;i&&t[i].sig[z]==p;i=t[i].f)
33                 t[i].sig[z]=y;
34         }
35     }
36     la=x;
37 }
38 void pre(){
39     int z;
40     for(int i=1;i;
41     for(int i=1;i1];
42     for(int i=tot;i;i--)b[cnt[t[i].len]--]=i;
43     for(int i=1;i1;
44     for(int i=tot;i;i--){
45         z=b[i];t[t[z].f].num+=t[z].num;
46     }
47     for(int i=1;i){
48         z=b[i];t[z].d[0]=t[z].f;
49         for(int j=1;j20;j++)
50             t[z].d[j]=t[t[z].d[j-1]].d[j-1];
51     }
52 }
53 void getit(int l,int r){
54     r=r/2;l=(l+1)/2;
55     if(rreturn;
56     int z=loc[r];
57     for(int i=20;i>-1;i--){
58         if(t[t[z].d[i]].len>=r-l+1){
59             z=t[z].d[i];        
60         }
61     }
62     long long sum=(long long)t[z].num*(r-l+1);
63     if(anssum;
64 }
65 void Manacher(){
66     for(int i=1;i2]=ch[i],ch1[i*2-1]=#;
67     ch1[siz*2+1]=#;ch1[siz*2+2]=\0;ch1[0]=$;
68     siz1=siz*2+1;
69     int mx=0,id=0;
70     for(int i=1;i){
71         if(mx>=i)p[i]=min(mx-i,p[id*2-i]);
72         else p[i]=0;
73         while(ch1[i+p[i]]==ch1[i-p[i]]){++p[i];getit(i-p[i]+1,i+p[i]-1);}
74         if(mxi;}
75     }
76 }
77 int main(){
78     memset(t,0,sizeof(t));
79     scanf("%s",ch+1);siz=strlen(ch+1);
80     for(int i=1;i1;add(ch[i]-a);}
81     pre();
82     Manacher();
83     printf("%lld\n",ans);
84     return 0;
85 }
View Code

 


评论


亲,登录后才可以留言!