【hdu4436/LA6387-str2int】sam处理不同子串

题意：给出n个数字，数字很长，用字符串读入，长度总和为10^5。求这n个字符串的所有子串（不重复）的和取模2012 。

例如字符串101，和就是1+10+101=112。

题解：

就是求不同的子串连成一个数。
sam的拓扑序真的很有用！按拓扑序可以保证能转移到当前x的节点都在之前被更新了。
每个节点x维护cnt表示root到x的方案数，sum表示以x为结尾的子串和。
找拓扑序有两种方法：1.拓扑排序一样bfs(O(n)) 2.按照step[x]排序（若x有个孩子是y，step[x]<=step[y](O(nlogn))
按照拓扑序for一遍，对于x，它的孩子y,cnt[y]+=cnt[x],sum[y]+=sum[x]*10+cnt[x]*(y所代表的数字)。
不能有前缀0----->根节点不走0孩子。
节点数开了10^6才AC。tell me why?!

打了两种求拓扑序的方法。

贴一下代码。

  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<cstring>
  4 #include<iostream>
  5 #include<queue>
  6 #include<ctime>
  7 #include<algorithm>
  8 using namespace std;
  9
 10 const int N=10*100010,Mod=2012;
 11 char s[N];
 12 int n,sl,cl,ans,tot,last;
 13 int son[N][30],pre[N],step[N],sum[N],in[N],cnt[N],c[N];
 14 bool vis[N];
 15 queue<int> Q;
 16
 17 int add_node(int x)
 18 {
 19     step[++tot]=x;
 20     return tot;
 21 }
 22
 23 void extend(int ch)
 24 {
 25     int p=last,np=add_node(step[last]+1);
 26     while(p && !son[p][ch]) son[p][ch]=np,in[np]++,p=pre[p];
 27     if(!p) pre[np]=1;
 28     else
 29     {
 30         int q=son[p][ch];
 31         if(step[q]==step[p]+1) pre[np]=q;
 32         else
 33         {
 34             int nq=add_node(step[p]+1);
 35             memcpy(son[nq],son[q],sizeof(son[q]));
 36             for(int i=0;i<=9;i++)
 37                 if(son[q][i]) in[son[q][i]]++;
 38             pre[nq]=pre[q];
 39             pre[np]=pre[q]=nq;
 40             while(son[p][ch]==q) in[q]--,in[nq]++,son[p][ch]=nq,p=pre[p];
 41         }
 42     }
 43     last=np;
 44 }
 45
 46 void find_tp_1()
 47 {
 48     for(int i=2;i<=tot;i++)
 49     {
 50         if(in[i]==0)
 51         {
 52             for(int j=0;j<=9;j++)
 53             {
 54                 int y=son[i][j];
 55                 if(!y) continue;
 56                 in[y]--;
 57             }
 58         }
 59     }
 60     while(!Q.empty()) Q.pop();
 61     Q.push(1);vis[1]=1;cl=0;
 62     while(!Q.empty())
 63     {
 64         int x=Q.front();vis[x]=0;c[++cl]=x;Q.pop();
 65         for(int i=0;i<=9;i++)
 66         {
 67             int y=son[x][i];
 68             if(!y) continue;
 69             in[y]--;
 70             if(!in[y] && !vis[y]) vis[y]=1,Q.push(y);
 71         }
 72     }
 73 }
 74
 75 bool cmp(int x,int y){return step[x]<step[y];}
 76 void find_tp_2()
 77 {
 78     cl=0;
 79     for(int i=1;i<=tot;i++)
 80     {
 81         if(in[i] || i==1) c[++cl]=i;
 82     }
 83     sort(c+1,c+1+cl,cmp);
 84 }
 85
 86 void solve()
 87 {
 88     cnt[1]=1;
 89     for(int i=1;i<=cl;i++)
 90     {
 91         int x=c[i];
 92         for(int j=0;j<=9;j++)
 93         {
 94             int y=son[x][j];
 95             if(!y || (x==1 && j==0)) continue;
 96             cnt[y]=(cnt[y]+cnt[x])%Mod;
 97             sum[y]=(sum[y]+(sum[x]*10)%Mod+(j*cnt[x])%Mod)%Mod;
 98         }
 99     }
100 }
101
102 int main()
103 {
104     freopen("a.in","r",stdin);
105     // freopen("me.out","w",stdout);
106     while(scanf("%d",&n)!=EOF)
107     {
108         memset(son,0,sizeof(son));
109         memset(step,0,sizeof(step));
110         memset(pre,0,sizeof(pre));
111         memset(in,0,sizeof(in));
112         memset(cnt,0,sizeof(cnt));
113         memset(sum,0,sizeof(sum));
114         tot=0;add_node(0);last=1;
115         for(int i=1;i<=n;i++)
116         {
117             scanf("%s",s+1);
118             sl=strlen(s+1);
119             last=1;
120             for(int j=1;j<=sl;j++)     extend(s[j]-‘0‘);
121         }
122         find_tp_1();//找拓扑序方法一 ： bfs=拓扑排序
123         // find_tp_2();//方法二：直接对step[x]进行排序
124         solve();
125         ans=0;
126         for(int i=1;i<=tot;i++) ans=(ans+sum[i])%Mod;
127         printf("%d\n",ans);
128         // cout<<"Time Used : "<<(double)clock()/CLOCKS_PER_SEC<<" s."<<endl;
129     }
130     return 0;
131 }