[luoguP2601] [ZJOI2009]对称的正方形（二维Hash + 二分 || Manacher）

传送门

 

很蒙蔽，不知道怎么搞。

网上看题解有说可以哈希+二分搞，也有的人说用Manacher搞，Manacher是什么鬼？以后再学。

 

对于这个题，可以从矩阵4个角hash一遍，然后枚举矩阵中的点，再二分半径。

但是得考虑边的长度为奇偶所带来的影响。

比如

1 1

1 1

这个边数为偶数的矩阵显然没法搞。

所以得在矩阵中插入0，

变成

0 0 0 0 0

0 1 0 1 0

0 0 0 0 0

0 1 0 1 0

0 0 0 0 0

具体操作就看代码好了。

然后只枚举 行 + 列 为偶数的点就行。

注意 用 unsigned long long 会超时和超空间，数据允许用 unsigned int

 

——代码

#include <cstdio>
#include <iostream>
#define UI unsigned int

const int MAXN = 2010, bs1 = 19260817, bs2 = 20011001;
int n, m, ans;
UI sum[4][MAXN][MAXN], base1[MAXN], base2[MAXN];

inline int read()
{
    int x = 0, f = 1;
    char ch = getchar();
    for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -1;
    for(; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + ch - '0';
    return x * f;
}

inline int min(int x, int y)
{
    return x < y ? x : y;
}

inline bool pd(int x, int y, int l)
{
    UI t, h;
    h = sum[0][x + l - 1][y + l - 1]
      -    sum[0][x - l][y + l - 1] * base1[l + l - 1]
      - sum[0][x + l - 1][y - l] * base2[l + l - 1]
      + sum[0][x - l][y - l] * base1[l + l - 1] * base2[l + l - 1];
    t = sum[1][x + l - 1][y - l + 1]
      - sum[1][x - l][y - l + 1] * base1[l + l - 1]
      - sum[1][x + l - 1][y + l] * base2[l + l - 1]
      + sum[1][x - l][y + l] * base1[l + l - 1] * base2[l + l - 1];
    if(h ^ t) return 0;
    t = sum[2][x - l + 1][y + l - 1]
      - sum[2][x + l][y + l - 1] * base1[l + l - 1]
      - sum[2][x - l + 1][y - l] * base2[l + l - 1]
      + sum[2][x + l][y - l] * base1[l + l - 1] * base2[l + l - 1];
    if(h ^ t) return 0;
    t = sum[3][x - l + 1][y - l + 1]
      - sum[3][x + l][y - l + 1] * base1[l + l - 1]
      - sum[3][x - l + 1][y + l] * base2[l + l - 1]
      + sum[3][x + l][y + l] * base1[l + l - 1] * base2[l + l - 1];
    if(h ^ t) return 0;
    return 1;
}

inline int work(int i, int j)
{
    int mid, s = 0, x = 1, y = min(min(i, n - i + 1), min(j, m - j + 1));//二分半径 
    while(x <= y)
    {
        mid = (x + y) >> 1;
        if(pd(i, j, mid)) s = mid, x = mid + 1;
        else y = mid - 1;
    }
    return s;
}

int main()
{
    int i, j, k, x;
    n = read();
    m = read();
    n = n << 1 | 1;
    m = m << 1 | 1;
    for(i = 2; i <= n; i += 2)
        for(j = 2; j <= m; j += 2)
        {
            x = read();
            for(k = 0; k < 4; k++) sum[k][i][j] = x;
        }
    base1[0] = base2[0] = 1;
    for(i = 1; i <= n; i++) base1[i] = base1[i - 1] * bs1;
    for(i = 1; i <= m; i++) base2[i] = base2[i - 1] * bs2;
    for(i = 1; i <= n; i++)
        for(j = 1; j <= m; j++)
            sum[0][i][j] += sum[0][i - 1][j] * bs1;
    for(i = 1; i <= n; i++)
        for(j = 1; j <= m; j++)
            sum[0][i][j] += sum[0][i][j - 1] * bs2;
    for(i = 1; i <= n; i++)
        for(j = m; j; j--)
            sum[1][i][j] += sum[1][i - 1][j] * bs1;
    for(i = 1; i <= n; i++)
        for(j = m; j; j--)
            sum[1][i][j] += sum[1][i][j + 1] * bs2;
    for(i = n; i; i--)
        for(j = 1; j <= m; j++)
            sum[2][i][j] += sum[2][i + 1][j] * bs1;
    for(i = n; i; i--)
        for(j = 1; j <= m; j++)
            sum[2][i][j] += sum[2][i][j - 1] * bs2;
    for(i = n; i; i--)
        for(j = m; j; j--)
            sum[3][i][j] += sum[3][i + 1][j] * bs1;
    for(i = n; i; i--)
        for(j = m; j; j--)
            sum[3][i][j] += sum[3][i][j + 1] * bs2;
    for(i = 1; i <= n; i++)
        for(j = 1; j <= m; j++)
            if((i ^ j ^ 1) & 1)
                ans += work(i, j) >> 1;
    printf("%d\n", ans);
    return 0;
}

 

Manacher的话，学完再搞吧。