bzoj2764 [JLOI2011]基因补全

题目

大概又是一道dp题吧，把碱基先处理一下，就成为匹配方案数了，就是简单dp。

f[i]表示前j个的方案数。

转移就是：

$f[j]=f[j]+f[j-1](A[i]==B[j])$

注意枚举顺序，m要从后面开始。

比较麻烦的是，要用高精度才能过去。

#include<bits/stdc++.h>
using namespace std;
struct bigint{
    int len,data[2005];
    bigint operator += (const bigint a)
    {
        int i;
        for(i=0;i<len||i<a.len||data[i];i++)
            data[i]+=a.data[i],data[i+1]+=data[i]/10,data[i]%=10;
        len=i;
        return *this;
    }
    void output()
    {
        int i;
        printf("%d",data[len-1]);
        for(i=len-2;i>=0;i--)printf("%d",data[i]);
        printf("\n");
    }
};
int n,m;
char A[2005],B[2005];
bigint f[2005];
inline bool judge(char a,char b)
{
    if(a=='A'&&b=='T')return true;
    if(a=='T'&&b=='A')return true;
    if(a=='C'&&b=='G')return true;
    if(a=='G'&&b=='C')return true;
    return false;
}
inline void write(bigint A)
{
    for(int i=A.len;i>=1;i--)putchar(A.data[i]+'0');
    putchar('\n');
}
int main()
{
    freopen("in.txt","r",stdin);
    scanf("%d%d\n",&n,&m);
    scanf("%s\n%s",A+1,B+1);
    for(int i=0;i<=m;i++)f[i].len=1;
    for(int i=0;i<=m;i++)f[i].data[0]=0;
    f[0].data[0]=1;
    for(int i=1;i<=n;i++)
        for(int j=m;j>=1;j--)
            if(judge(A[i],B[j]))f[j]+=f[j-1];
    f[m].output();
    return 0;
}