美しい整数の世界 (780ﾚｽ)
美しい整数の世界 http://rio2016.5ch.net/test/read.cgi/math/1686749616/
上
下
前次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
抽出解除
必死ﾁｪｯｶｰ(本家)
(べ)
自ID
ﾚｽ栞
あぼーん
ﾘﾛｰﾄﾞ規制
です｡10分ほどで解除するので､
他のﾌﾞﾗｳｻﾞ
へ避難してください｡
537: １３２人目の素数さん [sage] 2024/04/20(土) 22:42:50.60 ID:bVNPGaYh ■superPCM関数とは？ 奇数の数列2n-1から 合成数を取り除くアルゴリズム Product Combination Mod によって素数を1 合成数を0に振り分ける これはアナログをデジタルに変換する PCM(Pulse Coded Modulation)と 同じ発想 奇数の数列2n-1は乗積Πを掛けると その都度出力されてしまうので、 C(0,3-a)を使って一度だけ出力する Table[Product[(2n-1)^(C(0,3-a)) C(0,C(0,((n-a)^(2a-2)mod(2a-1)))),{a,3,30}],{n,50,232}] ◆aの範囲{a,3,30} 3は固定値、 終値の30は最大50まで設定できる これはnの初期値 しかし、aの終値は40や50に設定しても 30の時と精度に差は生じない http://rio2016.5ch.net/test/read.cgi/math/1686749616/537
538: １３２人目の素数さん [sage] 2024/04/20(土) 23:19:25.82 ID:bVNPGaYh ■合成数はどうやって取り除く？ 奇数の数列1,3,5,7,9,11,13,15,17,19… に対して 数列1,1,0,1,1,0,1,1,0,1,1,0…は a_n=n^2 mod3 数列1,1,1,1,0,1,1,1,1,0,1,1,1,1,0…は a_n=n^4 mod5 これを繰り返してゆくと、 Table[(C(0,n-1))+{(2n-1) {C(0,n-2)+((n+1)^2mod3)} {C(0,n-3)+((n-3)^4mod5)} {C(0,n-4)+((n-4)^6mod7)} {C(0,n-6)+((n-6)^10mod11)} {C(0,n-7)+((n-7)^12mod13)} {C(0,n-9)+((n-9)^16mod17)}},{n,1,180}] {n,1,180}の範囲で精度100％が得られる modの前後の数値を変数aとnで 置き換えると Table[Product[(2n-1)^(C(0,3-a)) C(0,C(0,((n-a)^(2a-2)mod(2a-1)))),{a,3,30}],{n,50,232}] 変数aとnを使うと乗積の計算が入るので 概ね200より大きな素数の判定となる http://rio2016.5ch.net/test/read.cgi/math/1686749616/538
539: １３２人目の素数さん [sage] 2024/04/20(土) 23:38:18.83 ID:bVNPGaYh ▲ ▼ http://rio2016.5ch.net/test/read.cgi/math/1686749616/539
540: １３２人目の素数さん [sage] 2024/04/20(土) 23:41:33.04 ID:bVNPGaYh ■superPCM関数とは？ 奇数の数列2n-1から 合成数を取り除くアルゴリズム PCM(Product Combination Mod) によって素数を1 合成数を0に振り分ける(量子化) これはアナログをデジタルに変換する PCM(Pulse Coded Modulation)と 同じ発想 奇数の数列2n-1は乗積Πを掛けると その都度出力されてしまうので、 C(0,3-a)を使って一度だけ出力する Table[Product[(2n-1)^(C(0,3-a)) C(0,C(0,((n-a)^(2a-2)mod(2a-1)))),{a,3,30}],{n,50,232}] ◆aの範囲{a,3,30} 3は固定値、 終値の30は最大50まで設定できる これはnの初期値 しかし、aの終値は40や50に設定しても 30の時と精度に差は生じない http://rio2016.5ch.net/test/read.cgi/math/1686749616/540
541: １３２人目の素数さん [sage] 2024/04/20(土) 23:52:16.46 ID:bVNPGaYh ■合成数はどうやって取り除く？ 奇数の数列1,3,5,7,9,11,13,15,17,19… に対して 数列1,1,0,1,1,0,1,1,0,1,1,0…は a_n=n^2 mod3 数列1,1,1,1,0,1,1,1,1,0,1,1,1,1,0…は a_n=n^4 mod5 これを繰り返してゆくと、 Table[(C(0,n-1))+{(2n-1) {C(0,n-2)+((n+1)^2mod3)} {C(0,n-3)+((n-3)^4mod5)} {C(0,n-4)+((n-4)^6mod7)} {C(0,n-6)+((n-6)^10mod11)} {C(0,n-7)+((n-7)^12mod13)} {C(0,n-9)+((n-9)^16mod17)}},{n,1,180}] {n,1,180}の範囲で精度100％が得られる +((n-5)^8mod9)と +((n-8)^14mod15)が抜けているが これらは1と0以外を出力するので、 0とのコンビネーションを二回かけて 1と0 だけにする さらに、 modの前後の数値を変数aとnで 置き換えると Table[Product[(2n-1)^(C(0,3-a)) C(0,C(0,((n-a)^(2a-2)mod(2a-1)))),{a,3,30}],{n,50,232}] 変数aとnを使うと乗積の計算が入るので 概ね100より大きな素数の判定となる http://rio2016.5ch.net/test/read.cgi/math/1686749616/541
ﾒﾓ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.024s