数论编程.docVIP

数 论 大 作 业 姓 名：叶明杰 任课教师：马华教授 导 师：齐小刚教授 学 号：1207122433 数论中一些基本问题的程序求解 下列算法的基本思想是课本第一章第一节定理七：设n1，若对所有的素数有,有,则n是素数。 判断一个数是不是素数的程序 #includestdio.h #includemath.h void main() { int i,n; printf(please input the number ); scanf(%d,n); for(i=2;i=sqrt(n);i++) /*判断n是不是素数的循环*/ if(n%i==0) break; /*n不是素数，跳出循环体*/ if (isqrt(n)) /*n是素数*/ printf(%d is a prime number\n,n); else printf(%d is not a prime number\n,n); } 运行结果： （2）求解n以内的所有素数 #includestdio.h #includemath.h void main() {int i,j,n; printf(请输入n的数值:); scanf(%d,n); printf(%d以内的所有素数为\n,n); for(i=2;in;i++) {for(j=2;j=sqrt(i);j++) /*判断小于n的数是不是素数*/ if(i%j==0) break; /*i不是素数，舍去*/ if(jsqrt(i)) /*i是素数，打印i*/ printf(%d ,i); } printf(\n); } 运行结果： （3）求解两个整数的最大公约数与最小公倍数 #includestdio.h void main() {int a,b,t,p,q,r; printf(please input two integers:); scanf(%d,%d,a,b); /*输入两个整数的数值*/ if(ab) t=a;a=b;b=t; /*当ab时，交换a,b的数值*/ p=a; q=b; r=p%q; /*r是a/b的余数*/ while(r) /*当余数不为0时，进行循环*/ {p=q; q=r; r=p%q; } printf(最大公约数为 %d\n,q); printf(最小公倍数为 %d\n,a*b/q); } 运行结果： （4）用爱拉托斯散筛选法求不超过n的素数 #includestdio.h #includemath.h void main() { int n,i,j,m,p,k,count=0,q=0,a[100]; printf(请输入n的数值:); scanf(%d,n); m=(int)sqrt((double)n); /*强制类型转换*/ printf(小于等于 %d 的素数为:,m); for(i=2;i=m;i++) { for(j=2;j=(int)sqrt((double)i);j++) if(i%j==0) break; /*i不是素数，舍去*/ if(j(int)sqrt((double)i)) /*i是素数*/ { printf(%d ,i); a[q]=i; /*把素数放入定义的数组中*/ q++; } } printf(\n); printf(小于等于 %d 的素数为:,n); for(p=0;pq;p++) { printf(%d ,a[p]); count++; /*计算sqrt(n)内的素数个数*/ } for(k=2;k=n;k++) { int flag=1; /*定义一个旗帜变量*/ for(p=0;pq;p++) { if(k%a[p]==0) /*删除sqrt(n)内素数的倍数*/ { flag=0; break; } } if(flag==1) /*旗帜变量为1，输出素数*/ { printf(%d ,k); count++; /*计算全部素数的个数*/ } } printf(\n); printf(小于等于 %d 的素数有 %d 个,n,count); printf(\n); } 运行结果： （5）合数的素分解 #includestdio.h #includemath.h void main() { int flag(int m); /*函数声明*/ int w=1,n,p,q; printf(请输入一个整数: ); scanf(%d,n); if(flag(n)==1) printf