Fortran 素因数分解

Fortranで素因数分解を実行したいのですが、最後まで分解できずに困っています。
たとえば
8 → 2, 2, 2
と出力したいのですが、実際には
8 → 2, 4
になってしまいます…！

問題のコードは以下の通りです。

Fortran
1program Prime_factorization
2
3	implicit none
4
5	integer n, m, i
6	integer remainder
7
8
9	do
10
11		write(*, *) ' 好きな数字を入力してください (if n <= 1 STOP) '
12		read(*, *) n
13		if (n <= 1 ) stop ' Please input the number that is bigger than 1 ! '
14		if (n > 10000) stop ' Is it bigger than 10,000 ! '
15
16		m = n
17
18		do i = 2, int(sqrt(dble(n)))
19
20
21			if (mod(m, i) == 0) then
22				write(*, *) i
23				m = m / i
24				cycle
25			else
26				m = m
27				write(*, *) 'a'
28			endif
29				
30		enddo
31		
32		write (*, *) m
33
34	enddo
35
36end program Prime_factorization

出力が以下の通りになります。

output
1  好きな数字を入力してください (if n <= 1 STOP) 
28
3           2
4           4
5

何卒ご教示よろしくお願いいたします…

行動規範の内容に同意します

回答1件

ベストアンサー

割り切れた場合に、cycleで次のiの値になってしまうのが原因でしょう。
同じiの値のままで、割り切れるかチェックするというように、ループを二重にする必要があります。

投稿2020/04/13 09:56

otn

総合スコア84557

BKS

2020/04/13 13:32

なるほど…気付きませんでした！以下のように改善しましたところ、問題なく出力することができました！ --------------------------------------------------------------------------------------------------- program Prime_factorization implicit none integer n, i, j integer remainder do write(*, *) ' 好きな数字を入力してください (if n <= 1 STOP) ' read(*, *) n if (n <= 1 ) stop ' Please input the number that is bigger than 1 ! ' if (n > 10000) stop ' Is it bigger than 10,000 ! ' i = 2 do if (i <= int(sqrt(dble(n))) ) then if (mod(n, i) == 0) then n = n / i write(*, *) i cycle else i = i + 1 cycle endif else exit endif enddo write(*, *) n enddo end program Prime_factorization --------------------------------------------------------------------------------------------------- ありがとうございました。

行動規範の内容に同意します