马拉车算法 最长回文串

1. const int maxn=1000010;  
2. char str[maxn];//原字符串  
3. char tmp[maxn<<1];//转换后的字符串  
4. int Len[maxn<<1];  
5. //转换原始串  
6. int INIT(char *st)  
7. {  
8.     int i,len=strlen(st);  
9.     tmp[0]='@';//字符串开头增加一个特殊字符，防止越界  
10.     for(i=1;i<=2*len;i+=2)  
11.     {  
12.         tmp[i]='#';  
13.         tmp[i+1]=st[i/2];  
14.     }  
15.     tmp[2*len+1]='#';  
16.     tmp[2*len+2]='$';//字符串结尾加一个字符，防止越界  
17.     tmp[2*len+3]=0;  
18.     return 2*len+1;//返回转换字符串的长度  
19. }  
20. //Manacher算法计算过程  
21. int MANACHER(char *st,int len)  
22. {  
23.      int mx=0,ans=0,po=0;//mx即为当前计算回文串最右边字符的最大值  
24.      for(int i=1;i<=len;i++)  
25.      {  
26.          if(mx>i)  
27.          Len[i]=min(mx-i,Len[2*po-i]);//在Len[j]和mx-i中取个小  
28.          else  
29.          Len[i]=1;//如果i>=mx，要从头开始匹配 
    
30. 
    
31.          while(st[i-Len[i]]==st[i+Len[i]])  
32.          Len[i]++; 
    
33. 
    
34.          if(Len[i]+i>mx)//若新计算的回文串右端点位置大于mx，要更新po和mx的值  
35.          {  
36.              mx=Len[i]+i;  
37.              po=i;  
38.          }  
39.          ans=max(ans,Len[i]);  
40.      }  
41.      return ans-1;//返回Len[i]中的最大值-1即为原串的最长回文子串额长度   
42.   }