前回（↓）の続き。 nazodane.hatenadiary.org前回のADDVAL>>SHIFTを2進法に直すと以下のようになる。 >>> [bin((0xaaaaaaaaaaaaaaab * ((1 << (i+1)) - 1) &0xffffffffffffffff) >> (i+1)) for i in range(63)] ['0b101010101010101010101010101010101010101…

前回（↓）の続き。 nazodane.hatenadiary.org前回のコードを適当にベクトライズするとこうなる。 // gcc rcollatz.c -O3 -march=haswell // clang rcollatz.c -O3 -mavx #include <stdint.h> #include <stdio.h> //#include <assert.h> void loop(uint64_t x, int l){ // assert(x&1); uin</assert.h></stdio.h></stdint.h>…

前回（↓）の続き。 nazodane.hatenadiary.orgコラッツ問題の計算で、最後の手前の数字は「3の倍数かつ1を足すと2の累乗数になる数」になるはず（偶数処理を遅延した場合はそれに2の累乗数を掛ける必要あるけど、今回は遅延しないで考えることとする）。1を足…

1億円欲しいからコラッツ予想に手を出してみる。science.srad.jp ja.wikipedia.orgまず偶数処理をまとめたコラッツ問題の計算プログラムを適当に書く。 // gcc collatz.c -march=haswell // ./a.out 13179928405231 // Copyright: Nazodane // License: GPLv…

今は別の界隈に居たりしますが、長年の病気の原因がやっとこさ分かって治る見込みが出来てきたのでここにも書いておこうかなと思いました。 ついでに昔の日記を読み返してみたけど当時のこと思い出して何とも言えない気分になってきてアレ(^^; 文体が安定し…