C语言 基因相关性

Description

为了获知基因序列在功能和结构上的相似性，经常需要将几条不同序列的DNA进行比对，以判断该比对的DNA是否具有相关性。

现比对两条长度相同的DNA序列。首先定义两条DNA序列相同位置的碱基为一个碱基对，如果一个碱基对中的两个碱基相同的话，则称为相同碱基对。接着计算相同碱基对占总碱基对数量的比例，如果该比例大于等于给定阈值时则判定该两条DNA序列是相关的，否则不相关。

Input

有三行，第一行是用来判定出两条DNA序列是否相关的阈值，随后2行是两条DNA序列（长度不大于500）。

Output

若两条DNA序列相关，则输出“yes”，否则输出“no”。

Sample Input

0.85
ATCGCCGTAAGTAACGGTTTTAAATAGGCC
ATCGCCGGAAGTAACGGTCTTAAATAGGCC

Sample Output

yes

Source Code

#include<stdio.h>
#include<string.h>
int main()
{
	double a;
	double ratio;	//计算相同碱基对占总碱基对数量的比例
	int count=0;	//统计相同碱基对的数量
	char b[500],c[500];
	scanf("%lf\n",&a);
	gets(b);
	gets(c);
	for(int i=0;i<strlen(b);i++)
	{
		if(b[i]==c[i])
			count++;
	}
	ratio=(count*1.0)/strlen(b);
	if(ratio>=a)
		printf("yes\n");
	else
		printf("no\n");
	return 0;
}

Computational Results

Hint

在scanf("%lf\n",&a);中，要加一换行符。
gets()函数：

• gets()函数用来从标准输入设备（键盘）读取字符串直至接受到换行符或EOF时停止结束，并将读取的结果存放在字符数组中，但换行符会被丢弃，然后在末尾添加’\0’字符。