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 牛的基因组学的主要内容,如果未能解决你的问题,请参考以下文章
Science综述 | 用单细胞基因组学将人类细胞表型匹配到基因型