2764: [JLOI2011]基因补全

题目链接

题目大意：给定一个长度为n的碱基序列S和一个长度为m的碱基序列T，现在希望向序列T里补一定的碱基使得序列S和序列T配对，求方案数

题解：$f[i][j]表示长串匹配到i，短串匹配到j的方案数$
$f[i][j]=f[i-1][j]，短串+一个字母$
$f[i][j]+=f[i-1][j-1]，长短串可以匹配$

大力高精度

我的收获：naive

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

#define rad 100000000
#define B 8

int n,m,p=1,l=0;

char s1[2005],s2[2005];

struct Bi{
    int s[105],l;
    Bi(){l=1;}
}f[2][2005];

Bi operator + (Bi a,Bi b){
    Bi c;memset(c.s,0,sizeof(c.s));
    c.l=a.l>b.l?a.l:b.l;
    for(int i=1;i<=c.l;i++){  
        c.s[i]+=a.s[i]+b.s[i];
        if(c.s[i]>=rad){c.s[i]-=rad;c.s[i+1]++;}
    }
    if(c.s[c.l+1]) c.l++;
    return c;
}

void print(Bi a){printf("%d",a.s[a.l]);for(int i=a.l-1;i;i--) printf("%08d",a.s[i]);}

bool judge(char a,char b){return (a=='A'&&b=='T')||(a=='G'&&b=='C')||(a=='C'&&b=='G')||(a=='T'&&b=='A');}

void work()
{
    f[0][0].s[1]=1,f[1][0].s[1]=1;
    for(int i=1;i<=n;i++,p^=1,l^=1)
        for(int j=1;j<=m;j++){
            f[p][j]=f[l][j];
            if(judge(s1[i],s2[j])) f[p][j]=f[p][j]+f[l][j-1];
        }
    print(f[l][m]);
}

void init(){
    scanf("%d%d%s%s",&n,&m,s1+1,s2+1);
}

int main()
{
    init();
    work();
    return 0;
}