6817. 【2020.10.07提高组模拟】DNA 序列

Description

Input

Output

Sample Input

Sample Input1
AAAAA
1


Sample Input2
ACTCACTC
4
 

Sample Output

Sample Output1
5


【样例 1 解释】
对于这段 DNA 序列，连续的 1 个碱基组成的碱基序列只有 A，共出现 5 次，所以答案为 5。 


Sample Output2
2

【样例 2 解释】
对于这段 DNA 序列，连续的 4 个碱基组成的碱基序列为：ACTC, CTCA, TCAC 与CACT。其中 ACTC 出现 2 次，其余均出现 1 次，所以出现最多的次数为 2，即为答案。 
 

Data Constraint

Solution

转成4进制标记数组即可。

Code 

#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define I int
#define ll long long
#define F(i,a,b) for(I i=a;i<=b;i++)
#define Fd(i,a,b) for(I i=a;i>=b;i--)
#define N 5000004
using namespace std;
I k,x,a[N],len,f[13],ans;
I bz[N<<1];
char c;
I main(){
	freopen("dna.in","r",stdin);
	freopen("dna.out","w",stdout);
	c=getchar();
	while(c<'A'||c>'T') c=getchar();
	while(c>='A'&&c<='T'){
		len++;
		a[len]=c=='A'?0:(c=='C'?1:(c=='G'?2:3));
		c=getchar();
	}
	f[0]=f[1]=1;
	F(i,2,10) f[i]=f[i-1]*4;
	k=x=0;
	while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9'){k=k*10+c-'0';c=getchar();}
	F(i,1,k-1) x=x*4+a[i];
	F(i,k,len){
		x=(x-f[k]*(x/f[k]))*4+a[i];
		ans=max(ans,++bz[x]);
	}
	printf("%d\n",ans);
	return 0;
}