- #1
- 25
- 0
..and tell what the last line they get in doing so.
Fortran:
program sieve_t
implicit none
integer*1 s(1000000), offset (10), sequence
integer i, j, n
n=0
sequence=4
do i=1, 10
offset(i)=0
enddo
offset(1)=2
offset(2)=4
offset(3)=2
do i=1, 1000000
s(i) = 1
enddo
do i=2, 1000000
if (s(i).eq.1) then
do j=2, (1000000/i)
s(i*j)=0
enddo
endif
enddo
do i=2, 1000000
if (s(i).eq.1) then
do j=1,sequence-1
if (s(i+offset(j)).ne.1) goto 10
enddo
n=n+1
write(*,*) n, i,i+2,i+6,i+8
endif
10 continue
enddo
end
Last edited by a moderator: