bzoj 1414 && bzoj 3705: [ZJOI2009]对称的正方形(二维Hash)
来源:互联网 发布:python解析网页源代码 编辑:程序博客网 时间:2024/06/11 05:50
1414: [ZJOI2009]对称的正方形
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 727 Solved: 345
[Submit][Status][Discuss]
Description
Orez很喜欢搜集一些神秘的数据,并经常把它们排成一个矩阵进行研究。最近,Orez又得到了一些数据,并已经把它们排成了一个n行m列的矩阵。通过观察,Orez发现这些数据蕴涵了一个奇特的数,就是矩阵中上下对称且左右对称的正方形子矩阵的个数。 Orez自然很想知道这个数是多少,可是矩阵太大,无法去数。只能请你编个程序来计算出这个数。
Input
文件的第一行为两个整数n和m。接下来n行每行包含m个正整数,表示Orez得到的矩阵。
Output
文件中仅包含一个整数answer,表示矩阵中有answer个上下左右对称的正方形子矩阵。
Sample Input
5 5
4 2 4 4 4
3 1 4 4 3
3 5 3 3 3
3 1 5 3 3
4 2 1 2 4
Sample Output
27
暴力枚举每个正方形的中点,然后二分+Hash判断当前最大的对称正方形
不过:
①正方形的边长可能是奇数或者偶数,所以要先像manacher那样处理再hash,最后答案/2
②维护四个方向的Hash,当且仅当四个方向Hash值相同,正方形对称
③复杂度O(n²logn),常数特别特别大,所以要优化下
#include<stdio.h>#include<algorithm>using namespace std;#define uint unsigned int#define x 731729523#define y 518758277uint a[4][2005][2005], px[2005] = {1}, py[2005] = {1};int Jud(int x1, int x2, int y1, int y2){uint X, Y, c1, c2, c3, c4;X = py[x2-x1+1], Y = px[y2-y1+1];c1 = a[0][x2][y2]-a[0][x2][y1-1]*Y-a[0][x1-1][y2]*X+a[0][x1-1][y1-1]*X*Y;c2 = a[1][x2][y1]-a[1][x2][y2+1]*Y-a[1][x1-1][y1]*X+a[1][x1-1][y2+1]*X*Y;if(c2!=c1)return 0;c3 = a[2][x1][y2]-a[2][x1][y1-1]*Y-a[2][x2+1][y2]*X+a[2][x2+1][y1-1]*X*Y;if(c2!=c3)return 0;c4 = a[3][x1][y1]-a[3][x1][y2+1]*Y-a[3][x2+1][y1]*X+a[3][x2+1][y2+1]*X*Y;if(c3!=c4)return 0;return 1;}int main(void){int n, m, i, j, ans, l, r, mid;scanf("%d%d", &n, &m);for(i=1;i<=n;i++){for(j=1;j<=m;j++){scanf("%u", &a[0][i*2][j*2]);a[1][i*2][j*2] = a[2][i*2][j*2] = a[3][i*2][j*2] = a[0][i*2][j*2];}}n = n*2+1;m = m*2+1;for(i=1;i<=2002;i++)px[i] = px[i-1]*x, py[i] = py[i-1]*y;for(i=1;i<=n;i++)for(j=1;j<=m;j++)a[0][i][j] += a[0][i][j-1]*x;for(i=1;i<=n;i++)for(j=1;j<=m;j++)a[0][i][j] += a[0][i-1][j]*y;for(i=1;i<=n;i++)for(j=m;j>=1;j--)a[1][i][j] += a[1][i][j+1]*x;for(i=1;i<=n;i++)for(j=m;j>=1;j--)a[1][i][j] += a[1][i-1][j]*y;for(i=n;i>=1;i--)for(j=1;j<=m;j++)a[2][i][j] += a[2][i][j-1]*x;for(i=n;i>=1;i--)for(j=1;j<=m;j++)a[2][i][j] += a[2][i+1][j]*y;for(i=n;i>=1;i--)for(j=m;j>=1;j--)a[3][i][j] += a[3][i][j+1]*x;for(i=n;i>=1;i--)for(j=m;j>=1;j--)a[3][i][j] += a[3][i+1][j]*y;ans = 0;for(i=1;i<=n;i++){for(j=1;j<=m;j++){if((i+j)%2==0){l = 1, r = min(min(i, n-i+1), min(j, m-j+1));while(l<r){mid = (l+r+1)/2;if(Jud(i-mid+1, i+mid-1, j-mid+1, j+mid-1))l = mid;elser = mid-1;}ans += l/2;}}}printf("%d\n", ans);return 0;}
阅读全文
1 0
- bzoj 1414 && bzoj 3705: [ZJOI2009]对称的正方形(二维Hash)
- BZOJ 1414 ZJOI2009 对称的正方形 Hash+二分
- [BZOJ]1414: [ZJOI2009]对称的正方形 二分+hash
- bzoj 1414: [ZJOI2009]对称的正方形
- bzoj 1414: [ZJOI2009]对称的正方形
- 1414: [ZJOI2009]对称的正方形 Hash+二分
- bzoj 1414 对称的正方形
- bzoj1414: [ZJOI2009]对称的正方形
- [BZOJ1414][ZJOI2009]对称的正方形(manacher+单调栈+二分)
- bzoj1414 [ZJOI2009]对称的正方形 && bzoj3705 对称的正方形
- [ZJOI2009对称的正方形]ST+Manacher
- bzoj 2351: [BeiJing2011]Matrix(二维Hash)
- BZOJ 1047 [HAOI2007]理想的正方形 二维ST表+压维
- bzoj 1433: [ZJOI2009]假期的宿舍(最大流)
- BZOJ [ZJOI2009]狼和羊的故事
- bzoj 1433: [ZJOI2009]假期的宿舍
- 【BZOJ 1433】 [ZJOI2009]假期的宿舍
- 【BZOJ 1433】 [ZJOI2009]假期的宿舍
- TEXMACS在ubuntu下的使用
- Blade:一款简洁优雅、微内核设计的Java Web框架
- 打造自己的Ubuntu渗透利器--msf、nessus安装
- so big
- 接口与类型信息
- bzoj 1414 && bzoj 3705: [ZJOI2009]对称的正方形(二维Hash)
- Android--hw_get_module解析
- 数字图像处理(1)
- 在linux中编译ACE静态库
- BZOJ-3209 (数论)
- 更换Eclipse黑色皮肤
- 数组计数某一个数的数量
- JS 对象与数组的变化监听
- 记一次混淆后引起的异常,被覆盖的方法未抛出Exception