Fortran程序设计第六章循环结构下.pptVIP

* * integer m, i, k read *,m k=sqrt(m); do i=2, k if (mod(m, I)= =0) exit if ( i=k+1) then print*,m,’ is a prime number’ else print *,m,’ is not a prime number’ endif enddo end 例4：歌德巴赫猜想（将一个100~200的偶数分解成两个素数之和） 程序： PROGRAM MAIN IMPLICIT NONE INTEGER::N,N1,N2,K,KEND LOGICAL PRIME N: DO N=100,200,2 N1=1 DO IF(N13)THEN N1=N1+1 ELSE N1=N1+2 ENDIF PRIME=.TRUE. KEND=SQRT(REAL(N1)) KK1: DO K=2,KEND IF(MOD(N1,K)==0)PRIME=.FALSE. ENDDO KK1 IF(.NOT.PRIME) CYCLE N2=N-N1;PRIME=.TRUE. KEND=SQRT(REAL(N2)) KK2: DO K=2,KEND IF(MOD(N2,K)==0) PRIME=.FALSE. ENDDO KK2 IF(.NOT.PRIME) THEN CYCLE ELSE EXIT ENDIF ENDDO WRITE (*,(I3,=,I3,+,I3)) N,N1,N2 ENDDO N END 例5：对于-5≤x≤11；-10≤y≤9；-6 ≤z ≤ 18，求方程：x3+y3+z3=3的全部解。 程序设计有三种基本的思维方法： 枚举法、归纳法、抽象法 枚举法：就是逐一列举出可能解的各个元素，并加以判断，直到求得所需要的解。常用在排列、组合、数据分类、信息检索、多解方程的求解上； 程序如下： Integer:: x, y, z do x= -5,11 do y=-10,9 do z=-6,18 if (x*x*x+y*y*y+z*z*z= =3) then print *,x,y,z end if End do End do End do end 使用枚举法，必须掌握两条原则： 确定搜索的范围（这个范围必须是有限的）； 选择枚举的策略（按一条怎样的路径来逐一枚举）； 例6：求整数a和b的最小公倍数。 分析： （1）如果 i 是a 和b的最小公倍数，则 i 必能被 a 和 b 整除，同时 i 必须是自然数，所以其取值范围是１～∞； （2）方法一：i 从1开始，依次增加１，直到第一个能被a和b整除为止，这个 i 就是a和b的最小公倍数； （3）方法二：令从i 从a开始，并使 i 每次增加a而不是增加1，这就保证了 i 总是a的倍数。故每次只要判断 i 能否被b整除就可以了。一旦判断成立，i 就是a和b的最小公倍数； integer a,b, i read*,a,b i=0 do while .true. I=I+1 if (mod(I,a)= =0) then if (mod(I,b)= =0) printf*, I; eixt endif endif enddo end 程序一： 例7：百钱百鸡问题。用100元钱买100只鸡，每只公鸡5元，每只母鸡3元，每3只小鸡1元，要求每种鸡至少买一只，且必须是整只的，问各种鸡各买多少只？ 分析： （1）这是一个组合问题，归根到底是求三元一次方程的一组解； （2）设 i，j，k分别表示公鸡、母鸡和小鸡的只数。为了确定i，j，k的取值范围，可以有不同方法。不同的方法，程序的计算量相差甚远； （3）方法一：i ：1~ 20；j：1 ~ 33；k：1~100； （4）方法二：i ：1~ 18；j：1 ~ 31；k：100－i－j； （5）方法三：由i +j+k=100及5i +3j+k/3=100得14i+8j=200，由此可得： i ：1~ 13；j：1 ~ 23；k=100－i－j