25/10/19 21:36:07.90 1trCfbwI.net
>>866
R
URLﾘﾝｸ(ideone.com)
868:デフォルトの名無しさん
25/10/22 21:34:32.83 vm0Iby1T.net
>>866
金額が大きい場合でも高速に求められるようにした。
R
URLﾘﾝｸ(ideone.com)
C++
URLﾘﾝｸ(ideone.com)
869:デフォルトの名無しさん
25/10/26 09:31:44.23 Y3+SSpql.net
お題というか、協力してほしい感じなんですが、素因数分解関数をHaskellで書いて色んな数を素因数分解して遊んでいたら確認したい事実に出くわしたので。
31 <- 素数
331 <- 素数
3331 <- 素数
33331 <- 素数
333331 <- 素数
と、３が５個並んで末尾が１の数字までは素数という事が分かりましたが、いかんせん、ノートだと力不足。
それにCとかで書き直したらもっと先まで行けるかも？という事で、この先、どこまで33...31が素数なのかを調べて欲しいのです。
協力お願いします<(_ _)>
一応、Haskellではこんなコードです。
factorization n = f primes n
where primes = 2:(sieve [3,5..])
where sieve (p:xs) = p:(sieve [x | x <- xs, x `mod` p /= 0])
f (p:ps) n | n <= p = [n]
f (p:ps) n | n `mod` p == 0 = p:f (p:ps) (n `div` p)
f (p:ps) n = f ps n
870:デフォルトの名無しさん
25/10/26 09:50:23.90 Y3+SSpql.net
あ、ただの素数判定でも良いです。
ちなみに、66..61の場合は6661までは素数ですが66661は素数じゃなくなりました。
なので、33..31もどこかで素数じゃなくなるのか？それともずっと素数になりそうなのか？って疑問が持ち上がりました。
871:デフォルトの名無しさん
25/10/26 10:13:01.91 XLS0tlS8.net
>>869
興を削いですまんが、「33...331は素数か」でググったら、AIが（あまり大きくない桁数で）答えを示してくれた…
872:デフォルトの名無しさん
25/10/26 10:19:15.38 0X7G2IAI.net
near-repdigit素数とかで研究されてるらしい
結果だけ知りたいなら↓がまとめてる
URLﾘﾝｸ(stdkmd.net)
873:869
25/10/26 10:21:11.68 XLS0tlS8.net
>>871
ちなみに、以下の思考経路だったのでプログラミング的な思考がゼロだった訳では無い。
1.多倍長整数で組むべきかな？
2.でも64bit整数の範囲で合成数だったら馬鹿馬鹿しいな
3.組む前にカンニングしちゃえ
4.>>871
874:デフォルトの名無しさん
25/10/26 17:08:32.87 Y3+SSpql.net
>>871-873
いえいえ、ああ、やっぱりずっと素数という訳にはいかないんですね…。
何か素数の秘密に触れるヒントか？と心躍ったけど、そんな訳なかったですね(´・ω・｀)
あやうく数学スレで鼻息荒く書き込むところでした。
ありがとうございました<(_ _)>
875:デフォルトの名無しさん
25/10/26 17:55:12.41 N6SeZsiy.net
29bitで収まる範囲内
333333331 = 17 × 19607843
これを求められなかったHaskellはすごく遅い？
876:デフォルトの名無しさん
25/10/26 23:06:07.50 Y3+SSpql.net
>>875
速いアルゴリズムに変えたら何とか1分ほどでその数字まで届きました。
（そもそも、>869 のは美しいとか短いとかの枕詞が付くコードですし。33...31に気付かなかったら4桁ぐらいが実用的ならおｋだったので）
改良版Haskell
pfactorization = f primes
where primes =2:3:5#primes
where n#x@(m:q:y)=[n|gcd m n<2]++(n+2)#last(x:[m*q:y|q^2-3<n])
というか、それより1桁少ない方が少し時間かかりますね。
19607843 < 33333331 なので、素数比較回数が多いのかと。
877:デフォルトの名無しさん
25/10/26 23:07:37.04 Y3+SSpql.net
あ、f の方を忘れた。
f のコードは変更なしです。
878:デフォルトの名無しさん
25/11/07 05:48:17.42 ckPLmv2U.net
>>870
cだとこんな感じでいいのかな？
3333333333333331くらいまで一瞬でできる
#include <stdio.h>
int main(void) {
long long d, n;
printf("Enter a number: ");
scanf("%lld", &n);
for (d = 2; d * d <= n; d++) {
if (n % d == 0)
break;
}
if (d * d <= n)
printf("%lld is divisible by %lld\n", n, d);
else
printf("%lld is prime.\n", n);
return 0;
}
879:デフォルトの名無しさん
25/11/25 16:00:02.31 COUmbZtB.net
お題：ulp()の実装
Pythonのmath.ulp()に相当する関数の実装
特殊値(非数や無限大や非正規化数)の考慮は任意
引数がゼロの場合はDBL_MINやDBL_TRUE_MINに相当する値を
返すこと(プラットホーム依存)
発展的なお題：
共用体やfrexp()やldexp()を使わずに実装
880:デフォルトの名無しさん
25/12/06 16:35:13.06 LjhUSqqq.net
お題：ページャーがある。範囲表示に使う配列を出力せよ。現在のページ番号はp，総ページ数はn、切り取る範囲はrとする。
例
p=1, n=10, r=5
[1,2,3,4,5]
p=3, n=10, r=5
[1,2,3,4,5]
p=5, n=10, r=5
[3,4,5,6,7]
p=8, n=10, r=5
[6,7,8,9,10]
p=10, n=10, r=5
[6,7,8,9,10]
881:デフォルトの名無しさん
25/12/06 19:35:04.35 cqDpuJik.net
>>880
R
URLﾘﾝｸ(ideone.com)
お題に指定がない仕様は適当に決めた。
882:デフォルトの名無しさん
25/12/30 21:09:13.42 NmUg6+Gs.net
>>880
Java
URLﾘﾝｸ(paiza.io)
883:デフォルトの名無しさん
26/03/02 03:36:41.07 e5Z4CqgwT
人手不足なんて嘘八百世界最悪の殺人テ□組織公明党強盜殺人の首魁蓄財3億円超の斎藤鉄夫國土破壞省らクソ公務員利権の弊害が表面化
してるた゛けの話，高額ナマポ集団カンコ一た゛のクソ航空テ口リストだのか゛力による－方的な現状変更によってクソ航空機倍増させて
閑静な住宅地から都心まで数珠つなぎで鉄道の30倍以上もの莫大な温室効果カ゛スまき散らして騷音まみれ静音が生命線の知的産業壊滅
システ厶障害に情報漏洩に誤送金にと連發させて後始末にまて゛莫大な税金て゛私腹を肥やして作為的に気候変動災害連發
わさ゛わざ税金て゛無駄な仕事作って日本中コンクリ一トまみれにしてさらに温暖化加速，世界的にも曰本近海か゛突出して海水温上昇してるし
何人で老人を支えてるかの図も大嘘て゛何人で税金泥棒公務員を支えてるかか゛正しい理解だか゛五輪や万博なんてのも邪悪な利権の最たるもの
諸悪の根源公務員を皆殺しにすれは゛公務員利権のネ夕が原因の人手不足なんてただちに解消して國民の生活か゛豊かになるのか゛現実
ます゛はクソ公務員を最低賃金に統‐して介護た゛の人手不足た゛のってところに出向させろ
航空機連絡先情報 noise.web.fc2.com
884:デフォルトの名無しさん
26/03/19 09:07:43.95 0MXzP+SS.net
URLﾘﾝｸ(youtu.be)
PC88SR FinalZoneと言うゲームのエンディング（動画の57分過ぎ）
横スクロールしながら背景の縦スクロールのプログラムの仕組みは？
コメント欄に作者が降臨してて、バグから生まれた効果と言ってるが詳細までは書かれていない
・約150秒間で横80バイト×縦100ライン相当を、70ラインほど下へスクロール
・各バイトを上から下へ左ローテイト処理しながら転送
・単なる下への転送ではなくある程度揺らぎで、でも全体が崩れるわけでもなく形を整えたままスクロール
ざっと2秒(120回処理）で800バイト動かして1ラインスクロールか？
1回(16ms）で8バイトずつくらいか
ローテイトを見せているので処理が終わりきるまで転送元と転送先をその間記録してるということ
ローテイト8回やって1つ終了となれば、1回（16ms）で8バイト×8回の64個を演算処理してる感じ
それとも単に1ラインスクロールしたあとに、全画面のバイトを左ローテイトしてるだけか？
それなら処理中のバイトを記録しておかなくてよい
885:デフォルトの名無しさん
26/03/19 09:27:44.56 0MXzP+SS.net
んー
1バイト転送してから自分をローテイトするだけだから、転送元と先はもう関係ないな
886:デフォルトの名無しさん
26/04/01 21:47:41.96 vg1bcFzn.net
>>880 octave
URLﾘﾝｸ(ideone.com)
・勝手に解釈変更。p=1, n=10, r=5 →[9,10,1,2,3]
887:デフォルトの名無しさん
26/04/04 17:26:39.83 sHUrEpm3.net
>>880 ocaml
URLﾘﾝｸ(ideone.com)