[Fortran] 纯文本查看 复制代码
MODULE my_Constants
implicit none
real(8),PARAMETER :: df=1
real(8),PARAMETER :: fmax=196078
real(8),PARAMETER :: f00=191693
real(8),PARAMETER :: f01=191694
real(8),PARAMETER :: f02=191696
end MODULE my_Constants
use my_Constants
implicit none
integer i,j,k,l,m,n,s,z,a,b,c,d
real(8),dimension(4385)::fijk,f
open(11,file='useful.dat')
do i=1,4385
f(i)=0.0
enddo
f(1)=f00
f(2)=f01
f(3)=f02
n=4
z=1
write(11,*)f(1),z
write(11,*)f(2),z
write(11,*)f(3),z
do 10 i=4,4385
m=n
f(m)=f00+i*df
do j=1,m
do k=1,m
do l=1,m
fijk(l)=f(j)+f(k)-f(l)
do 50 s=1,m
if((fijk(l).ne.f(s)).and.(fijk(l).le.fmax).and.(fijk(l).ge.f00)) then
write(11,*)f(m),z
n=m+1
goto 10
endif
50 continue
enddo
enddo
enddo
10 continue
end