首页 > 代码库 > BZOJ2423: [HAOI2010]最长公共子序列

BZOJ2423: [HAOI2010]最长公共子序列

题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2423

题解:方案数比较麻烦。写在注释里吧,

技术分享
  1 #include<cstdio>  2   3 #include<cstdlib>  4   5 #include<cmath>  6   7 #include<cstring>  8   9 #include<algorithm> 10  11 #include<iostream> 12  13 #include<vector> 14  15 #include<map> 16  17 #include<set> 18  19 #include<queue> 20  21 #include<string> 22  23 #define inf 1000000000 24  25 #define maxn 50000+5 26  27 #define maxm 200000+5 28  29 #define eps 1e-10 30  31 #define ll long long 32  33 #define pa pair<int,int> 34  35 #define for0(i,n) for(int i=0;i<=(n);i++) 36  37 #define for1(i,n) for(int i=1;i<=(n);i++) 38  39 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 40  41 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 42  43 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go) 44  45 #define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) 46  47 #define mod 100000000 48  49 using namespace std; 50  51 inline int read() 52  53 { 54  55     int x=0,f=1;char ch=getchar(); 56  57     while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();} 58  59     while(ch>=0&&ch<=9){x=10*x+ch-0;ch=getchar();} 60  61     return x*f; 62  63 } 64 int n,m,f[2][maxn],g[2][maxn]; 65 char a[maxn],b[maxn]; 66  67 int main() 68  69 { 70  71     freopen("input.txt","r",stdin); 72  73     freopen("output.txt","w",stdout); 74  75     scanf("%s",a+1);scanf("%s",b+1);n=strlen(a+1)-1;m=strlen(b+1)-1; 76     for0(i,m)g[0][i]=1;//注意初值的设立 77     int t; 78     for1(i,n) 79     { 80         t=i&1; 81         memset(f[t],0,sizeof(f[t]));//f[i][j]表示用a[1--i],b[1--j]能得到的最长公共子序列的长度 82         memset(g[t],0,sizeof(g[t]));//g[i][j]表示用a[1--i],b[1--j]能得到的最长公共子序列的方案数 83         g[t][0]=g[!t][0]=1;//注意初值的设立 84         for1(j,m) 85         if(a[i]==b[j]) 86         { 87             f[t][j]=f[!t][j-1]+1; 88             g[t][j]=g[!t][j-1]; 89             if(f[!t][j]==f[t][j])g[t][j]+=g[!t][j]; 90             if(f[t][j-1]==f[t][j])g[t][j]+=g[t][j-1]; 91                 //这三种情况都是互相不包含的,所以要加法原理 92             g[t][j]%=mod; 93         }else 94         { 95             f[t][j]=max(f[!t][j],f[t][j-1]); 96             if(f[!t][j]>f[t][j-1])g[t][j]=g[!t][j]; 97             else if(f[t][j-1]>f[!t][j])g[t][j]=g[t][j-1]; 98             else  99             {100                 g[t][j]=g[!t][j]+g[t][j-1];101                 if(f[!t][j-1]==f[t][j])g[t][j]-=g[!t][j-1];102                 //如果f[!t][j-1]==f[t][j]说明g[!t][j]和g[t][j-1]都计算了g[!t][j-1],所以算了两次,要去掉一次。103             }104             g[t][j]%=mod;105         }106     }107     cout<<f[t][m]<<endl<<(g[t][m]+mod)%mod<<endl;108 109     return 0;110 111 }  
View Code

 

BZOJ2423: [HAOI2010]最长公共子序列