[KMP][倍增求LCA]JZOJ 4669 弄提纲
分析
字符串前缀,我们可以想到KMP
然后我们就能想到,l的失配数组往前跳到第一个和r的适配数组往前跳到相同的,那个点就是它们最长公共前缀的点,然后能往前跳多少次就是公共前缀的个数
由于失配数组是一个数组,可以构成一个树的模型,就可以抽象成求LCA和LCA的深度了
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include <iostream> #include <cstdio> #include <cstring> #include <string> using namespace std; const int N=3e4+10; struct Edge { int u,v,nx; }g[N]; int list[N],cnt; string s; int m,len,f[N][20],dep[N]; void KMP() { int next=0; f[1][0]=0;g[++cnt]=(Edge){0,1,list[0]};list[0]=cnt; for (int i=1;i<len;i++) { while (next!=0&&s[next]!=s[i]) next=f[next][0]; if (s[next]==s[i]) next++; f[i+1][0]=next; g[++cnt]=(Edge){next,i+1,list[next]};list[next]=cnt; } } void DFS(int u) { for (int i=list[u];i;i=g[i].nx) dep[g[i].v]=dep[u]+1,DFS(g[i].v); } int LCA(int a,int b) { if (dep[a]<dep[b]) swap(a,b); for (int i=19;i>=0;i--) if (dep[f[a][i]]>=dep[b]) a=f[a][i]; if (a==b) return a; for (int i=19;i>=0;i--) if (f[a][i]!=f[b][i]) a=f[a][i],b=f[b][i]; return f[a][0]; } int main() { cin>>s; len=s.length(); KMP(); DFS(0); for (int k=1;k<20;k++) for (int i=0;i<=len;i++) f[i][k]=f[f[i][k-1]][k-1]; for (scanf("%d",&m);m;m--) { int l,r; scanf("%d%d",&l,&r); int lca=LCA(l,r); printf("%d %d\n",dep[lca],lca); } }
在日渐沉没的世界里,我发现了你。