【XSY2384】【GDOI2017】微信

致去年的我：这是道广义SAM模板题啊……

题意：

Description

Input

Output

HINT

$1\leq N\leq 20$，$1\leq Q\leq 10^5$，字符串总长$\leq 10^6$

题解：

题意就是求若干个trie的最长公共子串……

先把所有trie并起来建广义SAM，记录一下每个节点原来属于哪个trie，由于$N$很小，可以直接状压DP，按照parent树从上往下转移就好了……

预处理答案，询问可以$O(1)$处理；

但是本题有一个坑点是广义SAM必须要bfs建立以达到严格$O(n\times 字符集大小)$的时间复杂度，否则会被特殊构造的数据卡到85分。

代码：

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<queue>
 #define inf 2147483647
 #define eps 1e-9
 using namespace std;
 typedef long long ll;
 typedef double db;
 int n,qq,t,top,len,last,cnt=,tot=,rt=,rts[],s1[],s2[],son[][],fa[],mx[],s[],p[],ch[][],ss[],f[];
 char st[],nw[],qr[];
 queue<int>q;
 void ins(char *s,int len,int id){
     int nw=rts[len-];
     if(!ch[nw][s[len]-'a'])ch[nw][s[len]-'a']=++cnt;
     ss[ch[nw][s[len]-'a']]|=id;
     rts[len]=ch[nw][s[len]-'a'];
 }
 int extend(int p,int ch,int id){
     int np=++tot;
     mx[np]=mx[p]+;
     s[np]=id;
     for(;p&&!son[p][ch];p=fa[p])son[p][ch]=np;
     if(!p)fa[np]=rt;
     else{
         int q=son[p][ch];
         if(mx[q]==mx[p]+)fa[np]=q;
         else{
             int nq=++tot;
             s[nq]=id;
             mx[nq]=mx[p]+;
             memcpy(son[nq],son[q],sizeof(son[q]));
             fa[nq]=fa[q];
             fa[q]=fa[np]=nq;
             for(;p&&son[p][ch]==q;p=fa[p])son[p][ch]=nq;
         }
     }
     return last=np;
 }
 void build(){
     q.push(rt);
     p[rt]=;
     while(!q.empty()){
         int u=q.front();
         q.pop();
         for(int i=;i<;i++){
             int v=ch[u][i];
             if(v){
                 p[v]=extend(p[u],i,ss[v]);
                 q.push(v);
             }
         }
     }
     for(int i=;i<=tot;i++)s1[mx[i]]++;
     for(int i=;i<=tot;i++)s1[i]+=s1[i-];
     for(int i=;i<=tot;i++)s2[s1[mx[i]]--]=i;
     for(int i=tot;i;i--){
         s[fa[s2[i]]]|=s[s2[i]];
         f[s[s2[i]]]=max(f[s[s2[i]]],mx[s2[i]]);
     }
 }
 int main(){
     scanf("%d",&n);
     rts[]=;
     for(int i=;i<=n;i++){
         scanf("%s",st);
         len=strlen(st);
         top=;
         for(int j=;j<len;j++){
             if(st[j]=='<')top--;
             else{
                 nw[++top]=st[j];
                 ins(nw,top,<<(i-));
             }
         }
     }
     build();
     for(int i=(<<n)-;i>=;i--){
         for(int j=;j<n;j++){
             if((<<j)&i){
                 f[i^(<<j)]=max(f[i^(<<j)],f[i]);
             }
         }
     }
     scanf("%d",&qq);
     while(qq--){
         scanf("%s",qr);
         len=strlen(qr);
         t=;
         for(int i=len-;i>=;i--)t=t*+qr[i]-'';
         printf("%d\n",f[t]);
     }
     return ;
 }