《应用密码学》求给定n阶LFSR在指定q元域上的理想初始向量集合C

/*
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：