病毒感染检测问题

医学研究者最近发现了某些新病毒,通过对这些病毒的分析,得知他们的DNA序列都是环状的。现在研究者已收集了大量的病毒DNA和人的DNA数据,想快速检测出这些人是否感染了相应的病毒。为了方便研究,研究者将人的DNA和病毒DNA均表示成由一些字母组成的字符串序列,然后检测某种病毒DNA序列是否在患者的DNA序列中出现过,如果出现过,这此人感染了该病毒,否则没有感染。例如,假设病毒的DNA序列为baa,患者1的DNA序列为aaabbba,则感染。患者2的DNA序列为babbba,则未感染。(注意:人的DNA序列是线性的,而病毒的DNA序列是环状的)。

研究者将待检测的数据保存在一个文本文件中,文件格式和内容规定如下:文件有num+1行,第一行有一个整数num,表示有num个待检测的任务(num<=300)。接下来每行i(2<=i<=num+1)对应一个任务,每行有两个数据,用空格分隔,第一个数据表示病毒的DNA序列(长度<=6000),第二个数据表示人的DNA序列(长度<=10000)。

要求将检测结果输出到文件中,文件中包括num行,每行有三个数据,用空格分隔,前两个数据分别表示输入文件中对应病毒的DNA序列和人的DNA序列,如果该人感染了对应的病毒,该行第三个数据则为“YES”,否则为“NO”。

 

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
char virus[200];
char DNA[200];
//BF算法 
bool Index(char *S,char *T)
{
    int i=0,j=0;
    while(i<strlen(S) && j<strlen(T))
    {
        if(S[i]==T[j])
        {
            i++;j++;
        }
        else
        {
            i=i-j+1;
            j=0; 
        }
    }
    if(j>=strlen(T))return true;
    else return false;
}
//字符串旋转,把病毒第一个字符变为最后一个字符 
char *turn(char *T)
{
    char *temp=T;
    char tem[2]={temp[0]};
    strcpy(T,temp+1);
    strcat(T,tem);
    return T;
} 
//判断DNA是否感染 
bool judge(char *S,char *T)
{
    int len=strlen(T);
    for(int i=0;i<len;i++)
    {
        if(i!=0)                    //第一次病毒串不需要旋转 
            strcpy(T,turn(T));
        if(Index(S,T))return true;
    }
    return false;
}
int main()
{
    FILE *fp1,*fp2;
    //打开文件,没有则创建,打开失败退出程序 
    if((fp1=fopen("待检测数据-2018-10-10.txt","r"))==NULL)
    {
        printf("can't open the file\n");
        exit(0);
    }
    if((fp2=fopen("检测结果数据-2018-10-10.txt","w"))==NULL)
    {
        printf("can't open the file\n");
        exit(0);
    }
    int t;          //接受测试数据数目 
    fscanf(fp1,"%d",&t);
    fprintf(fp2,"%d\n",t);
    if(t<=0)                                //测试数据<0,结束 
        return 0;
    //t组测试数据 
    while(t--)
    {
        //从文件读取virus和DNA 
        fscanf(fp1,"%s",virus);
        fscanf(fp1,"%s",DNA);
        //virus和DNS输出到另一个文件中 
        fprintf(fp2,"%s\t",virus);
        fprintf(fp2,"%s\t",DNA);
        //判断是否被感染 
        if(judge(DNA,virus))
            fprintf(fp2,"YES",virus);
        else
            fprintf(fp2,"NO",virus);
        fprintf(fp2,"%\n",virus);
    }
    fclose(fp1);
	fclose(fp2);
	return 0;
}

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值