Toru3 @tetu@mathtod.online

こんな素数判定あったのか
https://en.wikipedia.org/wiki/Baillie%E2%80%93PSW_primality_test

ファンの掃除したら収まったので熱が原因っぽい

uop cache tag parity errorで死んだっぽい。調べてみたら同じ現象の人がいるみたい。

https://community.amd.com/thread/230743

係数それぞれに係数体の情報をもたせればきれいに実装できるが効率が悪い。
グローバル変数で係数体の情報を管理すると、複数の体で並列計算するときに困る。
クラスに係数体の情報を持たせるとクラスが巨大化する上演算子オーバーロードが出来なくてつらい。
そういう訳で多項式に一つ持たせる実装にすることが多い。
しかし、単純な多項式($0$, $1$, $x$など)はたまに型から生成したくなるので困る。(num::Zeroやnum::One実装するのに困った)

式テンプレートとautoの相性が悪い
右辺値入れて死んだ
Rustで書きたい

私は、数学で出てくるときはベクトル、プログラムで出てくるときはベクターと呼んでることが多い気がする。自分でも理由はよくわからないが。

面白そうな本を見つけたので買った
https://mathtod.online/media/Q8TRRK0DdBRz_x8f0YM

見つけた
https://ja.m.wikisource.org/wiki/%E5%88%9D%E7%AD%89%E6%95%B4%E6%95%B0%E8%AB%96%E8%AC%9B%E7%BE%A9/%E7%AC%AC1%E7%AB%A0/%E9%99%84%E8%A8%98_%E5%BE%AA%E7%92%B0%E5%B0%8F%E6%95%B0

循環小数の話、初等整数論講義に書いてあった気がする（うろおぼえ）

RAMに乗り切らない計算どうやっても辛そう。OS任せでスワップしまくると遅そうだし。RAM大量に積んだマシン複数台でMPIして処理するのが良さそうだけど...。

Typed Arenaを使うとうまく行くようです。>Rustでグラフ

https://qnighy.hatenablog.com/entry/2017/04/28/070000

素因数分解の一意性を用いて2が偶数乗と奇数乗で矛盾という証明法もあります。

Rustは(unsafeを使わない限り)スレッドセーフなのでグラフを書き換えながら処理するのは(私が知っている範囲では)unsafeを使わないと難しそうですね。

https://mathtod.online/@cmplstofB/2532004

個人的には latexmk をオススメします

$\require{mhchem}$$\ce{M_p(CO)_q}$
返信を押すと$\rm\TeX$のコマンドが読めるよ

出来た！

using Quaternions
using LinearAlgebra
function quat_lu(N)
A=fill(quat(BigInt(0)//1), N, N)
for i=1:N, j=1:N
A[i,j]=quat([rand(1:10)//rand(1:10) for i=1:4]...)
end
println(A)
println(lu(A))
end
quat_lu(4)

#Julia言語

今日はJuliaで四元数の計算をしようと思ったが、やり方が分からなかったのでとりあえず行列表現で計算した。juliaで四元数を扱う方法を知っている方がいたら教えて欲しい。

最大公約数を求めるアルゴリズムとしてはユークリッドの互除法が有名ですが、（二進の）計算機ではbinary GCDの方が平均的には速いはずです。

https://en.m.wikipedia.org/wiki/Binary_GCD_algorithm