BZOJ 1170 [Balkan2007]Cipher Hash
题目大意：给定一个二维矩阵，求出现次数最多的a*b的子矩阵
二维Hash，只要记住横纵的BASE不能相同就可以，爱怎么搞怎么搞
一开始写的自然溢出 结果OLE 以为是自然溢出被卡掉了于是写了双取模…… 结果还是OLE
最后发现尼玛这题读入坑爹……字符串里有空格不说，满满的不可见字符是咋回事……
记住不要用scanf读入……可以用gets，或者fread，注意要把一开始的回车过滤掉
getchar读进来全是错的 不知道怎么回事……
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 1010
#define BASE1 2591
#define BASE2 2593
#define MOD1 1000000007
#define MOD2 1000000009
using namespace std;
typedef long long ll;
int m,n,a,b,ans_time;
char s[M][M];
ll sum1[M][M],sum2[M][M],ans1,ans2,BASE1_a[3]={0,1,1},BASE2_b[3]={0,1,1};
namespace Hash_Table{
struct list{
ll hash_value1,hash_value2;
int times;
list *next;
list(ll _,ll __,list *___):
hash_value1(_),hash_value2(__),times(0),next(___) {}
}*head[1001001];
inline int& Hash(ll hash1,ll hash2)
{
list *temp;
int x=hash1*hash2%1001001;
for(temp=head[x];temp;temp=temp->next)
if(temp->hash_value1==hash1&&temp->hash_value2==hash2)
return temp->times;
return head[x]=new list(hash1,hash2,head[x]),head[x]->times;
}
}
void Output()
{
int i,j,k,l;
for(i=a;i<=m;i++)
for(j=b;j<=n;j++)
{
ll hash1=((sum1[i][j]
-sum1[i-a][j]*BASE1_a[1]
-sum1[i][j-b]*BASE2_b[1]
+sum1[i-a][j-b]*BASE1_a[1]%MOD1*BASE2_b[1]
)%MOD1+MOD1)%MOD1;
ll hash2=((sum2[i][j]
-sum2[i-a][j]*BASE1_a[2]
-sum2[i][j-b]*BASE2_b[2]
+sum2[i-a][j-b]*BASE1_a[2]%MOD2*BASE2_b[2]
)%MOD2+MOD2)%MOD2;
if(hash1==ans1&&hash2==ans2)
{
for(k=i-a+1;k<=i;k++,puts(""))
for(l=j-b+1;l<=j;l++)
putchar(s[k][l]);
return ;
}
}
}
void Find_Position()
{
int i,j;
for(i=a;i<=m;i++)
for(j=b;j<=n;j++)
{
ll hash1=((sum1[i][j]
-sum1[i-a][j]*BASE1_a[1]
-sum1[i][j-b]*BASE2_b[1]
+sum1[i-a][j-b]*BASE1_a[1]%MOD1*BASE2_b[1]
)%MOD1+MOD1)%MOD1;
ll hash2=((sum2[i][j]
-sum2[i-a][j]*BASE1_a[2]
-sum2[i][j-b]*BASE2_b[2]
+sum2[i-a][j-b]*BASE1_a[2]%MOD2*BASE2_b[2]
)%MOD2+MOD2)%MOD2;
if(hash1==ans1&&hash2==ans2)
printf("%d %d\n",i-a+1,j-b+1);
}
}
int main()
{
int i,j;
scanf("%d %d\n",&m,&n);
for(i=1;i<=m;i++)
{
fread(s[i]+1,n+1,1,stdin);
for(j=1;j<=n;j++)
sum1[i][j]=sum2[i][j]=s[i][j];
}
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
sum1[i][j]=(sum1[i][j]+sum1[i-1][j]*BASE1)%MOD1,
sum2[i][j]=(sum2[i][j]+sum2[i-1][j]*BASE1)%MOD2;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
sum1[i][j]=(sum1[i][j]+sum1[i][j-1]*BASE2)%MOD1,
sum2[i][j]=(sum2[i][j]+sum2[i][j-1]*BASE2)%MOD2;
cin>>a>>b;
cout<<a<<' '<<b<<endl;
for(i=1;i<=a;i++)
{
BASE1_a[1]*=BASE1;BASE1_a[1]%=MOD1;
BASE1_a[2]*=BASE1;BASE1_a[2]%=MOD2;
}
for(j=1;j<=b;j++)
{
BASE2_b[1]*=BASE2;BASE2_b[1]%=MOD1;
BASE2_b[2]*=BASE2;BASE2_b[2]%=MOD2;
}
for(i=a;i<=m;i++)
for(j=b;j<=n;j++)
{
ll hash1=((sum1[i][j]
-sum1[i-a][j]*BASE1_a[1]
-sum1[i][j-b]*BASE2_b[1]
+sum1[i-a][j-b]*BASE1_a[1]%MOD1*BASE2_b[1]
)%MOD1+MOD1)%MOD1;
ll hash2=((sum2[i][j]
-sum2[i-a][j]*BASE1_a[2]
-sum2[i][j-b]*BASE2_b[2]
+sum2[i-a][j-b]*BASE1_a[2]%MOD2*BASE2_b[2]
)%MOD2+MOD2)%MOD2;
int& val=Hash_Table::Hash(hash1,hash2);
if(++val>ans_time)
ans_time=val,ans1=hash1,ans2=hash2;
}
Output();
cout<<ans_time<<endl;
Find_Position();
return 0;
}