后缀自动机求endpos集大小

#include<bits/stdc++.h>
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define LNF 0x3f3f3f3f3f3f3f3f
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
#define pqueue priority_queue
#define NEW(a,b) memset(a,b,sizeof(a))
const double pi=4.0*atan(1.0);
const double e=exp(1.0);
const int maxn=4e5+;
typedef long long LL;
typedef unsigned long long ULL;
const LL mod=1e9+;
const ULL base=1e7+;
const int maxp=+;
using namespace std;
struct Suffix_Node{
    int ch[maxp],par,len,jg;
    LL dd[];
    void init(){
        NEW(ch,);
        par=len=;
    }
};
int num[maxn];
int la[maxn],ra[maxn];
int st[maxn];
int n,k;
ULL ans[maxn];
pair<ULL,int> w[maxn];
class Suffix_Automation{
private:
    Suffix_Node s[maxn];
    int cur,las,siz,crp,dx;
    vector<int> gg[maxn];
    struct node{
        int to,nxt;
    }g[maxn];
    int cnt,head[maxn],pk[maxn];
public:
    Suffix_Automation():las(),cur(),siz(),crp(){}
    const void init(){
        for(int i=;i<=siz;i++) {
            pk[i]=;s[i].init();
            for(int j=;j<=;j++) s[i].dd[j]=;
            gg[i].clear();
            ans[i]=;
            for(int j=;j<;j++){
                f[i][j]=;
            }
        }
        las=cur=siz=crp=;
        cnt=;
        head[]=head[]=-;
        dx=;
    }
    const int match(const char c)const{
        return s[crp].ch[c-'a'];
    }
    const void withdraw(const int len){
        while(crp!=&&s[s[crp].par].len>=len) crp=s[crp].par;
        if(crp==) crp=;
    }
    const void Transfer(const int len,const char c){
        crp=s[crp].ch[c-'a'];
        if(crp==) crp=;
        withdraw(len);
    }
    const void ex_tend(const char c){
        int x=c-'a';
        cur=++siz;
        head[siz]=-;
        s[cur].len =s[las].len+;
        s[cur].jg=;
        while(las!=&&!s[las].ch[x])
            s[las].ch[x]=cur,las=s[las].par;
        if(las==) s[cur].par=;
        else{
            int q,nq;
            q=s[las].ch[x];
            if(s[q].len==s[las].len+)
                s[cur].par=q;
            else{
                nq=++siz;
                head[siz]=-;
                s[nq]=s[q],s[nq].len=s[las].len+;
                s[nq].jg=;
                s[cur].par=s[q].par=nq;
                while(las!=&&s[las].ch[x]==q)
                    s[las].ch[x]=nq,las=s[las].par;
            }
        }
        las=cur;
    }
    void dfs(int u){
        if(f[u][]) return ;
        add(s[u].par,u);
        if(s[u].par!=){
            dfs(s[u].par);
        }
        f[u][]=s[u].par;
    }
    void solve(){
        for(int i=;i<=siz;i++){
            if(!f[i][]){
                dfs(i);
            }
        }
        for(int j=;j<;j++){
            for(int i=;i<=siz;i++){
                f[i][j]=f[f[i][j-]][j-];
            }
        }
    }
    void add(int x,int y){
        g[cnt].to=y;
        g[cnt].nxt=head[x];
        head[x]=cnt++;
    }
    int f[maxn][];
    int d[maxn];
    int fid(int u,int le){
        int y;
        for(int j=;j>=;j--){
            y=f[u][j];
            if(y!=&&s[y].len>=le){
                u=y;
            }
        }
        return u;
    }
    void dfs2(int u,int fa){
        d[u]=d[fa]+;
        la[u]=dx+;
        if(pk[u]!=){
            dx++;
            num[dx]=pk[u];
        }
        for(int i=head[u];i!=-;i=g[i].nxt){
            if(g[i].to!=fa){
                dfs2(g[i].to,u);
                s[u].jg+=s[g[i].to].jg;
            }
        }
        ra[u]=dx;
        gg[s[u].jg].push_back(u);
    }
    int get_dx(){
        return dx;
    }
    void ins(int x){
        pk[crp]=x;
        st[x]=crp;
    }
    void sss(){
        for(int i=n;i>=;i--){
            for(int j=;j<gg[w[i].se].size();j++){
                int y=gg[w[i].se][j];
                while(y)
                {
                    LL er=w[i].fi,flag=;
                    for(int j=;j>=;j--)
                    {
                        if(!(er>>j)) continue;
                        if(!s[y].dd[j]) {s[y].dd[j]=er;flag=;break;}
                        er^=s[y].dd[j];
                    }
                    if(!flag) break;
                    ans[y]+=w[i].fi;y=s[y].par;
                }
            }

        }
    }

}SAM;
char ss[maxn];
int main(){
    int t;
    scanf("%d",&t);
    int x,y,z;
    int ai;
    int gh;
    while(t--){
        SAM.init();
        scanf("%d",&n);
        scanf("%s",ss);
        for(int i=;i<=n;i++){
            scanf("%llu",&w[i].fi);
            w[i].se=i;
        }
        sort(w+,w++n);
        for(int i=;ss[i];i++){
            SAM.ex_tend(ss[i]);
        }
        SAM.solve();
        for(int i=;ss[i];i++){
            SAM.Transfer(i+,ss[i]);
            SAM.ins(i+);
        }
        SAM.dfs2(,);
        scanf("%d",&k);
        int g=SAM.get_dx();
        SAM.sss();
        while(k--){
            scanf("%d%d",&x,&y);
            x=y-x+;
            z=SAM.fid(st[y],x);
            printf("%llu\n",ans[z]);
        }
    }
}

http://acm.hdu.edu.cn/showproblem.php?pid=6694