B-M求线性移位寄存器

1、B-M算法求线性综合解的过程
2、假设a(11)=(00100011101)是二元域GF(2)上的一个长度为11的序列，试用B-M算法求其线性综合解。

#include<iostream>
using namespace std;
void B_M(int a,int nn)
{
    int i,j,k,n0;
    int c=1,m,l[15],d=0,d1,cc[15];
    cc[0]=c;
    for(i=0;i<nn;i++)//初始化
    {
        if(a>>i&1)
        {
            n0=i;
            d=d^((a>>i&1)<<i);
            l[i+1]=n0+1;
            break;
        }
        else
        {
            l[i+1]=0;
            cc[i+1]=1;
        }
    }
    c=c^(1<<n0+1);
    cc[n0+1]=c;
    m=n0-1;
    for(i=n0+1;i<nn;i++)
    {
        d1=a>>i&1;
        for(j=1;j<=l[i];j++)
            d1=d1^((cc[i]>>j&1)&(a>>i-j&1));
        d=d^(d1<<i);
        if(d1==0)
        {
            cc[i+1]=cc[i];
            l[i+1]=l[i];
        }
        else
        {
            if(l[i]<(i+1-l[i]))
                l[i+1]=i+1-l[i];
            else
                l[i+1]=l[i];
            if(l[i]>l[i-1])
                m=i-1;
            cc[i+1]=cc[i]^(cc[m]<<(i-m));
        }
    }
    cout<<"a("<<nn<<")=(";
    for(i=0;i<nn;i++)
        cout<<(a>>i&1);
    cout<<")的线性综合解为："<<endl<<"f(x)=1";
    for(j=1;j<=nn;j++)
        if(cc[nn]>>j&1)
            cout<<"+"<<"x^"<<j;
    cout<<endl;
}
int main()
{
    int a=0,nn;
    int b;
    cout<<"请输入序列长度n：";
    cin>>nn;
    cout<<"请输入序列a：";
    for(int i=0;i<nn;i++)
    {
        cin>>b;
        a=a^(b<<i);
    }
    B_M(a,nn);
    return 0;
}