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