bzoj1729 [Usaco2005 dec]Cow Patterns 牛的模式匹配 kmp+乱搞

题意：给定一个顺序序列表示奶牛的排列顺序，定义串a=b当且仅当a的每位的排名=b的每位的排名。（注意排名仅相对于当前区间）。

一开始想的是一边kmp一边维护排名，但是有点复杂，然后学习了一下更为简单的一种做法。
设rank[i][j]表示前i位数字j出现的次数，那么我们在kmp的时候就可以一边跑一边看排名是否相同，具体的话，是这样：
假设kmp到当前第i位，已经匹配了j位，那么我们要求i,j排名相同，那么怎么求排名呢。
首先算出在当前范围内小于s2[i]的个数，然后算出i位置上的s2[i]已经出现了几次，这就是i的排名，j同理，然后这些都可暴力算出来，S<=25完全不会超。
细节不懂看代码。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N = 1e5 + 5;
char s1[N], s2[N];
int rank1[N][30], rank2[N][30];
int next[N], ans[N];
int tot, n, m, s;
inline void get()
{
    int k = 0, i = 1;
    while (i <= m)
    {
        int x1 = 0, y1 = 0, x2 = 0, y2 = 0;
        fo(j, 1, s2[i] - 1)
            x1 += rank2[i][j] - rank2[i - k][j];
        y1 = rank2[i][s2[i]] - rank2[i - k][s2[i]];

        fo(j, 1, s2[k] - 1)
            x2 += rank2[k][j];
        y2 = rank2[k][s2[k]];
        if (k == 0 || (x1 == x2&&y1 == y2))i++, k++, next[i] = k;
        else k = next[k];
    }
}
inline void kmp()
{
    int i = 1, j = 1;
    while (i <= n)
    {
        int x1 = 0, x2 = 0, y1 = 0, y2 = 0;
        fo(k, 1, s1[i] - 1)
            x1 += rank1[i][k] - rank1[i - j][k];
        y1 = rank1[i][s1[i]] - rank1[i - j][s1[i]];
        fo(k, 1, s2[j] - 1)
            x2 += rank2[j][k];
        y2 = rank2[j][s2[j]];
        if (!j || (x1 == x2&&y1 == y2))j++, i++;
        else j = next[j];
        if (j == m + 1)
        {
            ans[++tot] = i - m;
            j = next[j];
        }
    }
}
int main()
{
    scanf("%d%d%d", &n, &m, &s);
    fo(i, 1, n)scanf("%d", &s1[i]);
    fo(i, 1, m)scanf("%d", &s2[i]);
    fo(i, 1, n)
    {
        fo(j, 1, s)rank1[i][j] = rank1[i - 1][j];
        rank1[i][s1[i]]++;
    }
    fo(i, 1, m)
    {
        fo(j, 1, s)rank2[i][j] = rank2[i - 1][j];
        rank2[i][s2[i]]++;
    }
    get();
    kmp();
    printf("%d\n", tot);
    fo(i, 1, tot)printf("%d\n", ans[i]);
    return 0;
}