C语言程序设计第5章(姜恒远著).pptVIP

主要内容 5.1 素数判断 5.2 最大公约数求解 5.3 穷举法求满足条件的一组解 5.4 级数近似计算 5.5 一元非线性方程求根 *5.6 定积分近似计算 5.1 素数判断 5.1.1 最简单素数判断算法 5.1.2 改进后的素数判断算法 5.1.2 改进后的素数判断算法 算法描述 如果n2且n不能被2～sqrt(n)之间的所有整数整除，n就是素数。 程序实现 r=sqrt(n); for(t=1,i=2; ti=r; i++) if(n%i==0) t=0; if(t)printf(“%d is prime”,n); 5.2 最大公约数求解 最大公约数 能同时被x和y整除的最大整数是整数x和y的最大公约数（GCD） 。 例如：30和45的公约数有3、5、15，其中15是30和45的最大公约数。 5.2 最大公约数求解 5.2.1 brute-force算法 5.2.2 欧几里德算法 5.2.1 brute-force算法 算法描述 (1) 设x（或y）是x、y的最大公约数z； (2) 判断x和y是否均能被z整除。若z不能同时整除x和y则z-1→z，重复S2步。否则，做下一步操作； (3) 输出（或返回）z。 5.2.1 brute-force算法 #include stdio.h int main(void) { int x=30,y=45,z; z=x; while(!(x%z==0y%z==0)) z--; printf(“%d”,z); return 0; } 5.2.2 欧几里德算法 算法描述 (1) 判断x除以y的余数r是否为0。若r为0则y是x、y的最大公约数，继续做下步操作；否则y→x，r→y重复做第(1)步。 (2) 输出（或返回）y 5.2.2 欧几里德算法 #include stdio.h int main(void) { long x=100000,y=100005, r; while((r=x%y)!=0){ x=y; y=r; } printf(“%ld”,y); return 0; } 5.3 穷举法求满足条件的一组解 例1：某人在纸上写了一个四位整数3a45（a代表一个数字）。已知这个四位整数被3除后的值是1115，请问这个四位整数是什么？ 算法描述 用a的所有可能取值（0?9）分别形成10个四位整数3045、3145、…、3945 依次判断这10个四位整数中的每一个被3除后的值是否1115，若是则输出。 5.3 穷举法求满足条件的一组解 #include stdio.h int main(void) { int a,b; for(a=0;a=9;a++){ b=3*1000+a*100+45; if(b/3==1115) printf(“%d”,b); } } 5.3 穷举法求满足条件的一组解 例2 中国古代数学著作《算经》中提出一个问题：公鸡每只5钱，母鸡每只3钱，小鸡1钱3只。若用100钱买100只鸡，可以有多少种买法？ 用x表示可以买到的公鸡数量、y表示可以买到的母鸡数量、z表示可以买到的小鸡数量。 则： x+y+z=100， 5x+3y+z/3=100 5.4 级数近似计算 一个函数可以近似地表示为一个收敛的幂级数。因此函数求值问题可以归结为级数项求和问题。例如 ： 用计算机计算无穷级数的值时只能计算有限个级数项的和。因此，用程序求得的级数值只能是在给定误差范围内的近似值。即：计算级数前n个项之和，当通项（第n项）的绝对值小于所给误差时停止级数项累加计算。 5.4 级数近似计算 5.4.1 简单方法 5.4.2 递推法 5.4.1 简单方法 程序实现 计算ex级数展开式的近似值，误差为键盘输入的eps。 5.4.2 递推法 算法描述 在计算级数的各项值时从前项结果推出后项结果 例如： 将级数： 表示为： 则： a1=1， an=an-1 x/(n-1） （n=2，3，4……） 5.4.2 递推法 采用递推法计算ex的级数展开式近似值，允许误差在eps范围内。 5.5 一元非线性方程求根 5.5.1 牛顿迭代法 5.5.2 二分法和弦截法 5.5.1 牛顿迭代法 牛顿迭代公式 xn+1= xn－f(xn)/f (xn) 求方程f(x)=0在x0附近的一个近似实根。 5.5.1 牛顿迭代法 (1)