[TJOI2017]DNA

题目描述

加里敦大学的生物研究所，发现了决定人喜不喜欢吃藕的基因序列S,有这个序列的碱基序列就会表现出喜欢吃藕的性状，但是研究人员发现对碱基序列S,任意修改其中不超过3个碱基，依然能够表现出吃藕的性状。现在研究人员想知道这个基因在DNA链S0上的位置。所以你需要统计在一个表现出吃藕性状的人的DNA序列S0上，有多少个连续子串可能是该基因，即有多少个S0的连续子串修改小于等于三个字母能够变成S。

输入输出格式

输入格式：

第一行有一个数T,表示有几组数据 每组数据第一行一个长度不超过10^5的碱基序列S0

每组数据第二行一个长度不超过10^5的吃藕基因序列S

输出格式：

共T行，第i行表示第i组数据中，在S0中有多少个与S等长的连续子串可能是表现吃藕性状的碱基序列

输入输出样例

输入样例#1：

1
ATCGCCCTA
CTTCA

输出样例#1：

2

说明

对于20%的数据，S0,S的长度不超过10^4

对于20%的数据，S0,S的长度不超过10^5，0<T<=10

为什么我用hash做只有40分啊， 到底哪里错了啊

----------------------------------------------------------------------------------

2018年3月2日22:32:03

更新， 终于AC了啊！！！！！检查了两天了....... 其实就是q数组初始化要到seed^100000， 我开小了才seed^10000，末尾少个零真的就是40分和100分的差距啊。。。。细心真的很重要！！！！ 

//2018年2月28日20:29:09
#include <iostream>
#include <cstdio>
#include <cstring> 
#include <algorithm>
using namespace std;

typedef unsigned long long ull;
const int N = 100005;
const int seed = 107;

int T;
char a[N], b[N];
ull Hash1[N], Hash2[N];
ull q[N];
int len1, len2;

ull get(int l, int r, ull *g){
    return g[r] - g[l-1] * q[r-l+1];
}



int getlcp(int x1, int x2){
    int l = 1, r = min(len1-x1+1, len2-x2+1), res = 0;
    while(l <= r){
        int mid = l+r >> 1;
        if(get(x1, x1+mid-1, Hash1) == get(x2, x2+mid-1, Hash2)) res=mid, l=mid+1;
        else r = mid-1;
    }
    return res;
}

bool check(int x){
    int y = 1, tmp;
    for(int i=1; i<=3; i++){
        tmp = getlcp(x, y);
        x += tmp + 1; y += tmp + 1;
        if(y > len2) return true;
    }
    
    tmp = getlcp(x , y); 
    x += tmp; y += tmp;
    if(y > len2) return true;
    else return false;

}

int main(){
    scanf("%d", &T);
    q[0] = 1;
    for(int i=1; i<=100000; i++)
        q[i] = (ull)q[i-1] * seed;
    
    while(T--){
        scanf("%s", a+1);
        scanf("%s", b+1);
        len1 = strlen(a+1); len2 = strlen(b+1);
        for(int i=1; i<=len1; i++)
            Hash1[i] = Hash1[i-1]*seed + a[i]-'A'+1;
        for(int i=1; i<=len2; i++)
            Hash2[i] = Hash2[i-1]*seed + b[i]-'A'+1;
        
        int ans = 0;
        for(int i=1; i+len2-1<=len1; i++){
            ans += check(i);
        }
        printf("%d
", ans);
    }

    return 0;
}

 

值得纪念一下