【bzoj2764】[JLOI2011]基因补全 dp+高精度

#include <cstdio>
#define mod 100000000
char s1[2010] , s2[2010];
struct data
{
    int len , num[100];
    data operator+=(const data a)
    {
        int i;
        for(i = 0 ; i < len || i < a.len || num[i] ; i ++ )
            num[i] += a.num[i] , num[i + 1] += num[i] / mod , num[i] %= mod;
        len = i;
        return *this;
    }
    void output()
    {
        int i;
        printf("%d" , num[len - 1]);
        for(i = len - 2 ; i >= 0 ; i -- )
            printf("%08d" , num[i]);
        printf("\n");
    }
}f[2010];
bool judge(char a , char b)
{
    return (a == ‘A‘ && b == ‘T‘) || (a == ‘G‘ && b == ‘C‘) || (a == ‘C‘ && b == ‘G‘) || (a == ‘T‘ && b == ‘A‘);
}
int main()
{
    int n , m , i , j;
    scanf("%d%d%s%s" , &n , &m , s1 + 1 , s2 + 1);
    f[0].len = f[0].num[0] = 1;
    for(i = 1 ; i <= n ; i ++ )
        for(j = m ; j ; j -- )
            if(judge(s1[i] , s2[j]))
                f[j] += f[j - 1];
    f[m].output();
    return 0;
}