7.26 机房欢乐赛 T1 无尽的矩阵 （hash + KMP）

1 无尽的矩阵(matrix.c/cpp/pas)

1.1 题目描述

从前有一个的小矩阵，矩阵的每个元素是一个字母(区分大小写)，突然有一天它发生了变异，覆盖了整个二维空间，即不停自我复制产生相同的矩阵然后无隙放置。现在二维空间已经被它占领了，但你只被告知了大小为R*C空间的内容(可能包含不完整的原矩阵)，为了将它恢复原状，你需要找到满足条件的面积最小的原矩阵。
奇怪的是，同时有 T 个二维空间发生了变异，你需要尽快解决这些变异。

1.2 输入格式

第一行为一个整数T，表示二维空间数目。
接下来T组数据。每组数据第一行包含两个数 R，C，表示你被告知的空间大小；接下来 R 行，每行包含 C 个字母，表示你被告知的空间内容。

1.3 输出格式

对于每一组数据输出一行，每行只包含一个数，表示最小的原矩阵面积。

1.4 样例输入

2
2 5
ABABA
ABABA
2 8
ABCDEFAB
AAAABAAA

1.5 样例输出

2
12

1.6 数据范围与约定

对于前20％的数据R<=20，C<=20；
对于前40％的数据R<=400，C<=100；
对于100％的数据R<=5000 ，C<=100，T<=50。

1.7 题解

将每一行hash 为一个数，对得到的新数组直接跑KMP 求最小循环节长度，列
同理。将两次求得的最小循环节长度相乘即为答案。这就是std 做法。

1.8 代码

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

char s[5010][150];
int t, n, m, ans1, ans2, nxt[5010], h[5010];

unsigned int ELFhash1(int a) {
    unsigned int hash = 0, x = 0, p = 0;
    while(p < m) {
        hash = (hash << 4) + s[a][p];
        if((x = hash & 0xf0000000)) {  
            hash ^= (x >> 24);
            hash &=~ x;
        }
        p ++; 
    }
    return (hash & 0x7fffffff);
}

unsigned int ELFhash2(int a) {
    unsigned int hash = 0, x = 0, p = 0;
    while(p < n) {
        hash = (hash << 4) + s[p][a];
        if((x = hash & 0xf0000000)) {  
            hash ^= (x >> 24);
            hash &=~ x;
        }
        p ++;
    }
    return (hash & 0x7fffffff);
}

int kmp(int x) {
    int i = 0, j = -1;
    nxt[0] = -1;
    while(i <= x) {
        if(j == -1 || h[i] == h[j]) {
            i ++, j ++;
            nxt[i] = j;
        }
        else j = nxt[j];
    }
    return x - nxt[x];
}

int main() {
    freopen("matrix.in", "r", stdin);
    freopen("matrix.out", "w", stdout);
    scanf("%d", &t);
    while (t --) {
        ans1 = 0, ans2 = 0;
        memset(s, 0, sizeof(s));
        scanf("%d %d", &n, &m);
        for(int i = 0; i < n; i ++)
            scanf("%s", s[i]), h[i] = ELFhash1(i);
        ans1 = kmp(n);
        for(int i = 0; i < m; i ++)
            h[i] = ELFhash2(i);
        ans2 = kmp(m);
        printf("%d\n", ans1 * ans2);
    }
    return 0;
}