Fortran学习笔记(5)

!循环结构程序设计
!do循环结构：用来实现循环次数确定的循环
do 循环变量v=初值e1,终值e2,步长e3
 循环体
enddo


!举例:
!1.编写程序输出1-999之间的所有"水仙花数(各位数字的立方和等于该数本身)"
integer a,b,c,n
do n=1,999
 a=n/100
 b=mod(n/10,10)
 c=mod(n,10)
 if(a**3+b**3+c**3==n) print*,n,'是水仙花数'
enddo
end


!do循环结构嵌套
!求1!+2!+...+n!
integer i,j,n,sum,np
print*,"请输入自然数值n:"
read*,n
sum=0
do i=1,n
 np=1
 do j=1,i
  np=np*j
 enddo
 sum=sum+np
enddo
print*,"自然数n的阶乘和是:",sum
end


!隐含do循环结构
print*,(w,v=e1,e2,e3)
read*,(w,v=e1,e2,e3)


!举例:
!1.打印9*9乘法表
integer i,j
do i=1,9
 print 10,(j,'*',i,'=',i*j,j=1,i)
 10 format(<i>(i1,a,i1,a,i2,2x))
enddo
end


!隐含do循环结构的嵌套
print*,((w,v1=e11,e12,e13),v2=e21.e22,e23)
read*,((w,v1=e11,e12,e13),v2=e21.e22,e23)


!do while循环结构
do while(逻辑表达式)
 循环体
enddo


!循环的流程控制:exit和cycle语句
!exit语句:直接跳出正在运行的循环，转到循环结构enddo后的下一个语句执行
!do和do while都可以使用
!cycle语句:略过循环体中在cycle后面的所有语句，直接跳回循环的开头进行下一次循环


!举例:
!1.打印出1-10内的数字，3和6不打印
integer i
do i=1,10
 if(i==3.or.i==6) cycle
 print*,i
enddo
end


!程序练习
!1.求[2,10000]内的守形数
!守形数:该数平方的低位数等于它本身
integer i,j,k,l,n,s
do n=2,10000
 s=n**2
 if(n/10==0)then
  if(mod(s,10)==n) print*,n
 else if(n/100==0)then
  if(mod(s,100)==n) print*,n
 else if(n/1000==0)then
  if(mod(s,1000)==n) print*,n
 else
  if(mod(s,10000)==n) print*,n
endif
enddo
end


!2.求2-100之间的素数
integer i,j,n,s
print*,"2-100间的素数有:"
print*,2
do n=3,100
 j=int(sqrt(1.0*n))
 do i=2,j
  s=mod(n,i)
  if(s==0) exit
 enddo
 if(s/=0) print*,n
enddo
end


!3.编写程序打印以下图形(n=5):
5 4 4 4 4 4 4 4 5
  4 3 3 3 3 3 4
    3 2 2 2 3
      2 1 2
        1
integer:: i,j,k,l=-1,n
read*,n
do k=n,1,-1
 l=l+1
 if(k>1)then
  print*,(' ',i=k,n+l),char(k+48),' ',(char(k-1+48),' ',j=1,k-1),&
  &(char(k-1+48),' ',j=2,k-1),char(k+48)
 else
  print*,(' ',i=1,n+l),char(k+48)
 endif
enddo
end