bzoj1141（矩阵hash）

一：编译超时，原因是数组开得过大！！

二：如果更改数组大小，在往小改的时候，注意有没有预处理的情况，防止越界

关于hash：base1，base2这样的取值，差不多取就行，根本不必纠结于此

关于题目：求有多少个向下左右对称的正方形子矩阵，

显然，这样的矩阵是  靶形 样式的圈绕在那里

其实就是，枚举位置，判断以这个位置为中心有多少层圈，其实就是二分多少层，矩阵hash判断左上右上左下右下是否相同？

这需要我们预处理从各个角开始的前缀矩阵hash，四个方向来作hash

因为有可能

2 2 

2 2  这样的矩阵，就是说要分单双两种情况，这样就是向manacher一样，相邻的之间加上0（不与矩阵中number重合的）的数

0 0 0 0 0

0 2 0 2 0

0 0 0 0 0

0 2 0 2 0

0 0 0 0 0

算是单双两种情况问题，通过灵活处理转化为一种情况

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
typedef unsigned long long ull;
const int base1=113;
const int base2=999983;
ull ul[2005][2005],ur[2005][2005],
	dl[2005][2005],dr[2005][2005];
ull b1[2020],b2[2020];
int n,m,mp[1005][1005],a[2020][2020];


inline int read()
{
	int ans,f=1;char ch;
	while ((ch=getchar())<'0'||ch>'9') if (ch=='-') f=-1;
	ans=ch-'0';
	while ((ch=getchar())>='0'&&ch<='9') ans=ans*10+ch-'0';
	return ans*f;
}


void init()
{
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) ul[i][j]=ul[i][j-1]*base1+a[i][j];
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) ul[i][j]=ul[i-1][j]*base2+ul[i][j];
	
	for (int i=n;i>=1;i--)
		for (int j=1;j<=m;j++) dl[i][j]=dl[i][j-1]*base1+a[i][j];
	for (int i=n;i>=1;i--)
		for (int j=1;j<=m;j++) dl[i][j]=dl[i+1][j]*base2+dl[i][j];
	
	for (int i=1;i<=n;i++)
		for (int j=m;j>=1;j--) ur[i][j]=ur[i][j+1]*base1+a[i][j];
	for (int i=1;i<=n;i++)
		for (int j=m;j>=1;j--) ur[i][j]=ur[i-1][j]*base2+ur[i][j];
		
	for (int i=n;i>=1;i--)
		for (int j=m;j>=1;j--) dr[i][j]=dr[i][j+1]*base1+a[i][j];
	for (int i=n;i>=1;i--)
		for (int j=m;j>=1;j--) dr[i][j]=dr[i+1][j]*base2+dr[i][j];	 			
}

int main()
{
	b1[0]=b2[0]=1;
	for (int i=1;i<=2005;i++) b1[i]=b1[i-1]*base1,b2[i]=b2[i-1]*base2;
	n=read();m=read();
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) mp[i][j]=read();
	int nn=0,mm=0;
	for (int i=1;i<=n*2+1;i++)
		for (int j=1;j<=m*2+1;j++) a[i][j]=0;
	for (int i=1;i<=n;i++)
	{
		nn+=2,mm=0;
		for (int j=1;j<=m;j++)
		mm+=2,
		a[nn][mm]=mp[i][j];
	}
	n=n*2+1,m=m*2+1;
	init();
	int ans=0;
	for (int i=1;i<=n;i++)


		for (int j=1;j<=m;j++) if (!((i+j)&1))
		{
			int l=1,r=min(min(i,n-i+1),min(j,m-j+1)),mid,tt=1,
				hul,hur,hdl,hdr;
			while (l<=r)
			{
				mid=(l+r)>>1;
				
				hul=
				ul[i][j]-
				ul[i][j-mid]*b1[mid]-
				ul[i-mid][j]*b2[mid]+
				ul[i-mid][j-mid]*b1[mid]*b2[mid];
				
				hur=
				ur[i][j]-
				ur[i][j+mid]*b1[mid]-
				ur[i-mid][j]*b2[mid]+
				ur[i-mid][j+mid]*b1[mid]*b2[mid];
				
				hdl=
				dl[i][j]-
				dl[i][j-mid]*b1[mid]-
				dl[i+mid][j]*b2[mid]+
				dl[i+mid][j-mid]*b1[mid]*b2[mid];
				
				hdr=
				dr[i][j]-
				dr[i][j+mid]*b1[mid]-
				dr[i+mid][j]*b2[mid]+
				dr[i+mid][j+mid]*b1[mid]*b2[mid];
				
				if ((hul==hur)&&(hur==hdr)&&(hdr==hdl)) 
				tt=mid,l=mid+1;
				else r=mid-1;
			}
			ans+=tt/2;
		}

		printf("%d",ans);
	return 0;
}