poj2185 Milking Grid

关于如何求一个(不完整)字符串的最小循环节：

n - nex[n - 1]即可。

若n % 循环节长度 == 0则表示这个字符串可由若干个循环节构成。

否则最后一个循环节有残缺。

回到本题：

首先行列分开。

一个直观的想法是求出最小循环节的lcm，但是诸君请看：

2 8

abDQRVab

aaaabaaa

然后怎么做呢？

有一种做法是枚举行的长度，然后列用KMP

然后，略加思索，我们行也可以用KMP啊......

然后双KMP搞定。

调了好久是因为把while写成了if......

 

#include <cstdio>

const int N = 10010, M = 100;

char s[N][M];
int nex[N], m, n;

inline bool equala(int a, int b) {
    for(int i = 0; i < m; i++) {
        if(s[a][i] != s[b][i]) {
            return 0;
        }
    }
    return 1;
}

inline bool equalb(int a, int b) {
    for(int i = 0; i < n; i++) {
        if(s[i][a] != s[i][b]) {
            return 0;
        }
    }
    return 1;
}

inline int geta() {
    nex[0] = 0;
    for(int i = 1, j = 0; i < n; i++) {
        while(j && !equala(i, j)) {
            j = nex[j - 1];
        }
        if(equala(i, j)) {
            j++;
        }
        nex[i] = j;
    }
    return n - nex[n - 1];
}

inline int getb() {
    nex[0] = 0;
    for(int i = 1, j = 0; i < m; i++) {
        while(j && !equalb(i, j)) {
            j = nex[j - 1];
        }
        if(equalb(i, j)) {
            j++;
        }
        nex[i] = j;
    }
    return m - nex[m - 1];
}

int main() {
    while(scanf("%d%d", &n, &m) != EOF) {
        for(int i = 0; i < n; i++) {
            scanf("%s", s[i]);
        }

        int c = geta();
        int d = getb();
        printf("%d\n", c * d);
    }
    return 0;
}