プログラミングのお題スレ Part22

**デフォルトの名無しさん** · 2024/02/20(火) 08:46:37.56

【㋮㋑㋣㋹㊀㋳】　チャールズ３世戴冠式に｀死神´
https://rio2016.5ch.net/test/read.cgi/kokusai/1690352002/l50

17 · 2024/02/20(火) 10:47:13.25

>>254
しらみ潰しって、どんなテストしたの？

**デフォルトの名無しさん** · 2024/02/20(火) 12:59:46.62

しらみ潰しとは例えば1から順番に見つかるまで全てを試していく最悪な方法を指す
今回の場合だと与えた数から順番に見つかるまで±1を続けて全てを試していって探すプログラムが該当する

9 · 2024/02/20(火) 17:18:07.59

「どんなテストしたの？」
って質問だよ

◆QZaw55cn4c · 2024/02/20(火) 20:48:17.02

>>262
私は >>248
だけれども、解法としてはそれしかないと思いますね

**デフォルトの名無しさん** · 2024/02/20(火) 22:00:55.76

>>249 と >>238 は
しらみ潰しではなく
きちんとプログラミングして算出しているようにみえますね

**デフォルトの名無しさん** · 2024/02/21(水) 13:54:29.89

>>264
私は解答は提出していないが、ざっくりと自分が思いついた方法

まず、以下のような操作を考える
A. 1234という入力に対して1234321を返す
B. 1234という入力に対して12344321を返す
ここで、xという入力に対してA,Bが返す数をA(x),B(x)と表すことにする

次に、与えられた数の桁数で場合分け
(1)与えられた数字の桁数が奇数の場合
例として5桁の数字を考える
N=a*10000+b*1000+c*100+d*10+e*1 (a~eは1桁の自然数, aは0でない)
が与えられたとき、
M=a*100+b*10+c*1
とすると、N=10000の場合を除いて、Nに最も近い回文数は
A(M), A(M+1), A(M-1)
の3つの候補に絞られる(厳密にはA(M)とNとの大小比較からA(M±1)の何れかは明らかに候補にならないので2つを考えれば良い)
N=10000の場合は9999と10001が答え

(2)与えられた数の桁数が偶数の場合
例として6桁の数を考える
(1)と同様に
N=a*100000+b*10000+c*1000+d*100+e*10+f*1
に対して
M=a*100+b*10+c*1
とすると、N=100000の場合を除いて
B(M), B(M+1), B(M-1)
のどれかがNに最も近い回文数（厳密には以下略）
N=100000の場合は99999と100001が答え

十分大きな数に対しては虱潰しに回文判定していくより速く求まる

**デフォルトの名無しさん** · 2024/02/21(水) 16:02:55.06

>>266
N=17
のときは？

**259** · 2024/02/21(水) 23:06:20.13

>>259
Iterable.generate(n).map(f)は単に
Iterable.generate(n, f)で良かったと判明

**デフォルトの名無しさん** · 2024/02/21(水) 23:42:23.78

>>234
>>249をC++で書き換え（入力値は64ビット整数の範囲内限定）
https://ideone.com/e1AM8A

元々はCで書き、4行目はなし、15行目と24行目はstrrev(s + i);だったが、Windowsのgccでは
コンパイルできたのにideoneではできなかったので、仕方なくC++にしてstd::reverseで代用した。

**デフォルトの名無しさん** · 2024/02/22(木) 00:34:50.13

>>234 lisp
>>266を参考に>>249(C#)を移植
https://ideone.com/CUPTas

**デフォルトの名無しさん** · 2024/02/22(木) 01:30:35.61

>>234
Rust

fn nearest_palindrome_numbers(n: usize) -> Vec<usize> {
　let mut dd = DecimalDigits::new(n);
　dd.palindrome_using_upper_half();
　let n1 = dd.to_number();
　match compare(n, n1) {
　　Equal => return vec![n],
　　Greater => dd.increment_upper_half(),
　　Less => dd.decrement_upper_half(),
　}
　if dd.is_most_upper_zero() {
　　return vec![n - 1, n + 1];
　}
　dd.palindrome_using_upper_half();
　let n2 = dd.to_number();
　match compare_absolute_diff((n, n1), (n, n2)) {
　　Less => return vec![n1],
　　Greater => return vec![n2],
　　Equal => return if n1 < n2 { vec![n1, n2] } else { vec![n2, n1] },
　}
}

**266** · 2024/02/22(木) 01:47:44.56

>>267
N=17=1*10+7*1のとき、Nは2桁(偶数桁)でM=1*1=1
B(M)=B(1)=11はNより小さいのでB(M-1)は考えなくてよい
B(M+1)=B(2)=22なので11,22が答えの候補
11より22のほうが17に近いので22が答え

ちょっとNに対するMの説明が足りてなかったけど言葉で上手く言い表せないすみません（上位半分以上かつ最小の桁数を抜き出す、的な）

**デフォルトの名無しさん** · 2024/02/22(木) 20:54:45.16

>>234 ruby
https://ideone.com/rJCYgT
・それっぽい三個の候補から選んでるだけ

**デフォルトの名無しさん** · 2024/02/22(木) 21:04:16.48

>>234
>>269の一部でC++の機能をどうせ使ってしまったので、この際、全部をC++流に変えたら
C流よりすっきり書けた。
https://ideone.com/38bo2E

**273** · 2024/02/22(木) 21:48:36.22

>>273
> [1000, [1001]]

誤：ps = [p.(s), p.(t.to_i.pred.abs.to_s + u), p.(t.succ + u)]
正：ps = [p.(s), p.(t.to_i.pred.abs.to_s + u), p.(t.succ + u), p.(s.to_i.pred.abs.to_s)]

とりあえず雑に修正してみたが？
(ﾉ∀`)ｱﾁｬｰ

17 · 2024/02/23(金) 18:10:28.85

>>262
>>245のKotlinのプログラムは何も考えてなくて本当に馬鹿正直に±1して一つ一つ検査する方式で作ったんだけど、それでもあなたのテストではダメということになったの?
まあ Int (符号付32bit整数) 使ってるからその限界超えたらダメではあるんだけど、そういう問題ではなく?

**デフォルトの名無しさん** · 2024/02/23(金) 18:58:34.05

>>276
それはしらみつぶしと言われる駄目プログラミングだよ
例えば求める解法が存在する方程式を解くのに値を±1しながら順に代入して試していくのと同じ

**デフォルトの名無しさん** · 2024/02/23(金) 21:39:34.55

>>277
あー。プログラムにバグがあってまともに答えが出ないっていうことではなく何の捻りもないプログラムだからダメっていう感想ね。それならわかる。
こちらもアルゴリズム思い浮かばないけどとりあえず作ってみただけだし。ダメというほどではないが良いとも思えないプログラムなので。

**273** · 2024/02/23(金) 23:06:54.56

>>234 ruby
https://ideone.com/E9VSE3
・273の[1000, [1001]]バグ修正版
・275とは違う方法で修正してみたがやっつけ感大

>>234 ruby 2.5.5
https://ideone.com/1zqSr1
・いわゆる（？）ジェネレータ版
・「終端を持たない範囲オブジェクト」はRuby 2.6.0から

**デフォルトの名無しさん** · 2024/02/24(土) 00:25:03.27

>>234
Ruby
https://paiza.io/projects/2G8vPQJQOBZecXPD5ZIFTQ?language=ruby

**279** · 2024/02/24(土) 13:21:16.94

>>234 ruby 2.5.5
https://ideone.com/04fxGM
・ジェネレータ版ちょっとアレンジ
・to_sしてto_iするのをやめた

◆QZaw55cn4c · 2024/02/24(土) 14:25:41.40

異なる自然数 a, b (a > b) における a^3 - b^3 を「a, b の三乗差」と呼ぶことにする。
異なる5通りの組(a, b) (c, d) ... (j, k) について三乗差がすべて相等しいとき
その組(a, b)...(j, k) および三乗差自体を求めよ
異なる6通りの組で三乗差が相等しい場合があるかも検討せよ

**デフォルトの名無しさん** · 2024/02/24(土) 16:47:38.19

>>282
[(1134, 357), (1155, 504), (1246, 805), (2115, 2004), (4746, 4725)]
a^3 - b^3 == 1412774811

17 · 2024/02/24(土) 16:52:08.36

数学、か・・・

**デフォルトの名無しさん** · 2024/02/24(土) 18:16:29.01

>>282 ruby
https://ideone.com/wBX9Rs
・そのまま版
・5秒じゃ㍉

**デフォルトの名無しさん** · 2024/02/24(土) 20:03:55.52

>>234
>>274を巨大整数対応にした。
https://ideone.com/Hg3kGJ

**デフォルトの名無しさん** · 2024/02/24(土) 22:30:55.28

>>282
aが5000以下限定で>>283の解はRで5秒以内に求められた。8～9行目は下三角行列の要素番号から
行番号と列番号を求めているだけで、>>282を解く特別なアルゴリズムというわけではない。
https://ideone.com/w4X3DC

**285** · 2024/02/25(日) 08:52:58.18

>>282 ruby 2.5.5
https://ideone.com/4GdMaO
・285の無駄なループ回数を削減
・でも5秒じゃ㍉

**デフォルトの名無しさん** · 2024/02/25(日) 21:06:57.77

>>282
>>287は下三角行列の要素番号から列番号を求めるのに2次方程式を解くのが分かりにくかったから、
各列の最下行の要素番号を計算しておいてそれを二分探索するように変更した。
https://ideone.com/Tefgkh

n = 25000で実行してみても、n = 5000ときの解のa, bを両方とも2, 3, 4, 5倍した解しか別に
見つからなかった。

**デフォルトの名無しさん** · 2024/02/26(月) 22:18:39.21

>>282
>>289と似た手順でC++
https://ideone.com/nFRsrK

Rではソートする前に重複値だけを抽出しているが、C++のunique関数はソート済みデータにしか
使えないので使っていない。

**288** · 2024/02/27(火) 21:45:30.42

>>282 c
https://ideone.com/eM18H1
・288の移植

**デフォルトの名無しさん** · 2024/02/27(火) 22:30:42.88

>>282
下三角行列の各列内の要素は昇順で既に並んでいるのに、>>290は下三角行列の全要素を
ソートして無駄なので、列のマージに変更（要するにマージソートを途中段階から開始）
したら少し速くなった。
https://ideone.com/EZSvB3

n = 10000でも5秒以内に終わった。
https://ideone.com/huQiBe

**288** · 2024/02/28(水) 22:00:05.18

>>282 c++
https://ideone.com/teo5Mm
・288の移植

>>282 rust
https://ideone.com/G94aN3
・288の移植

**デフォルトの名無しさん** · 2024/02/28(水) 23:21:09.16

>>282
>>292とは別の方法で>>290を高速化
https://ideone.com/QxjmyT

1³, 2³, 3³, …, 5000³をD = 5001で割った余りはすべて異なる値になるから、d = a³ − b³を
Dで割った余りはどれか1つの値に偏ることなく均等に分布する。dをDで割った余りによりdを
区分すれば、各区分に入る個数はどれも多すぎないのでソートに時間が余りかからない。
Dの値はconstexpr関数によりコンパイラに計算させている。n = 10000, 15000のときは
それぞれD = 10002, 15009になる。

**291** · 2024/02/29(木) 22:20:45.85

>>282 c
https://ideone.com/JnJpJW
・291から省メモリ化
　旧：unsigned int values[10];
　新：unsinged short values[4];

**デフォルトの名無しさん** · 2024/03/01(金) 22:22:26.10

>>282
C++
https://ideone.com/1c4s5I
>>294はa, bの二重ループ内でa³ − b³をD = 5001で割った余りrにより区分していたが、
rのループ内でa, bを変化させるように変更したら、2次元配列がなくなってすっきりした。
その結果、メモリ使用量が激減し、nが大きい場合でも実行できるようになった。

**デフォルトの名無しさん** · 2024/03/01(金) 22:23:18.03

>>296の続き
n = 1000000, m = 6で実行すると、12通りの解が見つかった。

[6つ組解]
424910390480793: (75978, 23919), (77385, 33768), (83482, 53935), (141705, 134268), (317982, 316575), (596001, 595602)
620174235433536: (86184, 27132), (87780, 38304), (90237, 48573), (94696, 61180), (160740, 152304), (360696, 359100)
1238805803151000: (107487, 14487), (108540, 34170), (110550, 48240), (119260, 77050), (454260, 452250), (851430, 850860)
1384074844012224: (112152, 29844), (125324, 83600), (130050, 93426), (159372, 138624), (224928, 215412), (357447, 353799)
1936290882196125: (127629, 52254), (133320, 75675), (149285, 111620), (228525, 215430), (246510, 235395), (290214, 282339)
4589726535576000: (170172, 69672), (177760, 100900), (185265, 120945), (304700, 287240), (328680, 313860), (386952, 376452)
4961393883468288: (172368, 54264), (175560, 76608), (180474, 97146), (189392, 122360), (321480, 304608), (721392, 718200)
11072598752097792: (224304, 59688), (250648, 167200), (260100, 186852), (318744, 277248), (449856, 430824), (714894, 707598)
36717812284608000: (340344, 139344), (355520, 201800), (370530, 241890), (609400, 574480), (657360, 627720), (773904, 752904)
52279853819295375: (382887, 156762), (399960, 227025), (447855, 334860), (685575, 646290), (739530, 706185), (870642, 847017)

[7つ組解]
15490327057569000: (249281, 6281), (255258, 104508), (266640, 151350), (298570, 223240), (457050, 430860), (493020, 470790), (580428, 564678)
123922616460552000: (498562, 12562), (510516, 209016), (533280, 302700), (555795, 362835), (597140, 446480), (914100, 861720), (986040, 941580)

6つ組解の(2, 7), (4, 8), (5, 10), (6, 9)番目は各括弧内で自然数比になっている。
6つ組解の5番目の2倍は7つ組解の1番目のうちの6組を構成している。

**デフォルトの名無しさん** · 2024/03/01(金) 22:24:42.81

>>282
C++
https://ideone.com/tM1cuo
>>296のrのループ内でa³ − b³をD2 = 5003で割った余りr2により区分し、それぞれの区分ごとに
解を探すようにしたら速くなった。ただし、nが大きい場合にはかえって遅くなる。

◆QZaw55cn4c · 2024/03/03(日) 19:08:39.58

>>297
出題者です。
すごいです。ありがとうございます。私の手元ではまだ6通り解、7通り解のひとつも入手できていないので、参考になりました
私のアルゴリズムは効率が悪いようですね

**デフォルトの名無しさん** · 2024/03/03(日) 22:19:54.92

>>282
C++
https://ideone.com/LEU7EV

>>298でnを大きくするにつれ>>296に対する高速化効果が薄れていくのは、ABをvectorでなく
配列にしたらある程度改善された。n = 5000のときの実行時間は>>296の半分以下になった。
ただし、n = 1000000まで大きくすると、296よりやっぱり遅くなる。

>>299
どんなプログラムを書いたのか見せて。

**デフォルトの名無しさん** · 2024/03/06(水) 22:35:52.23

>>282
C++
https://ideone.com/PG6UiY
>>300の実行時間を分析すると、最も時間が掛かっているのは46～と47行目だと判明した。
そこで配列ABの第1次元と第2次元を入れ替えてみると、n = 5000では変わらないが、
1万, 2万, 5万, 10万, 20万では35%前後高速になった。これは、改良前には第2次元の添字が
小さい要素に書き込みが集中しているため、改良後のように第1次元に入れ替えた方が
纏まったメモリ領域に書き込みが集中しキャッシュの効きが良くなるからだと考えられる。
一方、n = 100万で高速化しないのは、書き込み集中領域が大きすぎるからだろう。

https://ideone.com/6RzW0n
n = 100万の場合にはr2の値によってデータを多数の列へ振り分けるのをやめ、列を1つにして、
その内部でr2の値により2種類に区分し、それぞれの内部で2種類にさらに区分し、…と再帰的に
区分していけば（要するにクイックソートの変形版）、1つの配列内での要素のスワップだけで済み、
キャッシュの効きが改善されるとの予想通り、n = 100万で実行速度は>>296より25%速くなった。
(原理的には>>300より非効率なのでn = 5000では>>300より当然遅い)

**デフォルトの名無しさん** · 2024/03/08(金) 19:02:53.21

>>301
ハズレが多いから2passは効果ある？

**デフォルトの名無しさん** · 2024/03/09(土) 22:13:47.64

>>282
C++
https://ideone.com/xQD1W8
関数mainのループで配列A, B, Pに書き込まずdにだけ書き込むようにし、関数FindDuplicatesで
dの添字Pではなくdそのものをソートするように変えて、n = 1000000の場合に>>301より10%高速化。
関数PrintSolutionでa, bをmainでと同じ方法で再計算するのは非効率だが、PrintSolutionは僅か12回しか
呼ばれないため、全体の実行時間への影響は無視できる。

**デフォルトの名無しさん** · 2024/03/09(土) 22:47:01.30

お題

460円 580円 600円の3種類の商品があります
これらを組み合わせて合計10個買ったら5360円になりました
組み合わせを求めるプログラムを書いてください

ちなみに答えの一つは
・600円×2
・580円×4
・460円×4
だそうです

https://rio2016.5ch.net/test/read.cgi/cigaret/1706726196/56-57

**デフォルトの名無しさん** · 2024/03/09(土) 23:59:51.39

>>304
面倒なのでRで全探索
https://ideone.com/vrtYvk

**デフォルトの名無しさん** · 2024/03/10(日) 01:20:18.65

>>304
面倒なので全て460円を引くと
A=0円 B=120円 C=140円
10個で760円という問題

面倒なのでさらに20で割ると
A=0円 B=6 C=7円
10個で38円という問題

つまり唯一奇数のCは偶数個が確定
Cが6個以上だと42円以上でオーバーしてNG
Cが4個だと28円で残り10円をA,Bで作れないからNG
Cが2個だと14円で残り24円はBが4個で残り4個がA
Cが0個だと0円で残り38円をA,Bで作れないからNG
つまり解は(A,B,C)=(4,4,2)しかない

**デフォルトの名無しさん** · 2024/03/10(日) 11:20:30.42

>>306
すごすぎるだろ、日本の未来を頼む

**デフォルトの名無しさん** · 2024/03/10(日) 19:06:13.20

>>304
Rで全探索でなくちゃんと解くと
https://ideone.com/F44pCL

解が複数ある場合と全くない場合の例として、600円を540円と520円に変更したときの出力も載せた。

**デフォルトの名無しさん** · 2024/03/10(日) 20:08:55.20

2pass案は多少工夫したらかなり速い

n ␣␣m ␣296␣ ␣301-1 ␣301-2 ␣303␣ ␣2pass
5k␣␣5 ␣ 0.5s ␣ 0.1s ␣ 0.5s ␣ 0.4s ␣ 0.1s
25k ␣5 ␣12.7s ␣ 2.5s ␣13.9s ␣11.1s ␣ 1.7s
100k␣5 ␣3m52s ␣49.3s ␣4m13s ␣3m26s ␣38.9s
1M* ␣6 ␣8h23m ␣2h50m ␣8h51m ␣6h43m ␣1h11m
*n=100万は1万サンプルの部分ループ500k≦r＜510kから100倍

>>301の296と301-2の比較記述と違う傾向があるのはキャッシュ階層の違いだと思う
2passは301-1に近いけど1pass目でのランダムアクセスサイズを落としながらも
誤判定率を低く抑える(0.2%～2%)工夫をするのがお楽しみだと思う

**デフォルトの名無しさん** · 2024/03/14(木) 14:43:15.33

ｔ

**デフォルトの名無しさん** · 2024/03/27(水) 23:42:08.75

>>304

a = (600, 580, 460)
m = min(a)
h = set()

def buy(b, yen):
if yen < m: return
for i in range(0, len(a)):
v = a[i]
if yen >= v:
b[i] += 1
if yen == v:
h.add(str(b))
else:
buy(b, yen - v)
b[i] -= 1

buy([0, 0, 0], 5360)
for s in h: print(s)

**デフォルトの名無しさん** · 2024/03/27(水) 23:55:15.74

>>304
Haskell

[(a, b, c) | a <- [0..20], b <- [0..20], c <- [0..20], a * 460 + b * 580 + c * 600 == 5360]

output: [(0,2,7),(4,4,2)]

**デフォルトの名無しさん** · 2024/03/28(木) 00:00:41.99

合計10個という条件忘れてた。

[(a, b, c) | a <- [0..20], b <- [0..20], c <- [0..20], a + b + c == 10, a * 460 + b * 580 + c * 600 == 5360]

output: [(4,4,2)]

**デフォルトの名無しさん** · 2024/03/31(日) 11:57:53.31

大幅に手直しした
特に前回数値が一部出てこない状態になっていたので色々と手動で最適化した
新しいアイディアを思いつかない限りはシングルスレッドでの限界に近いと思う

n m 301-1 303 2pass 2pass'
5k 5 0.1s 0.4s 0.1s 0.1s
25k 5 2.5s 11.1s 2.3s* 1.7s
100k 5 49.3s 3m26s 38.9s 27.7s
1M* 6 2h50m 6h43m 1h11m 48m10s
2M* 6 17h06m 28h27m 5h47m 3h13m
Max* 6 35h51m 51h23m 11h09m 5h47m

*前回>>309 2pass n=25kの再計測値
*n=1Mは部分ループ500k<=r<510kから100倍
*n=2Mは部分ループ500k<=r<505kから400倍
*Max:=2642245は3乗がUINT64に収まる最大
*n=Maxは部分ループ500k<=r<500k+3785から2642245/3785倍

ヒント含みの数値がこちら

n D1 D2 D3 = 5000 5001 5003 5009
false_positive = 23 / 5001 = 0.46%
total_t_pass1 = 64.220 ms 2.568 ns/iter
total_t_pass2 = 0.044 ms 0.381 ns/iter
real 0m0.097s

**デフォルトの名無しさん** · 2024/03/31(日) 11:58:50.32

n D1 D2 D3 = 25000 25003 25005 25006
false_positive = 171 / 25003 = 0.68%
total_t_pass1 = 1654.681 ms 2.647 ns/iter
total_t_pass2 = 1.407 ms 0.329 ns/iter
real 0m1.709s

false_positive = 2211 / 100005 = 2.21%
total_t_pass1 = 27338.298 ms 2.734 ns/iter
total_t_pass2 = 78.402 ms 0.355 ns/iter
real 0m27.692s

n D1 D2 D3 = 1000000 1000002 1000009 1000015
false_positive = 18 / 10000 = 0.18%
total_t_pass1 = 28674.338 ms 2.867 ns/iter
total_t_pass2 = 5.642 ms 0.313 ns/iter
real 0m28.897s

n D1 D2 D3 = 2000000 2000003 2000013 2000015
false_positive = 13 / 5000 = 0.26%
total_t_pass1 = 28777.424 ms 2.878 ns/iter
total_t_pass2 = 8.620 ms 0.332 ns/iter
real 0m29.015s

n D1 D2 D3 = 2642245 2642246 2642253 2642258
false_positive = 315 / 3785 = 8.32%
total_t_pass1 = 29210.857 ms 2.921 ns/iter
total_t_pass2 = 336.864 ms 0.405 ns/iter
real 0m29.800s

**デフォルトの名無しさん** · 2024/03/31(日) 22:30:39.09

>>304
ぶっちゃけ、他の言語の人と同じっぽくないので心配なんだが…。
自分なりにHaskellで全探索じゃないバージョン書いてみた。

Haskell

[(a, b, c) | a <- [0..10], b <- [0..10 - a], c <- [0..10 - (a + b)], a * 460 + b * 580 + c * 600 == 5360, a + b + c == 10]

答えは同じ[(4,4,2)]。

**デフォルトの名無しさん** · 2024/04/01(月) 04:52:23.91

少し一般化して、N個の商品があり、i番目の商品はA_i円です
合計M個購入し、価格の合計がS円であるような購入の仕方を998244353で割った余りを求めてください
だとO(N M S)より小さい計算量で解けるのかな

**デフォルトの名無しさん** · 2024/04/01(月) 16:50:08.47

2個、4個、8個…みたいにメモ化すればMはlogMにできるかもしれんね
空間がlogM倍されそうだが

**デフォルトの名無しさん** · 2024/04/13(土) 11:43:17.27

ヘロンの公式を実装せよ

使用言語：C

17 · 2024/04/13(土) 16:57:10.76

>>319
https://paiza.io/projects/_ZdSzHtV9YdEzV-oOySQWQ

Wikipedia でヘロンの公式を調べてそのまま実装しただけで、ほとんど何も考えてない。

**デフォルトの名無しさん** · 2024/04/13(土) 23:01:22.75

>>319
https://ideone.com/YCi6qe
ヘロンが作ったもう1つの式である平方根を加算と除算の繰り返しで求める式も使用。
sqrt関数を呼び出すより実行形式ファイルサイズがほんの少しだけ小さくなる。

**デフォルトの名無しさん** · 2024/04/14(日) 00:59:32.83

>>321
無限ループにならない?
機械イプシロン(DBL_EPSILON)とか気になる

**デフォルトの名無しさん** · 2024/04/14(日) 18:34:21.49

解答例

#include <stdio.h>
#include <math.h>

void heron(double, double, double);

int main(void)
{
double a, b, c;
printf("３辺a, b, cを入力せよ ");
scanf("%lf,%lf,%lf", &a, &b, &c);

heron(a, b, c);
}

void heron(double x, double y, double z) // heronの定義
{
double s, t;

s = (x+y+z)*0.5;
t = s*(s-x)*(s-y)*(s-z);

printf("３角形の面積は S=%g\n", sqrt(t));

return;
}

**デフォルトの名無しさん** · 2024/04/14(日) 18:36:52.16

>>321 さすがですね

**デフォルトの名無しさん** · 2024/04/15(月) 21:01:04.41

>>322
p < 0 のとき（= 三角形を作れない場合）は浮動小数点数の特性に関係なく無限ループになる。
sqrt(p) と同様にNANを返すには、if (p < 0) return 0 / (p - p); を追加すれば良い。

p > 0 のときは無限ループにならないはず。以下が検証プログラム。
https://ideone.com/mzemEM

x = sqrt(p), y = p / x とすると、浮動小数点数の特性により x == y とならない場合は存在する。
このとき、xとyの仮数部を整数と見なした値（以降では「仮数整数」と呼ぶ）の差は1なので、
z = (x + y) / 2 はxとyのうち仮数整数が偶数の方に一致する。zを新たなxとして代入しyとzを
再計算すれば、今度はxの仮数整数が偶数なのでzはxに必ず一致し、>>321の収束判定条件が成立する。

具体例で見ると、p = 2 のときはxの仮数整数が奇数なので x != z となるが、zを新たなxとして代入し
再計算すれば x == z が成立する。桁上がりが起こる p = 3.9999999999999996 のときも、同様に
再計算で x == z が成立する。p = 3 のときはxの仮数整数が偶数なので x == z が成立し再計算は不要。

**デフォルトの名無しさん** · 2024/04/15(月) 22:06:46.39

>>325
解説ありがとう
俺には理解できないレベルだと分かりましたｗ
俺なら収束の自信が無くてDBL_EPSILONを使った判定と
ループ回数上限を組み合わせて実装しそうだ

**デフォルトの名無しさん** · 2024/04/17(水) 05:47:35.77

ヘロンの公式はそのままだと、数値計算での安定性が良くないらしいぞ
解決策は、Wikipediaの英語版の方に…
tps://en.wikipedia.org/wiki/Heron%27s_formula#Numerical_stability

**327** · 2024/04/17(水) 05:52:23.77

そしてこんなとこでもカハンせんせーの名前がが

**デフォルトの名無しさん** · 2024/04/17(水) 16:28:33.14

の長さ
この公式で計算される面積は、理論的には正しい値です。しかし、実際には、以下の理由で誤差が生じる可能性があります。

数値計算の誤差: 計算機で数値を扱う場合、有限桁しか扱えないため、丸め誤差が生じます。特に、辺の長さの値が大きく異なる三角形の場合、この誤差が顕著になります。
四捨五入誤差: 計算結果を小数点以下n桁まで表示する場合、n桁目以降の数字を切り捨てます。この四捨五入誤差も、面積の誤差に影響を与えます。

by Gemini

**デフォルトの名無しさん** · 2024/04/17(水) 23:38:33.35

>>327に載っている参考文献
　William M. Kahan, ‘Miscalculating Area and Angles of a Needle-like Triangle’
　http://www.cs.berkeley.edu/~wkahan/Triangle.pdf
のTable 1の問題がパソコン等でのC++プログラムでも再現されるか試してみた。
　https://ideone.com/r4toUS

Table 1とは違い、Accurate Δが概ね正確な場合にHeron's Δ'が大きく懸け離れた不正確な値に
なってしまうことはなく、ほぼ同じ値になり差はごく僅かしかない。Table 1のような不安定性は
Table 1の計算に使われたプログラマブル関数電卓に特有の問題で、パソコン等のプログラムでは
再現されない。(パソコン等のdoubleの方が精度が高いので当然と言えば当然だが)

一方、(a, b, c) = (5278.64055, 94721.35941, 99999.99996)の場合は、逆にHeron's Δ' = 0が
正確なのにAccurate Δ = 9.53674324543714が大きく懸け離れた不正確な値になってしまう
重大な欠点がある。これは、Accurate Δの式の根号内の第2因数c - (a - b)が正確には0なのに
3.63797880709171e-12と計算されてしまい、この誤差が他の因数との乗算により増幅されるから。
Heron's Δ'の式の根号内の第4因数s - cは0と計算されるので問題ない。

double向けの入力値(a, b, c) = (31622.77777777662, 0.000000000023, 31622.77777777661)を
作れば、Heron's Δ' = 2.30085990753844e-07, Accurate Δ = 3.20111707955507e-07となり、
相対差は確かに大きくなるが、200ビットで計算したほぼ正確な値3.27490470056059e-07から
見れば両方とも不正確だから、Accurate Δの利点はない。

だから、パソコン等のプログラムでは改良版の式を使う必要がないどころか使うべきではなく、
ヘロンの公式をそのまま使う方が良い。

**デフォルトの名無しさん** · 2024/04/18(木) 07:16:50.63

>(a, b, c) = (5278.64055, 94721.35941, 99999.99996)
>c - (a - b)が正確には0なのに3.63797880709171e-12と計算されてしまい

この例に限らず、たいていの場合a,b,cはdoubleでexactに格納されて無くて
この例では「c - (a - b)が正確には0」なのをチョイスしただけでは？

**デフォルトの名無しさん** · 2024/04/18(木) 07:30:10.21

パソロジカルな三角形をパラメトライズして面積を積分する検証はどう？
数式計算での正確な値
Heronで面積計算した時の数値積分
Accurateで面積計算した時の数値積分
を比べるのがフェアかなぁと

**デフォルトの名無しさん** · 2024/04/18(木) 07:34:09.77

> 200ビットで計算したほぼ正確な値3.27490470056059e-07
この例だけ見るとAccurate Δの方が優れているように見えるので

>>331の様なチェリーピックはどちらの計算式でも出来るので平均的に近似が近い方が精度的に優れているかと

**デフォルトの名無しさん** · 2024/04/18(木) 22:41:59.70

>>331
その通り。そして、(a, b, c) = (10000.1, 10000.2, 20000.3)とすれば、正しい面積は0なのに
Heron's Δ' = 2.69745899635295とAccurate Δ = 1.34872949817647は両方とも大間違いになる。
この場合のようにHeron's Δ'での問題がAccurate Δで改善されないだけでなく、>>331の引用の
場合のようにHeron's Δ'では結果的に問題ないのにAccurate Δでは新たな問題が生じてしまうのは、
参考文献の11ページで述べられた

　An algorithm stood convicted of numerical instability if it could be replaced by
　a new algorithm at least about as fast and accurate as the old for all data,
　and good for all data for which the old algorithm was bad.

　すべてのデータに対して旧アルゴリズムと少なくとも同じくらい高速かつ正確であり、
　かつ旧アルゴリズムが悪くなるすべてのデータに対して良くなる新アルゴリズムによって
　置き換えることができるとしたら、旧アルゴリズムは数値的に不安定と判定される。

という判定条件を満たさないから、Accurate Δは改良版としての適性を欠く。

>>333
その例では有効桁数がHeron's Δ'は0桁、Accurate Δは1桁しかなく、どちらの品質も絶対的に
劣悪で、それらの間の相対的な優劣に大した意味はない。

そもそも針のように異様に細長い三角形が重箱の隅をつつくような話で、普通はそんな場合は
想定しなくても良く、ヘロンの公式で充分。そこを敢えてつつくなら、ヘロンの公式だけでなく
改良式もぼろが出てしまうだけ。

**デフォルトの名無しさん** · 2024/04/18(木) 22:55:38.47

総合すると有効桁じゃなくて精度が2桁良いし実装上は大差ないから改良版を使う、と言う方が自然では？

**デフォルトの名無しさん** · 2024/05/01(水) 12:56:47.83

スレ落ちそうなのであげ

**デフォルトの名無しさん** · 2024/05/01(水) 15:39:17.16

>>336
嵐を呼び込むために・・・

**デフォルトの名無しさん** · 2024/05/01(水) 22:59:10.72

何でこんなに過疎化しちゃったのか。前に頻繁に出題していた人がいなくなったのか。

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

お題を出してみてください

**デフォルトの名無しさん** · 2024/05/02(木) 16:59:52.63

>>338
お題が出尽くしたってことはあるんじゃないか？
過去のお題拾ってきてそれを投稿すればいいぐらいまでスレが成熟してしまったのでは？

**デフォルトの名無しさん** · 2024/05/02(木) 17:21:22.07

PC買って、脱衣AIで遊びまくってる「
一日一回無料で使えるみたい「
https://mao.5ch.net/test/read.cgi/gymnastics/1322657462/98

17 · 2024/05/02(木) 18:44:04.16

>>340
なるほど。それをやるか。

17 · 2024/05/14(火) 05:34:03.62

じゃあ10年前のこのお題(URLを書くとNGになるようなので書かない)。

プログラミングのお題スレ Part4
115 ：デフォルトの名無しさん：2014/06/21(土) 18:36:45.72 ID:/fMJIWig.net

お題：文字列Aを1回以上繰り返した文字列Bが与えられたとき
文字列Aを求める。ただしAの候補が複数ある場合は最短のものとする。
例
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa -> a
123412312341231234123123412312341231234123 -> 1234123
oxoxoxoxoxoxoxoxxoxoxoxoxoxoxoxoxx -> oxoxoxoxoxoxoxoxx

**デフォルトの名無しさん** · 2024/05/14(火) 17:27:18.46

ttps://ideone.com/KrUq7e
Z-algorithm を使って O(|B|) で解いてみた

警備員[Lv.4][初] · 2024/05/14(火) 20:59:46.84

>>343
R
https://ideone.com/ITR83u

C
https://ideone.com/g2INyj

警備員[Lv.18] · 2024/05/23(木) 14:16:50.64

あれ？どんぐりの都合か？URL書いてあると書けなくなったような？

警備員[Lv.18] · 2024/05/23(木) 14:17:56.56

URLの先頭のhを抜いて書いてみよう。

>>343
Kotlin

こちらは普通に自作したやつ。
ttps://paiza.io/projects/OYy-A5rfKg7RqLzv-DKMIA

こちらは正規表現使ってとても小さくなったやつ。
ttps://paiza.io/projects/jgmtMRDhKfcjYfGAglEl3g

**デフォルトの名無しさん** · 2024/06/01(土) 10:16:34.91

お題:　コロン区切りの時分秒の時刻が与えられるので時分秒をそれぞれ掛け算した結果を表示せよ

例:
04:05:06
120

**デフォルトの名無しさん** · 2024/06/01(土) 11:08:12.83

お題: バイト列が与えられる。先頭から解析した場合にバイトが1だったら次の4バイトを読み込んで整数として出力し、バイトが2だったら次のバイトを0が来るまで読み込んで文字列として出力せよ

入力
1 1 0 0 0 2 65 66 67 0 1 128 0 0 0

出力
1ABC128

**デフォルトの名無しさん** · 2024/06/01(土) 12:57:50.32

知らんがな

**デフォルトの名無しさん** · 2024/06/01(土) 23:31:08.51

>>348
R
https://ideone.com/x4hIYG

>>349
C (データ識別子は1か2しかないものとし、整数のエンディアンは実行環境依存とする)
https://ideone.com/aJts6n

警備員[Lv.19] · 2024/06/02(日) 04:45:03.04

>>348
Perl

bash のコマンドラインから入力して実行(ワンライナー)

$ perl -ne 'if(/(\d+):(\d+):(\d+)/){print $1*$2*$3,"\n"}else{print"入力エラー\n"}'
1:2:3
6
3:4:5
60
04:05:06
120
$

警備員[Lv.20] · 2024/06/02(日) 05:19:02.45

>>349
Kotlin
https://paiza.io/projects/lk2q4-Ciny3H2NptrkAbcw

**デフォルトの名無しさん** · 2024/06/03(月) 13:25:02.42

>>348
Windows のPowershell 上で、Ruby の1-liner を使う

末尾の改行を削除して、: で分割して、
文字列を数値型に変換してから、全ての要素を掛ける。
%Q で、ダブルクォーテーションをエスケープする。つまり、split(":")

echo '01:2:09' | ruby -ne 'puts $_.chomp.split(%Q[:]).map(&:to_i).inject(:*)'

18

**デフォルトの名無しさん** · 2024/06/03(月) 15:28:01.42

>ruby -ne 'p eval($_.gsub(":","*"))'
1:2:9
18

**デフォルトの名無しさん** · 2024/06/07(金) 06:27:47.87

お題：引数sとnを取りシーザー暗号化を行う関数を作れ
sは平文、nはずらす文字数（負数可）、返り値は暗号化後の文字列
同様の関数で「Hello, World!」を暗号化し復号化せよ

**デフォルトの名無しさん** · 2024/06/07(金) 09:04:03.36

#! ruby

class String
def to_c(n)
if %r|^n|=~n
n=(n.sub(%r|^n|,"").to_i+26)%26
lb=("A".."Z").to_a.join
sb=("a".."z").to_a.join
la=lb[n..25]+lb[0..n-1]
sa=sb[n..25]+sb[0..n-1]
return self.tr(lb,la).tr(sb,sa)
else
return self
end
end
end

p "Hello,World!".to_c("n3") #=>"Khoor,Zruog!"
p "Hello,World!".to_c("n-5") #=>"Czggj,Rjmgy!"
p "Hello,World!".to_c("s") #=>"Hello,World!"
p "Khoor,Zruog!".to_c("n-3") #=>"Hello,World!"
p "Czggj,Rjmgy!".to_c("n5") #=>"Hello,World!"

**デフォルトの名無しさん** · 2024/06/07(金) 09:04:13.66

#! ruby

class String
def to_c(n)
if %r|^n|=~n
n=(n.sub(%r|^n|,"").to_i+26)%26
lb=("A".."Z").to_a.join
sb=("a".."z").to_a.join
la=lb[n..25]+lb[0..n-1]
sa=sb[n..25]+sb[0..n-1]
return self.tr(lb,la).tr(sb,sa)
else
return self
end
end
end

p "Hello,World!".to_c("n3") #=>"Khoor,Zruog!"
p "Hello,World!".to_c("n-5") #=>"Czggj,Rjmgy!"
p "Hello,World!".to_c("s") #=>"Hello,World!"
p "Khoor,Zruog!".to_c("n-3") #=>"Hello,World!"
p "Czggj,Rjmgy!".to_c("n5") #=>"Hello,World!"

**デフォルトの名無しさん** · 2024/06/07(金) 21:13:26.78

>>356
R
https://ideone.com/Qj6B3w

**デフォルトの名無しさん** · 2024/06/07(金) 23:24:01.54

Haskell
範囲外の数値は平文字をそのまま返すこととした。

import Data.Char

cearsar n |(-26) <= n && n <= 26 = map (f n)
where
f x = chr.(+x).ord
cearsar _ = cearsar 0

sample:
ghci> cearsar 50 "Hello, World!"
"Hello, World!"
ghci> cearsar 3 "Hello, World!"
"Khoor/#Zruog$"
ghci> cearsar (-3) "Khoor/#Zruog$"
"Hello, World!"

**デフォルトの名無しさん** · 2024/06/07(金) 23:28:58.71

TA = [ * ?\x20 .. ?\x7E ]
TS = TA.join
def caesar( s, n ) s.tr( TS, TA.rotate( n ).join ) end

s = "Hello, World!"
p caesar( s, 0 ) #=> "Hello, World!"
p caesar( s, 1 ) #=> "Ifmmp-!Xpsme\""
p caesar( s, -1 ) #=> "Gdkkn+~Vnqkc "
p caesar( s, 20240607 ) #=> "Jgnnq.\"Yqtnf#"
p caesar( caesar( s, 20240607 ), -20240607 ) #=> "Hello, World!"
p caesar( 'HAL9000', 1 ) #=> "IBM:111"

9 · 2024/06/11(火) 14:41:10.23

>>348 Perl5

($x = '04:05:06') =~ s/:/*/g;
print eval $x;

9 · 2024/06/13(木) 14:34:57.00

>>349
> 「バイトが1だったら次の4バイトを読み込んで整数として出力し、」
正直、意味がわからんかった

例で見ると
1 1 0 0 0 → 1
1 128 0 0 0 → 128
ということだが

1に続く4バイトを加算して出力するって意味だったのかいな

**デフォルトの名無しさん** · 2024/06/13(木) 14:48:43.48

>>363
整数が4バイト型という脳内補完だった

**デフォルトの名無しさん** · 2024/06/13(木) 14:54:19.57

>>363
出題者がエンディアンを知らなくて説明もなくリトル環境を前提にしてしまっている
エンディアンを知っている人たちは出題には書かれてないけど例よりリトル前提だと読み取ってこたえている

**デフォルトの名無しさん** · 2024/06/13(木) 14:59:03.73

>>365
リトルエンディアンはビットが逆って事じゃ無いぞ

10 00なら
00 10だぞ

**デフォルトの名無しさん** · 2024/06/13(木) 15:00:51.22

>>366
例>>363を見てないのか？

9 · 2024/06/13(木) 17:03:39.85

>>365
ああそういうことか「4バイトを読み込んで整数」と書いてあるのはそういう意味だったのか
ならわかるかも。
オレは4バイト一個一個が整数だと捉えて、それを「4バイトを読み込んで整数」とは何のこっちゃと？になってたわ

9 · 2024/06/13(木) 17:07:01.55

すまんね68系で育ったもんですぐ連想できなんだ

**デフォルトの名無しさん** · 2024/06/14(金) 21:10:23.18

>>367
例が間違えてるか説明が足りて無い

要は4バイトを読み込んでと説明してるが「一気に4バイト読み込む」とおかしくなる
1バイトずつ順に4バイトを読み込んでという説明なら例が腑に落ちる

**デフォルトの名無しさん** · 2024/06/14(金) 21:15:46.04

例を示しての出題だから

警備員[Lv.23] · 2024/06/15(土) 16:15:42.96

>>356
Kotlin

面倒なのでASCIIコード(0x20-0x7e)でしかシフトしないやつを作った。
まあでも Kotlin は Java 同様に内部でUnicodeで扱っているので平仮名とか漢字とか全然違う言語の文字とかも比較的楽に追加できると思う。

https://paiza.io/projects/5H9H1zSjDnVshGCf4JaQJg

**デフォルトの名無しさん** · 2024/06/19(水) 15:22:47.09

>>348 awk
https://ideone.com/dMvvZG

**デフォルトの名無しさん** · 2024/06/20(木) 17:43:48.64

お題：迷路生成を様々な言語で

例：
C
https://ideone.com/a527mc

**デフォルトの名無しさん** · 2024/06/20(木) 20:17:47.11

https://paiza.io/projects/a6ZS3co-gsEV1tN57VUctA

**デフォルトの名無しさん** · 2024/06/21(金) 02:17:42.15

迷路。やっつけで汚い。乱数自前

C
https://paiza.io/projects/5vyyygrG7exduzqDa-jBAA

Rust
https://paiza.io/projects/g1UmjFWydgmjAEZAgqS5xw

**デフォルトの名無しさん** · 2024/06/21(金) 05:01:41.33

乱数改良

C
https://paiza.io/projects/UWRNkhfbdscMIwR0bTvqHA

Rust
https://paiza.io/projects/g1UmjFWydgmjAEZAgqS5xw

**デフォルトの名無しさん** · 2024/07/26(金) 09:46:37.73

お題：配管サイズの「A呼称」「B呼称」の相互変換
tps://www.keyence.co.jp/ss/products/process/flowmeter/technique/size.jsp
などに掲載されている「A呼称」から「B呼称」への変換。及びその逆変換。
「A呼称」は整数値で(8とか)、「B呼称」は実数(0.125とか)または分母を8とする分数の分子(1とか)で表すものとする。
規格に存在しない場合(42Aとか)は考慮しなくてもよい。
表を引くだけなら簡単過ぎるので、ツェラーの公式の様な「技巧的」な解法を求む。

**デフォルトの名無しさん** · 2024/07/27(土) 06:14:08.30

REM 呼称変換.bat
chcp 65001

ruby -x "%~f0"
ruby tmp.rb
goto end

#~
#!ruby
#encoding: utf-8
Encoding.default_external='UTF-8'

require 'base64'

st=<<'EOS'
IyFydWJ5CiNlbmNvZGluZzogdXRmLTgKRW5jb2RpbmcuZGVmYXVsdF9leHRl
cm5hbD0nVVRGLTgnCgpjbGFzcyBPYmplY3QKICBkZWYgdG9fZmwKICAgIHJl
dHVybiBldmFsKHNlbGYuc3BsaXQoIisiKS5tYXB7fGV8CiAgICAgIGUuc3Bs
aXQoIi8iKS5tYXB7fGZ8IGYudG9fZi50b19zfS4KICAgICAgam9pbigiLyIp
fS5qb2luKCIrIikpCiAgZW5kCmVuZAoKY2xhc3MgTnVtZXJpYwogIEBAYT1b
Niw4LDEwLDE1LDIwLDI1LDMyLDQwLDUwLDY1LAogICAgICA4MCw5MCwxMDAs
MTI1LDE1MCwxNzUsMjAwLDIyNV0KICBAQGI9WyIxLzgiLCIxLzQiLCIzLzgi
LCIxLzIiLCIzLzQiLCIxIiwKICAgICAgIjErMS80IiwiMSsxLzIiLCIyIiwi
MisxLzIiLCIzIiwKICAgICAgIjMrMS8yIiwiNCIsIjUiLCI2IiwiNyIsIjgi
LCI5Il0KICBAQGM9QEBiLm1hcHt8ZXwgZS50b19mbH0KICBkZWYgYV90b19i
CiAgICBuPXNlbGYKICAgIGlmIG48MjUwCiAgICAgIGE9NgogICAgICBAQGEu
ZWFjaHt8ZXwgYT1lIGlmIGUvbjw9MS4wfQogICAgICBiPUBAYltAQGEuaW5k
ZXgoYSldCiAgICAgIGM9QEBjW0BAYS5pbmRleChhKV0KICAgICAgcmV0dXJu
ICIoYT0je2F9KSBiPSN7Yn09I3tjfSIKICAgIGVsc2UKICAgICAgYT0obi81
MCkudG9faSo1MAogICAgICBiPWEvMjUKICAgICAgcmV0dXJuICIoYT0je2F9

**デフォルトの名無しさん** · 2024/07/27(土) 06:15:03.76

KSBiPSN7Yn0iCiAgICBlbmQKICBlbmQKICBkZWYgYl90b19hCiAgICBuPXNl
bGYKICAgIGlmIG48MTAKICAgICAgYz0wLjEyNQogICAgICBAQGMuZWFjaHt8
ZXwgYz1lIGlmIGUvbjw9MS4wfQogICAgICBhPUBAYVtAQGMuaW5kZXgoYyld
CiAgICAgIGI9QEBiW0BAYy5pbmRleChjKV0KICAgICAgcmV0dXJuICIoYj0j
e2J9PSN7Y30pIGE9I3thfSIKICAgIGVsc2UKICAgICAgYj1uLnRvX2kKICAg
ICAgYT1iKjI1CiAgICAgIHJldHVybiAiKGI9I3tifSkgYT0je2F9IgogICAg
ZW5kCiAgZW5kCmVuZAoKd2hpbGUgMQogIGFyPVsi77yh5ZG856ew44GL44KJ
77yi5ZG856ewIiwi77yi5ZG856ew44GL44KJ77yh5ZG856ewIiwi57WC5LqG
Il0KICBtc2c9YXIubWFwLndpdGhfaW5kZXh7fGUsaXwiI3tpKzF9OiN7ZX0i
fS5qb2luKCJcbiIpKyJcbiIKICBwcmludCBtc2crImlucHV0IG51bWJlciA+
ICIKICBuPSRzdGRpbi5nZXRzLnRvX2kKICBicmVhayBpZiBuPT0zCiAgbmV4
dCBpZiBuPT0wCiAgd2hpbGUgMQogICAgcHJpbnQgIiN7YXJbbi0xXX06aW5w
dXQgdmFsdWUgPiAiCiAgICBtPSRzdGRpbi5nZXRzLmNob21wCiAgICBicmVh
ayBpZiBtPT0iIgogICAgbT1tLnRvX2ZsCiAgICBwdXRzIG0uYV90b19iIGlm
IG49PTEKICAgIHB1dHMgbS5iX3RvX2EgaWYgbj09MgogIGVuZAplbmQK
EOS

File.write("tmp.rb",Base64.decode64(st))
__END__
:end

**デフォルトの名無しさん** · 2024/07/30(火) 12:03:21.52

お題
C言語のトライグラフ(可能ならダイグラフも)と普通のテキストとの相互変換。

ファイルまたは標準入力から読んで変換した結果を標準出力に出力する。
面倒なら文字列変換する関数とそれをテストするメイン関数のみでも良い。

警備員[Lv.8] · 2024/08/03(土) 07:16:23.70

お題
Unicode 文字列を UTF-9 へ変換する。また UTF-9 を Unicode 文字列に変換する。
ただし1バイトが9bitではないコンピュータを使用する場合は9bit以上で扱いやすいbit数の変数(例えば 16bitの変数)を代用して下9bitのみを使用する等しても良い。
UTF-9 の仕様は RFC 4042 を見るか、または下記URLのページを参照。
https://www.wdic.org/w/WDIC/UTF-9%20%28RFC%29

警備員[Lv.9] · 2024/08/03(土) 17:38:42.19

>>381
Kotlin

Digraph, Trigraph の変換と逆変換

https://paiza.io/projects/FPUStm3O4tZZMYi1Cw2ruw
https://paiza.io/projects/4vD2Ux8Jd6EzLIynkzTXtw

Map と正規表現使った簡単な変換なので他の言語で作っても大差ないと思う。

**デフォルトの名無しさん** · 2024/08/03(土) 20:58:23.65

お題
整数の格子がある（伝われ）
任意の2点間を線で結ぶ（座標成分は実数）
格子と交わるところで線を分割せよ

**デフォルトの名無しさん** · 2024/08/03(土) 21:19:51.17

>>384
jwwの外部変形としてこんな感じで交点で分割させるの作ったことあるけどプログラミングスレののお題としては座標指定もなく漠然としすぎていていかがなものか。ちょっと伝わんなかった。

**デフォルトの名無しさん** · 2024/08/03(土) 21:36:58.95

ごめん、画像から察してほしい
遊びだから曖昧なところは適当で
https://i.imgur.com/YUKpLcI.png

**デフォルトの名無しさん** · 2024/08/03(土) 21:50:20.46

>>386
二直線の交点は連立方程式を解けば良くて行列を使えば簡単だけど、どういう回答が求められてるのかが分からず途方に暮れております

**デフォルトの名無しさん** · 2024/08/03(土) 22:08:55.42

>>386
出直し

**デフォルトの名無しさん** · 2024/08/03(土) 22:28:33.50

クイズと「プログラミングのお題」との決定的な違い
検証するための入力(問題)と出力(解答)例が明記されているかどうか

**デフォルトの名無しさん** · 2024/08/03(土) 22:29:35.16

テンプレに書いてなかったしいいやって思っちゃった
ここまで伝わらんもんなのか…

**デフォルトの名無しさん** · 2024/08/03(土) 22:35:56.24

具体例を計算するのはめんどくさいけど、例えばさっきの絵で言えば右上が始点で左下が終点なら始点側から順番に座標を出力するとかね

**デフォルトの名無しさん** · 2024/08/03(土) 22:41:57.69

入力例とその時の出力例を
データとして数種類用意すればOK
例えば>>234とか>>343を参考

**デフォルトの名無しさん** · 2024/08/03(土) 22:42:49.82

>>390
そういうことは伝える努力をしてから言おうな

**デフォルトの名無しさん** · 2024/08/03(土) 22:44:30.36

めんどくさい
そこまでして出そうと思わんから興味ないならスルーで

**デフォルトの名無しさん** · 2024/08/03(土) 22:52:49.60

この頭の悪さを見るに学校の課題か何かを丸投げしようとしてるのかな

**デフォルトの名無しさん** · 2024/08/03(土) 23:04:31.31

入出力例があると仕様が定まってプログラミングのお題が成立する
入出力例がないと宿題を手伝わせようとしているだけの可能性も否定できないね

**デフォルトの名無しさん** · 2024/08/04(日) 00:49:43.30

>>394
ドンマイ、出題の意図は分かった気がするので改めて出題
お題
１０mm方眼用紙、左下原点、入力された任意の直線と方眼線との交点を求めよ
（例）始点（3.5,5.2）終点（6.3,8.1）

**デフォルトの名無しさん** · 2024/08/04(日) 04:51:29.74

>>384
nodejs
https://www.mycompiler.io/view/8GupmXA8eQS

**デフォルトの名無しさん** · 2024/08/09(金) 23:38:00.57

お題：与えられた線分を指定された長さLで分割せよ。ただし最後に余る分は長さLでなくてもよい。
線分([始点,終点])=[[0,0],[100,100]]
L=50
↓
[[[0,0],[50,50]],[[50,50],[100,100]]]

警備員[Lv.7] · 2024/08/10(土) 04:24:35.61

>>399
その例おかしくない？それか問題がおかしいか。
長さがLということはx=yの線を分割する場合は座標はx, y 共にL/√2倍にならないか?

**デフォルトの名無しさん** · 2024/08/10(土) 08:05:08.91

ごめん、ミスった

**デフォルトの名無しさん** · 2024/08/10(土) 08:07:10.85

修正
線分([始点,終点])=[[0,0],[120,0]]
L=50
↓
[[[0,0],[50,0]],[[50,0],[100,0]],[[100,0],[120,0]]]

**デフォルトの名無しさん** · 2024/08/10(土) 08:54:23.54

REM 座標計算.bat
chcp 65001
ruby -x "%~f0"
goto end
#~
#!ruby
#encoding: utf-8
Encoding.default_external='UTF-8'
require "matrix"
while 1
print "input start point > " #0,0
p1=$stdin.gets.split(",").map{|e| e.to_f}
print "input end point > " #120,0
p2=$stdin.gets.split(",").map{|e| e.to_f}
print "input length > " #50
l=$stdin.gets.to_f
exit if l<=0
p1[1]=0 unless p1[1]
p2[1]=0 unless p2[1]
v1=Vector.elements(p1)
v2=Vector.elements(p2)
m=(v2-v1).r
n=(m/l).to_i
a=(0..n).map{|e| l*e}
b=a.map{|e| v1+(v2-v1)*(e/(v2-v1).r)}.map{|e| e.to_a}
b<<p2
p b.uniq
end
__END__
:end

**デフォルトの名無しさん** · 2024/08/20(火) 15:28:18.45

お題:今日、または任意の日付から、もういくつ寝るとお正月かを求めよ。（昼寝は除く）

**デフォルトの名無しさん** · 2024/08/20(火) 18:03:45.73

>>404 lisp
https://ideone.com/RhvMjC

**デフォルトの名無しさん** · 2024/08/21(水) 11:05:18.96

>>404
徹夜はどうしますか？

警備員[Lv.12] · 2024/08/22(木) 08:01:35.77

>>404
bash のコマンドラインから expr コマンドと date コマンドを利用して次のようにして計算した。
実際に日数が経過後に1月1日になっているかの確認もした。

$ expr $( date -d '2024/12/31' +%j ) - $( date +%j ) + 1
132
$ date -d 'now +132 days'
2025年 1月 1日水曜日 07:58:54 JST
$

警備員[Lv.12] · 2024/08/22(木) 08:04:39.02

>>407 をちょっと改良
今年に限定しないやりかた。

$ expr $( date -d $(date +%Y)'/12/31' +%j ) - $( date +%j ) + 1
132
$

**デフォルトの名無しさん** · 2024/08/24(土) 01:28:41.99

ロケール設定なしとかジョークかよ

**デフォルトの名無しさん** · 2024/08/24(土) 07:25:49.23

こまけぇこたぁいいんだよ（AA略）

**デフォルトの名無しさん** · 2024/08/24(土) 23:36:40.04

てかこのお題で私、期間中にヨーロッパから帰国しますんでってな特殊ケースまで折り込めって方が無理に思えるけど、ロケール設定した模範解答をどぞ、

**デフォルトの名無しさん** · 2024/08/24(土) 23:58:49.15

せっかくだから､旧正月も考慮して

**デフォルトの名無しさん** · 2024/09/01(日) 04:58:12.84

お題：複数の文字列が格納されている配列と配列の添字がある。標準入力からnextが入力されたら文字列を表示し添字を進め、prevが入力されたら文字列を表示して添字を戻す。添字が始点と終点に達したら始点は終点に、終点は始点に移動させる。動作は例に合うように実装せよ。

v = ["aaa","bbb","ccc"]
next
aaa
next
bbb
prev
aaa
next
bbb
next
ccc
next
aaa

**デフォルトの名無しさん** · 2024/09/01(日) 08:19:39.59

>>413 lua
https://ideone.com/OEEqAA

**デフォルトの名無しさん** · 2024/09/01(日) 11:51:46.92

>>413
REM nextprev.bat
ruby -x "%~f0"
exit
#~
#!ruby
#encoding: utf-8
i=0;v=["aaa","bbb","ccc"]
while 1
puts v[i]
print "input next or prev [n/p]> "
a=$stdin.gets.chomp
if ""==a
exit
elsif /^n/=~a
i=i+1
elsif /^p/=~a
i=i-1
end
i=(i+3)%3
end

**デフォルトの名無しさん** · 2024/09/03(火) 07:13:35.79

お題
引数aから引数bまでの数列を返す関数がある。この関数を数列にゾロ目が含まれていた場合に空の数列を返すように実装せよ

**デフォルトの名無しさん** · 2024/09/03(火) 07:56:47.18

>>416 lisp
https://ideone.com/WMC2aN

**デフォルトの名無しさん** · 2024/09/03(火) 09:22:31.13

>>417
windowsを使ってるんですが>>417を実行するにはどのような環境を整えれば良いでしょうか?

**デフォルトの名無しさん** · 2024/09/03(火) 10:00:07.45

417ではないが，clispかsbclをインストールすればいいよ．Windowsでもできるし，WSLならさらに簡単にできる

**デフォルトの名無しさん** · 2024/09/03(火) 10:17:03.41

>>419
ありがとうございます
そうするとemacsが使えるようになるのでしょうか?
それともプログラミング言語としてのlisp環境が整うだけでしょうか?

**デフォルトの名無しさん** · 2024/09/03(火) 13:27:22.48

>>418
Windowsは不断使いしていないんで開発環境とかemacsとかはわからないけど、動かすだけなら
https://www.sbcl.org/platform-table.html からWindows版のSBCLインストーラーを落としてインストール

コマンドプロンプトなどのコマンド行で
sbcl --script ファイル名とすると「ファイル名」のファイルをスクリプトとして実行します

--script を付けない場合REPLが起動します(Pythonでいうところの対話モード)
これの終了は (quit) と入力してEnterキー

**デフォルトの名無しさん** · 2024/09/03(火) 13:38:16.60

ちょっと試すだけならウェブサービスを使う手も

Common Lispを動かせるオンラインコンパイラーサイト
https://rextester.com/l/common_lisp_online_compiler
https://www.jdoodle.com/execute-clisp-online
https://wandbox.org/

これらはideoneとは異なり "Share" とか "Save" とかをクリックしない限り保存・公開されない(はず)

**デフォルトの名無しさん** · 2024/09/03(火) 22:13:24.20

レス頂いた方、ありがとうございます
とりあえず>>422辺りで試してみます

**デフォルトの名無しさん** · 2024/09/04(水) 03:34:25.27

REM zorome.bat
ruby -x "%~f0"
exit
#~
#!ruby
while 1
print "input a > "
a=$stdin.gets.chomp
exit if a==""
a=a.to_i
print "input b > "
b=$stdin.gets.to_i
ar=(a..b).to_a
z=0
ar.each{|e|
n=e.to_s.gsub(/^-/,"")
m=n.split("").uniq
if n.size>1 && m.size==1
z=1
break
end
}
ar=[] if z==1
p ar
end
__END__
a=0;b=10 => [0,1,2,3,4,5,6,7,8,9,10]
a=0;b=11 => []

**デフォルトの名無しさん** · 2024/09/04(水) 06:34:37.52

>>416
fn is_zorome(n: i32) -> bool {
let s = n.abs().to_string();
let b = s.as_bytes();
if b.len() > 1 {
b[1..].iter().all(|&d| d == b[0])
} else {
false
}
}

fn odai(a: i32, b: i32) -> Vec<i32> {
if (a..=b).any(|n| is_zorome(n)) {
vec![]
} else {
(a..=b).collect()
}
}

fn main() {
assert_eq!(odai(1, 10), [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]);
assert_eq!(odai(1, 11), []);
}

警備員[Lv.7] · 2024/09/29(日) 14:12:38.24

>>416
Kotlin
(URL投稿禁止になってるようなので先頭の https:／／は抜きで)。
paiza.io/projects/nPK7U1uDtt4b7O2saECP0A

[入力]タブの方で開始、終了の値を入れている。とりあえず 10 40 を入れてある。
一桁の数を指定した場合は全てゾロ目扱いになるので出力は空のリストになる。

**デフォルトの名無しさん** · 2024/10/08(火) 17:17:22.09

お題：文字列を中央寄せする関数を作れ。
関数は文字列、全体の横幅（整数）、余白を埋める文字、を引数に取る。横幅分の文字列を引数の文字で作成し、その中央に引数の文字列を埋め込み戻り値として返せ。

**デフォルトの名無しさん** · 2024/10/08(火) 19:52:28.65

py
text.center

警備員[Lv.12] · 2024/10/09(水) 00:07:23.87

>>427
文字幅についてはどう考えれば良い？
エンコーディングが UTF-8 や UTF-16 だとバイト数や要素数が画面の表示幅と一致しないわけだが。
一致するような Shift-JIS だったとしても端末のフォント設定次第で表示幅が変わる。

**デフォルトの名無しさん** · 2024/10/09(水) 00:38:16.24

EUCにも半角2バイト文字なんてあったりする

**デフォルトの名無しさん** · 2024/10/09(水) 00:56:06.70

asciiの前提では
全角判定関数は切り分けるべきだし

**デフォルトの名無しさん** · 2024/10/09(水) 01:18:55.72

#ドラクエ３が超好きだったので、お題にしてみる。ただしロト３のファミコン版のドラクエな！ｗ
"こちらの職業は勇者、戦士、空手家、僧侶、柔道家の５人がかりｗ
勇者から柔道家までの5人の攻撃回数が異なります。"まずは変数と確率の問題。
勇者＝you戦士＝sen,空手家＝kara,僧侶=sou
まず勇者は1~２回攻撃、戦士は1~３回攻撃、空手は2~３回攻撃、柔道は足払いで相手を２０％で麻痺させることができる。（１ターン何もできない）１０％では超麻痺。（２ターン何も敵は動かない）
５％でクリティカルヒットの足払いで３ターン敵は何もできない。
柔道は敵に攻撃を加えられない分、足払いで相手を痙攣させる。
（勇者の１回分の攻撃は100〜200、２０％の確率でクリティカルヒットの攻撃になることがある。
（戦士の攻撃は２００〜３００、２０％の確率で攻撃がクリティカルヒット１回になることがある。
（空手家）空手家の攻撃は３５０〜５００、ただし２０％の確率で攻撃がクリティカルヒットとなるときがある。youはぞろ目（１１１とか２２２とか９９９まで）の場合はロトの剣を取得する。
攻撃回数無限で全てクリティカルヒットになります。
クリティカルヒットは攻撃力X２です。さてプログラムをやっていくか。"
＊ただし、ロトの剣は勇者しか使えない。
”まずギルガメッシュっていう名前のenemy
ラスボスの名前はenemyとする。ボスのHPは5万。
ボスは基本的に何もしません、できません。これやるとさすがに面倒なのでｗ
基本的にボスは攻撃しません。（面倒だからｗ）
一方的になかまあつめてボコすｗ”
＊ロトの剣で攻撃したら攻撃力が1.5倍あがります。＊
＊＊癒しの穂。食べたら防御力と賢さがあがります＊＊
ロトの剣の出現条件は3並びで一つ食べることができます。
例えば１１１攻撃とか２２２攻撃とか３３３攻撃とか〜９９９まで。
次の発言からコードを書いていく。
ニーズあるか？

**デフォルトの名無しさん** · 2024/10/09(水) 01:22:59.53

プログラミングする以前に読むのが面倒くさい

**デフォルトの名無しさん** · 2024/10/09(水) 01:43:57.73

しまった。初期化してなかったなｗ変数のところ。
あとsouの回復もホイミ、ベホイミ、ベホマラーの３種類だけ。
その３種類がぞろ目の時（１１１，２２２など）
そしたらベホマズンを使える。ほぼ運ゲーですｗ
人修羅みたいに取り巻きがいます。それはバラモスエビルと神龍です。
ある一定時期を過ぎると取り巻きが一斉にHPを連続してあげまくるので
取り巻き先に倒さないと厳しいです。取り巻き＝tori種類１の取り巻き＝tori0
種類２の取り巻き＝tori1　5人だとなかなか倒せないから６人目の魔法剣士＝maken
だれか
これをいれよう。
さてif文を完成させよう！if zoma =1
print("ゾーマが現れた”）
import random
random.random()x = random.uniform(0.1,0.9)
print(x)
if x >200 :
print( 'enemyに' +200のダメージを与えた)
elif x > 0.5 :
print( 'enemyに１００＋のダメージを与えた' )
elif x > 0.2 :
print( '相手が躱した' )
else :
print( "まだゾーマは死んでません")
つかれたーｗ
もっと簡単なやつにしまーすｗ

**デフォルトの名無しさん** · 2024/10/09(水) 01:45:35.37

頑張って読んだけどこのままでは多分スレチ、ずっと条件を簡単にして出題し直してみる。

お題
勇者（HP1000）対魔王（HP5000）、先鋒勇者、サイコロ（目は１から６）を振って出た目だけ相手にダメージを与えられ、交互に攻撃しHP0となったら負け。魔王のサイコロは１個だけど勇者のサイコロはN個（N=1～10）の時の勇者の勝率を求めよ。

**デフォルトの名無しさん** · 2024/10/09(水) 02:10:34.29

# RPGの攻撃状況を再現

# モンスターと戦っている。
# 1から10のサイコロをふって、
# 6未満の場合、サイコロの目だけダメージを与えたと表示。
# 6以上の場合、クリティカルヒットとして、100のダメージを与えたと表示。
# 0の場合、防御されたとして、０ダメージを表示。

import random
hit = random.randint(0,10)
if 1 <= hit < 6:
print("スライムに、" + str(hit) + "のダメージを与えた！")
elif hit >= 6:
print("クリティカルヒットダメージ１００を与えた")
else:
print("ダメージは、" + str(hit) + "防御された")
#楽しすぎだなｗこれはｗ

**デフォルトの名無しさん** · 2024/10/09(水) 18:59:48.61

>>429
asciiで。

**デフォルトの名無しさん** · 2024/10/11(金) 05:53:27.29

お題：文字列の中で1つ以上連続している空白をすべて空白4つに置換せよ。ただし正規表現は使ってはならない。

例：
"abc def ghi" -> "abc def ghi"

**438** · 2024/10/11(金) 05:54:39.00

5chは空白を1つに置換するらしく例は出せない。脳内保管せよ。

**デフォルトの名無しさん** · 2024/10/11(金) 20:53:22.06

>>438
perl -lpE'$_=join$"x4,split$"'
(文字列を入れる)

**デフォルトの名無しさん** · 2024/10/11(金) 20:58:31.10

>>440
splitのあとの$"は不要だった。
どちらにしてもタブも置換してまうのでアウトかも。

**デフォルトの名無しさん** · 2024/10/11(金) 22:15:53.82

>>438
R
https://ideone.com/LO7vyR

**デフォルトの名無しさん** · 2024/10/12(土) 01:34:10.91

>>435
こういう場合の確率の計算はよく分からなかったので実際にそれぞれ一万回ずつ戦わせてみた。
サイコロ1個で100000戦0勝
サイコロ2個で100000戦0勝
サイコロ3個で100000戦0勝
サイコロ4個で100000戦0勝
サイコロ5個で100000戦47171勝
サイコロ6個で100000戦100000勝
サイコロ7個で100000戦100000勝
サイコロ8個で100000戦100000勝
サイコロ9個で100000戦100000勝
サイコロ10個で100000戦100000勝

**デフォルトの名無しさん** · 2024/10/12(土) 01:39:02.69

>>443
#!ruby
#encoding: utf-8
$stdout=open("kekka.txt","w")
i=100000 #対戦回数
k=0 #勇者の勝ち
n=1 #勇者のサイコロの数
hpm=5000 #魔王のＨＰ
hpy=1000 #勇者のＨＰ

**デフォルトの名無しさん** · 2024/10/12(土) 01:39:40.49

>>444
while n<=10
j,t=0,0
while j<i
hpy-=rand(1..6)
if hpy<=0
j+=1
t,hpy,hpm=0,1000,5000
next
end
m=0
while m<n
hpm-=rand(1..6)
m+=1
end
if hpm<=0
k+=1
j+=1
t,hpy,hpm=0,1000,5000
next
end
t+=1
end
puts "サイコロ#{n}個で#{i}戦#{k}勝"
k,hpy,hpm=0,1000,5000
n+=1
end

警備員[Lv.12] · 2024/11/05(火) 05:27:21.14

>>427
Kotlin

スペースが何個なのか分り易くするために '*' をスペースの代わりに入れて呼び出している。

https://paiza.io/projects/xSP0nP203HXTwqyxlQD0Gg

**デフォルトの名無しさん** · 2024/11/13(水) 22:01:15.62

>>435
R

f=\(N){
u=1000
m=5000
while(u>0 & m>0){
m=m-sum(sample(6,N,replace=TRUE))
u=u-sample(6,1)
}
u>m
}
sim=\(N) mean(replicate(1e3,f(N)))
p=sapply(1:10,sim)
plot(p)
data.frame(N=1:10,p=p)

N p
1 1 0.000
2 2 0.000
3 3 0.000
4 4 0.000
5 5 0.498
6 6 1.000
7 7 1.000
8 8 1.000
9 9 1.000
10 10 1.000

**デフォルトの名無しさん** · 2024/11/13(水) 22:22:31.42

>>435
Wolfram

f[N_] :=(
u=1000;
m=5000;
Until[u<=0 || m<=0, m-=Total@RandomChoice[Range[6],N];u-=Total@RandomChoice[Range[6],1]];
Boole[u >= m]
)
Table[Mean@Table[f[x],10^3],{x,1,10}] // N

{0., 0., 0., 0., 0.511, 1., 1., 1., 1., 1.}