ICode9

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

ybtoj·子正方形【Hash】

2021-05-21 21:01:18  阅读:190  来源: 互联网

标签:Hash 正方形 int ybtoj ll h1 h2 55 xr


子正方形

Description–

高效进阶「字符串算法」第2章 Hash 和 Hash 表课堂过关 例题5


代码–

#include <iostream>
#include <cstdio>
#define ull long long

using namespace std;

int n, ans, a[55][55], b[55][55];
ull t1[55], t2[55], h1[55][55], h2[55][55];

bool xyn(int x, int y, int xr, int yr, int ll) //求这两个矩阵的哈希值
{
	ull aa, bb;
	aa = h1[x][y] - h1[x - ll][y] * t2[ll] - h1[x][y - ll] * t1[ll] + h1[x - ll][y - ll] * t1[ll] * t2[ll];
	bb = h2[xr][yr] - h2[xr - ll][yr] * t2[ll] - h2[xr][yr - ll] * t1[ll] + h2[xr - ll][yr - ll] * t1[ll] * t2[ll];
	return aa == bb;
}

void hll(int x, int y, int xr, int yr) //二分
{
	int l = 0, r = min(min(x, y), min(xr, yr));
	while (l <= r)
	{
		int mid = (l + r + 1) >> 1;
		if (xyn(x, y, xr, yr, mid))
		  l = mid + 1, ans = max(ans, mid);
		else r = mid - 1;
	}
}

int main()
{
	scanf("%d", &n);
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= n; ++j)
	    scanf("%d", &a[i][j]);
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= n; ++j)
	    scanf("%d", &b[i][j]);
	t1[0] = t2[0] = 1;
	for (int i = 1; i <= n; ++i)
	  t1[i] = t1[i - 1] * 131, t2[i] = t2[i - 1] * 61;
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= n; ++j)
	  {
	  	  h1[i][j] = h1[i][j - 1] * 131 + a[i][j];
	  	  h2[i][j] = h2[i][j - 1] * 131 + b[i][j];
	  }
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= n; ++j)
	  {
	  	  h1[i][j] += h1[i - 1][j] * 61;
	  	  h2[i][j] += h2[i - 1][j] * 61;
	  }
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= n; ++j)
	    for (int k = 1; k <= n; ++k)
	      for (int l = 1; l <= n; ++l)
	         if (a[i][j] == b[k][l])
			   hll(i, j, k, l);
	printf("%d", ans);
	
	return 0;
}

标签:Hash,正方形,int,ybtoj,ll,h1,h2,55,xr
来源: https://blog.csdn.net/qq_43654542/article/details/116533956

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

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

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

ICode9版权所有