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;
}