C++相談室 part160

**デフォルトの名無しさん** · 2022/04/15(金) 08:39:32.42

前スレ
C++相談室 part158
https://mevius.5ch.net/test/read.cgi/tech/1645239402/

**デフォルトの名無しさん** · 2022/05/04(水) 11:15:08.92

>>599
10進にして全桁表示可能だろうかと思って
10進での答えが知りたいわけだ。

**デフォルトの名無しさん** · 2022/05/04(水) 11:16:25.61

>>600
可能だと思うよ
自分は詳しくないけどそういう大きい桁に関するライブラリの書き込みをここでよく見るので

**デフォルトの名無しさん** · 2022/05/04(水) 11:19:13.87

でもその前にどのくらいのメモリが必要そうかの見積もりはざっくりやったほうが良さそうね

**デフォルトの名無しさん** · 2022/05/04(水) 11:20:53.78

ざっくり2GBなのか
なんか現実的ではない気がしてきた
少なくとも家庭用のPCなんかでは無理そうかも

**デフォルトの名無しさん** · 2022/05/04(水) 11:48:39.54

>>603
ありがとう。
まだ時間がかかりそうだな。
ちなみに128倍精度浮動小数点数の最大値の計算の途中で出てくる計算。

**デフォルトの名無しさん** · 2022/05/04(水) 14:50:41.22

128倍精度ってなんですか（・∀・）？
倍精度で表せる最大値（(2^(1023+1)) - epsilon）の2^128倍は
(2^(130944+1)) - epsilon
であって
2^17179869183
よりｶﾅｰﾘ小さげ

**デフォルトの名無しさん** · 2022/05/04(水) 15:26:18.04

ｺﾞﾒｽ落ち着いて考えたら倍精度で表せる最大値は
　((2^(1023+1)) - epsilon*(2^1023)　=　(2^(1023+1)) - (2^970)
で、その128倍の精度だとしたら
　((2^(130944+1)) - (128倍精度のepsilon)*2^130944) = ((2^(130944+1)) - (2^122771)
ここで
　128倍精度のepsilon = 2^(-8173)
と仮定ェ、
（∵全体が倍精度の128倍=1024 byte = 8192 bitで指数部が8 bit増加で19 bit、符号ビット1 bit
　やったら仮数部は8192-19-1 = 8172 bitで、最上位の1を省略する
　ケチ表現なので実質2^(-8173)が仮数部で表せる最小値すわなち128倍精度のepsilon

**デフォルトの名無しさん** · 2022/05/04(水) 15:36:55.26

倍精度浮動小数点数が64ビット浮動小数点数なので、
128倍精度浮動小数点数は4096ビット浮動小数点数。

IEEE754-2008の推奨による仕様での最大値の計算をしていたんです。

**デフォルトの名無しさん** · 2022/05/04(水) 15:38:07.79

k bit浮動小数点数
https://irem.univ-reunion.fr/IMG/pdf/ieee-754-2008.pdf#page=25

**デフォルトの名無しさん** · 2022/05/04(水) 15:39:34.72

ちなみに128倍精度の最大値は
(2^((2^(35-1))-1))*(2-(2^-4060))
で計算可能かと。

**デフォルトの名無しさん** · 2022/05/04(水) 15:59:04.91

倍精度浮動小数点数の場合。

>(2^((2^(11-1))-1))*(2-(2^-52))
=
+ 1797 69313
48623 15708 14527 42373 17043 56798 07056 75258 44996 59891
74768 03157 26078 00285 38760 58955 86327 66878 17154 04589
53514 38246 42343 21326 88946 41827 68467 54670 35375 16986
04991 05765 51282 07624 54900 90389 32894 40758 68508 45513
39423 04583 23690 32229 48165 80855 93321 23348 27479 78262
04144 72316 87381 77180 91929 98812 50404 02618 41248 58368.

309桁の値。

となるので、計算式は合っているはずだが根拠を失念してしまったｗ
ネットのどこかに書いてあったはず。その式を使っただけ。
128倍精度の場合、35が指数部(符号含む)で4060が仮数部(符号含まない)

**デフォルトの名無しさん** · 2022/05/04(水) 17:03:38.52

>>597
2GiBのメモリをゼロクリアしてからMSBを1にするだけ
簡単すぎて話にもならん

**デフォルトの名無しさん** · 2022/05/04(水) 17:27:38.78

10進数でデータが欲しいって書いてあるのが見えないかｗ

**デフォルトの名無しさん** · 2022/05/04(水) 17:31:26.47

C++でやることにこだわりないならPythonが手っ取り早いぞ
あれのintはデフォルトで多倍長だから

>>> f=open(r"C:\pow2.txt", "w")
>>> print(2**17179869183, file=f)

**デフォルトの名無しさん** · 2022/05/04(水) 17:31:51.73

10進数でデータの意味が分からんが一桁づつ印字すれば良いだけだろ

**デフォルトの名無しさん** · 2022/05/04(水) 17:32:35.63

>>613
ありがとう。
参考になります。

**デフォルトの名無しさん** · 2022/05/04(水) 17:57:22.19

スレチだけど
https://paiza.io/projects/
のpython3で
print(2**17179869183)
をやってみたところエラーになった。
print(2**171798)
では答えが出る。
メモリかな？

**デフォルトの名無しさん** · 2022/05/04(水) 19:39:52.57

http://codepad.org/ でも print(2**171798) までやね
print(2**1717986) だと time out
print(2**171798691) だと MemoryError

**デフォルトの名無しさん** · 2022/05/05(木) 00:09:45.17

下記は一部ですので不完全です、型をソートしています。
index_sequenceで展開が面倒なんですよね、何かうまい方法ありませんか

template
<
　template<typename> class Orderable,
　typename　　　　 ListT
>
class SortTypes;

template
<
　template<typename>　class Orderable,
　template<typename...> class List,
　typename...　　　　 Ts
>
class SortTypes<Orderable,List<Ts...>>
{
　template<typename SeqT> struct Impl;

　template<std::size_t... Is>
　struct Impl<std::index_sequence<Is...>>
　{
　　static constexpr std::array<std::size_t,sizeof...(Ts)> indexes =
　　　[]{
　　　　　using Ord = std::common_type_t<decltype(Orderable<Ts>::value)...>;
　　　　　struct Pair {Ord;std::size_t ndx;};
　　　　　Pair arr[] = {{Orderable<Ts>::value,Is}...};
　　　　　insert_sort(arr,0,sizeof...(Ts),[](auto x,auto y){
　　　　　　return x.ord < y.ord;
　　　　　});
　　　　　return std::array<std::size_t,sizeof...(Ts)>{ arr[Is].ndx...};
　　　 }();
　　
　　using type = List<std::tuple_element<indexes[Is],TypelList<Ts...>>...>;
　　　
　};
};

**デフォルトの名無しさん** · 2022/05/05(木) 02:32:03.93

フリーランスに立ちはだかる「常駐」の壁。慣例を打ち壊し、
“テレワーク”案件3割→8割へと成長を遂げた「クラウドテック」の軌跡

リモートワーク求人専門サイト「プロリモート」がリニューアルオープン、
業務委託契約の求職者と企業をマッチング

1/3以上が採用につながる高マッチング率、リモートワーク×エンジニア・デザイナー専門の
人材紹介サービス「ReworkerAgent」正式リリース場所からも時間からも自由な働き方を実現！

『ReWorks（リワークス）』リモートワーク特化型転職サイトとして 3月5日リニューアル

副業・兼業マッチングサービス「クラウドリンクス」登録者数2万人突破
中小企業で進む副業人材の採用、96％が継続採用を希望

茨城県日立市、県外からの「テレワーク移住者」に最大151万円の助成金

長野市、市内に移転・事業所設置し、移住することで最大550万円の支援金を支給

フリーランスが活用できる「最大1,000～3,000万円・補助率50%～75%」の
『ものづくり・商業・サービス補助金』とは？概要や条件を解説

**デフォルトの名無しさん** · 2022/05/05(木) 03:10:15.14

>>597
やったったが愚直に10で割って下の桁から求めるだと手元のPCでは200年ぐらいかかりそう……
https://ideone.com/6BG7hf

つなみにこれは2^0、2^8、2^64、2^1000には一瞬で正解する
[302] 100.000000 % (38/38) --- 0 sec
2^1000 = 1.071508607186267320948425049060001810561404811705533607443750388 x 10^301

**デフォルトの名無しさん** · 2022/05/05(木) 03:10:34.57

10より大きい10^n（100万とか）で割ることにしたら多少は改善するかもしれないが、
最終的に10^nのn桁を10進数表示することを繰り返すことを考えると10で割るトータルの回数に違いが出ないから改善しなさげ……orz

むしろ2進数で計算して10進数にするのではなく最初から10^n進数表現（ただし10^n < UINT32_MAXとか）で2^Nの計算の方を実質10進数でやった方が速そう
計算の手間がO(N^2)なのは変わらないが、割り算無しで済むのは大きいい（多分、

**デフォルトの名無しさん** · 2022/05/05(木) 03:19:38.09

2^Nを10^n進数表現（ただし10^n < UINT32_MAXとか）で計算するのはうまくやれば（10^nをUINT_MAXとかよりやや小さめにしたら）
掛け算によるキャリーの伝搬頻度を減らせるから副次的にさらなる高速化の余地が生まれ得る、

**デフォルトの名無しさん** · 2022/05/05(木) 06:04:52.05

繰り返し二乗法で大きい素数のmodを計算すれば良いだけ
整数の商がいくらかもすぐ計算できるので、簡単に復元できる
復元した値をメモリに持つかどうかは別の話
はい終わり

これしきのことも分からない奴ら頭ヤバいよ
言語仕様の話ばっかり大好きでアルゴリズムの知識が微塵もないボンクラ
競プロキッズのこと馬鹿にできないね

**デフォルトの名無しさん** · 2022/05/05(木) 06:34:28.46

素直に10進多倍長で計算すればいいだけ
今のPCならすぐ

**デフォルトの名無しさん** · 2022/05/05(木) 06:58:12.94

>>623-624
すごいね
じゃあ 2^17179869183 の解をここに書いてみて
あとそのプログラムコードも

**デフォルトの名無しさん** · 2022/05/05(木) 08:43:37.74

自演ちゃうねん
繰り返し二乗法の適用を思いつかなかったのは太くの至りで宦官の至りやが
2^N計算を10^n進数ベースでやるときに適用して高速化するやり方しか思いつかん
>大きい素数のmodを計算すれば良いだけ
とは何ですか（・∀・）？

**デフォルトの名無しさん** · 2022/05/05(木) 10:32:36.82

興味本位で全桁表示したいだけなら5171655947桁あることを知って諦めてほしい
/と%による単純な基数変換を使うにしても、/も%も桁数Nに対してO(NlogN)程度かかるらしく、それだけかかってようやく1桁求まるレベルなんで、まあ無理よ

概算したいだけなら高校生でも知ってる対数を使う方法で4.637182185+e5171655947くらいと分かる

**デフォルトの名無しさん** · 2022/05/05(木) 10:34:08.63

10000007とかで割れって言ってるんだろ
これじゃ余分な7の誤差蓄積しまくって上位めちゃくちゃになるし、影響ないほど巨大な10^n近似素数使っちゃったら結局計算量は大差ない
競プロキッズの考えそうなこと

**デフォルトの名無しさん** · 2022/05/05(木) 10:39:30.08

あごめん桁数は5171655948が正しいね

**デフォルトの名無しさん** · 2022/05/05(木) 10:41:25.83

>>628
7の誤差ってなにけ？
整数なのに誤差とかあるん？

**デフォルトの名無しさん** · 2022/05/05(木) 10:44:36.09

>>627
10進での桁数は
17179869183*log2/log10
https://ja.wolframalpha.com/input?i=17179869183*log2%2Flog10
で計算できることは知っているが、
先頭の値はどうやって出すの？

**デフォルトの名無しさん** · 2022/05/05(木) 11:11:25.53

ていうか対数知らんのだろうね

**デフォルトの名無しさん** · 2022/05/05(木) 11:31:29.78

>>613
自動的に昇進するだけで、デフォルトで多倍長って表現は違和感。
2**-1は0.5だけどデフォルトで浮動小数点ってわけじゃないよね

**デフォルトの名無しさん** · 2022/05/05(木) 11:36:06.02

>>631
小数部だけをまた10の肩に乗っければいい

**デフォルトの名無しさん** · 2022/05/05(木) 11:47:06.17

整数で返すべき関数の実装あるんか？

**デフォルトの名無しさん** · 2022/05/05(木) 12:32:26.31

GMPでやってみたけどoverflowで落ちるな

**デフォルトの名無しさん** · 2022/05/05(木) 12:52:41.09

皆さんすみません>>597です。

2^17179869183
は
とりあえず
10^5171655945.66625415874401993568579902587321942922330151106009...
で勘弁してください。

実数が分かれば文字数で桁数も分かりますので、
実数が知りたかったんですが、ここまでしか分かりません。
お手数おかけしました。
多倍長電卓で計算できないんでプログラムならできるのかな？と思ったわけです。
現時点では計算時間とメモリを食いそうですね。

**デフォルトの名無しさん** · 2022/05/05(木) 17:24:23.79

ざっくりとした大きさや上下数桁はわかっても全体を丸々把握するのは厳しそう

**デフォルトの名無しさん** · 2022/05/05(木) 18:10:40.93

https://ideone.com/O1rs45
誰か助けてー。
ウエブのレクチャー見て書いたんだけど、一件通ったと思ったらほかのテストが通らない。
デコーダが悪いのだろうか？というか、解釈狂ってるかなぁ？？マジ解らん。

アドバイスクレクレ。

**デフォルトの名無しさん** · 2022/05/05(木) 18:12:32.63

まじ　わか　らん

**蟻人間** ◆T6xkBnTXz7B0 · 2022/05/05(木) 18:17:21.57

>>639
std::uniform_int_distribution<int> ui('A', 'z');

大文字小文字？

**デフォルトの名無しさん** · 2022/05/05(木) 19:01:19.69

Aのほうが数字的に小さいのですよ。

**デフォルトの名無しさん** · 2022/05/05(木) 19:06:48.31

>>639
V.push_back({ { D[i] }, i }); は V.push_back({ D[i], i }); じゃないの？

**デフォルトの名無しさん** · 2022/05/05(木) 19:40:11.06

>>639
プログラムの意図がわからんから投げたわ

1. Data型の<1>要素の意味はなんだよ。typedefした行のコメントに書いといてほしい
2. 関数名を流し見る限り、ソート・エンコード・デコードの3つの仕事がある。中途半端に仕事が混ざってない？

**デフォルトの名無しさん** · 2022/05/05(木) 19:50:34.56

>>644
ブロックソートアルゴリズムの実装をしているだろう。
とすると size_t はソートバリエーションのインデックスを表していると思われる。

**デフォルトの名無しさん** · 2022/05/05(木) 21:14:55.46

Wikipediaのブロックソートを実装しようとしているのなら、
あちこち違ってるんで、修正版を貼っておく。
https://ideone.com/q008E8

**デフォルトの名無しさん** · 2022/05/05(木) 22:10:21.31

int64_tとsize_tを暗黙の型変換を期待して混ぜて計算するとはまる？

(int64t - int64t) + sizet >= int64t
の比較をしたら2つ目のint64tが大きくなっても(とは言っても全ての値は10～-10くらい)boolの値はずっと1でした。
比較の前まで切り出してprintfすると正しい値が出るのですが、比較式に組み込むとおかしくなります。
sizetをint64tにキャストすると正しく動きました。

暗黙の型変換はsizetをint64tにしてくれないのでしょうか

環境はWSL2のg++です。

教えて下さい

**はちみつ餃子** ◆8X2XSCHEME · 2022/05/05(木) 22:59:30.02

>>647
色々と面倒な細かいルールはあるんだが、大まかには

・二項演算子は原則としては大きい側の型に揃えてから計算する (符号は符号無しが優先。結果の型もそれ)
・ size_t の大きさは処理系定義 (符号無しであることは保証される)
・ signed な型が unsigned に変換されると変換後の型の最大値より 1 大きい値で mod をとった値になる

というのがこの場合にかかわる重要なルールだと思う。
size_t は処理系定義だからその環境でどうだかわからんが uint64_t 相当だと仮定するぞ。
int64_t - int64_t の型は int64_t で、 int64_t + size_t は uint64_t になる。
つまり int64_t-int64_t の値が負数の場合でも常に正の値に修正される。
しかも値として常に 10～-10 くらいということは、
unsigned な値への型変換規則によってかなり巨大な値になってしまう。
たぶんこれが意図しない挙動の原因だと思う。

**はちみつ餃子** ◆8X2XSCHEME · 2022/05/05(木) 23:03:22.90

>>648
符号無しが優先って書いちゃったけど、これは大きさが同等なら符号無しが優先ってことね。
大きさに差があるなら大きい方の符号が採用される。
(大きさって言ってるけど実際には順位 (rank) の序列が定められている。まあ大きさの順なので実質同じやろ。)

**デフォルトの名無しさん** · 2022/05/05(木) 23:06:38.57

size_t が非負の32 bit なら、int64_t の範囲内だから、OK だけど、
size_t が非負の64bitなら、int64_tの範囲外で無理

int64_tの非負の範囲は、63bitしかないから

**はちみつ餃子** ◆8X2XSCHEME · 2022/05/05(木) 23:20:02.75

例としてこういうのがあるとわかりやすいかな。

#include <iostream>
#include <cstdint>

int main(void){
auto a=std::int64_t(6)-std::int64_t(8)+std::uint64_t(1);
std::cout << a << std::endl;
}

結果として 18446744073709551615 が表示される。

**デフォルトの名無しさん** · 2022/05/05(木) 23:39:48.09

というわけで高速化すた、
https://ideone.com/TWKDT6

動きを見る限り、2^17179869183が計算が終わるまで手元のPCで8時間ぐらいの見込み

（検算）
次のぐらいなら一瞬で出る（小数点の位置が手抜きなので指数の解釈に注意
2^1000 = 0000001 0715086 0718626 7320948 4250490 6000181 0561404 8117055 3360744 3750388. x 10^298
2^171798 = 0000002 2448900 7901275 1700257 1990086 5699838 8910190 8350859 7305786 1643657. x 10^51713

2^171798は>>620のでは4秒かかかっていたやつが100 msぐらいで出る

**647** · 2022/05/05(木) 23:52:14.24

皆様ありがとう。
たいへん納得できました

**デフォルトの名無しさん** · 2022/05/06(金) 07:02:57.27

ｺﾞﾒｽ8時間と言う見積は指数の桁を一つ少なかった2^1717986918で見積もってたorz
やったら最新状況は↓こんな感じでこの後4^9×5815 = 1,524,367,360秒 = 48.304年かかることに、。n_
Calculating 2^17179869183 ...
n=17179869183, mbuf sz=2 --- 0 sec
n=8589934591, mbuf sz=2 --- 0 sec
n=4294967295, mbuf sz=2 --- 0 sec
...（中略）...
n=8191, mbuf sz=180374 --- 24 sec
n=4095, mbuf sz=360748 --- 96 sec x 4
n=2047, mbuf sz=721494 --- 385 sec x 4.10
n=1023, mbuf sz=1442986 --- 1467 sec x 3.81
n=511, mbuf sz=2885970 --- 5815 sec x 3.96

**デフォルトの名無しさん** · 2022/05/06(金) 10:53:36.44

生きてるうちに終わりそうでよかった

**デフォルトの名無しさん** · 2022/05/06(金) 12:17:44.01

python超初心者スレでひな形を教えてもらってpythonでプログラミングしたところ、

2^17179869183
は
4.637182185095045843877145312775985050353215513064979818029199207764811280740258379690028942685133666E+5171655945

になるようです。
ご参考まで。精度等未確認ですが。

**デフォルトの名無しさん** · 2022/05/06(金) 12:40:01.31

>>656
ちなみにこの計算は、
128倍精度浮動小数点数の最大値の計算で出てくる式の掛け算で区切った左側の値です。

128倍精度浮動小数点数 (2^((2^(35-1))-1))*(2-(2^-4060)) 5171655946桁の値有効桁数1222桁

(2^((2^(35-1))-1))*(2-(2^-4060))
は
9.274364370190091687754290625551970100706431026129959636058398415529622561480516759380057885370267332E+5171655945
という値になりました。

精度等確認できる方はお願いします。
仮数部の最終桁の値は丸めになっているようなので+1されている可能性があります。

**デフォルトの名無しさん** · 2022/05/06(金) 13:36:52.82

>>646
修正感謝。自分のコードに転写してみたらちゃんと動いた。
これちゃんとブロックソートか解らないけど、可逆ではあるので圧縮してみるね。
本当にありがとう。

**デフォルトの名無しさん** · 2022/05/06(金) 15:08:33.81

https://ideone.com/xxiusG
詐欺れなーい。LZWの1ワードが2バイトだと、若干損位で収束してしまうな。

何でだー。圧縮理論のはずだろ～～。

**デフォルトの名無しさん** · 2022/05/06(金) 15:20:41.30

延々同じようなことやってる低学歴さんたちって>>623見えんの？

**デフォルトの名無しさん** · 2022/05/06(金) 15:25:15.07

C++ のプログラムでサスペンド/復帰をお手軽に実装するための定番の方法ってないですかね？
USR1 を受け取ったらメモリの内容を丸々ファイルにダンプして、復帰するときはそれを読み込んで続きから実行、くらい簡単な仕組みが実現てきたら嬉しいです

**デフォルトの名無しさん** · 2022/05/06(金) 15:25:29.57

>>659
真の乱数データならどのような圧縮アルゴリズムでも縮まんよ

**デフォルトの名無しさん** · 2022/05/06(金) 15:28:31.87

>>662
噂によるとLZWはエントロピー圧縮では無いので出来ても良いかと思ったんだけど、駄目か。
根拠は辞書捨てるからだけどね。

とりあえず完成。

**デフォルトの名無しさん** · 2022/05/06(金) 16:06:26.56

>>661
WM_POWERBROADCAST ggrks

**デフォルトの名無しさん** · 2022/05/06(金) 17:03:00.97

圧縮方式がなんだろうと可逆圧縮はデータ長が短くなるケースがあるならデータ長が長くなるケースが必ずある

65536通りの表記がある2byteデータで、1通りが1byteに圧縮できるとしよう。(0x0000→0x00)
残りの65535通りから256通りの表現を奪うので(0x00から始まるパターンが作れなくなる)、3byte長で表現しなきゃいけないやつが出てくる。

このへんはトレードオフだが、テキストデータとか一定の傾向があるデータのデータ長を短くできりゃいいだろ？
ランダムデータとか圧縮済みデータを圧縮するなんて考えない

**デフォルトの名無しさん** · 2022/05/06(金) 17:26:03.51

>>665
まぁ、そうですね。
LZWの性質上、一回出たモノを辞書から洗って一個追加して辞書に返す。という事をやるのです。
エントロピー圧縮では無いので乱数を詐欺れたら無敵だったのですけど、現実は厳しい。

**はちみつ餃子** ◆8X2XSCHEME · 2022/05/06(金) 18:18:33.94

データ圧縮時にあり得る最大伸長率をまとめたサイトをどこかで見たような記憶があったんだけど
あらためて探すと見つからない……。幻だったんだろうか。

**デフォルトの名無しさん** · 2022/05/06(金) 18:30:36.95

>>667
ZipBombとかやべーやつもいるから安全に展開できるかの仕様の確認目的かな？

正攻法での最大圧縮率とかだと話が変わってきそうだが

**デフォルトの名無しさん** · 2022/05/06(金) 20:00:41.83

>>660
繰り返し二乗法は分かったのだが
modを取る大きい素数の選び方が分からないので教えてほしい

**デフォルトの名無しさん** · 2022/05/06(金) 20:02:30.43

>>664
Linuxではどうすればいいの？

**デフォルトの名無しさん** · 2022/05/06(金) 20:05:10.12

内部クラスを外部クラスが継承するって不可能？
struct A : A::B { struct B {}; };という特殊な事やりたいんやけど前方宣言とかできないっぽいし無理やろか

**デフォルトの名無しさん** · 2022/05/06(金) 20:10:36.04

>>670
SIGなんちゃら系だろ
知らんけど

**デフォルトの名無しさん** · 2022/05/06(金) 20:36:05.94

>>671
無理だろ

**デフォルトの名無しさん** · 2022/05/06(金) 22:11:28.52

>>623
無理だろ

**デフォルトの名無しさん** · 2022/05/06(金) 23:29:10.82

>>597
30分くらいで計算できた

**デフォルトの名無しさん** · 2022/05/06(金) 23:34:15.82

計算結果いる？
5GBくらいのテキストファイル

**デフォルトの名無しさん** · 2022/05/07(土) 00:00:32.86

ああもう48年後か…
おまいらまだ生きてるか？

**デフォルトの名無しさん** · 2022/05/07(土) 00:32:32.99

>>676
先頭100桁とE+の値(文字数でもいいです)と最後の20桁ください。

**デフォルトの名無しさん** · 2022/05/07(土) 01:43:20.82

1ｹﾞｯﾄズザー！！！１１１！！１１１！c⌒っﾟДﾟ)っ
https://ideone.com/R3szJK

2^17179869183 is between
0000000 0000046 3718218 5095045 8438771 4531277 5985050 3531670 5893980 9439277. x 10^5171655888
0000000 0000046 3718218 5095045 8438771 4531277 5985050 3531670 5894606 3902528. x 10^5171655888

100 msで計算できるようにしたったわ！

>>675
暗算はスレチ

**デフォルトの名無しさん** · 2022/05/07(土) 01:50:30.88

つなみに基数変換アルゴリズムにはトーナメント方式というものがあって並列計算に持ち込めるらし
星の数ほどプロセサーが使える場合は有効かもしれん
富岳を使える機会があったらやってみたい

**デフォルトの名無しさん** · 2022/05/07(土) 02:01:09.56

>>678 >>679
それだけなら>>657の128倍精度のライブラリで計算外すればすぐでは？

**デフォルトの名無しさん** · 2022/05/07(土) 02:05:54.40

なにおうヽ(｀Д´#)ﾉﾑｷｰ!
2^17179869183 is between
0000000 0000046 3718218 5095045 8438771 4531277 5985050 3532155 1306497 9818029
1992077 6481128 0740258 3796900 2894268 5133665 8876785 8487588 5591832 7186419
0576212 7656864 8464720 4285493 3923042 9654230 7916817 0776591 7647233 1128072
2760420 2266531 9507091 8968753 5716084 2962215 4728982 6664323 6555234 1274646
4896521 9429451 2975764 8380822 6363423 2746025 0005054 5131536 2217832 5023004
1605005 9971605 5861808 1342288 9537948 8099475 2493069 5439539 2974540 8604555
3827776 8910650 1779698 0664637 7405672 7472854 5855966 0105172 1933481 1967719
5700701 3333182 0937700 3166605 4871435 3734677 0010217 6905552 3397265 0416277
0094988 0613651 2390001 0317736 9989814 4998699 5188652 7509316 0899781 9619323
5715923 0653762 3370836 8150405 6226421 7362228 9333270 6654384 1686901 0642565. x 10^5171655258
0000000 0000046 3718218 5095045 8438771 4531277 5985050 3532155 1306497 9818029
1992077 6481128 0740258 3796900 2894268 5133665 8876785 8487588 5591832 7186419
0576212 7656864 8464720 4285493 3923042 9654230 7916817 0776591 7647233 1128072
2760420 2266531 9507091 8968753 5716084 2962215 4728982 6664323 6555234 1274646
4896521 9429451 2975764 8380822 6363423 2746025 0005054 5131536 2217832 5023004
1605005 9971605 5861808 1342288 9537948 8099475 2493069 5439539 2974540 8604555
3827776 8910650 1779698 0664637 7405672 7472854 5855966 0105172 1933481 1967719
5700701 3333182 0937700 3166605 4871435 3734677 0010217 6905552 3397265 0416277
0094988 0613651 2390001 0317736 9989814 4998699 5188652 7509316 0899781 9619323
5715923 0653762 3370836 8150405 6226421 7362228 9333270 6654384 1687053 4062589. x 10^5171655258

**デフォルトの名無しさん** · 2022/05/07(土) 02:06:56.87

>>682
筆算法じゃ全桁計算は無理

**デフォルトの名無しさん** · 2022/05/07(土) 02:12:56.32

>>683
48年ぐらい待ったらできうる

**デフォルトの名無しさん** · 2022/05/07(土) 02:17:43.49

2^17179869183 全桁

https://5.gigafile.nu/0514-n79efacc97fda1fc1ecce0eeb21225ab7
(ギガファイル便テキストファイル)

**デフォルトの名無しさん** · 2022/05/07(土) 02:20:59.45

なぜコードでなく計算結果なのか

**デフォルトの名無しさん** · 2022/05/07(土) 03:43:06.05

>>661,664,670,672
すみません。Linuxを想定していました

**デフォルトの名無しさん** · 2022/05/07(土) 04:06:12.29

>>678
素人なので全部は計算出来ないですが
最後の20桁ならすぐに計算できました
88085753133198737408

**デフォルトの名無しさん** · 2022/05/07(土) 04:20:35.96

なんかCRIUとかいうのでできるらしいぞ

**デフォルトの名無しさん** · 2022/05/07(土) 05:41:08.75

乗算は量子演算使わないと万年かかる

**デフォルトの名無しさん** · 2022/05/07(土) 07:36:33.84

>>685
全桁データだからって正確とはかぎらんのかーい！

**デフォルトの名無しさん** · 2022/05/07(土) 07:40:26.77

合ってると思うけど
2通りで計算して同じだったし

**デフォルトの名無しさん** · 2022/05/07(土) 08:15:43.50

先頭数十桁と最後の数十桁は合っていそう。
文字数は5171655817で合ってるの？

**デフォルトの名無しさん** · 2022/05/07(土) 08:21:32.36

>>637
はウルフラムでの17179869183*Log2/log10の計算結果。
>>656
はpythonで打ったプログラムで出たもの。

桁数はE+5171655945で間違いないと思ったが。

**デフォルトの名無しさん** · 2022/05/07(土) 08:49:15.35

2^17179869183

●計算結果
https://5.gigafile.nu/0514-n79efacc97fda1fc1ecce0eeb21225ab7
(5171655946バイト)

※実行時間: 約20分

●コード (C言語)
https://ideone.com/xSaa6x

※メモリが16GBくらい必要です

**デフォルトの名無しさん** · 2022/05/07(土) 09:26:24.51

>>693
文字数チェックに使ったGigaTextViewerにバグがあったらしい。
バイト数で5171655946で見たら合ってるな。
バグないみたい。ご苦労様。

あとは128倍精度浮動小数点数の最大値とか計算できると良いな。
相当凄腕なら64倍以上のn倍精度精度浮動小数点数の最大値を0.06秒くらいで出せると良いけどな。

**デフォルトの名無しさん** · 2022/05/07(土) 10:01:38.30

(2^((2^(35-1))-1))*(2-(2^-4060))

これの計算なら
>>695のコードを少しいじれば出来るから
やってみな

**デフォルトの名無しさん** · 2022/05/07(土) 10:05:53.71

ともかくどうもありがとう。

**デフォルトの名無しさん** · 2022/05/07(土) 10:27:32.54

128倍精度付近の計算は5GB程度と分かったが、256倍精度なら80GB程度とかになるのか？
全部のデータを吐くのには無理があるな。