【bzoj2423】最长公共子序列[HAOI2010](dp)
題目傳送門(mén):http://www.lydsy.com/JudgeOnline/problem.php?id=2423
題目大意:求兩個(gè)字符串的最長(zhǎng)公共子序列長(zhǎng)度和最長(zhǎng)公共子序列個(gè)數(shù)。
這道題的話,對(duì)于神犇來(lái)說(shuō),肯定是一眼看出狀態(tài)轉(zhuǎn)移方程的。而我這個(gè)蒟蒻,看了幾篇博客之后才看懂。。。
第一問(wèn)模板lcs,大家肯定都會(huì),就是設(shè)f[i][j]為A串跑到第i位,B串跑到第j為時(shí)的最長(zhǎng)公共子序列長(zhǎng)度,然后就有:
f[i][j]=f[i-1][j-1]+1 (a[i]==b[j])
?=max(f[i-1][j],f[i][j-1]) (a[i]!=b[j]) (原諒我不會(huì)編輯公式)
我就解釋一下第二問(wèn)的方程。先設(shè)g[i][j]為A串跑到第i位,B串跑到第j為時(shí)的最長(zhǎng)公共子序列個(gè)數(shù),方程就是這樣:
g[i][j]=g[i-1][j-1]
+g[i-1][j] (f[i][j]==f[i-1][j])
+g[i][j-1] (f[i][j]==f[i][j-1])
(a[i]==b[j])
=g[i-1][j] (f[i][j]==f[i-1][j])
+g[i][j-1] (f[i][j]==f[i][j-1])
-g[i-1][j-1] (f[i][j]==f[i-1][j-1])
(a[i]!=b[j])
這里當(dāng)a[i]和b[j]相同時(shí),g[i-1][j],g[i-1][j-1],g[i][j-1]這三個(gè)的最長(zhǎng)公共子序列不會(huì)重復(fù),因?yàn)檫@里的g[i-1][j-1]實(shí)際上還要在末尾添加上a[i](或b[j]),因此這些lcs全都是以a[i],b[j]結(jié)尾的,而g[i-1][j]不包含以a[i]結(jié)尾的lcs,g[i][j-1]不包含以b[j]結(jié)尾的lcs,因此這三類(lèi)lcs不會(huì)重復(fù),可以直接相加。
當(dāng)a[i]與b[j]不同時(shí),最后當(dāng)f[i][j]==f[i-1][j-1]時(shí)要減去g[i-1][j-1]就是因?yàn)檫@時(shí)g[i-1][j-1]被分別包含在g[i-1][j]和g[i][j-1]中,算了兩次,要把重復(fù)的減掉。
于是就可以愉快地AC這道題了。
還有,一定要用滾動(dòng)數(shù)組,不然爆!空!間!
丑代碼:
#include<cstdio> #include<cstring> #include<cmath> #include<algorithm> using namespace std; const int mod=100000000; int f[2][5010],g[2][5010]; char a[5010],b[5010]; int main() {int i,j,n,m;scanf("%s%s",a,b);n=strlen(a)-1; m=strlen(b)-1;for(i=0;i<=m;i++)g[0][i]=1; g[1][0]=1;for(i=1;i<=n;i++)for(j=1;j<=m;j++)if(a[i-1]==b[j-1]){f[i&1][j]=f[i&1^1][j-1]+1;g[i&1][j]=(g[i&1^1][j-1]+(f[i&1^1][j]==f[i&1][j])*g[i&1^1][j]+(f[i&1][j-1]==f[i&1][j])*g[i&1][j-1])%mod;}else{if(f[i&1^1][j]>f[i&1][j-1])f[i&1][j]=f[i&1^1][j];else f[i&1][j]=f[i&1][j-1];g[i&1][j]=((f[i&1^1][j]==f[i&1][j])*g[i&1^1][j]+(f[i&1][j-1]==f[i&1][j])*g[i&1][j-1]-(f[i&1][j]==f[i&1^1][j-1])*g[i&1^1][j-1]+mod)%mod;}printf("%d\n%d",f[n&1][m],g[n&1][m]); } View Code?
轉(zhuǎn)載于:https://www.cnblogs.com/quzhizhou/p/6978170.html
總結(jié)
以上是生活随笔為你收集整理的【bzoj2423】最长公共子序列[HAOI2010](dp)的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: AC日记——双栈排序 洛谷 P1155
- 下一篇: AC日记——[SDOI2010]大陆争霸