BZOJ 2434 最长公共子序列

int f[2][N],g[2][N];
char s[N],t[N];
int n,m;




void add(int &x,int y)
{
    x+=y;
    x%=mod;
}


int main()
{
    RD(s); n=strlen(s)-1; 
    RD(t); m=strlen(t)-1; 
    int i,j;
    int pre=0,cur=1;
    for(i=0;i<=m;i++) g[0][i]=1;
    g[pre][0]=1;
    FOR1(i,n) 
    {
        clr(f[cur],0);
        clr(g[cur],0);
        g[cur][0]=1;
        FOR1(j,m)
        {
            if(s[i-1]==t[j-1])
            {
                f[cur][j]=f[pre][j-1]+1;
                g[cur][j]=g[pre][j-1];
                if(f[cur][j-1]==f[cur][j]) add(g[cur][j],g[cur][j-1]);
                if(f[pre][j]==f[cur][j]) add(g[cur][j],g[pre][j]);
            }
            else
            {
                f[cur][j]=max(f[cur][j-1],f[pre][j]);
                if(f[cur][j-1]==f[cur][j]) add(g[cur][j],g[cur][j-1]);
                if(f[pre][j]==f[cur][j]) add(g[cur][j],g[pre][j]);
                if(f[pre][j-1]==f[cur][j]) add(g[cur][j],-g[pre][j-1]);
            }
        }
        swap(pre,cur);
    }
    if(g[pre][m]<0) g[pre][m]+=mod;
    PR(f[pre][m]);
    PR(g[pre][m]);
}