bzoj2764 基因补全

标签：

Description

在生物课中我们学过，碱基组成了DNA（脱氧核糖核酸），他们分别可以用大写字母A,C,T,G表示，其中A总与T配对，C总与G配对。两个碱基序列能相互匹配，当且仅当它们等长，并且任意相同位置的碱基都是能相互配对的。例如ACGTC能且仅能与TGCAG配对。一个相对短的碱基序列能通过往该序列中任意位置补足碱基来与一个相对长的碱基序列配对。补全碱基的位置、数量不同，都将视为不同的补全方案。现在有两串碱基序列S和T，分别有n和m个碱基(n>=m)，问一共有多少种补全方案。

Input

数据包括三行。

第一行有两个整数n，m，表示碱基序列的长度。

第二行包含n个字符，表示碱基序列S。

第三行包含m个字符，表示碱基序列T。

两个碱基序列的字符种类只有A,C,G,T这4个大写字母。

Output

答案只包含一行，表示补全方案的个数。

dp+高精度

把第二个串换为与其配对的串

若s1[i]==s2[j]则 f[i][j]=f[i-1][j]+f[i-1][j-1]

否则 f[i][j]=f[i-1][j]

使用滚动数组以节约空间

#include<cstdio>
#include<cstring>
int n,m;
char s1[2002],s2[2002];
char ch[256];
int f[2002][70];
int main(){
    scanf("%d%d",&n,&m);
    scanf("%s%s",s1+1,s2+1);
    ch[‘A‘]=‘T‘;ch[‘T‘]=‘A‘;ch[‘G‘]=‘C‘;ch[‘C‘]=‘G‘;
    for(int i=1;i<=m;i++)s2[i]=ch[s2[i]];
    f[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=m;j;j--){
            if(s1[i]==s2[j]){
                int*w=f[j];
                for(int x=0;x<70;x++)w[x]+=f[j-1][x];
                for(int x=1;x<70;x++){
                    if(w[x-1]>=1000000000)
                    ++w[x],w[x-1]-=1000000000;
                }
            }
        }
    }
    bool d=0;
    for(int i=69;i>=0;i--){
        if(d)printf("%09d",f[m][i]);
        if(f[m][i]&&!d)printf("%d",f[m][i]),d=1;
    }
    if(!d)puts("0");
    return 0;
}

bzoj2764 基因补全

标签：

原文地址：http://www.cnblogs.com/ccz181078/p/5136168.html