您现在的位置:计算机技术学习网 > 技术中心 > 编程开发 > C >

梅西迭代算法的实现

来源:网络收集 责任编辑:栏目编辑 发表时间:2013-07-02 05:30 点击:
/*
  Name:  梅西迭代算法的实现   
  Copyright: 2003(C)
  Author:  徐岩柏
  Date: 31-10-03 16:09
  Description: 该问题是线性移位寄存器的综合问题提出的,给定一个N长的
               二元序列,如何求出产生这一序列的级数最小的线性移位寄存
               器,即最短的线性移位寄存器
*/



#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <fstream>
using namespace std;

int main(int argc, char *argv[])
{
  typedef vector<bool> MyArray; file://定义自己的数据类
  MyArray a;   file://状态数组 ,用来保存给定的m序列?
  vector<int>l;   file://线性移位寄存器的级数数组
  MyArray temp;
  vector<MyArray> f; file://为线性移位寄存器对应的多项式
  int N =0;
  ifstream infile("input.txt"); file://从文件中读取数据
  ofstream outfile("output.txt"); file://把结果写入文件中。
  if(!infile)
  {
    cout << "Aata file input.txt is not ready! Please check it"<<endl;
    exit (-1);
  }
  string strTemp;
  infile>>strTemp; file://读入N
 // cout << strTemp <<endl;
  N = atoi(strTemp.c_str());
  infile>>strTemp; file://读入序列
  for( int i = 0 ;i<N; ++i)
  {
     a.push_back(strTemp[i]=='1'); file://把序列保存
  }
 file://cout << strTemp <<endl;
 int n=0;
 l.push_back(0);
 temp.push_back(1); file://形成f0
 f.push_back(temp); file://把f0加到f数组中。
 do{
     int dn =0;
     for (int i = 0;i<=l[n];++i)
     {

       dn += f[n][i]*a[n-i];   file://计算dn
     }
     dn = dn % 2; file://取余
  //   cout <<"d"<<n<<"= "<<dn <<endl;
     int isum = 0;
     for (int i =0;i<=n ;++i)
     {
          isum +=l[i]; file://判定是不是所有的ln都等于0就是把所有的ln加起来看是不是零
     }
     if( dn == 0) file://如果dn =0
     {
          f.push_back(f[n]); file://fn+1 = fn
          l.push_back(l[n]); file://ln+1 = ln
     }
     else if (!isum) file://所有的ln都是零
     {
          temp.clear();
          for(int i =0;i<n+2;++i)
          {
               temp.push_back(0);
          }
          temp[0] = 1;
          temp[n+1] = 1;
          f.push_back(temp);
          l.push_back(n+1);
     }
     else  file://lm<lm+1 = ...=ln
     {
 

    相关新闻>>

      发表评论
      请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
      用户名: 验证码:点击我更换图片
      最新评论 更多>>
      网站首页 - 友情链接 - 网站地图 - TAG标签 - RSS订阅 - 内容搜索
      Copyright © 2008-2015 计算机技术学习交流网. 版权所有

      豫ICP备11007008号-1