您的位置:首页 > 其它

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

2017-07-12 10:50 344 查看
题目大意:

给出一个字符串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
;
int zf[128];
int bit[11][11][4]
;
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;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: 
相关文章推荐