円周率を計算するプログラムを作りました。
lang
1#include <stdio.h> 2 3int main(void) 4{ 5 double pi,i,a=1; 6 while(i!=1000000){ 7 pi=pi+(1/a); 8 a=a+2; 9 pi=pi-(1/a); 10 a=a+2; 11 i++; 12 } 13 printf("%f\n",pi*4); 14 return 0; 15} 16
で、実行結果が、
3.141592
になります。
これより小数点の桁を下げることってできないでしょうか?
気になる質問をクリップする
クリップした質問は、後からいつでもMYページで確認できます。
またクリップした質問に回答があった際、通知やメールを受け取ることができます。
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
回答5件
0
ベストアンサー
pi
とi
が初期化されていない等の問題はありますが、そこをきちんと対応すれば、1000000
の数を増やしてprintf
で%.15f
などと表示桁数を増やせば、理論上doubleの精度の限界(通常は約15桁)までは計算して表示できます。ただし、桁が上がる度に処理時間が10倍増えますので、現実的ではありません。
※ 1000000
の数を増やさないと8桁ぐらいまでの精度しかありませんので、表示だけ増やしても無意味です。
C
1#include <stdio.h> 2 3int main(void) 4{ 5 double pi = 0; 6 for (long long i = 0; i < 10000000000; i++) { 7 pi += 1.0 / (4 * i + 1); 8 pi -= 1.0 / (4 * i + 3); 9 } 10 printf("%.32f\n", pi * 4); 11 return 0; 12}
出力結果(上が結果、下が本当の円周率)
3.14159265353839645129596647166181 3.14159265358979323846264338327950288...
上は10桁目までは出せますが、私の手元では86.15秒もかかります。11桁目を出すには約861秒(約15分)、12桁目は8615秒(約2時間半)かかるでしょう。現実的な時間で桁数を増やすにはもっと収束が速いアルゴリズムを用いる必要があります。
投稿2016/07/31 07:27
編集2016/07/31 07:44総合スコア21735
0
math.hをインクルードして以下の式で計算すればいいと思います。
atan2( (double)1, (double)1) * 4
表示の際には
printf("%lf\n", double型変数名)
で表示しましょう。
投稿2016/07/31 07:20
編集2016/07/31 07:21総合スコア4830
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2016/07/31 07:47
0
これ本当に動きます?
###追記
これで 9 桁まで出ます。
pi を初期化し、i を long にして while を for に変え、printf の書式を変えました。
精度はもう少し高くできますが、いずれにしろ double の限界が来るので、精度が重要であればアルゴリズムを根本的に見直す必要があると思います。
π を求めるアルゴリズムはググればいくらでも出てくると思いますので、探してみてください。
C
1#include <stdio.h> 2#include <limits.h> 3 4int main(void) 5{ 6 double pi=0, a=1; 7 for (long i = 0; i < LONG_MAX; i++) { 8 pi += 1/a; 9 a += 2; 10 pi -= 1/a; 11 a += 2; 12 } 13 printf("%.9lf\n",pi*4); 14 return 0; 15}
###訂正
C
1printf("%.9lf\n",pi*4);
上記でなく下記が正しいそうです。
C
1printf("%.9f\n",pi*4);
投稿2016/07/31 06:39
編集2016/07/31 07:29総合スコア28660
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2016/07/31 06:52
2016/07/31 06:57
2016/07/31 06:59
2016/07/31 07:21
2016/07/31 07:23
2016/07/31 07:27
0
printf のフォーマットを変更してはいかが?
投稿2016/07/31 06:46
総合スコア610
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
あなたの回答
tips
太字
斜体
打ち消し線
見出し
引用テキストの挿入
コードの挿入
リンクの挿入
リストの挿入
番号リストの挿入
表の挿入
水平線の挿入
プレビュー
質問の解決につながる回答をしましょう。 サンプルコードなど、より具体的な説明があると質問者の理解の助けになります。 また、読む側のことを考えた、分かりやすい文章を心がけましょう。
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。