codeforces-827C-树状数组,同余定理

题目大意:给定一段DNA序列,和q个query,每次可以改变一个碱基,或者查询一个区间里面有多少个碱基与给定碱基序列(可以无限重复)相同;

题目解析:首先因为匹配序列可以无限循环,所以可以根据位置%长度来确定一般性,因为如果一个碱基的位置mod (len)与匹配序列一个碱基的位置mod(len)一样,那么这个无限延长的碱基序列一定能匹配到第一个碱基,所以我们可以BIT来维护前缀和,根据碱基种类,匹配序列长度,pos%匹配序列长度来分类;

AC代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
using namespace std;
const int maxn=100010;
int dp[4][11][11][maxn],has[200];
void init()
{
	memset(dp,0,sizeof(dp));
	has['A']=0;
	has['T']=1;
	has['C']=2;
	has['G']=3; 
}
int lowbit(int k)
{
	return k&(-k);
}
void add(int pos,int type,int val)
{
	for(int i=pos;i<maxn;i+=lowbit(i))
	{
		for(int j=1;j<=10;j++)
		{
			dp[type][j][pos%j][i]+=val;
		}
	}
}
int query(int type,int len,int mod,int x)
{
	int ans=0;
	for(int i=x;i>0;i-=lowbit(i))
	{
		ans+=dp[type][len][mod][i];
	}
	return ans;
}
int sum(int l,int r,int mod,int type,int len)
{
	return query(type,len,mod,r)-query(type,len,mod,l-1);
}
char s[maxn];
int main()
{
	init();
	scanf("%s",&s[1]);	
	int l=strlen(s+1);
    for (int i = 1; i <=l; ++i) {  
        add(i, has[s[i]], 1);  
    }  
    char ch[111];  
    int q,a,b,c; scanf("%d", &q);  
    while (q--) {  
        scanf("%d", &a);  
        if (a == 1) {  
            scanf("%d %s", &b, ch);  
            add(b, has[s[b]], -1);  
            add(b, has[ch[0]], 1);  
            s[b] = ch[0];  
        }  
        else  
        {  
            scanf("%d %d %s", &b, &c, ch);  
            int len = strlen(ch);  
            int ans = 0;  
            for (int i = 0; i < len; ++i) {  
                ans += sum(b, c, (i + b) % len, has[ch[i]], len);

            }  
            printf("%d\n", ans);  
        }  
    }  
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值