/*
ID:53101304
PROG:"求q元有限域上的n级LFSR的理想初始向量C1--Cn"
LANG: C
TIME:2012/10/22
本程序只适用于q<10并且n<10的情况
*/
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<time.h>
/*计算矩阵周期 参数为初始友矩阵 n为级数 q为有限域元素个数
返回值为矩阵的周期
*/
int matrix_multi(int a[10][10],int n,int q)
{
int p[10][10];
int e[10][10];
int ans=(int)pow(q,n)+0.05;
int an =2;
//初始化e为矩阵Ap P为零矩阵
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
e[i][j]=a[i][j];
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
p[i][j]=0;
//an为矩阵乘法运行的次数
while(an<ans)
{
//矩阵乘法运算 最终值存在p中
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
for(int k=0;k<n;k++)
{
p[i][j]+=e[i][k]*a[k][j];
p[i][j]=p[i][j]%q;
}
//判断得到的p矩阵是否为单位阵 是返回an 否继续乘自己
bool temp = false;
bool temp1= true;
bool temp2 =true;
for(int i=0; (i<n)&&(temp1&&temp2);i++)
for(int j=0;j<n;j++)
{
if(i==j)
{
if(p[i][j]!=1) {temp1=false; break;}
}
else
{
if(!temp1) break;
else{
if(p[i][j]!=0){temp2=false;break;}
}
}
}
temp = temp1&&temp2;
if(temp) return an;
else an++;
//重置p为零矩阵 乘法得到的矩阵用e保存
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
e[i][j]=p[i][j];
p[i][j]=0;
}
}
return an-1;
}
/*
进行+1运算 q为进制(有限域中元素个数) 数组c储存位数值 n为位数即阶数
返回+1后MOD q的值 不同位上的数依旧储存在c中
如果超出所给位数n 返回0
*/
int add_1(int q,int *c,int n)
{
if(c[n-1]>q-1) return 0;
else
{
c[0]++;
for(int i=0;i<n-1;i++)
{
if(c[i]>q-1)
{
c[i]=c[i]%q;
c[i+1]++;
}
}
return 1;
}
}
int main()
{
int n=0;
int q=0;
printf("请输入LFSR的级数n:");
scanf("%d",&n);
printf("请输入有限域Fq的元素个数q:");
scanf("%d",&q);
//c保存初始向量的值
//Ap保存友矩阵
int c[10];
int Ap[10][10];
//初始化c
for(int i=0;i<n-1;i++)
{
c[i]=0;
}
c[n-1]=1;
//一种输出格式 无他
for(int i=1;i<=n;i++)
{
printf("c");
printf("%d ",i);
}
printf("\n");
//符合条件的友矩阵周期ans=q的n次方-1
int ans = (int)pow(q,n)-1+0.05;
while((c[n-1]<q))
{
//初始化Ap
for(int i=0;i<n-1;i++)
for(int j=0;j<n;j++)
{
if(j==i+1) Ap[i][j]=1;
else Ap[i][j]=0;
}
for(int j=0;j<n;j++)
Ap[n-1][j]=c[n-j-1];
//调用matrix_multi
int mm=matrix_multi(Ap,n,q);
//得到的友矩阵周期符合要求
if(mm== ans)
{
for(int i=0;i<n;i++)
printf("%d ",c[i]);
printf("\n");
}
//初始向量+1 再进行判断
int num = add_1(q,c,n);
if(!num) break;
}
//printf("Time use = %.2lf\n",(double)clock()/CLOCKS_PER_SEC);
return 0;
}
输出结果1:
输出结果2: