Codeforces Round #423 (Div. 1, rated, based on VK Cup Finals) C DNA Evolution

题目大意:

给出一个字符串S(|S|<=100,000),字符集ATCG,M次操作(M<=100,000)。1操作是修改一个位置的字符,2操作是询问给定的l,r之间有多少个字符可以和给定的字符串匹配上。匹配方式如下:

先将E的开头和S中的第l个对齐,然后一个一个匹配。如果到了E的结尾,那么将E转到开头,继续匹配。如果S中的位置超过r,则退出。

一道比较简单的C题(然而我并没有打这场)

不难发现,由于字符集很小,我们可以分开考虑每个字符。而询问的位置,显然是只和l和|E|有关。那么我们可以再分开考虑|E|中的每一位。这里可以建立10*10*4个树状数组来解决,线段树会炸空间。那么我们就能在O(logn)时间解决修改,O(logn)时间解决询问。虽然要乘上常数10,但是树状数组常数很小,可以轻松跑过。

附代码:

#include<bits/stdc++.h>
#define N 100100
using namespace std;
char s[N];
int zf[128];
int bit[11][11][4][N];
int n,m;
inline int lowbit(int k){return k&-k;}
inline void add(int a,int b,int c,int d,int t){
	while(d<=n){
	//cout<<a<<" "<<b<<" "<<c<<" "<<d<<" "<<t<<"\n";
		bit[a][b][c][d]+=t;
		d+=lowbit(d);
	}
}
inline int que(int a,int b,int c,int d){
	int t=0;
	//cout<<a<<" "<<b<<" "<<c<<" "<<d<<" ";
	while(d){
		t+=bit[a][b][c][d];
		d-=lowbit(d);
	}
	//cout<<t<<"\n";
	return t;
}
int main(){
	zf['A']=0;zf['C']=1;zf['T']=2,zf['G']=3;
	scanf("%s%d",s+1,&m);
	n=strlen(s+1);
	for(int i=1;i<=n;i++)
		for(int l=1;l<=10;l++)
			add(i%l,l,zf[s[i]],i,1);
	while(m--){
		int a,l,r;char ch[12];
		scanf("%d",&a);
		if(a==1){
			scanf("%d%s",&l,ch);
			for(int i=1;i<=10;i++) add(l%i,i,zf[s[l]],l,-1),add(l%i,i,zf[ch[0]],l,1);
			s[l]=ch[0];
		}
		else{
			scanf("%d%d%s",&l,&r,ch+1);
			int ans=0,len=strlen(ch+1);
			for(int i=1;i<=len;i++){
				ans+=que((i+l-1)%len,len,zf[ch[i]],r)-que((i+l-1)%len,len,zf[ch[i]],l-1);
			}
			printf("%d\n",ans);
		}
	}
	return 0;
}


评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值