Description
字符串是oi界常考的问题。现在给定你n个字符串，询问每个字符串有多少子串（不包括空串）是所有n个字符串中
至少k个字符串的子串（注意包括本身）。
Input
第一行两个整数n，k。
接下来n行每行一个字符串。
n,k,l<=100000
Output
输出一行n个整数，第i个整数表示第i个字符串的答案。
Sample Input
3 1
abc
a
ab
Sample Output
6 1 3
思路
假设不管k的限制，答案应该就是在后缀自动机上面跑到每个点在parent树上到根路径上贡献的和
每个点的贡献就是(p->maxl-p->prt->maxl)
然后这个维护直接(o(n))递推就可以了
然后维护每个子串在哪些子串中出现过也很简单
虽然很暴力
直接建出广义后缀自动机之后再把每个串在后缀自动机上跑一遍
跑到的每个节点只需要更新它所有的prt节点就可以了
为了避免算重记录下当前节点上一次是被哪个标记访问的就可以了
#include<bits/stdc++.h>
using namespace std;
#define N 100010
#define fu(a,b,c) for(int a=b;a<=c;++a)
const int CHARSET_SIZE=26;
int n,k;
struct Sam{
Sam *ch[CHARSET_SIZE],*prt;
int maxl,val,pre;
Sam(int maxl=0):maxl(maxl){}
};
struct Suffix_Automaton{
Sam pool[N<<1],*cur,*root;
int buc[N];vector<Sam*> topo;
Suffix_Automaton(){cur=pool;root=new (cur++)Sam;}
Sam *extend(int c,Sam *last){
Sam *u=new (cur++)Sam(last->maxl+1),*v=last;
for(;v&&!v->ch[c];v=v->prt)v->ch[c]=u;
if(!v)u->prt=root;
else if(v->ch[c]->maxl==v->maxl+1){
u->prt=v->ch[c];
}else{
Sam *n=new (cur++)Sam(v->maxl+1),*o=v->ch[c];
copy(o->ch,o->ch+CHARSET_SIZE,n->ch);
n->prt=o->prt;
o->prt=u->prt=n;
for(;v&&v->ch[c]==o;v=v->prt)v->ch[c]=n;
}
return u;
}
void init(){
for(Sam *p=pool+1;p!=cur;p++){
if((signed)p->val<k)p->val=0;
else p->val=p->maxl-p->prt->maxl;
}
int maxv=0;
for(Sam *p=pool;p!=cur;p++){
maxv=max(maxv,p->maxl);
buc[p->maxl]++;
}
fu(i,1,maxv)buc[i]+=buc[i-1];
topo.resize(cur-pool);
for(Sam *p=pool;p!=cur;p++)topo[--buc[p->maxl]]=p;
fu(p,1,(signed)topo.size()-1)
topo[p]->val+=topo[p]->prt->val;
}
}sam;
struct Node{
Node *ch[CHARSET_SIZE];
};
struct Trie{
Node pool[N],*cur,*root;
Trie(){cur=pool;root=new (cur++)Node;}
Node *insert(Node *last,int c){
if(!last->ch[c])last->ch[c]=new (cur++)Node;
return last->ch[c];
}
void insert(char *s,int id){
int len=strlen(s);
Node *last=root;
fu(i,0,len-1)last=insert(last,s[i]-'a');
}
#define pi pair<Node*,Sam*>
void bfs(){
queue<pi> q;
q.push(pi(root,sam.root));
while(!q.empty()){
pi u=q.front();q.pop();
fu(i,0,25)if(u.first->ch[i]){
Sam *tmp=sam.extend(i,u.second);
q.push(pi(u.first->ch[i],tmp));
}
}
}
}trie;
char s[N];
int bg[N],en[N],c[N];
int main(){
freopen("input.txt","r",stdin);
scanf("%d%d",&n,&k);
int tot=0;
fu(i,1,n){
scanf("%s",s);
int len=strlen(s);
trie.insert(s,i);
bg[i]=tot+1;
fu(j,0,len-1)c[++tot]=s[j]-'a';
en[i]=tot;
}
trie.bfs();
fu(i,1,n){
Sam *now=sam.root;
fu(j,bg[i],en[i]){
now=now->ch[c[j]];
Sam *tmp=now;
for(;tmp!=sam.root&&tmp->pre!=i;tmp=tmp->prt)
tmp->pre=i,tmp->val++;
}
}
sam.init();
fu(i,1,n){
int ans=0;
Sam *now=sam.root;
fu(j,bg[i],en[i]){
for(;now&&!now->ch[c[j]];now=now->prt);
if(now){
now=now->ch[c[j]];
ans+=now->val;
}else now=sam.root;
}
printf("%d ",ans);
}
return 0;
}