[bzoj3670] [Noi2014]动物园

　　初二那时候这题成功0分。

　　然而现在对于复杂度还是有点乱。。。

　　暴力做法的话。。显然就是使劲在fail树上跳，看fail链上有多少个<=长度/2的节点。

　　如果对于fail树上的节点，记录一下根到它的节点数的话。。每次查询的时候我们只要找到刚好<=长度/2的那个分界点就行了。

　　对于一条fail链，分界点显然是随着深度而单调的。。。所以每次从fail链最末的节点开始，把整条链都扒出来，利用分界点的单调性求解。

　　复杂度的话。。并不是严格O(n)...要看这些fail链有多少节点会被重复计算。

　　比方说。。字符串 aaaaaa.......如果后面有很多个节点的fail都指向那一坨'a'的末端的话。。。

　　似乎就会有很多重复计算，但实际上最末的那个'a'一定是分界点。

　　因为fail的大小的增长是线性的...每往后一位，fail的大小最多+1

　　所以大概复杂度就不用担心了吧= =

速度极慢，应该是比较丑的。。

#include<cstdio>
#include<iostream>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=1000233;
const ll modd=1000000007;
char s[maxn];
int fail[maxn],sm[maxn],st[maxn],top;
int i,j,k,n,m,l;
ll an,ans[maxn];
bool u[maxn];

inline void getfail(int n){
    int i,p=0;fail[1]=0;
    for(i=2;i<=n;i++){
        while(p&&s[p+1]!=s[i])p=fail[p];
        p+=s[p+1]==s[i],fail[i]=p;
    }
}
int main(){
    int T;scanf("%d",&T);
    for(int ii=1;ii<=T;ii++){
        
        scanf("%s",s+1);
    //    for(i=1;i<=1000000;i++)if(i&2)s[i]='b';else s[i]='a';
        s[0]='%';n=strlen(s)-1;
        if(ii>1)memset(sm,0,(n+1)<<2),memset(u,0,n+1);
        getfail(n);
        for(i=1;i<=n;i++)sm[i]=sm[fail[i]]+1;
        
        for(i=n;i;i--)if(!u[i])
            if((fail[i]<<1)<=i)u[i]=1,ans[i]=sm[fail[i]];
            else{
                for(st[top=1]=j=i;fail[j];st[++top]=(j=fail[j]));
                l=top+1;st[top+1]=0;
                for(j=top;j;j--){
                    while((st[l-1]<<1)<=st[j])l--;
                    ans[st[j]]=sm[st[l]];
                    u[st[j]]=1;
                }
            }
        for(i=an=1;i<=n;i++){
            an*=ans[i]+1;
            if(an>=modd)an%=modd;
        //    printf("  %d %lld
",i,ans[i]);
        }
        printf("%lld
",an);
    }
    return 0;
}