HDU 6046 hash（鸽巢原理+unordered_map）

Description
给出一个hash函数，用该hash函数可以得到一个1e6*1e6的01矩阵，现在要在这个大矩阵中找到一个1e3*1e3的小矩阵，问该小矩阵的左上角元素在大矩阵中所处位置
Input
第一行一整数T表示用例组数，每组用例输入一个1e3*1e3的01矩阵(T<=3)
Output
输出小矩阵的左上角元素在大矩阵中所处位置
Sample Input
1
0000011100
0000110011
0111111100
0011110010
0110101010
1001001001
0100111110
1111001010
0011101110
1100110100
Sample Output
Case #1 :123456 234567
Solution
从大矩阵左上角开始每隔1000-16行和1000-16列选出若干8*8小矩阵，那么由鸽巢原理至少有一个这种小矩阵会出现在所给的1000*1000的矩阵中，把这些8*8的小矩阵看作一个64位数的二进制表示，用unordered_map存下这些值与其位置之间的关系，然后枚举所给1000*1000矩阵中所有8*8矩阵，只要找到一个即固定了这个矩阵在大矩阵中的位置
Code

#include<cstdio>#include<map>#include<unordered_map>using namespace std;typedef long long ll;typedef pair<int,int>P;const int maxn=1001;unordered_map<ll,P>M;char s[maxn][maxn];inline unsigned sfr(unsigned h, unsigned x){    return h>>x;}int f(ll i, ll j){    ll w=i*1000000ll+j;    int h=0;    for(int k=0;k<5;k++)    {        h+=(int)((w>>(8*k))&255);        h+=(h<<10);        h^=sfr(h,6);    }    h+= h<<3;    h^=sfr(h,11);    h+=h<<15;    return sfr(h,27)&1;}int Case=1,T;int main(){    scanf("%d",&T);    while(T--)    {        for(int i=1;i<=1000;i++)scanf("%s",s[i]+1);        M.clear();        for(int i=1;i<=992;i++)            for(int j=1;j<=992;j++)            {                ll temp=0;                for(int k=0;k<8;k++)                    for(int l=0;l<8;l++)                    {                        temp<<=1;                        if(s[i+k][j+l]=='1')temp|=1;                    }                M[temp]=P(i,j);            }        int gg=1;        for(int i=1;i<=1e6&&gg;i+=984)            for(int j=1;j<=1e6&&gg;j+=984)            {                ll temp=0;                for(int k=0;k<8;k++)                    for(int l=0;l<8;l++)                    {                        temp<<=1;                        if(f(i+k,j+l))temp|=1;                    }                if(M.find(temp)!=M.end())                {                    P ans=M[temp];                    gg=0;                    printf("Case #%d :%d %d\n",Case++,i-ans.first+1,j-ans.second+1);                }            }    }}