b-m算法常用来研究某一的01比特串的线性复杂度,换句话说通过BM算法可以得到01序列的最小线性多项式。
在介绍b-m算法之前,我先列出基本的公式
从公式就可以看出b-m算法是一个迭代算法,它的基本思想是利用前n-1个01序列所求出的最小线性多项式来求前n个01 序列的最小线性多项式,初始值容易明确,其重点在于整个算法如何如何迭代。
首先要明确的是公式中的加法并不是‘或’ 的意思,而是‘异或’运算,参数c,x和a之间是‘与’运算,x的n次方是指,第n号线性移位寄存器的值。
迭代有三种情况:
第一种情况比较特殊,就是前n个01序列全为0的情况,这时的L0,L1,...,Ln全为0,很容易得到fn(x)=1;
第二种情况就是当dn==0时, 这时应当认为fn-1(x)是可以推到出前n个01 序列的,即fn-1(x)==fn(x)
第三种情况则是dn==1,这时fn-1(x)不适用前n个01 序列,如何进行调整?先把问题放在这里,我们先讨论一下L,L是指fn(x)所能取到的最大的次数,当然由于fn(x) 可能会退化,fn(x)的次数可能会小于L。明白了这一点,就很容易知道L是单调不减的,同时它增长的时候必然是dn==1的。同样是两种情况,一种是前n-1个01序列全为0,第n个序列为1的情况,易得Ln=n,当然这是一种特殊情况; 同样Ln可以看成生成前n个序列所需要的最大寄存器数,那么这第二种情况就可以看成当fn(x)不适用前n+1个01 序列时,比较ln与n+1-ln,如果n+1-ln比较大,那就说明原来的那ln个寄存器已经不够用了,这时需要n+1-ln个寄存器。 ok,当我们在回来看m的取值时,就很容易的理解,m就是记录最新的一次L的更新,值得注意的是fm所求出的dn==1。回到最开始的问题,我们需要考虑两个方面,一个是如何把dn调整到0,另一个是如何不会破化前面已经修正的n-1序列。方法就是fn-1 + x^(n-m)*fm
个人能力有限,无法给出严格的数学证明,只能做一些直观的理解,如果有错误的地方,欢迎指正
附上c语言代码,采用c99标准,在命令行下运行,用了几个01序列验证,尚未发现bug
typedef unsigned int bin; //将接受的01字符串转化为一个无符号整形 bin strToBinary(char* str) { bin result=0; int len=strlen(str); for(int i=0;i<len;i++) { if(str[i]=='1') { result=(result<<1)+1; } else { result<<=1; } } return result; } //利用list序列,cn参数,以及ln+1长度求出dn(dn为0或1) bool getbn(bin list,bin cn,int seek) { bool bn=((cn>>seek) & (list>>seek)) & 1; while((--seek)>-1) { bn^=((cn>>seek) & (list>>seek)) & 1; } return bn; } //通过命令行接收01序列 int main(int argc,char* argv[]) { //简单处理得到的参数 if(2!=argc) { perror("error:parameter is not correct\n"); exit(22); } else if(strlen(argv[1])>32) { perror("error:parameter is too big\n"); exit(22); } //序列的长度 int N=strlen(argv[1]); //list为格式化的序列 bin list=strToBinary(argv[1]); int n=0,ln=0; bool bn; int m; //初始化fm,每一位分别是cln,cln-1,...,c1,c0 bin fm=0; //初始化f(x),每一位分别是cln,cln-1,...c1,c0 bin f=1; //判断l0,l1,...,ln是否全为0 bool flagIs0=true; //临时保存fn-1 bin fn_1; do { bn=getbn(list>>(N-n-1),f,ln); if(bn) { if(flagIs0) { m=n; fm=f; f=(1<<(n+1)) | 1; ln=n+1; } else { fn_1=f; f^=fm<<(n-m); if(ln<n+1-ln) { m=n; fm=fn_1; ln=n+1-ln; } } if(0!=ln) flagIs0=false; } n++; }while(n<N); //打印f printf("f(x)=1"); for(int i=1;i<N;i++) { if((f>>i)&1) { printf("+x^%d",i); } } printf("\n"); return 0; }