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

  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 }