素数判定与程序设计基础.pdfVIP

素数判定

§素数:整数p1,只有1和p自身能整除p,p为素数

如:2,3,5,7…

§合数:大于1且不是素数的数

如:4,6,8,9…

定理：如果a是合数,则a必有小于等于a的

真因子。

证：如果a是合数，则a可表示成abc,其中

1ba,1ca.显然,b和c中必有一个小于等

于a，否则,bc(a)2a,。

素数判定算法

/**********************************************

功能：判定x是否是素数，x≥2

返回值：返回1表示是素数，返回0则不是素数

***********************************************/

intis_prime(intx)

{

intm;

if(x2)return1;

if(x%20)return0;/*偶数不是素数*/

m=sqrt(x);

for(i=3;i=m;i++)

if(x%i0)/*x能被i整除，x不是素数*/

return0;

return1;/*循环结束还没有返回，说明x是素数*/

}

【例3-17】编一程序打印出2至99之间的所有素数。

分析：本例是上面素数判定算法的一个简单应用

#includestdio.h

intis_prime(intx);

voidmain()

{

for(i=2;i100;i=i+1)

if(is_prime(i))

printf(%3d,i);

}

/*上机调试时，将is_prime函数定义放在这儿*/

运行结果

235711131719232931374143475359616771

7379838997

素数设计程序求解

p27

1.梅森数:(2-1),p为素数,如:2-1=3,2-1=127

2.孪生素数:两个素数的差值是2,如2-99之间的孪生素

数

3.可逆素数:一个素数将其各位数字的顺序倒过来构成的反序

数也是素数,如:1009,1021

4.回文素数:一个素数从左到右和从右向读的结果相同且是素

数,如:101,131,151,181,191

5.歌德巴赫猜想:任何大于4的偶数都可以表示成两个奇素数的

和,如:1978=5

练习直接思路

穷举法

编程求解两个数a，b的最大公约数

如，(8,12)=4

编程求解两个数a，b的最小公倍数

如，LCM(20,15)=60

求最大公约数的算法--辗转相除法

一个非常著名的古老算法是用于求两个整数的最大公约数

的欧几里德算法。欧几里德算法也称为辗转相除法。

递推：

(a,b)=(b,a%b)

这个的含义是a与b的最大公约数等于b与a%b

的最大公约数。一直递推下去，直到(rn，0)，

此时rn即为所要求的最大公约数。

计算91和52的最大公约数，求解过程如下：

（1）mod(91,52)=39

（2）mod(52,39)=13

（3）mod(39,13)=0

所以，13就是91和52的最大公约数。

自然语言表示的欧几里德算法如下：

输入：两个正