BZOJ2423/HAOI2010 最长公共子序列

最新推荐文章于 2020-06-06 11:43:00 发布

sherlock_zhuang

最新推荐文章于 2020-06-06 11:43:00 发布

阅读量624

点赞数 1

分类专栏： my-code 文章标签： dp 最长公共子序列

本文链接：https://blog.csdn.net/sherlock_zhuang/article/details/50811307

版权

my-code 专栏收录该内容

28 篇文章 0 订阅

订阅专栏

思路：
非常经典的问题就不讲了。后面求方案的时候注意不要忘记讨论情况。
还有要滚动一维。


#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
using namespace std;
const int imax=5000+9;
const int mod=1e8;
char s1[imax],s2[imax];
int n,m; 
int f[2][imax],g[2][imax];

void iread()
{
    scanf(" %s",s1+1);
    scanf(" %s",s2+1);
    n=strlen(s1+1); m=strlen(s2+1); 
    n--; m--;
}

void iwork()
{
    for(int i=0;i<=m;i++) g[0][i]=1; 
    g[1][0]=1;

    int now,pre;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) 
        {
            now=i&1; pre=now^1;
            if(s1[i]==s2[j])
            {
                f[now][j]=f[pre][j-1]+1;
                g[now][j]=g[pre][j-1];
                if(f[now][j]==f[pre][j])   g[now][j]+=g[pre][j];
                if(f[now][j]==f[now][j-1]) g[now][j]+=g[now][j-1];

                g[now][j]=(g[now][j]+mod)%mod;  
            }
            else
            {
                f[now][j]=max(f[pre][j],f[now][j-1]);
                if(f[now][j-1]!=f[pre][j])
                {
                    if(f[now][j-1]>f[pre][j]) g[now][j]=g[now][j-1];
                    else g[now][j]=g[pre][j];   
                }
                else g[now][j]=(g[pre][j]+g[now][j-1]);
                if(f[pre][j-1]==f[now][j]) g[now][j]-=g[pre][j-1];

                g[now][j]=(g[now][j]+mod)%mod;
            }   
        }
    printf("%d\n",f[now][m]);
    printf("%d\n",g[now][m]);
}

int main()
{
    iread();
    iwork();
    return 0;
}