【bzoj4084】【sdoi2015】双旋转字符串

• 题解

  • 首先题中说了$n>=m$;
  • 分成的循环串左右两边为本质相同的单循环串循环串，分别长为$l = frac{n + m}{2} $;
  • 所以$S$串的前$l$位为双循环串的一半$S1$，后一半为$S2$;
  • 倍长$S1$找到$S2$在$S1$中出现的位置，把单循环串匹配位置后的剩下的$m$位哈希存下来；
  • 读入每个$T$串查找即可
  • 注意有可能$S2$在$S1$中会重复出现，记录$last$防止重复计算；

  • #include<bits/stdc++.h>
    #define ull unsigned long long 
    #define base 1234567891
    #define il inline
    #define rg register 
    using namespace std;
    const int N=8000010,sz=1e7;
    int n,m,k,l,t1,t2,o,nt[sz],hd[sz],w[sz],nxt[N<<1],lst[sz];
    char s[N<<1],t[N<<1];
    ull h[N<<1],v[sz],pw,ans;
    il char gc(){
        static char*p1,*p2,buf[1000000];
        if(p1==p2)p2=(p1=buf)+fread(buf,1,1000000,stdin);
        return(p1==p2)?EOF:*p1++;
    }
    il int rd(){
        int x=0; char c=gc();
        while(c<'0'||c>'9')c=gc();
        while(c>='0'&&c<='9')x=x*10+c-'0',c=gc();
        return x;
    }
    il char gt(){char c=gc();while(!isalpha(c))c=gc();return c;}
    void solve(){
        for(rg int i=1;i<=n;i++){
            for(rg int j=1;j<=l;j++)h[j]=h[j-1]*base+(s[j]=gt());
            for(rg int j=1;j<=l;j++)h[j+l]=h[j+l-1]*base+s[j];
            if(t1)t[1]=gt();
            for(rg int j=2,tj=0;j<=t1;nxt[++j]=tj){
                t[j]=gt();
                while(tj&&t[j]!=t[tj+1])tj=nxt[tj];
                if(t[j]==t[tj+1])tj++;
            }
            for(rg int j=1,tj=0;j<=l+max(0,t1-1);j++){
                while(tj&&s[j]!=t[tj+1])tj=nxt[tj];
                if(s[j]==t[tj+1])tj++;
                if(!t1||tj==t1){
                    ull x = h[j+t2] - h[j]*pw;
                    int fg=0;
                    for(int r=hd[x%sz];r;r=nt[r])if(v[r]==x){
                        fg=1;
                        if(lst[r]!=i)w[r]++,lst[r]=i;
                        break;
                    }
                    if(!fg)nt[++o]=hd[x%sz],v[hd[x%sz]=o]=x,w[o]=1,lst[o]=i;
                    tj=nxt[tj];
                }
            }
        } 
        t1-=l;
        for(rg int i=1;i<=m;i++){
            ull x=0;for(rg int j=1;j<=t2;j++)x=x*base+gt();
            for(int r=hd[x%sz];r;r=nt[r])if(v[r]==x){ans+=w[r];break;}    
        }
    }
    int main(){
        #ifndef ONLINE_JUDGE
        freopen("bzoj4084.in","r",stdin);
        freopen("bzoj4084.out","w",stdout);
        #endif
        n=rd();m=rd();t1=rd();t2=rd();l=(t1+t2)>>1;
        t1-=l;for(int i=pw=1;i<=t2;i++)pw=pw*base;
        solve();
        printf("%llu
    ",ans); 
        return 0;
    }