POJ 2185

题目链接：

　　http://poj.org/problem?id=2185

题目描述：

　　给定一个R * C的矩阵，在该矩阵中寻找一个子矩阵，要求子矩阵复制多次能够包含原矩阵。

解题思路：

　　如果在字符串中找到一个最小子串使其复制多次包含原字符串，那么我们可以很容易想到通过kmp算法的next数组来得出该子串长度，即len - next[len]。

　　那么回到这道题，要在矩阵中找到合法的最小子矩阵，其实跟找最小子串的原理是一样的，只需要对矩阵的行和列分别看成一个字符串用kmp算法过一遍，然后分别得到一个最小行数，该行数复制多遍可以包含原矩阵，以及一个最小列数，该列数复制多遍同样可以包含原矩阵，那么题目要求的子矩阵宽和高就这样得到了，最后相乘一下输出就搞定啦。

代码：

#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <map>
#include <vector>
#include <algorithm>
using namespace std;

#define ll long long
#define INF 0x3f3f3f3f3f
#define MAX 

char rectangle[10010][100];
int fail[10010];
int R, C;

/*
判断矩阵行与行之间是否相等
 */
bool cmp1(int j, int k)
{
    for(int i = 0; i < C; ++i) {
        if(rectangle[j][i] != rectangle[k][i])
            return 0;
    }
    return 1;
}

/*
判断矩阵列与列之间是否相等
 */
bool cmp2(int j, int k)
{
    for(int i = 0; i < R; ++i) {
        if(rectangle[i][j] != rectangle[i][k])
            return 0;
    }
    return 1;
}

int main()
{
    //freopen("debug\in.txt", "r", stdin);
    //freopen("CON", "w", stdout);
    int i, j, k;
    scanf("%d%d", &R, &C);
    for(i = 0; i < R; ++i) {
        scanf("%s", rectangle[i]);
    }
    fail[0] = k = -1;
    j = 0;
    //kmp获取next数组的超简洁版写法
    while(j < R) {
        while(k != -1 && !cmp1(j, k)) k = fail[k];
        fail[++j] = ++k;
    }
    int ans1 = R - fail[R]; //获得子矩阵高
    fail[0] = k = -1;
    j = 0;
    while(j < C) {
        while(k != -1 && !cmp2(j, k)) k = fail[k];
        fail[++j] = ++k;
    }
    int ans2 = C - fail[C]; //获得子矩阵宽
    printf("%d
", ans1 * ans2);
    return 0；
}