b-m算法

b-m算法原理及实现

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

  #include <stdio.h>
  #include <string.h>
  #include <stdlib.h>
  #include <stdbool.h>
  ​
  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;
  }