Skip to content

Instantly share code, notes, and snippets.

@komasaru
Created April 20, 2017 05:04
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save komasaru/11a2513001796fafd92de6f677732ec7 to your computer and use it in GitHub Desktop.
Save komasaru/11a2513001796fafd92de6f677732ec7 to your computer and use it in GitHub Desktop.
Fortran code to compute prime numbers.
!****************************************************
! 素数一覧
! : 入力値以下の素数を全て出力する
!
! date name version
! 2017.04.20 mk-mode.com 1.00 新規作成
!
! Copyright(C) 2017 mk-mode.com All Rights Reserved.
!****************************************************
!
program prime_numbers
implicit none
integer i, j, n
logical is_prime
print *, "自然数を入力してください N:"
read *, n
do i = 2, n
is_prime = .true.
do j = 2, int(sqrt(dble(i)))
if (mod(i,j) == 0) then
is_prime = .false. ! 割り切れるので素数ではない
exit
end if
end do
if (is_prime) print *, i ! もし素数ならば出力
end do
end program
@taymazalizadeh
Copy link

taymazalizadeh commented Feb 22, 2022

What change do we have to make in the code if we want to distinguish prime numbers from the numbers we enter?

program prime number
implicit none
integer::n,s,i
read*,n
s=0
do i=1,(n/2)+1
if (mod(n,i)==0)then
s=s+i
end if
end do
if (s==n)then
print*,"yes"
else
print*,"no"
end if
end program prime number

@komasaru
Copy link
Author

Thanks for your comment.
I'm sorry.
I'm not sure what you want to do with this code.
(Is this the code to determine if it is a prime number?)

@taymazalizadeh
Copy link

yes

@komasaru
Copy link
Author

If you ask me what you should do,
I'm very sorry, but I can only say that you should refer to this 'prime_numbers.f95' code.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment