bzoj2780 [Spoj]8093 Sevenk Love Oimaster（SAM+dfs序+BIT）

这个问题和bzoj3881基本一模一样。
3881是拿ACAM+树链的并写的，先把询问串建出ACAM,然后拿文本串去更新答案。

这道题呢是拿SAM+dfs序+离线+BIT做的。先把文本串建出SAM，在每个点上挂一个vector记录这个节点所属的文本串。然后对于每一个询问串，先跑出它在SAM上的代表节点p，那么他的答案就是p的子树中不同颜色的个数，用dfs序转化为序列就是求区间内不同数的个数。可以离线+BIT解决。
复杂度都是$O(nlogn)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n=0,m,rt,last,son[N<<1][26],par[N<<1],mx[N<<1],in[N<<1],out[N<<1],dfn=0,pos[N<<1];
vector<int>col[N<<1],Son[N<<1];int ans[60010],tot=0,lst[10010],c[N<<1];
char s[360010];
struct quer{
    int l,r,id;
    friend bool operator<(quer a,quer b){return a.r<b.r;}
}qq[60010];
inline void ins(int ch,int id){
    int p=last;
    if(son[p][ch]){
        int q=son[p][ch];
        if(mx[q]==mx[p]+1){last=q;col[q].push_back(id);return;}
        int nq=++n;mx[nq]=mx[p]+1;memcpy(son[nq],son[q],sizeof(son[q]));
        par[nq]=par[q];par[q]=nq;last=nq;col[nq].push_back(id);
        for(;p&&son[p][ch]==q;p=par[p]) son[p][ch]=nq;return;
    }int np=++n;mx[np]=mx[p]+1;col[np].push_back(id);last=np;
    for(;p&&!son[p][ch];p=par[p]) son[p][ch]=np;
    if(!p){par[np]=rt;return;}
    int q=son[p][ch];
    if(mx[q]==mx[p]+1){par[np]=q;return;}
    int nq=++n;mx[nq]=mx[p]+1;memcpy(son[nq],son[q],sizeof(son[q]));
    par[nq]=par[q];par[q]=par[np]=nq;
    for(;p&&son[p][ch]==q;p=par[p]) son[p][ch]=nq;
}
void dfs(int x){
    in[x]=++dfn;pos[dfn]=x;
    for(int i=0;i<Son[x].size();++i) dfs(Son[x][i]);
    out[x]=dfn;
}
inline void gao(int id){
    m=strlen(s+1);int p=rt;
    for(int i=1;i<=m;++i){
        int ch=s[i]-'a';
        if(!son[p][ch]) return;p=son[p][ch];
    }qq[++tot].id=id;qq[tot].l=in[p];qq[tot].r=out[p];
}
inline void add(int x,int val){for(;x<=n;x+=x&-x) c[x]+=val;}
inline int ask(int x){int res=0;for(;x;x-=x&-x) res+=c[x];return res;}
int main(){
//  freopen("a.in","r",stdin);
    int nn=read(),mm=read();last=rt=++n;
    for(int i=1;i<=nn;++i){
        scanf("%s",s+1);m=strlen(s+1);last=rt;
        for(int j=1;j<=m;++j) ins(s[j]-'a',i);
    }for(int i=2;i<=n;++i) Son[par[i]].push_back(i);dfs(1);
    for(int i=1;i<=mm;++i){
        scanf("%s",s+1);gao(i);
    }sort(qq+1,qq+tot+1);int r=1;
    for(int i=1;i<=tot;++i){
        while(r<=qq[i].r){
            for(int j=0;j<col[pos[r]].size();++j){
                int x=col[pos[r]][j];add(r+1,-1);
                add(lst[x]+1,1);lst[x]=r;
            }++r;
        }ans[qq[i].id]=ask(qq[i].l);
    }for(int i=1;i<=mm;++i) printf("%d\n",ans[i]);
    return 0;
}