5000万までの素数を求めるプログラムなのですが、私の作ったプログラムは実行時間26秒くらいかかります。
先生が言うには10秒台が出るとのことですが、私は頑張っても時間を短くすることができません。
下に私の作ったプログラムを載せますので短くする方法を教えて下さい。
integer table(2:50000000),pno(50000000),cnt,m,i,j
m=sqrt(50000000.)
do 10 i=2,m
do 10 j=i*2,50000000,i
table(j)=1
10continue
do 20 i=2,50000000
if(table(i).eq.0)then
cnt=cnt+1
end if
20continue
write(6,610)cnt
610format('sosu no goukei =',i8)
end
No.1ベストアンサー
- 回答日時:
ちょっといじってみました。
integer table(2:50000000) / 499999999*0 /
integer cnt / 49999999 /
integer m
integer i, j
m=int(sqrt(50000000.))+1
do 20 i=2,m
if(i .gt. 2 .and. mod(i,2) .eq. 0 ) go to 20
if(i .gt. 3 .and. mod(i,3) .eq. 0 ) go to 20
if(i .gt. 5 .and. mod(i,5) .eq. 0 ) go to 20
if(i .gt. 7 .and. mod(i,5) .eq. 0 ) go to 20
do 10 j=i*2,50000000,i
if(table(j) .eq. 0) then
table(j)=1
cnt=cnt-1
endif
10 continue
20 continue
m=2を判定すれば偶数はすべてチェック完了でしょ。4や6のループは不要で、これだけで計算量は半減です。ついでに3や5、7位まではやっておきましょう。
then~endif構文を使うと煩雑なのでgotoにしてあります。
と思っていたら、もっとスマートなのが
do 20 i=2,m
(table(i) .eq. 0) then
do 10 j=i*2,50000000,i
if(table(j) .eq. 0) then
table(j)=1
cnt=cnt-1
endif
10 continue
endif
20 continue
こちらはどうしてうまくいくか考えてみてね
あと、まず変数はdata文や宣言文で初期化しましょう。
FORTRANの言語仕様では宣言だけした変数の値は不定で、処理系に依存します。実際は0クリアされていることが多いですが、こういうのに頼ってはいけません。
4や6のループは不要というのでわかりました。
私はk=i-(i/2)*2を使って、2以外の偶数を省くというやり方でやりました。
そしたら13秒とかなり早くなりました。
ありがとうございました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Visual Basic(VBA) ファイル全てを .xlsm に変更したところ、プログラムが途中で落ちてしまっています 17 2022/12/07 12:03
- Visual Basic(VBA) Dateserialで データ抽出 2 2022/06/26 21:07
- Visual Basic(VBA) VBA初心者です。 VBAで行単位で条件付き書式の色をカウントしたいです。 大量のデータがあるExc 3 2022/06/08 10:02
- Visual Basic(VBA) VBA初心者です。 VBAで行単位で条件付き書式の色をカウントしたいです。 大量のデータがあるExc 3 2022/06/08 10:00
- Visual Basic(VBA) 配列の勉強をしています。使用する変数の意味、検索条件の書き方が難しいです。 2 2022/09/15 14:06
- Visual Basic(VBA) vba GetAsyncKeyState関数について 1 2023/08/24 12:08
- Excel(エクセル) VBAで組み合わせ算出やCOUNTIFSの処理を高速化したいです。 4 2022/04/07 02:38
- Visual Basic(VBA) 日付の計算方法 2 2022/07/01 23:35
- Visual Basic(VBA) VBAが止まります。 2 2022/09/02 14:02
- Visual Basic(VBA) vbaで色付けされたセルを除外したいのですが 5 2023/01/28 17:22
関連するカテゴリからQ&Aを探す
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
フローチャートの菱形が狭い。。。
-
for文のフローチャート
-
初心者のフローチャート
-
フローチャート
-
フローチャート以外の設計方法
-
フローチャート
-
fortran errorについて
-
フローチャートの演算記号
-
二分法(FORTRAN)
-
0除算して、落ちるプログラムと...
-
正しい五十音順について
-
65536は2の何乗なのでしょうか?
-
あるプログラムのコマンドライ...
-
C++ で、「)」が必要 というエ...
-
VBAで仕様書は書きますか?
-
Excelで4096点以上のFFTの方法
-
VBAにてメール作成した際、一部...
-
変数の値が勝手に変化する原因
-
Excelに埋め込んだVBAのプログ...
-
Vba 実数および実数タイプの変...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
フローチャートで 変数に代入す...
-
フローチャートの菱形が狭い。。。
-
fortran errorについて
-
フローチャートで。
-
二分法(FORTRAN)
-
連立1次方程式を計算するFortra...
-
フローチャート以外の設計方法
-
for文のフローチャート
-
TeXでフローチャート
-
フローチャートのループ
-
Fortranの素数のプログラム
-
学校でフローチャートって教わ...
-
フローチャート
-
fortran
-
フローチャート(本当に困って...
-
フローチャート
-
Rの質問です。 x<-rnorm(n=100,...
-
配列 x に入っているデータの最...
-
フローチャートこれのフローチ...
-
fortran go to 文
おすすめ情報