Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include<bits/stdc++.h>
- using namespace std;
- typedef long long ll;
- typedef unsigned long long ull;
- typedef pair<ll,ll>pll;
- typedef pair<ll,pair<ll,ll>>plll;
- #define fastread() (ios_base:: sync_with_stdio(false),cin.tie(NULL));
- #define vll(v) v.begin(),v.end()
- #define all(x) x.rbegin(),x.rend()
- #define min3(a, b, c) min(a, min(b, c))
- #define max3(a, b, c) max(a, max(b, c))
- #define F first
- #define S second
- #define in freopen("input.txt", "r", stdin)
- #define out freopen("output.txt", "w", stdout)
- #define minheap int,vector<int>,greater<int>
- #define pb push_back
- #define eb emplace_back
- #define ischar(x) (('a' <= x && x <= 'z') || ('A' <= x && x <= 'Z'))
- #define isvowel(ch) ((ch=='a'||ch=='e'||ch=='i'||ch=='o'||ch=='u')||(ch=='A'|| ch=='E' || ch=='I'|| ch=='O'|| ch=='U'))
- #define bug cout<<"BUG"<<endl;
- const int Max = 2e5 + 10;
- const int Mod = 1e9 + 7;
- const double PI =3.141592653589793238463;
- bool compare(const string &a, const string &b)
- {
- for(ll i=0; i<a.size(); i++)
- {
- if(a[i]!=b[i] && (a[i]=='$' || b[i]=='$'))
- {
- if(a[i]=='$')return 1;
- else if(b[i]=='$')return 0;
- }
- else if(a[i]<b[i])return 1;
- else if(a[i]>b[i])return 0;
- }
- return 1;
- }
- ll lcm(ll a,ll b)
- {
- if(a==0 || b==0)return 0;
- return a/__gcd(a,b)*b;
- }
- ll ara[5][1000];
- ll number(char ch)
- {
- if(ch=='A')return 1;
- else if(ch=='C')return 2;
- else if(ch=='G')return 3;
- else return 4;
- }
- int main()
- {
- fastread();
- ll i,j,n,m,p,sum=0,k,t,cnt=0,q,l,r;
- bool flag=false;
- cin>>n;
- m=0;
- string str;
- for(i=0; i<n; i++)
- {
- cin>>str;
- m=str.size();
- for(j=0; j<m; j++)
- {
- ara[number(str[j])][j]++;
- }
- }
- ll ans=0;
- string temp="ACGT",Ans;
- for(i=0; i<m; i++)
- {
- ll mx=max(ara[1][i],max3(ara[2][i],ara[3][i],ara[4][i]));
- ans+=(n-mx);
- for(j=0; j<4; j++)
- {
- if(mx==ara[j+1][i])
- {
- Ans+=temp[j];
- break;
- }
- }
- }
- cout<<"score "<<ans<<"\nConsensus(Motifs) : "<<Ans<<endl;
- }
- /*
- 10
- TCGGGGGTTTTT
- CCGGTGACTTAC
- ACGGGGATTTTC
- TTGGGGACTTTT
- AAGGGGACTTCC
- TTGGGGACTTCC
- TCGGGGATTCAT
- TCGGGGATTCCT
- TAGGGGAACTAC
- TCGGGTATAACC
- */
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement