[NOIP2017模拟]Bovine Genomics

题目描述
农夫约翰最近发现他的有些奶牛会画画，刚刚学完小学生物的他立刻开始研究起了这些奶牛的 DNA。
具体的，约翰有 N 头会画画的奶牛，N 头不会画画的奶牛，这些奶牛的 DNA 可以用 M 个字符来描述，每个字符是 A、C、G、T中的一个
比如说以下是 N=3,M=8 的一个示例：
位置: 1 2 3 4 5 6 7 8（哎呀对歪了）

会画画的牛 1: A A T C C C A T
会画画的牛 2: A C T T G C A A
会画画的牛 3: G G T C G C A A

不会画画的牛 1: A C T C C C A G
不会画画的牛 2: A C T C G C A T
不会画画的牛 3: A C T T C C A T

约翰发现，他只需要看[2,5]这个区间的 DNA就可以知道一头牛会不会画画（也就是说，没有一头会画画的牛的这个区间的 DNA是和某一头不会画画的牛的这个区间的 DNA相同的），他相信一头牛是否会画画是由这个区间的 DNA决定的，他称这样的区间为“决定区间”。
现在给你所有牛的 DNA，求最短的决定区间的长度。

输入格式
第一行两个数 N,M（N<=500；M<=500）。
接下来N 行，每行 M个字符，描述 N 头会画画的牛的DNA。
接下来N 行，每行 M个字符，描述 N 头不会画画的牛的 DNA。

输出格式
输出最短的决定区间的长度。

样例数据
输入
3 8
AATCCCAT
ACTTGCAA
GGTCGCAA
ACTCCCAG
ACTCGCAT
ACTTCCAT
输出
4

备注
【数据范围】
40% 的数据：N，M<=100。
100% 的数据：N，M<=500。

分析：数据很小，二分枚举长度，然后在每一只牛的DNA中从头开始截取这么长，把会画画的丢进哈希表，用不会画画的去匹配，如果满足没有一个相同的就返回true，继续二分找更短的；不然依次平移一个字符再次计算，如果到最后都没有满足false，就返回，继续二分找更长的。（hash差不多忘光了，回头写个hash总结吧）

代码：

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<string>#include<ctime>#include<cmath>#include<algorithm>#include<cctype>#include<iomanip>#include<queue>#include<set>using namespace std;int getint(){    int sum=0,f=1;    char ch;    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());    if(ch=='-')    {        f=-1;        ch=getchar();    }    for(;ch>='0'&&ch<='9';ch=getchar())        sum=(sum<<3)+(sum<<1)+ch-48;    return sum*f;}const int N=505;const long long p=31,H=23333;//听说进制用31进制最佳（听说而已），但是23333是质数这是真的！int n,m,ans,tot,hash[H],nxt[N];unsigned long long mi[N],num[N<<1][N],val[N*N];char s[N<<1][N];void prework(){    mi[0]=1;    for(int i=1;i<=m+1;++i)        mi[i]=mi[i-1]*p;    for(int i=1;i<=n*2;++i)        for(int j=1;j<=m;++j)            num[i][j]=num[i][j-1]*p+s[i][j]-'A'+1;}void insert(unsigned long long x){    int y=x%H;    for(int v=hash[y];v;v=nxt[v])        if(val[v]==x)            return;    nxt[++tot]=hash[y];    hash[y]=tot;    val[tot]=x;}bool find(unsigned long long x){    long long y=x%H;    for(int v=hash[y];v;v=nxt[v])        if(val[v]==x)            return true;    return false;}bool check(int len){    unsigned long long x;    for(int i=1;i<=m-len+1;++i)//在DNA中依次平移    {        memset(hash,0,sizeof(hash));        tot=0;        int bz=1;        for(int j=1;j<=n;++j)//把会画画的整个一段丢进hash表        {            x=num[j][i+len-1]-num[j][i-1]*mi[len];            insert(x);        }        for(int j=n+1;j<=2*n;++j)//用不会画画的这一段去匹配        {            x=num[j][i+len-1]-num[j][i-1]*mi[len];            if(find(x))            {                bz=0;                break;            }        }        if(bz)            return true;    }    return false;}int main(){    freopen("cownomics.in","r",stdin);    freopen("cownomics.out","w",stdout);    int len;    unsigned long long x;    n=getint();m=getint();    for(int i=1;i<=n*2;++i)        scanf("%s",s[i]+1);    prework();    int l=1,r=m;    while(l<=r)//二分查找    {        int mid=l+r>>1;        if(check(mid))            r=mid-1;        else            l=mid+1;    }    cout<<l<<endl;    return 0;}

本题结