码迷,mamicode.com
首页 > 其他好文 > 详细

C. DNA Alignment 数学公式推导 Codeforces Round #295 (Div. 2)

时间:2015-03-03 09:59:49      阅读:141      评论:0      收藏:0      [点我收藏+]

标签:c. dna alignment   数学公式推导   cf round #295 div. 2   

C. DNA Alignment
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output

Vasya became interested in bioinformatics. He‘s going to write an article about similar cyclic DNA sequences, so he invented a new method for determining the similarity of cyclic sequences.

Let‘s assume that strings s and t have the same length n, then the function h(s,?t) is defined as the number of positions in which the respective symbols of s and t are the same. Function h(s,?t) can be used to define the function of Vasya distance ρ(s,?t):

技术分享
where 技术分享 is obtained from string s, by applying left circular shift i times. For example,
ρ("AGC",?"CGT")?=?
h("AGC",?"CGT")?+?h("AGC",?"GTC")?+?h("AGC",?"TCG")?+?
h("GCA",?"CGT")?+?h("GCA",?"GTC")?+?h("GCA",?"TCG")?+?
h("CAG",?"CGT")?+?h("CAG",?"GTC")?+?h("CAG",?"TCG")?=?
1?+?1?+?0?+?0?+?1?+?1?+?1?+?0?+?1?=?6

Vasya found a string s of length n on the Internet. Now he wants to count how many strings t there are such that the Vasya distance from the string s attains maximum possible value. Formally speaking, t must satisfy the equation: 技术分享.

Vasya could not try all possible strings to find an answer, so he needs your help. As the answer may be very large, count the number of such strings modulo 109?+?7.

Input

The first line of the input contains a single integer n (1?≤?n?≤?105).

The second line of the input contains a single string of length n, consisting of characters "ACGT".

Output

Print a single number — the answer modulo 109?+?7.

Sample test(s)
input
1
C
output
1
input
2
AG
output
4
input
3
TTT
output
1
Note

Please note that if for two distinct strings t1 and t2 values ρ(s,?t1) и ρ(s,?t2) are maximum among all possible t, then both strings must be taken into account in the answer even if one of them can be obtained by a circular shift of another one.

In the first sample, there is ρ("C",?"C")?=?1, for the remaining strings t of length 1 the value of ρ(s,?t) is 0.

In the second sample, ρ("AG",?"AG")?=?ρ("AG",?"GA")?=?ρ("AG",?"AA")?=?ρ("AG",?"GG")?=?4.

In the third sample, ρ("TTT",?"TTT")?=?27


题意:给一个由A,T,C,G组成的字符串s,求使ρ(s,t)最大的t的个数。易知,t应该由字符串s中个数最多的字母组成(总数最多的字母只有一个的时候很好理解,如果有多个字母的个数相同,那么就可以将这些个数相同的字母看成一个字母,这样是等价的),求出字母个数最多的种类数x,根据排列组合,答案就是x^n

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#include <map>
#include <stack>
#include <vector>
#include <set>
#include <queue>
#pragma comment (linker,"/STACK:102400000,102400000")
#define maxn 100005
#define MAXN 2005
#define mod 1000000007
#define INF 0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-6
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
#define FRE(i,a,b)  for(i = a; i <= b; i++)
#define FRL(i,a,b)  for(i = a; i < b; i++)
#define mem(t, v)   memset ((t) , v, sizeof(t))
#define sf(n)       scanf("%d", &n)
#define sff(a,b)    scanf("%d %d", &a, &b)
#define sfff(a,b,c) scanf("%d %d %d", &a, &b, &c)
#define pf          printf
#define DBG         pf("Hi\n")
typedef long long ll;
using namespace std;

__int64 n;
char str[maxn];
__int64 num[30];

__int64 pow_m(__int64 a,__int64 n)
{
    __int64 ret=1;
    __int64 tmp=a%mod;
    while (n)
    {
        if (n&1) ret=(ret*tmp)%mod;
        tmp=tmp*tmp%mod;
        n>>=1;
    }
    return ret;
}

int main()
{
    while (~scanf("%I64d",&n))
    {
        scanf("%s",str);
        mem(num,0);
        for (__int64 i=0;i<n;i++)
            num[ str[i]-'A' ]++;
        __int64 maxx=0;
        for (__int64 i=0;i<26;i++)
            maxx=max(maxx,num[i]);
        __int64 m=0;
        for (__int64 i=0;i<26;i++)
            if (num[i]==maxx) m++;
        printf("%I64d\n",pow_m(m,n));
    }
    return 0;
}



C. DNA Alignment 数学公式推导 Codeforces Round #295 (Div. 2)

标签:c. dna alignment   数学公式推导   cf round #295 div. 2   

原文地址:http://blog.csdn.net/u014422052/article/details/44023967

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!