AHCPC2021本科生组 题I 病毒


一、题目

地球即将毁灭,科学家们派遣了一支队伍前往银河系外进行调查,寻找未来人类能居住的行星。幸运的是,在一个星球上,科学家们搜集到了许多 RNA 片段,科学家们把这些 RNA 片段运送往地球的实验基地中进行研究。经过几周的研究,他们发现这些 RNA 片段中存在着许多未知的病毒!
每个 RNA 片段都由 ACT 三种碱基组成,经过研究发现,如果一个 RNA 能够翻译出一个 AACCTT 无限循环的蛋白质序列,那么这个 RNA 就存在病毒!RNA 翻译蛋白质时,首先会随机选择一个位置的碱基,将该碱基翻译为蛋白质,然后继续随机选择相邻位置的碱基进行翻译,直到结束。现在给出一个 RNA,你需要判断这个RNA 是否存在病毒。

输入:
第一行两个整数N和 M,N 表示 RNA 中所含碱基总数,M 表示相邻碱基的对数。接下来一行,给出一个长 N 的字符串,表示每个位置碱基的种类。接下来 M 行, 每行两个整数 ai,bi,表示位置 ai 的碱基和位置 bi 的碱基相邻。
1 <= N, M <= 105
1 <= ai, bi <= N
保证字符串中只含 ACT。

输出:
如果存在病毒,输出“Yes”,否则输出“No”


提示:以下是本篇文章正文内容,下面代码仅供参考。

二、解题思路

首先想到使用邻接矩阵存储碱基的相邻关系,然后通过调用递归函数实现病毒特征的循环序列查找,找到返回Yes,否则返回No。

具体解法

(1)创建全局变量N、M,创建一维动态数组a[]储存每个位置碱基的种类,建立动态二维数组b[],通过邻接矩阵储存相邻关系, char c[8] = “AACCTTA”,存储题目给出的AACCTT无限循环的蛋白质序列,并且以上对象作用域均为全局;
(2)第一行输入两个整数 N 和 M (以空格隔开),N 表示 RNA 中所含碱基总数,M 表示相邻碱基的对数,通过cin读入;根据N的大小初始化动态二维数组b[],并赋值为0;
(3)cin不带走分隔符,故用cin.get()清除缓冲区里的分隔符,避免cin.getline()读空;输入一个长 N 的字符串,表示每个位置碱基的种类(保证字符串中只含 ACT),读取到a[];
(4)接下来 M 行每行输入两个整数 ai,bi(以空格隔开),表示位置 ai 的碱基和位置 bi 的碱基相邻,将相邻关系存储到邻接矩阵b[],赋值为1(b[ai-1][bi-1] = 1 b[bi-1][ai-1] = 1);
(5)for循环找到a[]里的每一个‘A’,对应循环序列的第一个值‘A’,调用递归函数find(int i, int t)实现循环序列查找(i是当前在a[]中的位置;t在c[]中的位置,表示查找进度),当t=7时说明循环序列已找到,向上返回true,输出结果“Yes”,若循环到最后一遍结束前仍未返回true,输出“No”。

三、源代码

代码如下(示例):

#include<iostream>
using namespace std;

int N, M; //全局变量
char* a = new char[100]; //创建动态数组储存每个位置碱基的种类
int** b = new int* [100]; //建立动态二维数组,通过邻接矩阵储存相邻关系
char c[8] = "AACCTTA"; //AACCTT无限循环的蛋白质序列

bool find(int i, int t) { //递归函数实现循环序列查找(i是当前在a[]中的位置;t在c[]中的位置,表示查找进度)
	bool x = false;
	if (t == 7) { //说明c[6]已经找到,有循环序列,向上返回
		return true;
	}
	for (int j = 0; j < N; j++) {
		if (a[j] == c[t] && b[i][j] == 1) {
			x=find(j, t + 1); //逐层向上返回
			return x; 
		}
	}  
}

int main() {
	cout << "第一行输入两个整数 N 和 M (以空格隔开),N 表示 RNA 中所含碱基总数,M 表示相邻碱基的对数"<<endl;
	cin >> N;
	cin >> M;

	for (int i = 0; i < N; i++) { //初始化动态二维数组
		b[i] = new int[N];
		for (int j = 0; j < N; j++) { b[i][j] = 0; }
	}

	cin.get(); //cin不带走分隔符,故清除缓冲区里的分隔符,避免getline读空
	cout << "输入一个长 N 的字符串,表示每个位置碱基的种类(保证字符串中只含 ACT)" << endl;
	cin.getline(a,N+1); //读取输入到a[]
	
	int ai, bi;
	cout << "接下来 M 行每行输入两个整数 ai,bi(以空格隔开),表示位置 ai 的碱基和位置 bi 的碱基相邻" << endl;
	for (int i = 0; i < M; i++) { //将相邻关系存储到邻接矩阵b[]
		cin >> ai; cin >> bi;
		b[ai-1][bi-1] = 1;
		b[bi-1][ai-1] = 1;
	}

	for (int i = 0; i < N; i++) { //找到a[]里的每一个‘A’,对应循环序列的第一个值‘A’,然后调用递归函数
		if (a[i] == 'A') {
			if (find(i, 1) == true) {
				cout << "Yes";
				break;
			}
		}
		if (i = N - 1) { //若循环到最后一遍结束前仍未返回true,输出“No”
			cout << "No";
		}
	}
}

四、运行样例

(1)
在这里插入图片描述
(2)
在这里插入图片描述
(3)
在这里插入图片描述
(4)
在这里插入图片描述

五、总结

本篇文章是作者第一次发布文章,文章肯定存在一定的错误和疏漏,代码和算法也一定不是最优的,欢迎大家评论交流,批评指正。

  • 4
    点赞
  • 14
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值