算法竞赛进阶指南：0x18：奶牛矩阵

每天早上，农夫约翰的奶牛们被挤奶的时候，都会站成一个 R 行 C 列的方阵。

现在在每个奶牛的身上标注表示其品种的大写字母，则所有奶牛共同构成了一个 R 行 C 列的字符矩阵。

现在给定由所有奶牛构成的矩阵，求它的最小覆盖子矩阵的面积是多少。

如果一个子矩阵无限复制扩张之后得到的矩阵能包含原来的矩阵，则称该子矩阵为覆盖子矩阵。

输入格式

第 1 行：输入两个用空格隔开的整数，R 和 C。

第 2..R+1 行：描绘由奶牛构成的 R 行 C 列的矩阵，每行 C 个字符，字符之间没有空格。

输出格式

输出最小覆盖子矩阵的面积。（每个字符的面积为 1）

数据范围

1≤R≤10000,
1≤C≤75

输入样例：

2 5
ABABA
ABABA

输出样例：

2

提示

样例中给出的矩阵的最小覆盖子矩阵为 AB，面积为 2。

 

思路：暴力枚举+kmp

求宽度：暴力枚举每个长度，测试是否匹配，由此可求出最小循环长度

求高度：kmp，根据求出的宽度将每一行变为一个字符串(0~宽度-1)，

将这一行字符串作为一个字母进行行的kmp匹配，求出行的最小循环节，可求出面积

 

#include<iostream>
#include<cstring>
using namespace std;

const int N=10010,M=80;

char s[N][M];

bool st[M];//st[i]表示长度为i的循环节，值为true表示不可用，为false表示可用

int r,c,ne[N];//ne表示kmp的next数组

int main()
{
    cin>>r>>c;
    
    for(int i=1;i<=r;++i)//对于每一列，暴力枚举每个长度
    {
        scanf("%s",s[i]);
        for(int j=1;j<c;++j)//枚举长度
        {
            if(st[j]) continue;//已经不可用直接跳过
            
            for(int k=j;k<c;k+=j)//k表示以j为循环长度，每个循环的下标
                for(int u=0;u<j&&u+k<c;++u)//u表示初始循环，u+k表示该匹配的循环
                    if(s[i][u]!=s[i][u+k]) st[j]=1;
        }
    }
    
    int weight;//求出最小长度
    for(int i=1;i<=c;++i)
        if(!st[i]) {
            weight=i;
            break;
        }
    
    for(int i=1;i<=r;++i)   s[i][weight]='\0';//将每一行循环外变为字符串
    
    for(int i=2,j=0;i<=r;++i){//行的kmp匹配
        while(j&&strcmp(s[i],s[j+1])) j=ne[j];
        if(!strcmp(s[i],s[j+1])) j++;
        ne[i]=j;
    }
    
    int height=r-ne[r];//最小循环节就是高度
    
    cout<<height*weight;
}