题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6138
题意:给n个串,每次询问x号串和y号串的最长公共子串的长度,这个子串必须是n个串中某个串的前缀
解法1:AC自动机。做法是把n个串建成AC自动机,前缀树中每个节点都当做结尾节点,val赋为trie树深度,然后把x串丢进自动机里,把匹配到的前缀节点染个色,再把y串丢进去,遇到同样颜色的前缀节点就更新一下答案。
#include <bits/stdc++.h> using namespace std; const int N = 1e5+10; const int M = 5e5+10; const int S = 26; struct AcAutomata{ int root,sz; int nxt[M][S],fail[M],val[M],col[N]; int newnode(){ val[sz] = col[sz] = 0; memset(nxt[sz], -1, sizeof(nxt[sz])); return sz++; } void init(){ memset(val, 0, sizeof(val)); sz = 0; root = newnode(); } void insert(char *s){ int u=root; int len=strlen(s); for(int i=0; i<len; i++){ int id=s[i]-‘a‘; if(nxt[u][id]==-1) nxt[u][id]=newnode(); val[nxt[u][id]]=val[u]+1; u=nxt[u][id]; } } void build(){ queue <int> q; fail[root] = root; for(int i=0; i<S; i++){ int &v = nxt[root][i]; if(~v){ fail[v] = root; q.push(v); } else{ v = root; } } while(q.size()){ int u = q.front(); q.pop(); for(int i = 0; i < S; i++){ int &v = nxt[u][i]; if(~v){ fail[v] = nxt[fail[u]][i]; q.push(v); } else{ v = nxt[fail[u]][i]; } } } } void update(char *s, int x){ int len = strlen(s); int u=root; for(int i=0; i<len; i++){ int id=s[i]-‘a‘; u=nxt[u][id]; int tmp=u; while(tmp){ col[tmp]=x; tmp=fail[tmp]; } } } int query(char *s, int x){ int len = strlen(s); int u = root; int ans = 0; for(int i=0; i<len; i++){ int id=s[i]-‘a‘; u=nxt[u][id]; int tmp=u; while(tmp){ if(col[tmp]==x) ans=max(ans, val[tmp]); tmp=fail[tmp]; } } return ans; } }ZXY; char s[N]; int pos[N]; int main() { int T,n,q; scanf("%d", &T); while(T--) { scanf("%d", &n); ZXY.init(); int d=1; for(int i=1; i<=n; i++){ pos[i]=d; scanf("%s", s+d); ZXY.insert(s+d); int len=strlen(s+d); d+=len+1; } ZXY.build(); scanf("%d", &q); int id=1; while(q--) { int x, y; scanf("%d%d",&x,&y); ZXY.update(s+pos[x],id); int ans = ZXY.query(s+pos[y],id); ++id; printf("%d\n", ans); } } return 0; }
解法2:KMP,直接枚举n个串做KMP。。
#include<bits/stdc++.h> using namespace std; const int maxn = 1e5 + 5; vector<int> Next[maxn]; string str[maxn]; void getnext(string &s, vector<int> &nxt) { int len = s.size(); nxt.resize(len); nxt[0] = -1; int i, j = -1; for(i = 1; i < len; i++) { while(j >= 0 && s[j + 1] != s[i]) j = nxt[j]; if(s[j + 1] == s[i]) j++; nxt[i] = j; } } int getMax(string &s, int strid) { int len = s.size(); int i, j = -1; int ret = 0; for(i = 0; i < len; i++) { while(j >= 0 && str[strid][j + 1] != s[i]) j = Next[strid][j]; if(str[strid][j + 1] == s[i]) j++; ret = max(ret, j + 1); } return ret; } char buf[maxn]; int main() { int T; scanf("%d", &T); while(T--) { int n; scanf("%d", &n); for(int i = 1; i <= n; i++) { scanf("%s", buf); str[i] = buf; getnext(str[i], Next[i]); } int q; scanf("%d", &q); while(q--) { int x, y; scanf("%d %d", &x, &y); int ans = 0; for(int i = 1; i <= n; i++) { int u = getMax(str[x], i); int v = getMax(str[y], i); ans = max(ans, min(u, v)); } printf("%d\n", ans); } } return 0; }
时间: 2024-10-13 09:24:14