CodeForces 520C DNA Alignment

标签：codeforces   贪心   

题意：

一段DNA序列（10^5长度）  定义h函数为两序列相同碱基个数  p函数为分别移动两个DNA序列后所有可能的h函数之和  问使p最大的序列有多少个

思路：

根据p函数的定义  我们发现p这个函数其实就是A序列每个碱基和B序列每个碱基比较再乘一个n

因此可以贪心构造B序列  即每次新加一个碱基必定是A序列中出现次数最多的碱基

那么最后的答案就是A序列中出现次数最多的碱基种类数的n次方

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<cstdlib>
#include<ctime>
#include<cmath>
using namespace std;
typedef long long LL;
#define N 100010
#define mod 1000000007

int n;
char s[N];
int f[10];

template<class T>
inline void RD(T &ret){
    char c;
    ret = 0;
    while ((c = getchar()) < '0' || c > '9');
    while (c >= '0' && c <= '9')
        ret = ret * 10 + (c - '0'), c = getchar();
}

LL quickpow(LL m , int n){
    LL ans = 1;
    while(n) {
        if(n&1) ans = (ans * m) % mod;
        n = n >> 1;
        m = (m * m) % mod;
    }
    return ans;
}

int main(){
    RD(n);
    scanf("%s",s);
    for(int i=0;i<n;i++){
        if(s[i]=='A') f[1]++;
        else if(s[i]=='G') f[2]++;
        else if(s[i]=='C') f[3]++;
        else if(s[i]=='T') f[4]++;
    }
    int m;
    m=max(max(f[1],f[2]),max(f[3],f[4]));
    int x=0;
    for(int i=1;i<=4;i++){
        if(f[i]==m) x++;
    }
    printf("%I64d\n",quickpow(x,n));
    return 0;
}

CodeForces 520C DNA Alignment

标签：codeforces   贪心   

原文地址：http://blog.csdn.net/houserabbit/article/details/44025505