计算代数与计算数论.doc

大素数的生成#include #include //大素数的生成int euild(int a,int b);main(){int n,i,t,b;printf("请输入一个奇数和安全系数(用逗号隔开)：");scanf("%d,%d",&n,&t);i=1; b=rand()%(n-3)+2;if(euild(b,n)==1){ for(;i<=t;i++) { b=rand()%(n-3)+2; if(euild(b^(n-1),n)!=1) break; }}if(i>t)printf("该奇数很可能是一个素数！\n");else{printf("该奇数是不一个素数！\n它对%d不成立",b);}}int euild(int a,int b){ int r[3],s[3],t[3],z,q; if(a0 i=i+1; n(i)=mod(N,2); N=fix(N/2);endif n(1)==0 a(1)=1;else a(1)=mod(b,m);endb(2)=mod(b*b,m);for j=2:1:i if n(j)==0 a(j)=a(j-1); else a(j)=mod(a(j-1)*b(j),m); end b(j+1)=mod(b(j)*b(j),m);endanswer=a(j) 求两个数的最大公因数#include main(){ int r[3],s[3],t[3],a,b,z,q; printf("请输入两个整数："); scanf("%d,%d",&a,&b); if(a0 i=i+1; n(i)=mod(N,2); N=fix(N/2);endif n(1)==0 a(1)=1;else a(1)=mod(b,m);endb(2)=mod(b*b,m);for j=2:1:i if n(j)==0 a(j)=a(j-1); else a(j)=mod(a(j-1)*b(j),m); end b(j+1)=mod(b(j)*b(j),m);end y=a(j);end。