bzoj1567[JSOI2008]Blue Mary的战役地图二分+矩阵hash

二分查找与矩阵Hash解决Blue Mary的战役地图问题

最新推荐文章于 2019-02-06 10:33:32 发布

EMber _

最新推荐文章于 2019-02-06 10:33:32 发布

阅读量264

点赞数

CC 4.0 BY-SA版权

分类专栏： bzoj hash 二分

本文链接：https://blog.youkuaiyun.com/qq_35866453/article/details/78439474

bzoj 同时被 3 个专栏收录

452 篇文章

订阅专栏

二分

28 篇文章

订阅专栏

hash

10 篇文章

订阅专栏

该博客介绍了如何利用二分查找和矩阵Hash技术解决JSOI2008竞赛中名为'Blue Mary的战役地图'的问题。主要内容包括计算两个矩阵的Hash值，并通过二分搜索找到矩阵间最大重合的正方形边长。文章强调了选择合适模数的挑战，建议使用自然溢出处理。

题意：给出两个矩阵，求最大重合。
做出两个矩阵的hash以后二分最大正方形的边长，基础hash。
模数不好设，用自然溢出就好。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int mo=2333333;
int n;
ll hash1[55][55],hash2[55][55];
ll bin[55],tmp[100005];
int mp1[55][55],mp2[55][55];
inline int calc1(int x,int y,int l)
{
    ll v=1;
    fo(i,x,x+l-1)v=(v*1789+hash1[i][y+l-1]-hash1[i][y-1]*bin[l]);
    return v;
}
inline int calc2(int x,int y,int l)
{
    ll v=1;
    fo(i,x,x+l-1)v=(v*1789+hash2[i][y+l-1]-hash2[i][y-1]*bin[l]);
    return v;
}
bool pd(int x)
{
    int tot=0;
    fo(i,1,n-x+1)
    {
        fo(j,1,n-x+1)
        tmp[++tot]=calc1(i,j,x);
    }
    sort(tmp+1,tmp+tot+1);
    fo(i,1,n-x+1)
    {
        fo(j,1,n-x+1)
        {
            ll v=calc2(i,j,x);
            ll pos=*lower_bound(tmp+1,tmp+tot+1,v);
            if (pos==v)return 1;
        }
    }
    return 0;
} 
int main()
{
    bin[0]=1;
    fo(i,1,50)bin[i]=bin[i-1]*1321;
    scanf("%d",&n);
    fo(i,1,n)
    {
        hash1[i][0]=1;
        fo(j,1,n)
        {
            scanf("%d",&mp1[i][j]);
            hash1[i][j]=(1ll*hash1[i][j-1]*1321+mp1[i][j]);
        }
    }
    fo(i,1,n)
    {
        hash2[i][0]=1;
        fo(j,1,n)
        {
            scanf("%d",&mp2[i][j]);
            hash2[i][j]=(1ll*hash2[i][j-1]*1321+mp2[i][j]);
        }
    }
    int l=1,r=n,ans=0;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (pd(mid))ans=mid,l=mid+1;
        else r=mid-1;
    }
    printf("%d\n",ans);
}