题目大意:求正方形回文子矩阵数量(即左右对称、上下对称的正方形子矩阵)
正解是Manacher……但是Hash+二分是能卡过去的0.0 我太丧病了0.0
首先为了避免边长奇偶性带来的WT要把矩阵扩大二倍 然后样例就变成了这样:
00000000000
04020404040
00000000000
03010404030
00000000000
03050303030
00000000000
03010503030
00000000000
04020102040
00000000000
把这个矩阵从四个角各哈希一遍 对于每个点二分答案 验证时将四个哈希值全都取出来对比即可
然后只枚举i+j为偶数的点 得到的边长除以2就是以这个点为中心的正方形回文子矩阵数量
记得二维哈希的时候横竖的BASE值不能相同 然后这个用unsigned long long是超时的 数据没有特殊构造 所以用unsigned int就能卡过去了
回头学学Manacher……我真是太丧病了
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 2020
#define BASE1 999911657
#define BASE2 999911659
using namespace std;
typedef unsigned int ll;
ll sum[4][M][M];
int m,n;
ll ans;
ll power1[M],power2[M];
bool Judge(int x,int y,int l)
{
ll hash,temp;
hash=sum[0][x+l-1][y+l-1]
-sum[0][x-l][y+l-1]*power1[l+l-1]
-sum[0][x+l-1][y-l]*power2[l+l-1]
+sum[0][x-l][y-l]*power1[l+l-1]*power2[l+l-1];
temp=sum[1][x+l-1][y-l+1]
-sum[1][x-l][y-l+1]*power1[l+l-1]
-sum[1][x+l-1][y+l]*power2[l+l-1]
+sum[1][x-l][y+l]*power1[l+l-1]*power2[l+l-1];
if(temp!=hash)
return false;
temp=sum[2][x-l+1][y+l-1]
-sum[2][x-l+1][y-l]*power2[l+l-1]
-sum[2][x+l][y+l-1]*power1[l+l-1]
+sum[2][x+l][y-l]*power1[l+l-1]*power2[l+l-1];
if(temp!=hash)
return false;
temp=sum[3][x-l+1][y-l+1]
-sum[3][x-l+1][y+l]*power2[l+l-1]
-sum[3][x+l][y-l+1]*power1[l+l-1]
+sum[3][x+l][y+l]*power1[l+l-1]*power2[l+l-1];
if(temp!=hash)
return false;
return true;
}
ll Bisection(int x,int y)
{
int l=1,r=min(min(x,m-x+1),min(y,n-y+1));
while(l+1<r)
{
int mid=l+r>>1;
if( Judge(x,y,mid) )
l=mid;
else
r=mid;
}
if( Judge(x,y,r) )
return r;
return l;
}
int main()
{
#ifdef PoPoQQQ
freopen("1414.in","r",stdin);
freopen("1414.out","w",stdout);
#endif
int i,j,k,x;
cin>>m>>n;
m=m<<1|1;n=n<<1|1;
for(i=2;i<=m;i+=2)
for(j=2;j<=n;j+=2)
{
scanf("%d",&x);
for(k=0;k<4;k++)
sum[k][i][j]=x;
}
power1[0]=power2[0]=1;
for(i=1;i<M;i++)
power1[i]=power1[i-1]*BASE1,
power2[i]=power2[i-1]*BASE2;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
sum[0][i][j]+=sum[0][i-1][j]*BASE1;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
sum[0][i][j]+=sum[0][i][j-1]*BASE2;
for(i=1;i<=m;i++)
for(j=n;j;j--)
sum[1][i][j]+=sum[1][i-1][j]*BASE1;
for(i=1;i<=m;i++)
for(j=n;j;j--)
sum[1][i][j]+=sum[1][i][j+1]*BASE2;
for(i=m;i;i--)
for(j=1;j<=n;j++)
sum[2][i][j]+=sum[2][i+1][j]*BASE1;
for(i=m;i;i--)
for(j=1;j<=n;j++)
sum[2][i][j]+=sum[2][i][j-1]*BASE2;
for(i=m;i;i--)
for(j=n;j;j--)
sum[3][i][j]+=sum[3][i+1][j]*BASE1;
for(i=m;i;i--)
for(j=n;j;j--)
sum[3][i][j]+=sum[3][i][j+1]*BASE2;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
if( (i^j^1)&1 )
ans+=Bisection(i,j)>>1;
cout<<ans<<endl;
}