C言語なら俺に聞け 146

2018/04/30(月) 04:47:37.50

C言語の話題のみ取り扱います C++の話題はC++スレへ
質問には最低限の情報(ソース/コンパイラ/OS)を付ける
数行で収まらないソースは以下を適当に使ってURLを晒す
https://paiza.io/
https://ideone.com/
http://codepad.org/

C11
http://www.open-std.org/jtc1/sc22/WG14/www/docs/n1570.pdf

C99
http://www.open-std.org/jtc1/sc22/WG14/www/docs/n1256.pdf
http://kikakurui.com/x3/X3010-2003-01.html

C FAQ 日本語訳
http://www.kouno.jp/home/c_faq/

JPCERT C コーディングスタンダード
https://www.jpcert.or.jp/sc-rules/

前
C言語なら俺に聞け 144
https://mevius.5ch.net/test/read.cgi/tech/1514025223/

次スレを立てる時は本文の1行目に以下を追加して下さい
!extend:on:vvvvv:1000:512
-
※前スレ
C言語なら俺に聞け 145
http://mevius.5ch.net/test/read.cgi/tech/1519046038/

VIPQ2_EXTDAT: default:vvvvv:1000:512:----: EXT was configured

2018/08/07(火) 12:50:36.34

誤差www
そりゃアホが作れば誤差が問題になるだろうねえwww

2018/08/07(火) 12:51:27.81

アホは素直に加算を繰り返せば良いよ

2018/08/07(火) 17:18:17.21

あれ？今日は半角にしないんだ。

2018/08/07(火) 19:11:24.43

おれは半角野郎じゃない
半角野郎はアホだから無理

◆QZaw55cn4c (ﾜｯﾁｮｲ 7f60-0Uuo) · 2018/08/07(火) 19:21:21.20

>>802
√5 の係数を整数で持てばいいのでは？
F_n = （φ^n - (-φ)^(-n)) / √5, φ=(1 + √5)/2

2018/08/07(火) 20:18:00.72

富士山麓オウム啼く・・・あれ？怖いなこれ

2018/08/07(火) 22:59:36.12

https://ideone.com/jzecIF

やはり
整数を足すにかぎるわ

2018/08/07(火) 23:08:45.73

IEEEの倍精度（64bit）で15～16桁程度の精度
つまりdouble、>>809の精度が概ねその程度

64bit整数なら19桁はいける

2018/08/07(火) 23:43:22.81

相変わらずアホだね
視野が狭い

2018/08/07(火) 23:47:20.78

ルートとｎ乗求める処理オーダーってどんなんじゃろ

６４bitまでならルート使ったほうが早いだろうけど
結局桁数が大きくなったらルートとってｎ乗してって結構処理重そう

じつは全部足し算したほうが早かったりして

2018/08/08(水) 05:09:15.02

重いのはべき乗のほうだね

2018/08/08(水) 07:15:16.24

nが大きい場合の計算オーダー

足し算の繰り返しの場合 : n^2
一般項を普通に計算した場合 : n (log n)^2

2018/08/08(水) 09:01:38.18

足し算の繰り返し：O(n)
一般項：O(1)
じゃないの？

2018/08/08(水) 12:44:41.36

加算、乗算、expなどの計算がO(1)で出来る範囲ならそうだね

2018/08/08(水) 12:50:59.40

>>816
まずはO(1)の定義を調べようか…

2018/08/08(水) 12:52:04.40

お前より詳しいと思うよ

2018/08/08(水) 15:07:25.41

お前のほうが可愛いよ

2018/08/08(水) 16:09:12.75

ありがと

2018/08/08(水) 18:58:28.93

フィボナッチの私の回答

nが比較的小さな数までである場合、テーブルを持つのが最速であることは間違いない

じゃあテーブルを持たない場合、

固定有効精度であれば、以下で求めるのが計算量のオーダーが1
(int)(exp(.4812118250596034475*n-.8047189562170501873)+0.5)

floatやdoubleではなく、4倍精度や100倍精度であってもオーダーは1
上の2個の定数は事前に計算しておけばいいし、
毎回計算してもnとは無関係
nが非常に小さい数でなければ単純に加算するよりは上記の方が速いでしょう
floatやdoubleであれば1行で済みます

2018/08/08(水) 19:03:59.97

じゃあ多倍長演算を用いて任意のnに対して正確に求める場合を考える

以下は基礎知識
n番目のフィボナッチ数の桁数はnに比例する
乗算の計算オーダーは、桁数をnとすると n log(n)
ルートの計算オーダーは n log(n)
n乗は、オーダーlog(n)回の乗算で求まる

以上より、一般項の式の通りまじめに計算しても
オーダー n (log n)^2 で計算できる

2018/08/08(水) 19:10:00.90

工夫すると、計算オーダーを n log n にすることが出来る

行列の形でフィボナッチ数列の漸化式を記述すると
|0 1|^n |0|
|1 1| . . |1|
でn番目とn+1番目のフィボナッチ数が求まることがわかる

この式を用いれば
計算オーダー n log n で正確に n 番目のフィボナッチ数を求めることができる

◆QZaw55cn4c (ﾜｯﾁｮｲ 7f60-0Uuo) · 2018/08/08(水) 19:13:46.21

>>810,812,814
漸化式 f(n+2)=f(n+1)+f(n) から一般のf(n) を求めるオーダーはΟ(n)
一般式 >>807 から求めるオーダーはΟ(log n)
だと思う
https://mevius.5ch.net/test/read.cgi/tech/1434079972/49
c を定数として、c^n を求めるにはΟ(log n) でいける

2018/08/08(水) 19:30:35.51

元は>>775の初心者の質問で
計算オーダーとか全く関係なかったんだけど
半角君があまりに無知な癖に自信過剰なので
つい書いてしまった

2018/08/08(水) 19:34:04.93

>>824
それは固定有効精度の場合

固定有効精度であればオーダーは1で計算出来る
>>821参照

2018/08/08(水) 19:36:02.18

固定精度であれば
四則演算もexpも丸めも全て
値によらず固定時間で計算出来る

2018/08/08(水) 19:41:56.80

てことで、
半角君に３勝ですね

コーディング
アルゴリズム
数学

全て勝ってしまった

まあ学歴も勝ってると思うので４勝か

2018/08/08(水) 19:42:03.28

>>824
桁数分のnかからんの？

◆QZaw55cn4c (ﾜｯﾁｮｲ 7f60-0Uuo) · 2018/08/08(水) 20:30:18.66

>>829
>>826 の指摘のとおり、固定有効精度だったら >>829 になるようです
>>821 は、ちょっと理解が及びません…

2018/08/08(水) 21:38:16.32

頭悪いバカが半角に勝ったつもりでいる
半角はこのスレでは間違いなく天才
他のやつらがお話にならないぐらい頭悪いからな
典型的な頭悪いやつが半角に勝とうと必死になってる

2018/08/08(水) 23:31:41.87

ここまで一般式も行列もコードなし多倍長も

まだはじまってもいない

2018/08/08(水) 23:37:23.86

ISO C99に準拠して、なおかつある程度開発が盛んなOSSを教えてください。
ソースコードが綺麗なことで有名ならなお嬉しいです。

◆QZaw55cn4c (ﾜｯﾁｮｲ 7f60-0Uuo) · 2018/08/08(水) 23:41:36.55

>>832
多桁長は結構手間、というか C/C++ でやるのは実は大変だと思います、キャリーフラグが使えない…

2018/08/09(木) 08:27:30.86

素直にライブラリを使えば良いよ

加減乗算しか使わないから
速度を求めなければ自作するのも大した手間では無いけど

2018/08/09(木) 08:39:55.95

>>833
勉強のためなのかな
もしそうならFreeBSDのコマンドのソースコードが良いと思う
短いものが多いからおすすめ
Linuxのものと比べると読みやすい

FreeBSDソースコードリポジトリ
https://svnweb.freebsd.org/base/release/11.2.0/
"bin"が名前に含まれるフォルダにコマンドのソースがある

2018/08/09(木) 10:23:18.39

GNUのツールチェインもいいんじゃないか
ストールマンのコードとか見れるよ

2018/08/09(木) 10:24:09.98

>>836
ありがとうございます。
> 勉強のため
その通りです。
プログラム言語自体がほぼ初めてなのですが
やはり教科書や仕様書ばかり見て実際を見ないのでは身に付かない
と思いまして、OSSのソースコードであれば
コーディング規約が守られていたり
優秀なアルゴリズムが使われていたりするかなと予測して
OSSのソースコードリーディングをすることにしました。
ところが私が勉強の対象にしているのはISO/IEC 9899:1999（C99）であるのに対し、
ほとんどの有名なOSSはC89/90に準拠し、その古い制約に縛られていてC99の勉強の妨げになりそうでした。
というわけでC99に準拠した有名なOSSを探しているのです。

**デフォルトの名無しさん** (JP 0Heb-5Ep3) · 2018/08/09(木) 10:43:20.34

C89とC99の差を何も見ずにここに書けるなら別だが、
そうでなければ勉強の妨げにはならんから安心しろ

2018/08/09(木) 10:54:41.59

多倍長ならGMPがあるやん

2018/08/09(木) 12:48:19.15

>>838
musl libc

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/09(木) 21:25:05.58

>>832
＞一般式
＞多倍長
https://mevius.5ch.net/test/read.cgi/tech/1434079972/50
遅くて残念

2018/08/09(木) 22:23:11.20

もうつくったんかはええｗ

順次加算とスピード同じくらい？

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/09(木) 22:27:53.84

>>843
よくわかりません
多桁長計算に手間取りすぎていて、順次加算だろうが、一般式だろうが、差が出ないのです、私の多桁長計算は亀のように遅いのです…

2018/08/09(木) 23:30:35.12

コードを見てないけど

多倍長は乗算がキモ
オーダーn log nで計算する方法を考えよう
フーリエ変換を使うのが普通

2018/08/09(木) 23:45:29.15

最終的に10進数で表記するのが目的であれば
10進数のまま演算した方が良い
2進10進変換の方がフィボナッチより重い

>>823の方法を使う場合
多倍長ライブラリの乗算を使うより
フーリエ変換のライブラリを使う方が速い
それは式の関係上多倍長乗算を使うと
同じ事を複数回行うために無駄が多いから

同じく>>823の方法を使う場合
桁数の多い最後の方の計算ほど時間がかかる
最後の方の桁数の多い乗算を減らすのが高速化のにつながる

2018/08/10(金) 01:21:54.58

>>841
ありがとうございます。
そもそもmuslを初めて聞きました（恥ずかしながら）
ちょうど，Linuxシステムも勉強したいと思っていたので，これを読もうかなと思います。

2018/08/10(金) 21:12:30.49

とりあえずこの前ｲﾝｽﾄｰﾙした仮想環境のlinuxにgmplibをｲﾝｽﾄｰﾙした
相変わらず低学歴知恵遅れどもは頭悪いﾃｷﾄｰなことばっかりいってるわ

2018/08/10(金) 23:14:02.99

俺なんかGoogleのAIライブラリのTensolFlowインストールしたぞ！

インストールしたぞ

2018/08/10(金) 23:55:30.59

>>823の方針でフィボナッチ数を計算するコードを書いてみました
F(1億) が10.5秒
F(10億) が137.1秒
Cで450行くらいのコード
計算時間より秀丸エディタで結果を開く方が時間がかかってました

F(10) : time = 0.000013 / err = 0.000000
F(100) : time = 0.000017 / err = 0.000000
F(1000) : time = 0.000021 / err = 0.000000
F(10000) : time = 0.000111 / err = 0.000000
F(100000) : time = 0.001233 / err = 0.000001
F(1000000) : time = 0.019993 / err = 0.000004
F(10000000) : time = 0.251717 / err = 0.000019
F(100000000) : time = 10.507523 / err = 0.000076
F(1000000000) : time = 137.143945 / err = 0.000305

2018/08/11(土) 00:06:54.68

10000進数多倍長
超単純なFFT
演算は乗算と加算のみ
誤差の感じから100000進数でも大丈夫そうですね

計算式は基本以下を多倍長にしただけ
多少の無駄は除いてますが

----
uint64_t f(uint64_t n){
n++;
uint64_t a = 1;
uint64_t b = 0;
uint64_t t;
for (int i = 0 ; i < 64 ; i++){
t = b * b;
b = 2 * a * b + t;
a = a * a + t;
if (n & 0x8000000000000000){
t = b;
b = a + b;
a = t;
}
n += n;
}
return a;
}

2018/08/11(土) 00:16:27.47

コードをアップしようと思ったけど
ideoneだとうまく動かないみたい

2018/08/11(土) 10:03:32.46

乙
じゃ、多倍長使って一般項の公式で算出したケースを晒します

スレチと知りつつ C++ で　boost/multiprecision バックエンドは gmp
https://wandbox.org/permlink/QqxaauQeHTjJRrYX

# fib(0) .. fib(1000) + fib(10億)
real 0m0.004s
user 0m0.004s
sys 0m0.000s

ま、そりゃ速いですわな……
科学的表記で出力すると速いけど　整数表記しようとするとfib(10億)で2分ぐらい
文字列化でかなり時間を食ってます

2018/08/11(土) 10:40:58.88

正確に全桁計算してその時間？
じゃないよね？

2018/08/11(土) 10:44:26.74

有効精度何桁で計算した場合の時間？

2018/08/11(土) 10:50:42.30

有効数字が少ない時は>>821の式の方が速いよ

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/11(土) 11:00:32.72

>>821　は桁数が非常に大きくなるとどうなるか、よくわからないんですよ、ちょっと試してみますね

2018/08/11(土) 11:02:48.70

2個の定数は求めたい有効精度+αの精度で求めておけばOK

2018/08/11(土) 11:04:23.75

>>853の式の、
nに依存せずに計算できる所だけ先に計算しただけなので

2018/08/11(土) 11:11:54.17

expの計算の内部ではlog 2を使うし
2進10進変換にはlog_10 (2)を使うので
doubleより高精度で求めるなら
この辺も先に計算しておくことになるでしょう

>>850ははじめから10進での計算なので
10進文字列に変換するのは一瞬

2018/08/11(土) 11:13:07.14

書き忘れましたが
>>850は10進数で全桁正確に求めた場合の時間です

2018/08/11(土) 11:19:04.68

>>855
>>853は仮数部のビット数を任意長に指定してあるから無限精度ですね
出力で桁落ちしてますが　正確に全桁計算してるはずです

2018/08/11(土) 11:22:04.43

>>862
どこで有効桁数を指定してます？
無理数の数値計算なので指定しないと低精度で計算しそうですが

2018/08/11(土) 11:25:21.63

f(10億)の下位100桁はいくつになります？
私は今出先で私の結果は夜貼ります

2018/08/11(土) 11:29:35.08

>>863
https://www.boost.org/doc/libs/1_68_0/libs/multiprecision/doc/html/boost_multiprecision/tut/floats/gmp_float.html
gmp_floatが　number<gmp_float<0>> の typedef なので
"variable precision by setting the template argument to zero"で無限精度指定になってると理解してます
正確に出力させると桁が膨大で検証しにくくてしょうがないな……　これ

2018/08/11(土) 11:39:50.69

とりあえずかわいそうなぐらい頭悪いﾔﾂしかいないのは分かった

一旦、多倍長演算向けに3つの方法を評価する
ちなみにgmpの関数にﾌｨﾎﾞﾅｯﾁの関数がついてる
きっとこの速度にすら届かないと考えられる（まだ動かしてない）

↓多倍長演算使ってない3つの方法の簡単なｺｰﾄﾞがｺﾚ
https://ideone.com/vhpLPV
※　ｵﾏｹでｱﾎが書いたｺｰﾄﾞ(>>851)も入ってる
※　ｵﾚの適切なありがたい注釈がついてる

1.ひたすら足し算

2.一般項
　多倍長演算をするまえに適切な精度を設定しないといけない
　どれぐらいの精度にすればいいかがまだ未解決　※ とりあえず2回計算することでいけるような気がしないでもない

3.再帰階乗演算
　https://www.ics.uci.edu/~eppstein/161/960109.html
　探した中でｺｲﾂが一番いい感じがする
　> This is a recursive algorithm, so as usual we get a recurrence relation defining time,
　> just by writing down the time spent in a call to matpow (O(1)) plus the time in each recursive call
　> (only one recursive call, with argument n/2). So the recurrence is
　> time(n) = O(1) + time(n / 2)

2018/08/11(土) 11:48:25.12

>>864
手計算だと下位3桁は875
合ってます？

2018/08/11(土) 11:52:02.92

>>865
√5を無限精度で求めるなんて不可能ですよ

2018/08/11(土) 11:55:08.80

>>866
>>851はそこのリンクのAlgorithm 5の無駄を省いたものですよ

2018/08/11(土) 11:58:25.35

桁数は>>821で瞬時に出ますね
>>850もまず桁数を求めて必要なメモリを確保してから計算してます

2018/08/11(土) 12:00:16.70

↓ｱﾙｺﾞﾘｽﾞﾑ5でｵﾚのｻﾝﾌﾟﾙではちゃんと動いてるからな
https://ideone.com/vhpLPV

↓ｵﾏｴのは動かしても適切な結果にならない
https://ideone.com/vhpLPV

ｵﾏｴのはお話にならないぐらいﾍﾟｹ
↓ﾍﾟｹの理由もちゃんとｺﾒﾝﾄに書いてある
https://ideone.com/vhpLPV

2018/08/11(土) 12:04:23.64

劣化コピーしたらそりゃ正しく動かないよ
nは64bitにしないと

2018/08/11(土) 12:06:06.40

>>868
あー、そりゃそうですね…

>>867
968 ずれてますね

2018/08/11(土) 12:06:47.07

ｵﾏｴは基本的にｵﾂﾑに問題がある
病院へいったほうがいい

2018/08/11(土) 12:07:58.68

nが32bitなら0x80000000にしないと

多倍長で動かすためのテストコードだから
桁数が小さい時の計算回数は全く無視
ループが進んでからやっと演算がはじまる

2018/08/11(土) 12:09:49.79

64回しかﾙｰﾌﾟしてない

　Σ n
　n→64

はいくつになる

2018/08/11(土) 12:14:07.21

n <<= 1;
って書いてあげた方がよかったかな
どうせコンパイラが足し算に直すだろうし

2018/08/11(土) 16:36:24.18

1.と3.は簡単にできた
https://ideone.com/cmAQzO

gmplibのﾌｨﾎﾞﾅｯﾁは n=10,000,000 を
0.7秒ぐらいで処理できる
ｵﾚが作ったのは4秒ぐらいかかる

m,10000000,4.008508,3.989280
g,10000000,0.724954,0.717120

全然かなわない
少しやってみて面倒そうだから後回しにしてた一般項やってみるか

2018/08/11(土) 18:32:01.22

>>853
同じライブラリを使って>>851のアルゴリズムで計算すると速いかもしれません
a, b, t だけ多倍長にすればいいです

2018/08/11(土) 18:34:26.69

ループの最後の64回目は
aだけを計算すればいいので
かなり省略出来ます

2018/08/11(土) 19:41:53.58

https://ideone.com/2by35t

精度評価用のﾛｸﾞ取得処理

32bit版GMPの仮数部は初期値が64bitで32bitずつ増えていく模様
ちなみに64bit精度で一瞬で求まるﾌｨﾎﾞﾅｯﾁ数のｻｲｽﾞ（bit）予測値はほぼ誤差なし
※ 9999件まで確認したが誤差は発生してない
※ 保険で、1dword(32bit - 最小単位) 分、余分に足しとけば、この予測値で間違いがおきることはまずないと考えられる

ﾌｨﾎﾞﾅｯﾁ数のｻｲｽﾞ（bit）予測値から必要な精度が予測できる
※ 9999件まで確認したが
※ ﾌｨﾎﾞﾅｯﾁ数のｻｲｽﾞ（bit）の1dword(32bit - 最小単位) 分より大きい誤差は発生してない

↓ｺﾚがその結果
http://fast-uploader.com/file/7089539563485/
http://fast-uploader.com/file/7089539604357/

数学的に証明したわけではないが、きっとｺﾚで間違いなくいける
つまり予測値に64bit足しとけば間違いがおきることはない

2018/08/11(土) 19:44:24.09

で、>>881の結果に基づいて
一般項で処理するｺｰﾄﾞを書いた

　https://ideone.com/QKTrLi
　一般項で処理

やってみたが
一般項で処理なんかするとともかく遅い

6,942,482 bitsの一般項の計算で
お話にならないぐらいものすごい時間がかかる

calculation 6942482bits
f,10000000,35.082393,34.855636
g,10000000,0.722054,0.720584

つまり、結論としてﾌｨﾎﾞﾅｯﾁ数を求めるなら
GMPに用意されてる関数を使うのが一番

再帰階乗演算使う方がはるかにﾏｼ
一般項で求めるのはｳﾝｺ

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/11(土) 19:56:34.45

>>882
gmp の生Cインターフェースを叩くとは、勇気がありますね
私は、C++ インターフェース gmp_class でやることはあっても gmp を C インターフェースで叩くことはありません
あれは、あとから見てなにがなんだかわからなくなってしまう…

a + b√5, a, b ∈Z で表される数は、それ同士を足しても掛けても、その結果はやはり a + b√5 の形になるのだから（そして結果は整数になることがわかっているのだったらなおさら）、
√5 を開かずにその係数を有理数の範囲で計算するようにすれば(>>807)、ずいぶんと違うと思います >>824, >>842

2018/08/11(土) 21:43:12.30

F(1千万)の計算時間
gmp_fib_uiは0.72秒
私のは0.25秒

私が適当に作ったヤツの方が勝っちゃいましたね

コードは以下
https://ideone.com/hklTK1

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/11(土) 21:47:29.90

>>884
＞F(1000000000) = 0
となっていますが、ちゃんと動いているのでしょうか？

2018/08/11(土) 21:50:41.32

ローカル環境では動いてますが、
>>852に書いた通りideoneだと動きません

2018/08/11(土) 21:58:41.55

floatとdoubleの値の範囲を教えてください。

-0.000034 < float < 0.000034
-0.00000000016 < ｌong < 0.0000000016

みたいな感じで教えてください＞＜

2018/08/11(土) 22:22:37.58

>>886
ライブラリがないとかならエラーになるし
そんなことってあるん？
何が原因だ

◆QZaw55cn4c (ﾜｯﾁｮｲ ea60-Qb5F) · 2018/08/11(土) 22:25:26.31

>>886
私も手元で確かめました、確かにこれは爆速ですね…私の古いPCでもf(1000万)までならば一瞬で答えが出ます。

2018/08/11(土) 22:26:14.75

C言語でゲームがアクションゲームが作りたくてプログラミングはじめた
実際にプログラミングって楽しいと感じてる
でも入門書通りソートアルゴリズムを勉強してて、基本交換法やら基本選択法は楽にフローチャートもコードも書けたのに、そのあとにやったクイックソートとか単純挿入法では面白いぐらい苦戦した（どっちももうできるようになったけど）
これって、俺のアルゴリズムを考える能力が乏しいってことなのかな
それともはじめはみんなこんなもの？

2018/08/11(土) 22:26:40.35

×ゲームがアクションが
◯アクションゲームが

2018/08/11(土) 22:35:15.61

再帰とポインタで躓くと昔から伝わっている
特に再帰は禅でもやって悟りを得ないと理解できない

2018/08/11(土) 22:35:57.56

だれでも最初は赤ちゃん歩き。私も最初はじゃんけんのプログラミングが全くできなかった。

2018/08/11(土) 22:42:41.92

>>884のバグを修正しました
https://ideone.com/4a3zE8

F(1千万)のideoneでの計算時間は0.31s

2018/08/11(土) 23:13:57.75

n log nより遅いように見えるのは
サイズが大きいとキャッシュから外れるため

フーリエ変換はメモリのランダムアクセスが頻発するので
キャッシュに入るかどうかで大きく時間が変わる

適当に並び変えながら変換することで
キャッシュヒット率を上げたり
メインメモリに収まりきらない巨大サイズをHDDを使いながら変換することも出来るのだが、
今回はシンプルさ重視の為見送り

世の中にある高速FFTライブラリに差し替えれば
もっと高速化する

2018/08/11(土) 23:25:51.63

>>887
環境依存

float, double がそれぞれIEEE754のbinary32とbinary64準拠の場合はおおよそ

-3.4028234663852885981*10^38 ≦ float ≦ 3.4028234663852885981*10^38
-1.7976931348623157081*10^308 ≦ double ≦ 1.7976931348623157081*10^308

正確には以下の範囲
float : ±(2^24-1)*2^104
double : ±(2^53-1)*2^971

これに追加して±∞, NaN などがある

組み込みなどではdoubleも32bitであることも多い

2018/08/11(土) 23:33:36.42

C言語スレだから ^ は排他的論理和の意味で使うべきでしたかね

-340282346638528859811704183484516925440 ≦ binary32 ≦ 340282346638528859811704183484516925440

-1797693134862315708145274237317043567980705675258449965989174768
0315726078002853876058955863276687817154045895351438246423432132
6889464182768467546703537516986049910576551282076245490090389328
9440758685084551339423045832369032229481658085593321233482747978
26204144723168738177180919299881250404026184124858368
≦ binary64 ≦
1797693134862315708145274237317043567980705675258449965989174768
0315726078002853876058955863276687817154045895351438246423432132
6889464182768467546703537516986049910576551282076245490090389328
9440758685084551339423045832369032229481658085593321233482747978
26204144723168738177180919299881250404026184124858368

2018/08/11(土) 23:41:42.36

>>888
単純なバグでした

size_t が64bit前提の部分があったのと、
非常に恥ずかしいですが以下のような動作不定なコードがありました
a->data[i++] = b->data[i] + carry;
私のローカル環境だとたまたま動いていたようです

intが16bitだとまだヤバそうですが、
16bit環境で動かすことはないと思うので気にしないことにします

2018/08/11(土) 23:47:51.76

>>896
ありがとうございました

例えば、986.0033っていう数字は
floatではないのでしょうか？

2018/08/11(土) 23:56:09.61

>>899
986.0033 はdoubleのリテラル
986.0033f はfloatのリテラル

2018/08/12(日) 00:23:34.17

>>898
ideoneだと100万ちょっと下でRuntimeErrorで途中でとまってるぽいね

2018/08/12(日) 07:44:42.21

蛇とJSとアセンブラ少しかじったんだけどCの効率的な勉強法ある？
主に蛇のチューニングと、これから勉強する予定なんだけど組み込み用マイコンのプログラムを書くために使う

小型のドローンもどきにのせたマイコンに、学習済みパターンをのせて動作制御したい