c++基因相关性题目为了获知基因序列在功能和结构上的相似性,经常需要将几条不同序列的DNA进行比对,以此来判断该比对的DNA是否具有相关性.现假设比对两条长度相同的DNA序列,当其中对应位

Question

c++基因相关性题目
为了获知基因序列在功能和结构上的相似性,经常需要将几条不同序列的DNA进行比对,以此来判断该比对的DNA是否具有相关性.
现假设比对两条长度相同的DNA序列,当其中对应位上的相同碱基的比率大于某给定值时则断定该两条DNA序列是相关的,否则不相关.
关于输入
有三行,第一行是用来分界出该两条DNA序列是否相关的比率,随后2行是两条DNA序列（长度不大于500）.
关于输出
若两条DNA序列相关,则输出“yes”；若不相关,则输出“no”；若输入的序列中有A、T、C、G以外的字母出现时或两条序列长度不相等时,则输出“error”.
例子输入
0.85
ATCGCCGTAAGTAACGGTTTTAAATAGGCC
ATCGCCGGAAGTAACGGTCTTAAATAGGCC
例子输出
yes
#include
#include
int main(){
\x09int i;
\x09double n;
\x09int g=0;
\x09char a[500],b[500];
\x09scanf("%d",&n);
\x09scanf("%s%s",&a,&b);
\x09if(strlen(a)!=strlen(b)){
\x09\x09printf("error");
\x09}
\x09else{
\x09\x09for(i=0;a[i]!='\0';i++){
\x09\x09\x09if(a[i]=b[i]){
\x09\x09\x09g++;}
\x09\x09}
\x09\x09if(g/(strlen(a)*1.0)>n){
\x09\x09printf("yes");
\x09\x09}
\x09\x09else{printf("no");}
\x09}
\x09return 0;
}
结果都是输出“error”,是不是长度比较的时候错了?还有“若输入的序列中有A、T、C、G以外的字母出现时”这部分要怎么写.

Answer 1

#include 
#include 
int main(){
\x09int i;
\x09float n;
\x09int g=0;
\x09char a[500], b[500];
\x09
\x09memset(a, 0, 500); /*初始化a b 数组,避免内存中其他程序遗留的垃圾数据 */
\x09memset(b, 0, 500);
\x09
\x09scanf("%f",&n);
\x09scanf("%s", a);
\x09scanf("%s", b);
\x09
\x09if(strlen(a)!=strlen(b)){
\x09\x09printf("error\n");
\x09}
\x09else{
\x09\x09for(i=0; a[i]!='\0'; i++){
\x09\x09\x09if(!strchr("ATGC", a[i]) ){
\x09\x09\x09    printf("error\n");
\x09\x09\x09    return 1;
\x09\x09\x09} else if( a[i] == b[i] ) {
\x09\x09\x09    g++;
\x09\x09\x09}
\x09\x09}
\x09\x09if(g / (strlen(a) * 1.0) > n){
\x09\x09\x09printf("yes\n");
\x09\x09}
\x09\x09else{printf("no\n");}
\x09}
\x09return 0;
}\x1a

程序出错的位置

数组名就是指针,因此在scanf 函数中 不能用&a, &b,而应该用a, b,这是错误的关键；

if( a[i] == b[i])时,用=代替了==是一处错误,而且没有达到设计目标,即非ATGC打印error.

将n设置成float类型, scanf 输入n 的值用%f,而不是%d.