BZOJ 1567 Blue Mary的战役地图(二维hash+二分)

题意: 求两个矩形最大公共子正方形。(n<=50)

范围这么小可以枚举子正方形的边长。那么可以对这个矩形进行二维hash，就可以在O(1)的时候求出任意子矩形的hash值。然后判断这些正方形的hash值有没有相同的

部分就行了。可以用二分来判断。

需要注意的是行和列乘的hash种子值需要不同的质数，否则可能出现冲突。

时间复杂度O(n^3logn).

# include <cstdio>
# include <cstring>
# include <cstdlib>
# include <iostream>
# include <vector>
# include <queue>
# include <stack>
# include <map>
# include <set>
# include <cmath>
# include <algorithm>
using namespace std;
# define lowbit(x) ((x)&(-x))
# define pi acos(-1.0)
# define eps 1e-
# define MOD
# define INF
# define mem(a,b) memset(a,b,sizeof(a))
# define FOR(i,a,n) for(int i=a; i<=n; ++i)
# define FO(i,a,n) for(int i=a; i<n; ++i)
# define bug puts("H");
# define lch p<<,l,mid
# define rch p<<|,mid+,r
# define mp make_pair
# define pb push_back
typedef pair<int,int> PII;
typedef vector<int> VI;
# pragma comment(linker, "/STACK:1024000000,1024000000")
typedef long long LL;
int Scan() {
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
void Out(int a) {
    if(a<) {putchar('-'); a=-a;}
    if(a>=) Out(a/);
    putchar(a%+'');
}
const int N=;
//Code begin...
 
int hash1[N][N], hash2[N][N], P[N*N], Q[N*N], n;
VI a;
VI::iterator it;
 
bool check(int ans){
    FOR(i,ans,n) FOR(j,ans,n) {
        int tmp=hash2[i][j]-hash2[i-ans][j]*Q[ans]-hash2[i][j-ans]*P[ans]+hash2[i-ans][j-ans]*P[ans]*Q[ans];
        it=lower_bound(a.begin(),a.end(),tmp);
        if (it==a.end()||*it!=tmp) continue;
        return true;
    }
    return false;
}
int main ()
{
    scanf("%d",&n);
    P[]=; FOR(i,,) P[i]=P[i-]*;
    Q[]=; FOR(i,,) Q[i]=Q[i-]*;
    FOR(i,,n) FOR(j,,n) scanf("%d",&hash1[i][j]), hash1[i][j]+=hash1[i][j-]*P[];
    FOR(i,,n) FOR(j,,n) hash1[i][j]+=hash1[i-][j]*Q[];
    FOR(i,,n) FOR(j,,n) scanf("%d",&hash2[i][j]), hash2[i][j]+=hash2[i][j-]*P[];
    FOR(i,,n) FOR(j,,n) hash2[i][j]+=hash2[i-][j]*Q[];
    int ans;
    for (ans=n; ans>=; --ans) {
        a.clear();
        FOR(i,ans,n) FOR(j,ans,n) a.pb(hash1[i][j]-hash1[i-ans][j]*Q[ans]-hash1[i][j-ans]*P[ans]+hash1[i-ans][j-ans]*P[ans]*Q[ans]);
        sort(a.begin(),a.end());
        if (check(ans)) break;
    }
    printf("%d\n",ans);
    return ;
}