[BZOJ3277]串

bzoj

题意

给出\(n\)个串，求每个串中有多少个子串在至少\(k\)个串中出现过。

sol

广义后缀自动机。
先建好广义SAM，也就是把这\(n\)个串全都插到SAM里面去，那么可以维护每个状态在哪些串里面出现过。
可以给每个状态维护一个\(set\)，然后向上合并就是直接并起来。
这里我写基数排序结果鬼掉了，写成建边然后dfs的形式才能过。难道广义后缀自动机不能写基数排序？还求大佬们能教教我。
然后因为是自下而上取并，所以往上跳时集合大小一定是单调不降的。所以在对每一个串在SAM上跑匹配的时候，要一路跳到\(sz[now]\ge k\)的状态，然后计算答案。

code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
using namespace std;
const int N = 2e5+5;
int n,k,last=1,tot=1,tr[N][26],fa[N],len[N],t[N],a[N],size[N];
string s[N];set<int>S[N];
void extend(int c,int id)
{
    int v=last,u=++tot;last=u;
    len[u]=len[v]+1;S[u].insert(id);
    while (v&&!tr[v][c]) tr[v][c]=u,v=fa[v];
    if (!v) fa[u]=1;
    else{
        int x=tr[v][c];
        if (len[x]==len[v]+1) fa[u]=x;x
        else{
            int y=++tot;
            memcpy(tr[y],tr[x],sizeof(tr[y]));
            fa[y]=fa[x];fa[x]=fa[u]=y;len[y]=len[v]+1;
            while (v&&tr[v][c]==x) tr[v][c]=y,v=fa[v];
        }
    }
}
int to[N],nxt[N],head[N],cnt;
void link(int u,int v){to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;}
void dfs(int u)
{
    for (int e=head[u];e;e=nxt[e])
    {
        dfs(to[e]);
        for (set<int>::iterator it=S[to[e]].begin();it!=S[to[e]].end();++it)
            S[u].insert(*it);
    }
    size[u]=S[u].size();
}
int main()
{
    ios::sync_with_stdio(false);
    cin>>n>>k;
    for (int i=1;i<=n;++i)
    {
        cin>>s[i];last=1;
        for (int j=0,l=s[i].length();j<l;++j) extend(s[i][j]-'a',i);
    }
    for (int i=1;i<=tot;++i) if (fa[i]) link(fa[i],i);
    dfs(1);
    if (k>n) {for (int i=1;i<=n;++i) cout<<"0 ";return 0;}
    for (int i=1;i<=n;++i)
    {
        long long ans=0;
        for (int j=0,l=s[i].length(),now=1;j<l;++j)
        {
            now=tr[now][s[i][j]-'a'];
            while (size[now]<k) now=fa[now];
            ans+=len[now];
        }
        cout<<ans<<' ';
    }
    puts("");return 0;
}

转载于:https://www.cnblogs.com/zhoushuyu/p/8660511.html