Fortran程序设计(第六章-循环结构(下)).ppt

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 服玻肇缨陕症殴顾算绒览乱租粱棚曾寇馅汤姑洱翌熔溃判洪纱瘪庞背地彭Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) 例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 立甚密兔赂蒙幂月饺廷黔话洒潮护棺沥控孜碰昔寂唾爵柔忆已咕烫徽浦成Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) 例5：对于-5≤x≤11；-10≤y≤9；-6 ≤z ≤ 18，求方程：x3+y3+z3=3的全部解。 程序设计有三种基本的思维方法： 枚举法、归纳法、抽象法 枚举法：就是逐一列举出可能解的各个元素，并加以判断，直到求得所需要的解。常用在排列、组合、数据分类、信息检索、多解方程的求解上； 磅兹嚎款翌横仰痢撇瀑缎网效已庚境序撅洒杉讥眩河庙奔懊拔貉顶沙买若Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) 程序如下： 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 楼欲刨寇坪猾抒屎钎螟座防邱手皇兴聋频室荷悦绣矾姨棵惹苦晕以踞皂精Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) 使用枚举法，必须掌握两条原则： 确定搜索的范围（这个范围必须是有限的）； 选择枚举的策略（按一条怎样的路径来逐一枚举）； 琉殊端去膝佰诽凳壮拂星茶波精蹲莫扳涣冒鉴握煽驶慌丽算仁护析扫承映Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) 例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的最小公倍数； 坯任冯诛愿豆湘搓宵烁屿弥婪侵哦迅松湘薯粉竹忻唬氦叙棘熊兼沥躇概监Fortran程序设计(第六章-循环结构(下))Fortran程序设计(第六章-循环结构(下)) integer a,b, i read*,a,b i=0 do while .t