hihoCoder #1457 : 后缀自动机四·重复旋律7（后缀自动机 + 拓扑排序)

Posted 2020-11-02 日拱一卒功不唐捐

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了hihoCoder #1457 : 后缀自动机四·重复旋律7（后缀自动机 + 拓扑排序)相关的知识，希望对你有一定的参考价值。

http://hihocoder.com/problemset/problem/1457

val[i] 表示状态i所表示的所有字符串的十进制之和

ans= ∑ val[i]
在后缀自动机上，从起始状态走任意一条路径到达任意一个状态，这条路径上的字符就是到达的状态的字符串之一

所以利用拓扑排序，记录从起始状态走到这个状态的路径数，即这个状态的字符串个数 sum

若后缀自动机上有边u-->v，加的是数字m，sum[v]+=sum[u],val[v]+=val[u]*10+sum[u]*m

至于多个串，中间加特殊字符，拓扑排序的时候不走特殊字符边即可

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

const int mod=1e9+7;

#define N 2000001

int tot=1,ch[N<<2][11];
int fa[N<<1],len[N<<1];
int dep[N<<1];
int last=1,p,q,np,nq;

char s[N];

int v[N<<1];
int sa[N<<1];

int sum[N<<1],val[N<<1];

void extend(int c)
{
    len[np=++tot]=len[last]+1;
    dep[np]=dep[last]+1;
    for(p=last;p && !ch[p][c];p=fa[p]) ch[p][c]=np;
    if(!p) fa[np]=1;
    else
    {
        q=ch[p][c];
        if(len[q]==len[p]+1) fa[np]=q;
        else
        {
            nq=++tot;
            fa[nq]=fa[q];
            memcpy(ch[nq],ch[q],sizeof(ch[nq]));
            fa[q]=fa[np]=nq;
            len[nq]=len[p]+1;
            dep[nq]=dep[p]+1;
            for(;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
        }
    }
    last=np;
}

void ADD(int &x,int y)
{
    x+=y;
    x-=x>=mod ? mod : 0;
}    

int main()
{
    int n;
    scanf("%d",&n);
    int L;
    while(n--)
    {
        scanf("%s",s+1);
        L=strlen(s+1);
        for(int i=1;i<=L;++i) extend(s[i]-‘0‘);
        if(n) extend(10);
    }
    for(int i=1;i<=tot;++i) v[dep[i]]++;
    for(int i=1;i<=tot;++i) v[i]+=v[i-1];
    for(int i=1;i<=tot;++i) sa[v[dep[i]]--]=i;
    sum[1]=1;
    int x;
    for(int i=1;i<tot;++i) 
    {
        x=sa[i];
        for(int j=0;j<10;++j)
            if(ch[x][j]) 
            {
                ADD(sum[ch[x][j]],sum[x]);
                ADD(val[ch[x][j]],(1LL*val[x]*10+1LL*j*sum[x])%mod);
            }
    }
    int ans=0;
    for(int i=1;i<=tot;++i) ADD(ans,val[i]);
    printf("%d",ans);
}