日常训练 20170619 被删除量子化粉碎未出生的回文矩阵

题意：
  一个 $X \times Y$ 的矩阵当且仅当它关于水平轴对称时才被称为回文矩阵。例如：
  a b a
  b a b
  现在给出一个 $N \times N\ (N \le 3000)$ 的母矩阵，对于所有的 $i\ (1\le i\le N)$ ,你要求出宽为 i <script type="math/tex" id="MathJax-Element-700">i</script> 的最大回文子矩阵，即最大化回文子矩阵的高。
题解：
  首先回文的话就先对每一行跑一遍manacher，然后竖着要算每种数字最长有多长，这个可以用单调栈，栈里放向后延伸的元素，记下这个元素最早可以在的时间即可，当遇到小于等于栈顶的元素就表示这个长度告一段落，弹栈更新。

#include<bits/stdc++.h>
const int N = 3050;
int n, mx, mxi, r[N * 2], a[N * 2][N], f[N][20], g[N][20], b[N], t[N], top, ans[N];
char ch[N], c[N * 2];
int main() {
    freopen("symmetric.in", "r", stdin);
    freopen("symmetric.out", "w", stdout);
    scanf("%d", &n);
    for (int k = 1; k <= n; k++) {
        scanf("%s", ch);
        for (int i = 0; i < n; i++)
            c[i * 2 + 1] = ch[i],
            c[i * 2 + 2] = '*';
        c[0] = '@'; c[n * 2] = '#';
        mxi = 0; mx = 0;
        for (int i = 1; i < n * 2; i++) {
            if (mx >= i) r[i] = std::min(r[mxi * 2 - i], (mxi * 2 - i - (mxi - r[mxi]))); else r[i] = 1;
            while (c[i + r[i]] == c[i - r[i]]) r[i]++;
            if (i + r[i] - 1 > mx) mx = i + r[i] - 1, mxi = i;
        }
        for (int i = 1; i < n * 2; i += 2) if (!(r[i] & 1)) r[i]--;
        for (int i = 2; i < n * 2; i += 2) if ((r[i] & 1)) r[i]--;
        for (int i = 1; i < n * 2; i++) a[i][k] = r[i];
    }
    for (int k = 1; k < n * 2; k++) {
        top = 0;
        for (int i = 1; i <= n; i++) {
            int idx = i;
            while (top && a[k][i] <= b[top])
                ans[b[top]] = std::max(ans[b[top]], i - t[top]), idx = t[top], top--;
            b[++top] = a[k][i], t[top] = idx;
        }
        for (int i = 1; i <= top; i++)
            ans[b[i]] = std::max(ans[b[i]], n - t[i] + 1);
    }
    for (int i = n - 2; i > 0; i--)
        ans[i] = std::max(ans[i], ans[i + 2]);
    for (int i = 1; i <= n; i++)
        printf("%d\n", ans[i]);
    return 0;
}