Manacher--雾窗寒对遥天暮，暮天遥对寒窗雾

POJ 3974: Palindrome

题意：

最长回文子串的长度...

分析：

Manacher板子题...

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
using namespace std;
//眉眼如初，岁月如故 

const int maxn=2000000+5;

int cas,len,p[maxn];

char s[maxn],str[maxn];

inline void prework(void){
    int i=0;
    for(i=0;s[i];i++)
        str[i*2+1]='#',str[(i+1)*2]=s[i];
    len=i*2+1;str[0]='$';str[len]=str[len+1]='#';
}

inline void manacher(void){
    int id,mx=0,ans=0;
    for(int i=1;i<=len;i++){
        p[i]=i<mx?min(mx-i,p[id*2-i]):1;
        while(str[i+p[i]]==str[i-p[i]])
            p[i]++;
        if(p[i]+i>mx)
            mx=p[i]+i,id=i;
        if(ans<p[i]-1)
            ans=p[i]-1;
    }
    printf("%d
",ans);
}

signed main(void){cas=0;
    while(scanf("%s",s)&&s[0]!='E'){
        printf("Case %d: ",++cas);
        prework();manacher();
    }
    return 0;
}//Cap ou pas cap. Pas cap. 

HDU 3948: The Number of Palindromes

题意：

求原串中本质不同的子串个数...

分析：

写完发现大家都用的后缀数组...然而我不会...只能用Manacher+Hash水...

 本质不同的回文串个数级别是O(n)的...所以我们可以暴力扫描回文串...我们枚举回文串的中心位置，然后判断最长回文串，如果没有出现过则判断其子串，否则直接break...

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
#define int long long
#define II int
using namespace std;
//眉眼如初，岁月如故 

const II maxn=200000+5,Mod=23333333;

int l,cas,len,p[maxn],po[maxn],Hash[maxn]; 

char s[maxn],str[maxn];

struct M{
    int hd[Mod+10],nxt[maxn],w[2][maxn],cnt;
    
    inline void clear(void){
        memset(hd,-1,sizeof(hd));cnt=0;
    }
    
    inline void insert(int x,int l,int r){
        w[0][cnt]=l,w[1][cnt]=r,nxt[cnt]=hd[x],hd[x]=cnt++;
    }
    
    inline bool find(int x,int l,int r){//cout<<"find"<<endl;cout<<l<<" "<<r<<endl;
        for(int i=hd[x];i!=-1;i=nxt[i])
            if(r-l+1==w[1][i]-w[0][i]+1){
                int j;
                for(j=1;j<=r-l+1;j++)
                    if(str[w[0][i]+j-1]!=str[w[0][i]+j-1])
                        break;
                if(j==r-l+2)//{
                    return true;//cout<<"finish find"<<endl;}
            }//cout<<"finish find"<<endl;
        return false;
    }
    
}mp;

inline void prework(void){
    int i=0;l=strlen(s);
    for(i=0;s[i];i++)
        str[i*2+1]='#',str[(i+1)*2]=s[i];
    len=i*2+1;str[0]='$';str[len]=str[len+1]='#';
    for(int i=1;i<=len;i++){
        if(str[i]>='a'&&str[i]<='z')
            Hash[i]=(Hash[i-1]*27%Mod+str[i]-'a'+1)%Mod;
        else
            Hash[i]=(Hash[i-1]*27%Mod+27)%Mod;
    }
}

inline void manacher(void){
    int id,mx=0;
    for(int i=1;i<=len;i++){
        p[i]=i<mx?min(mx-i,p[(id<<1)-i]):1;
        while(str[p[i]+i]==str[i-p[i]])
            p[i]++;
        if(i+p[i]>mx)
            mx=i+p[i],id=i;
    }
}

inline void find(void){
    int ans=0;
    for(int i=1;i<=len;i++){//cout<<i<<" "<<p[i]<<endl;
        if(!mp.find((Hash[i+p[i]-1]-Hash[i-1]*po[p[i]]%Mod+Mod)%Mod,i,i+p[i]-1)){
            ans++,mp.insert((Hash[i+p[i]-1]-Hash[i-1]*po[p[i]]%Mod+Mod)%Mod,i,i+p[i]-1);
            int lala=p[i]-2;
            while(lala>=1&&!mp.find((Hash[i+lala-1]-Hash[i-1]*po[lala]%Mod+Mod)%Mod,i,i+lala-1))
                ans++,mp.insert((Hash[i+lala-1]-Hash[i-1]*po[lala]%Mod+Mod)%Mod,i,i+lala-1),lala-=2;
        }
    }
    printf("%lld
",ans-1);
}

signed main(void){
    scanf("%lld",&cas);po[0]=1;Hash[0]=0;
    for(int i=1;i<=100000;i++)
        po[i]=po[i-1]*27%Mod;
    for(int i=1;i<=cas;i++){mp.clear();
        scanf("%s",s);printf("Case #%lld: ",i);
        prework();manacher();find();
    }
    return 0;
}//Cap ou pas cap. Pas cap. 

---

By NeighThorn