[TJOI2017] DNA 解题报告 （hash+二分）

题目链接：https://www.luogu.org/problemnew/show/P3763

题目大意：

给定原串S0，询问S0有多少个子串和给定串S相差不到3个字母

题解：

我们枚举S0的子串，问题转化为如何高效的判断两个串是否相差不到三个字母

考虑到数据范围，发现只能有log的时间余地

自然想到二分

solve每次找到第一个不同的位置并且返回，连续solve 4次，若S在这期间被处理完了，那么说明两个串相差不到3个字母

当然还有一些小细节

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
typedef unsigned long long ull;
const int N=1e5+;
int lena,lenb;
ull pin[N],haa[N],hab[N];
char a[N],b[N];
void H()
{
    haa[]=hab[]=;
    for (int i=;i<=lena;i++)
    {
        haa[i]=haa[i-]*+a[i]-'a';
    }
    for (int i=;i<=lenb;i++)
    {
        hab[i]=hab[i-]*+b[i]-'a';
    }
}
int solve(int x,int y)
{
    int l=,r=lenb;
    while (l<r)
    {
        int mid=l+r>>;
        if ((haa[x+mid-]-haa[x-]*pin[mid])==(hab[y+mid-]-hab[y-]*pin[mid]))
        {
            l=mid+;
        }
        else r=mid;
    }
    if (haa[x+l-]-haa[x-]*pin[l]!=hab[y+l-]-hab[y-]*pin[l])
    {
        return l;
    }
    else return l+;
}
int check(int x)
{
    int now=,k;
    for (int i=;i<=;i++)
    {
        k=solve(x,now);
        //printf("%d %d %d\n",x,now,k);
        x+=k;now+=k;
        if (now>lenb) return ;
    }
    k=solve(x,now);
    x+=k-;now+=k-;
    if (now>=lenb) return ;
    return ;
}
int main()
{
    pin[]=;
    for (int i=;i<=N;i++) pin[i]=pin[i-]*;
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%s",a+);
        scanf("%s",b+);
        lena=strlen(a+),lenb=strlen(b+);
        H();
        //printf("%d\n",solve(2,3));
        int re=;
        for (int i=;i<=lena-lenb+;i++)
        {
            if (check(i)) re++;
        }
        printf("%d\n",re);
    }
    return ;
}