求离散对数
#include
#define MAX 100
void main()
{ int i,j, num,sum=0, m[MAX]={0},mm[MAX]={0},mm1[MAX]={0},mmm=1,a,n,q,p,b;
int theone(int ,int ,int,int,int);
int poww(int ,int);
printf("***********************************************************\n");
printf("a^cmodp q(为分解得到的素因子) n（为素因子的幂）\n");
printf("测试 p=41 a=6 b=29 求得离散对数c=7\n");
printf("p=37 a=2 b=29 求得离散对数c=21\n");
printf("***********************************************************\n");
printf("分成几个素因子,请输入素因子的个数：\n");
scanf("%d",&num);
for(i=0;i
{
printf("b^cmodp;a^cmodp abp q n输入五个数值 以次是：\n");
scanf("%d%d%d%d%d",&a,&b,&p,&q,&n);
mm[i]=(p-1)/poww(q,n);
for(j=0;j
if(j*mm[i]%poww(q,n)==1)
mm1[i]=j;
m[i]=theone(a,b,p,q,n);
}
for(i=0;i
{ mmm*=mm[i];
sum=sum+mm[i]*mm1[i]*m[i];
}
printf("%d\n",sum%mmm);
}
int theone(int a,int bb,int p,int q,int n)//求出每个方程的解
{
int i=0,j=0,k,k1=0,ani,c,aa[MAX]={0},b[MAX]={0};
int mod(int,int,int);
int poww(int ,int);
b[0]=bb;
c=(p-1)/q;
if(mod(b[0],c,p)==mod(a,0,p))
{ aa[0]=0;
}else {
for(i=0;mod(b[0],c,p)!=mod(a,c*i,p);)
{
i++;
aa[0]=i;
}
}
if(n>1)
{
for(k=1;k
if(a*k%p==1)
ani=k;
for(j=0;j
{
b[j+1]=b[j]*poww(poww(ani,aa[j]),poww(q,j))%p;
if(mod(b[1+j],c/poww(q,j+1),p)==mod(a,0,p))
{ aa[j+1]=0;
} else{
for(i=0;mod(b[j+1],c/poww(q,j+1),p)!=mod(a,c*i,p);)
{
i++;
aa[j+1]=i;
}
}
}
}
for(i=0;i
k1=k1+aa[i]*poww(q,i);
return(k1);
}
int mod(int x,int r,int n)//模乘法
{
int i=0,temp,a[30],first=1;
for(;i<30;i++)
{temp=r%2;
a[i]=temp;
r=r/2;
if(r==0)
break;
}
for(;i>=0;i--)
{
first=first*first %n;
if(a[i]==1)
first=(first*x)%n;
}
return(first);
}
int poww(int x,int y)//指数函数
{
int i,sum=1;
if(y==0)
return(1);
else{
for(i=0;i
sum=sum*x;
return(sum);
}