Luogu4081 USACO17DEC Standing Out from the Herd(广义后缀自动机)

  建出广义SAM,通过parent树对每个节点求出其是否仅被一个子串包含及被哪个包含。

  写了无数个sam板子题一点意思都没啊

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1000000010
#define N 200010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n,m,son[N][26],fail[N],len[N],f[N],q[N],tmp[N],cnt=1,last;
ll ans[N];
char s[N];
int ins(int c,int p,int i)
{
	int u;
    if (son[p][c])
    {
        int q=son[p][c];
        if (len[p]+1==len[q]) u=q;
        else
        {
            int y=++cnt;
            len[y]=len[p]+1;
            memcpy(son[y],son[q],sizeof(son[q]));
            fail[y]=fail[q],fail[q]=y;
            while (son[p][c]==q) son[p][c]=y,p=fail[p];
            u=y;
        }
    }
    else
    {
        int x=++cnt;len[x]=len[p]+1;
        while (!son[p][c]&&p) son[p][c]=x,p=fail[p];
        if (!p) fail[x]=1;
        else
        {
            int q=son[p][c];
            if (len[p]+1==len[q]) fail[x]=q;
            else
            {
                int y=++cnt;
                len[y]=len[p]+1;
                memcpy(son[y],son[q],sizeof(son[q]));
                fail[y]=fail[q],fail[x]=fail[q]=y;
                while (son[p][c]==q) son[p][c]=y,p=fail[p];
            }
        }
        u=x;
    }
    if (f[u]==-1) f[u]=i;else if (f[u]!=i) f[u]=0;
    return u;
}
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
#endif
	n=read();memset(f,255,sizeof(f));
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		m=strlen(s+1);last=1;
		for (int j=1;j<=m;j++) last=ins(s[j]-'a',last,i);
	}
	for (int i=1;i<=cnt;i++) tmp[len[i]]++;
	for (int i=1;i<=cnt;i++) tmp[i]+=tmp[i-1];
	for (int i=1;i<=cnt;i++) q[tmp[len[i]]--]=i;
	for (int i=cnt;i>=1;i--)
	{
		int x=q[i];
		if (f[x]>0) ans[f[x]]+=len[x]-len[fail[x]];
		if (f[fail[x]]==-1) f[fail[x]]=f[x];
		else if (f[fail[x]]!=f[x]) f[fail[x]]=0;
	}
	for (int i=1;i<=n;i++) printf("%lld\n",ans[i]);
	return 0;
	//NOTICE LONG LONG!!!!!
}

  

posted @ 2019-05-08 15:31  Gloid  阅读(226)  评论(0编辑  收藏  举报