7.12T1序列

1．序列

【问题描述】 Hzy 得到了一个字符串，这个字符串只有’A’,’G’,’C’,’T’这四种字符，她发现这个 序列中连续 k 个字符可以形成一种新的字符序列，她称这种序列为 Hzy 序列，她现在想知 道在所有的 Hzy 序列中，出现最多的一种的出现次数。

【输入格式】 输入文件名为 uaena.in。 输入文件的第一行为一个字符串，保证合法。 输入文件的第二行为一个正整数 k。

【输出格式】 输出文件名为 uaena.out。 输出文件只有一个数，即所求答案。

【输入输出样例 1】

input

AAAAA

1

output

5

【样例解释 1】 对于这段字符串中，连续的 k 个字符组成的 Hzy 序列只有 A，共出现 5 次，所以答案 为 5。

【输入输出样例 2】

ueana.in

ACTCACTC

4

ueana.out

2

【样例解释 1】 对于这段字符串中，连续的 k 个字符组成的 Hzy 序列有 ACTC，CTCA，TCAC，CACT 其中 ACTC 共出现 2 次，其余只出现了 1 次， 所以答案为 2。

【数据规模与规定】

n<=5 × 10^6  k<=10

sol: 哈希一下，O(n）水过

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
inline ll read()
{
    ll s=0;
    bool f=0;
    char ch=' ';
    while(!isdigit(ch))
    {
        f|=(ch=='-'); ch=getchar();
    }
    while(isdigit(ch))
    {
        s=(s<<3)+(s<<1)+(ch^48); ch=getchar();
    }
    return (f)?(-s):(s);
}
#define R(x) x=read()
inline void write(ll x)
{
    if(x<0)
    {
        putchar('-'); x=-x;
    }
    if(x<10)
    {
        putchar(x+'0'); return;
    }
    write(x/10);
    putchar((x%10)+'0');
    return;
}
#define W(x) write(x),putchar(' ')
#define Wl(x) write(x),putchar('
')
const int N=5000005;
const ll Base=4,Mod=19260817;
int n,m,a[N];
char S[N];
inline int Check(char ch)
{
    if(ch=='A') return 0;
    if(ch=='T') return 1;
    if(ch=='C') return 2;
    if(ch=='G') return 3;
}
ll Seed[N];
inline void Pre(int n)
{
    int i;
    Seed[0]=1ll; for(i=1;i<=n;i++) Seed[i]=(ll)(Seed[i-1]*4%Mod);
}
struct Hash
{
    ll Hash[N];
    inline void Make(int n,int *a)
    {
        int i; Hash[0]=0; for(i=1;i<=n;i++) Hash[i]=(Hash[i-1]*4%Mod+a[i])%Mod;
    }
    inline ll Ask(int l,int r)
    {
        return (ll)(Hash[r]+Mod-Hash[l-1]*Seed[r-l+1]%Mod)%Mod;
    }
    inline ll Get(int l,int r)
    {
        return Ask(l,r)*Seed[n-r]%Mod;
    }
}Ha;
int Sum[3000005];
int main()
{
    freopen("ueana.in","r",stdin);
    freopen("ueana.out","w",stdout);
    int i,ans=0;
    scanf("%s",S+1); n=strlen(S+1); R(m);
    for(i=1;i<=n;i++) a[i]=Check(S[i]);
    Pre(n); Ha.Make(n,a);
    for(i=m;i<=n;i++)
    {
//        cout<<i-m+1<<' '<<i<<' '<<Ha.Ask(i-m+1,i)<<' '<<Ha.Get(i-m+1,i)<<endl;
        Sum[Ha.Ask(i-m+1,i)]++;
    }
    for(i=0;i<=Seed[m];i++) ans=max(ans,Sum[i]);
    Wl(ans);
    return 0;
}
/*
input
AAAAA
1
output
5

input
ACTCACTC
4
output
2
*/