ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

《ybtoj高效进阶》第二部分第二章例题5 子正方形

2021-03-24 20:29:19  阅读:176  来源: 互联网

标签:进阶 int ybtoj mn mid yy hb 例题 1001


题目大意

给2个正方形,求最大公共子正方形的边长。
2个正方形边长<=50

思路

显然可以通过二维hash+二分答案+暴力枚举达到目的
code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
using namespace std;
int a[1001][1001],b[1001][1001],ans;
int n,m;
unsigned long long ha[1001][1001],hb[1001][1001],p[1001],p2[1001];//数组开大又不会死
bool check(int x,int y,int l)
{
	if (x>n||y>m||x<l||y<l) return 0;
	unsigned long long aa,bb;
	aa=ha[x][y]-ha[x][y-l]*p[l]-ha[x-l][y]*p2[l]+ha[x-l][y-l]*p[l]*p2[l];
	for (int xx=l;xx<=n;xx++)
	{
		for (int yy=l;yy<=m;yy++)
		{
			if (xx>n||yy>m) continue;
			bb=hb[xx][yy]-hb[xx][yy-l]*p[l]-hb[xx-l][yy]*p2[l]+hb[xx-l][yy-l]*p[l]*p2[l];
			if (bb==aa) return 1;
		}
	}
	return 0;
}
int main()
{
	cin>>n;
	m=n;
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=m;j++)
		{
			cin>>a[i][j];
		}
	}
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=m;j++)
		{
			cin>>b[i][j];
		}
	}
	p2[0]=p[0]=1;
	for (int i=1;i<=1000;i++) p[i]=97*p[i-1],p2[i]=13331*p2[i-1];
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=m;j++)
		{
			ha[i][j]=ha[i][j-1]*97+a[i][j];
			hb[i][j]=hb[i][j-1]*97+b[i][j];
		}
	}
	for (int i=2;i<=n;i++)
	{
		for (int j=1;j<=m;j++)
		{
			ha[i][j]=ha[i-1][j]*13331+ha[i][j];
			hb[i][j]=hb[i-1][j]*13331+hb[i][j];
		}
	}
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=m;j++)
		{
			int l=1,r=n+m,mn=0;
			while (l<=r)
			{
				int mid=(l+r)>>1;
				if (check(i+mid,j+mid,mid*2+1)) l=mid+1,mn=max(mn,mid*2+1);
				else r=mid-1;
			}
			ans=max(ans,mn);
			l=1,r=n+m,mn=0;
			while (l<=r)
			{
				int mid=(l+r)>>1;
				if (check(i+mid,j+mid,mid*2)) l=mid+1,mn=max(mn,mid*2);
				else r=mid-1;
			}
			ans=max(ans,mn);
		}
	}
	cout<<ans;
    return 0;
}

标签:进阶,int,ybtoj,mn,mid,yy,hb,例题,1001
来源: https://blog.csdn.net/weixin_49843717/article/details/115186539

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有