【YBTOJ】子正方形

26 篇文章 0 订阅
15 篇文章 0 订阅
这篇博客探讨了如何通过四重循环枚举右下角坐标,再结合二分查找算法来判断是否存在特定大小的子正方形。内容主要集中在算法实现上。
摘要由CSDN通过智能技术生成

在这里插入图片描述

思路:

直接四重循环枚举右下角,然后二分长度hash判断

c o d e code code

#include<iostream>
#include<cstdio>

using namespace std;

unsigned long long mod1=1000000007ull, mod2=1000000009ull;
int n;
unsigned long long hasha[100][100], hashb[100][100];
unsigned long long p1[100], p2[100];

bool check(int mid, int x, int y, int xx, int yy)
{
	if(x-mid<0||y-mid<0||xx-mid<0||yy-mid<0)
		return 0;
	unsigned long long sa=hasha[x][y]-hasha[x-mid][y]*p2[mid]-hasha[x][y-mid]*p1[mid]+hasha[x-mid][y-mid]*p1[mid]*p2[mid];
	unsigned long long sb=hashb[xx][yy]-hashb[xx-mid][yy]*p2[mid]-hashb[xx][yy-mid]*p1[mid]+hashb[xx-mid][yy-mid]*p1[mid]*p2[mid];
	return sa==sb;
}
int main()
{
	scanf("%d", &n);
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			scanf("%d", &hasha[i][j]);
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			scanf("%d", &hashb[i][j]);
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			hasha[i][j]+=hasha[i][j-1]*mod1;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			hasha[i][j]+=hasha[i-1][j]*mod2;
	p1[0]=p2[0]=1ull;
	for(int i=1; i<=n; i++)
		p1[i]=p1[i-1]*mod1, p2[i]=p2[i-1]*mod2;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			hashb[i][j]+=hashb[i][j-1]*mod1;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			hashb[i][j]+=hashb[i-1][j]*mod2;
	int ans=0;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			for(int i1=1; i1<=n; i1++)
				for(int j1=1; j1<=n; j1++)
				{
					int l=0, r=n;
					while(l<=r)
					{
						int mid=l+r>>1;
						if(check(mid, i, j, i1, j1))
						{
							ans=max(ans, mid);
							l=mid+1;
						}
						else
							r=mid-1;
					}
				}
	printf("%d", ans);	
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值