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

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

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

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 50000+5

#define maxm 200000+5

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)

#define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)

#define mod 100000000

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int n,m,f[2][maxn],g[2][maxn];
char a[maxn],b[maxn];

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    scanf("%s",a+1);scanf("%s",b+1);n=strlen(a+1)-1;m=strlen(b+1)-1;
    for0(i,m)g[0][i]=1;//注意初值的设立
    int t;
    for1(i,n)
    {
        t=i&1;
        memset(f[t],0,sizeof(f[t]));//f[i][j]表示用a[1--i],b[1--j]能得到的最长公共子序列的长度
        memset(g[t],0,sizeof(g[t]));//g[i][j]表示用a[1--i],b[1--j]能得到的最长公共子序列的方案数
        g[t][0]=g[!t][0]=1;//注意初值的设立
        for1(j,m)
        if(a[i]==b[j])
        {
            f[t][j]=f[!t][j-1]+1;
            g[t][j]=g[!t][j-1];
            if(f[!t][j]==f[t][j])g[t][j]+=g[!t][j];
            if(f[t][j-1]==f[t][j])g[t][j]+=g[t][j-1];
                //这三种情况都是互相不包含的，所以要加法原理
            g[t][j]%=mod;
        }else
        {
            f[t][j]=max(f[!t][j],f[t][j-1]);
            if(f[!t][j]>f[t][j-1])g[t][j]=g[!t][j];
            else if(f[t][j-1]>f[!t][j])g[t][j]=g[t][j-1];
            else 
            {
                g[t][j]=g[!t][j]+g[t][j-1];
                if(f[!t][j-1]==f[t][j])g[t][j]-=g[!t][j-1];
                //如果f[!t][j-1]==f[t][j]说明g[!t][j]和g[t][j-1]都计算了g[!t][j-1]，所以算了两次，要去掉一次。
            }
            g[t][j]%=mod;
        }
    }
    cout<<f[t][m]<<endl<<(g[t][m]+mod)%mod<<endl;

    return 0;

}