SXOJ844 牛的基因组

最新推荐文章于 2021-06-15 07:00:00 发布

Legend4ever

最新推荐文章于 2021-06-15 07:00:00 发布

阅读量725

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/Legend4ever/article/details/78318720

本文介绍了一种关于奶牛基因匹配的问题解决方法，通过对基因片段进行数值化转换并采用三位4进制数组进行优化，实现了对可能基因组合的有效查找。文章分享了从输入处理到算法实现的全过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://119.23.35.203/oj/contest/194/problems/2

题目的解述其实有点混乱，实际上就是在前n行有污点的奶牛任选三个序号，

对应到后n行没有污点的奶牛的该序号上看是否相同，如果没有一只无污点的奶牛的基因片段与之匹配，

可能的方案数就要增加1.

该题目数据比较水，暴力枚举的话也有人能够拿到80分，但是如果数据较大的时候，就需要用到数据优化来减少时间复杂度。

因为基因片段只有A,C,G,T四种可能性，所以可以将每只牛的基因片段转化为一个4进制的数串，

分别用0,1,2,3对应A,C,G,T。

再开一个三位的4进制数组作为判断该基因组合在无污点奶牛对应位置上是否有出现就可以了。

在写代码的过程中，由于对字符串输入的极度不熟练，在这上面卡了很长时间。。因为不想用cin这个写法

一度写出了scanf("%c",&wd[i][j])导致读入数据都没完整读进来，最后开了char的二维数组，

for(int i=0;i<n;i++)
	scanf("%s",&wd[i]);

使用以上写法成功读入了全部是数据实在是惭愧

通过这道题目，暴露出了自己对输入都还未能熟练掌握，还是任重而道远啊。

下面是满分通过的代码：

#include<cstdio>
#include<string>
#include<iostream> 
using namespace std;
char wd[510][55],no[510][55];
int WD[510][55],NO[510][55],n,m,o,p,q;
int jud[64];
bool judge(int x,int y,int z)    //x,y,z为三个基因片段对应的下标位置
{
	bool flag=true;
	for(int i=0;i<n;i++) //对于前n行有污点的奶牛,设x,y,z位置的基因片段组合为形成污点的可能方案,并压缩成三位4进制数标记下来
		jud[WD[i][x]*16+WD[i][y]*4+WD[i][z]]=1;
	for(int i=0;i<n;i++)
		if(jud[NO[i][x]*16+NO[i][y]*4+NO[i][z]]==1)  //判断无污点奶牛的对应位置上的基因是否相同，同则为false
			flag=false;
	for(int i=0;i<n;i++)
		jud[WD[i][x]*16+WD[i][y]*4+WD[i][z]]=0;      //恢复标记
	return flag;
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=0;i<n;i++)
	{
		scanf("%s",&wd[i]);
	}
	for(int i=0;i<n;i++)
	{
		scanf("%s",&no[i]);
	}
	for(int i=0;i<n;i++)                           //优化奶牛基因，把A,C,G,T分别对应为0,1,2,3，优化时间复杂度
		for(int j=0;j<m;j++)
		{
			if(wd[i][j]=='A') WD[i][j]=0;
			if(wd[i][j]=='C') WD[i][j]=1;
			if(wd[i][j]=='G') WD[i][j]=2;
			if(wd[i][j]=='T') WD[i][j]=3;
		}
	for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
		{
			if(no[i][j]=='A') NO[i][j]=0;
			if(no[i][j]=='C') NO[i][j]=1;
			if(no[i][j]=='G') NO[i][j]=2;
			if(no[i][j]=='T') NO[i][j]=3;
		}
	int ans=0;
	for(o=0;o<m;o++)
		for(p=o+1;p<m;p++)
			for(q=p+1;q<m;q++)       //枚举所有三个基因组合，并判断是否为可能方案
				if(judge(o,p,q)==true)
					ans++;
	printf("%d",ans);
}

记录自己走过的OI之路，希望日后自己能越来越强！