BM算法总结

BM算法，可以求一个数列的最短递推式。
采用增量法，依次考虑每个数：
若在这个位置上正确，则忽略；
否则，类似拉格朗日插值法，找一个满足在前面位置都为0，这个位置上不为0的递推式，进行修补。
每当我们遇到一个这样的位置时，我们都可以得到一个这样的递推式：用目前的递推式，在0位置（即这个位置上），增加一个-1。相当于积累了一个“经验”。
这样，我们在所有的“经验”中，找到一个，并在它的前面补0，相当于移位。
若有多个，取补0后长度最小的。
时间复杂度：\(O(n(n+m))\)。（n为递推式长度，m为输入长度）。
代码：

int BM(int sz[10010],int n)
{
    int m=0;
    for(int i=1;i<=n;i++)
    {
        cd[i]=inf;
        int he=(md-sz[i])%md;
        for(int j=1;j<=m;j++)
            he=(he+1ll*ans[j]*sz[i-j])%md;
        if(he==0)
            continue;
        int ny=ksm(he,md-2);
        sb[i][0]=(md-ny)%md;
        for(int j=1;j<=m;j++)
            sb[i][j]=1ll*ans[j]*ny%md;
        cd[i]=m;int mi=inf,wz;
        for(int j=0;j<i;j++)
        {
            int t=cd[j]+(i-j);
            if(t<=mi)
                mi=t,wz=j;
        }
        for(int j=0;j<=cd[wz];j++)
            ans[i-wz+j]=(ans[i-wz+j]-1ll*sb[wz][j]*he%md+md)%md;
        if(mi>m)m=mi;
    }
    return m;
}