Poj3294 Life Forms

题目传送门

哈哈哈广义SAM真的好简(du)单(liu)啊

到时候讲课可以拿来祸害众生，Yeah!

好了开始讲题解，我们将所有字符串加入广义SAM里面

对每一个节点维护一个bitset表示它在哪些主串中出现过，让后标记上传就用或运算就好了

因为题目要求输出方案，加上一个dfs就可以了，复杂度O(L*n/64)

（其实题目本质就是parent树上每个节点为一个集合，求一个子树的并的大小而已，金牌之王zjt给出了一种更优秀的做法，每次将具有相同元素的节点按照dfs序排好，让后在相邻两个节点的lca上面做上标记-1，标记节点操作可以在构建广义SAM的时候就可以完成，求lca可以用树剖，这样复杂度就是O(L lg L)，不过鉴于此题n<=100，上面bitset做法会更快）

当然，这道题也是SA的经典题，不过还是广义SAM比较好看~

#include<bitset>#include<stdio.h>#include<string.h>#include<algorithm>#define N 350010using namespace std;char c[N];bitset<110> sz[N];int s[N][26],mx[N],f[N],ans=0;int n,m,lst=1,cnt=1,v[N],r[N];inline int extend(int c){int p=lst,np=lst=++cnt,q,nq;mx[np]=mx[p]+1;for(;p&&!s[p][c];p=f[p]) s[p][c]=np;if(!p) return f[np]=1;q=s[p][c];if(mx[q]==mx[p]+1) f[np]=q;else{nq=++cnt;mx[nq]=mx[p]+1;f[nq]=f[q]; f[q]=f[np]=nq;memcpy(s[nq],s[q],26<<2);for(;p&&s[p][c]==q;p=f[p]) s[p][c]=nq;}}inline void Ex_extend(int c){int p=lst,q=s[p][c],nq;if(q){if(mx[p]+1==mx[q]) lst=q;else{lst=nq=++cnt;mx[nq]=mx[p]+1;f[nq]=f[q]; f[q]=nq;memcpy(s[nq],s[q],26<<2);for(;p&&s[p][c]==q;p=f[p]) s[p][c]=nq;}} else extend(c);sz[lst][m]=1;}inline void dfs(int x,int d){if(d!=mx[x] || d>ans) return;if(mx[x]==ans && sz[x].count()>m){ puts(c); return; }for(int i=0;i<26;++i)if(s[x][i]){ c[d]=i+'a'; dfs(s[x][i],d+1); c[d]=0; }}inline void init(){memset(s,0,sizeof s); memset(f,0,sizeof f);memset(v,0,sizeof v);ans=0;for(;cnt;--cnt) sz[cnt]=0;cnt=1;}int _18520(){ init();scanf("%d",&n); if(!n) return 0;for(m=0;m<n;++m){scanf("%s",c); lst=1;for(int i=0,l=strlen(c);i<l;++i) Ex_extend(c[i]-'a');} m=n>>1; memset(c,0,sizeof c);for(int i=1;i<=cnt;++i) ++v[mx[i]];for(int i=1;i<=cnt;++i) v[i]+=v[i-1];for(int i=cnt;i;--i) r[v[mx[i]]--]=i;for(int p,i=cnt;i;--i){p=r[i]; sz[f[p]]|=sz[p];if(sz[p].count()>m) ans=max(ans,mx[p]);}if(!ans) puts("?"); else dfs(1,0); puts(""); return 1;} int main(){ while(_18520()); }