DNA Alignment CodeForces - 520C (思维+快速幂)

最新推荐文章于 2024-07-10 08:38:52 发布

原创最新推荐文章于 2024-07-10 08:38:52 发布 · 1k 阅读

2 ·

CC 4.0 BY-SA版权

思维技巧同时被 2 个专栏收录

181 篇文章

订阅专栏

数论

128 篇文章

订阅专栏

本文介绍了一种新的DNA序列相似性评估方法，并提供了一个具体的算法实现案例。该方法通过计算两个周期性DNA序列之间的相似度来确定它们的接近程度。

DNA Alignment

CodeForces - 520C

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):

$\text{[math]}$ where $\text{[math]}$ 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: $\text{[math]}$ .

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 10⁹ + 7.

Input

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

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 10⁹ + 7.

Examples

Input

1
C

Output

Input

2
AG

Output

Input

3
TTT

Output

Note

Please note that if for two distinct strings t₁ and t₂ values ρ(s, t₁) и ρ(s, t₂) 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

题意：有一个h（s，t）函数，这个函数的函数值是字符串s和t中字符相同的位置的个数。另一个函数

$\text{[math]}$

其实就是二重循环内循环字符串s不变，字符串t第一个元素每次移到最后一个，然后求h（），外循环s每次首字符移到最后然后再进行内循环，求这些h（）函数总和。就如题目中给的例子

ρ("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

题目让求有多少种t使得这个p（）最大

思路：这道题一开始就被题目吓到了，一直到比赛结束也没有想到思路。搜题解发现这是一个找规律的题目，题目让求有多少种t使得这个p（）最大，首先我们的任务就是要知道在知道s的情况下p（）的最大值是多少

观察发现它的最大值只取决于s字符串中某个出现次数最多的字符个数，为什么这么说呢？

举个例子AAACGT这个字符串我们发现次数最多的是A，那么对于

s：A A A C G T 最大的情况一定是

t：A A A A A A 因为每一个A都可以匹配三次，那么也就是在一次内循环中（t变换时）

每个A可以匹配三次那么一次内循环h（）的总和应该是3*n，假如说字符串中有别的字符，那么这个字符在t循环变换的过程中这个字符可以匹配的总次数等于这个字符在s串中的个数，一定小于出现最多的那个次数，比如有C那么t串变化一个循环，C最多可以匹配一次，小于3，那么总的也一定小于3*n；再加上外循环，其实无所谓每次都是3*n，所以是最大值3*n*n。

另外一个例子加入最多出现次数的字符有多个是什么情况呢

s：A A A C C C G T 其实可以这样想，因为这么一个字符串中字符最大出现次数是3

即使有多个出现3的，但是我们t字符串中的字符能和s匹配的字符的最大次数还是3次，不会因为有多个最大出现次数而改变，而会改变的就是我们有了更多的选择方案

为了让p（）函数最大，我们就会希望每个字符的匹配次数都是最大值，也就是都是都是s串中字符出现的最大次数，那么我们只要在其中最大出现次数的字符中任选一个即可

s : A A A C C C G T

t：A C C A C A C C

t : C C A C A A A A

t： _ _ _ _ _ _ _ _

……

所以这样每个字符的再一次变换循环中最大匹配次数都是s中字符出现次数的最大值

所以说到这里估计大家就可以看到大体思路了，每个位置我们又几个最大出现次数的字符就有几种选择x，一共n个位置，总共t可以满足最大值的种类就有x^n中

最终就转化成了求快速幂，哎感叹自己思维辣鸡啊

code：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int mod = 1e9+7;
const int maxn = 1e5+10;
int n;
char s[maxn];
int vis[4];
ll q_pow(ll a, ll b){
    ll ans = 1LL;
    while(b){
        if(b & 1)
            ans = ans * a % mod;
        b >>= 1;
        a = a * a % mod;
    }
    return ans;
}
int main(){
    memset(vis,0,sizeof(vis));
    scanf("%d",&n);
    scanf("%s",s);
    for(int i = 0; i < n; i++){
        if(s[i] == 'A') vis[0]++;
        if(s[i] == 'C') vis[1]++;
        if(s[i] == 'G') vis[2]++;
        if(s[i] == 'T') vis[3]++;
    }
    int maxcnt = 0;
    int cnt = 0;
    for(int i = 0; i < 4; i++){
        maxcnt = max(maxcnt,vis[i]);
    }
    for(int i = 0; i < 4; i++){
        if(vis[i] == maxcnt) cnt++;
    }
    printf("%lld\n",q_pow((ll)cnt,(ll)n));
    return 0;
}