acwing 1776 牛的基因组学

Posted 霜序0.2℃

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了acwing 1776 牛的基因组学相关的知识,希望对你有一定的参考价值。

题目

农夫约翰拥有 NN 头有斑点的牛和 NN 头没有斑点的牛。

他刚刚完成了牛遗传学课程,他确信奶牛上的斑点是由牛基因组突变引起的。

农夫约翰花了大价钱对他的奶牛的基因组进行了测序。

每个基因组都是一个由四个字符 A,C,G,TA,C,G,T 构成的长度为 MM 的字符串。

当他统计得到的奶牛的基因组序列时,可以得到一个如下所示的表:(此时,N=3N=3)

位置 :    1 2 3 4 5 6 7 ... M

斑点牛 1: A A T C C C A ... T
斑点牛 2: G A T T G C A ... A
斑点牛 3: G G T C G C A ... A

普通牛 1: A C T C C C A ... G
普通牛 2: A C T C G C A ... T
普通牛 3: A C T T C C A ... T

通过仔细观察该表,他发现通过位置 22 的字符足以判断奶牛是否存在斑点。

也就是说,仅通过查看这个位置上的字符,农夫约翰就可以判断他的哪些奶牛有斑点,哪些没有斑点。(在这里,AA 和 GG 表示有斑点,CC 表示无斑点,TT 无关紧要,因为没有任何一头牛的第二个位置上的字符是 TT)

位置 11 上的字符不足以判断奶牛是否存在斑点,因为 AA 既可以表示有斑点也可以表示无斑点。

给定约翰的奶牛的基因组序列列表,请你计算可以单独用来判断奶牛是否存在斑点的位置的数量。

输入格式

第一行包含两个整数 NN 和 MM。

接下来 NN 行,每行包含一个长度为 MM 的字符串,用来描述斑点牛的基因组序列。

再接下来 NN 行,每行包含一个长度为 MM 的字符串,用来描述普通牛的基因组序列。

输出格式

输出可以单独用来判断奶牛是否存在斑点的位置的数量。

数据范围

1≤N,M≤1001≤N,M≤100

输入样例:

3 8
AATCCCAT
GATTGCAA
GGTCGCAA
ACTCCCAG
ACTCGCAT
ACTTCCAT

输出样例:

1

解释与代码

遍历斑点牛基因的每一列,如果普通牛的每一列出现与斑点牛相同的基因就结束循环,如果没有出现相同的ans+1

#include <bits/stdc++.h>
#define pb          push_back
#define ppb         pop_back
#define lbnd        lower_bound
#define ubnd        upper_bound
#define endl        '\\n'
#define all(a)      (a).begin(),(a).end()
#define what_is(x)  cerr << #x << " is " << x << endl;
#define ini(a)      memset(a,0,sizeof(a))
#define case        ll T;read(T);for(ll Q=1;Q<=T;Q++)
#define lowbit(x)   x&(-x)
#define pr          printf
#define sc          scanf
#define TIE \\
    cin.tie(0);cout.tie(0);\\
    ios::sync_with_stdio(false);

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double PI    = acos(-1.0);
const int    INF   = 0x3f3f3f3f;
const ll     LLINF = 0x3f3f3f3f3f3f3f3f;
const int    maxn  = 5e5+10;
const ll     N     = 5;

char ss[999][999], ss2[999][999];

void solve()
	int n, m, ans = 0, flag = 0;
	map<char,int> mci;
	cin>>n>>m;
	for (int i=0; i<n; i++) 
		cin>>ss[i];
	
	for (int i=0; i<n; i++) 
		cin>>ss2[i];
	
	for (int i=0; i<m; i++) 
		flag = 1;
		mci.clear();
		for (int j=0; j<n; j++) 
			mci[ss[j][i]]++;
		
		for (int j=0; j<n; j++) 
			if (mci[ss2[j][i]] != 0) 
				flag = 0;
                break;
			
		
		if (flag) ans++;
	
	cout<<ans<<endl;


int main()

//	TIE;

	solve();
//    casesolve();
//    casecout<<"Case "<<Q<<":"<<endl;solve();

以上是关于acwing 1776 牛的基因组学的主要内容,如果未能解决你的问题,请参考以下文章

功能基因组学的预测网络模型识别炎症性肠病的调节因子

2023-05-27 量学基础 十八个涨停基因(转)

Science综述 | 用单细胞基因组学将人类细胞表型匹配到基因型

全程回顾丨肿瘤基因组学研究进展及公共数据库数据挖掘利用研讨会

国家基因组科学数据中心工作好吗

《抗癌前线:癌症研究的前沿进展》:肿瘤基因组学博士的肿瘤科普,三星推荐