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

**デフォルトの名無しさん** · 2019/11/17(日) 09:00:22.10

プログラミングのお題スレです。

【出題と回答例】
1 名前：デフォルトの名無しさん
　　お題：お題本文

2 名前：デフォルトの名無しさん
　　>>1 使用言語
　　回答本文
　　結果がある場合はそれも

【ソースコードが長くなったら】 (オンラインでコードを実行できる)
https://ideone.com/
http://codepad.org/
http://compileonline.com/
http://rextester.com/runcode
https://runnable.com/
https://code.hackerearth.com/
http://melpon.org/wandbox
https://paiza.io/

宿題は宿題スレがあるのでそちらへ。

※前スレ
プログラミングのお題スレ Part15
http://mevius.5ch.net/test/read.cgi/tech/1564310397/

**デフォルトの名無しさん** · 2019/12/09(月) 21:20:18.06

int prime[] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,91};

**デフォルトの名無しさん** · 2019/12/09(月) 21:22:25.12

91じゃなかった

◆QZaw55cn4c · 2019/12/09(月) 21:29:40.49

>>197
https://mevius.5ch.net/test/read.cgi/tech/1434079972/
2625個、でしょうか？

◆QZaw55cn4c · 2019/12/09(月) 21:35:00.48

やりなおします
>>197
https://mevius.5ch.net/test/read.cgi/tech/1434079972/62
2625個

**デフォルトの名無しさん** · 2019/12/09(月) 22:42:10.29

ふるいで5000以下の素数を求めて
2重ループで列挙

が速いかな

**デフォルトの名無しさん** · 2019/12/09(月) 22:48:18.91

>>201
正解。Rでエラトステネスの篩を使って求めるプログラムを一応貼っておく。
https://ideone.com/anqKGW

**デフォルトの名無しさん** · 2019/12/10(火) 06:55:28.37

>>197 J

(smoutput~ 2=#@q:)@>>:i.10000

4
6
9
10
...
9991
9993
9995
9997
9998

**デフォルトの名無しさん** · 2019/12/10(火) 07:06:32.14

問題の一区分である素数判定、並びに範囲内の素数列挙するコード
https://ideone.com/fGT115
~/bin/is_prime.exe 2 10000
とやれば1万までの素数が列挙され、
~/bin/is_prime.exe 2017
と1つ引数与えればそれだけ判定
引数無いとURLの用にOFする限界付近まで全部

1万までの素数出して、それパイプで処理したら楽かなと思った

**デフォルトの名無しさん** · 2019/12/10(火) 08:49:47.69

5000までの素数で十分だって言ってるのに

**デフォルトの名無しさん** · 2019/12/10(火) 09:32:39.88

>>197
お題：このお題の回答を論理式で表すとどうなるでしょうか。

**デフォルトの名無しさん** · 2019/12/10(火) 09:53:12.04

>>207
それはプログラミングのお題なのか？

**デフォルトの名無しさん** · 2019/12/10(火) 13:07:12.61

また数(ry

**デフォルトの名無しさん** · 2019/12/10(火) 13:10:37.69

出題者はいろんな言語の表記方法を知りたいだけか？
数学やアルゴリズム的には全然おもしろくないのばかり

**デフォルトの名無しさん** · 2019/12/10(火) 15:09:39.15

>>197
https://ideone.com/wVVsOF
C++。>>205 の素数判定パクりました。楽すぎ。
と、思ったらこれ、俺の回答間違ってる。

**デフォルトの名無しさん** · 2019/12/10(火) 15:15:40.70

>>197
https://ideone.com/1Vax7p

C++。こうかいな。

**デフォルトの名無しさん** · 2019/12/10(火) 16:37:43.07

>>210
お前が面白い問題出せばいいじゃん
たぶん誰も解かないだろうけどwww

**デフォルトの名無しさん** · 2019/12/10(火) 17:22:00.22

過去スレからお題引っ張りたいんだけど、有料会員じゃないからむりぽ

**デフォルトの名無しさん** · 2019/12/10(火) 17:45:48.68

2chscとかいうのが無料サルベージに向いていると聞いたことがある

**デフォルトの名無しさん** · 2019/12/10(火) 17:47:09.91

>>214
普通のブラウザで見ても出ないんだっけ？

**デフォルトの名無しさん** · 2019/12/10(火) 18:57:11.15

>>213
問題を出して人に解かせるのはあまり好きじゃない
解く方が好き

**デフォルトの名無しさん** · 2019/12/10(火) 19:01:43.67

過去の良問があればおしえろください

**デフォルトの名無しさん** · 2019/12/10(火) 19:30:15.36

>>216
見れた( ﾟДﾟ)

**デフォルトの名無しさん** · 2019/12/10(火) 20:07:19.37

>>219
そうか。普通のブラウザだとエロ広告が激しく付くからそれで過去スレ見せる料金なんとかしてるのかもね。

>>217
そんなあなたにとっておきのお題をひとつ。

お題: 面白いお題を作れ。

**デフォルトの名無しさん** · 2019/12/10(火) 20:25:27.92

じゃあ四角形を全部違う大きさの円で埋める

**デフォルトの名無しさん** · 2019/12/10(火) 21:05:07.05

>>207
答えは高々有限個の整数でしかないんだから論理式にはならない

**デフォルトの名無しさん** · 2019/12/10(火) 22:54:03.68

>>210
まあ、そんな所だね。このスレは競技プログラミングじゃないから、アルゴリズムや
パフォーマンスの追求よりは、各自が使う言語で楽な書き方ができるのを披露する方が多い。
>>203も可変長ベクトルへの再代入の繰り返しという非効率なことをやっているが、
自前のforループ不要で簡潔に書けるし、篩い落とす操作を忠実に表してもいる。
昔と違ってこれでも実用な速度で動くので、色々な書き方ができるようになった。

C#, Julia, PowerShellでも類似の書き方ができる（>>206の通り素数は5000までにした）。
C#　https://ideone.com/YWc3f8
Julia　https://ideone.com/8NdOe3
PowerShell　https://ideone.com/DRoOMU
R　https://ideone.com/34lnNQ

C#のstatic void Main()の中だけを抜き出したC#スクリプトも作って、それぞれの実行時間を
PowerShellの

measure-command {semiprime.exe}
measure-command {csi semiprime.csx}
measure-command {julia semiprime.jl}
measure-command {semiprime.ps1}
measure-command {rscript semiprime.R}

で計測すると、C#を1倍として、R 3.40倍、C#スクリプト 10.1倍、Julia 16.6倍、
PowerShell 199倍となった。PowerShellだけは延々と待たされ実用的でない。
他の4つはすぐに終わるが、右へ行くほどややもたつく。Rは遅いと言われているが、
スクリプト言語の中ではかなり速い。

**デフォルトの名無しさん** · 2019/12/10(火) 22:54:55.70

>>217
じゃあ、これ解いてみる?

整数x, y, z, kに関する方程式x^3 + y^3 + z^3 = kの解を、k = 1から100までの場合について求めよ。
https://engineer.fabcross.jp/archeive/190921_charity-engine.html

**デフォルトの名無しさん** · 2019/12/10(火) 23:25:26.63

k = 64, z = 4の時
任意の整数を+-反転した組が x,yの解であり、その個数は無限

**デフォルトの名無しさん** · 2019/12/10(火) 23:27:52.24

なぜ
k=1, z=1
じゃない？

**デフォルトの名無しさん** · 2019/12/10(火) 23:47:24.77

あとから追加されそうな条件の
仮に全部0以上の整数とした時に
5*5*5>125,100>4*4*4のメモ代り

**デフォルトの名無しさん** · 2019/12/10(火) 23:59:19.26

k=1から100のどれかに対してじゃなくて、
k=1から100のそれぞれすべてに対して求めるんじゃろ…

**デフォルトの名無しさん** · 2019/12/11(水) 00:11:02.49

https://ideone.com/WP5jkx

10000000以下で0.04秒
C++は速い！

**デフォルトの名無しさん** · 2019/12/11(水) 00:13:04.91

C#の10000以下より速い！

**デフォルトの名無しさん** · 2019/12/11(水) 00:30:48.20

>>224
k=1の時からいきなり難しいなあ

**デフォルトの名無しさん** · 2019/12/11(水) 00:40:32.37

>>224
ん？
解を全て求めるのではなく
各kに対して1個解を求めればいいの？

**デフォルトの名無しさん** · 2019/12/11(水) 09:11:06.45

>>197 seq, factor, awk

seq 10000 | factor | awk 'NF == 3'

**デフォルトの名無しさん** · 2019/12/11(水) 09:24:12.51

>>233
awkのとここれどういう意味？わたし女騎士だけど教えて！

**デフォルトの名無しさん** · 2019/12/11(水) 10:31:05.38

>>234
女騎士？

まあいいや。NFが3になる行だけ出力するんだよ。NFはフィールド数ね。
区切り文字がデフォルトのままだと空白文字で区切った時の個数。例えば行に a b c って入ってたら 3 になる。

**デフォルトの名無しさん** · 2019/12/11(水) 12:05:48.97

>>197
Kotlin
https://paiza.io/projects/6zHoXJXw1KuNLWZViVDYug

**デフォルトの名無しさん** · 2019/12/11(水) 13:38:08.03

>>233
素因数分解しちゃえばいいのか
サイコー

**デフォルトの名無しさん** · 2019/12/11(水) 13:42:56.98

>>233
やってみたら
--- Data stack:
って出力が10000行並ぶだけなんだけど…
なんかオプションいる？

**デフォルトの名無しさん** · 2019/12/11(水) 14:03:27.78

普通にできたけど
何のシェル使ってるの？

**デフォルトの名無しさん** · 2019/12/11(水) 14:35:31.47

bash。macで。

**デフォルトの名無しさん** · 2019/12/11(水) 14:38:38.09

俺もbashもだけど
seq 100くらいなら動くの？

**デフォルトの名無しさん** · 2019/12/11(水) 15:36:13.32

$ seq 3 | factor で止めてawk飛ばすと以下の出力です。
Factor 0.98 x86.64 (1886, heads/master-211d69561a, Jul 2 2018 17:46:19)
[Clang (GCC 4.2.1 Compatible Apple LLVM 7.3.0 (clang-703.0.29))] on macosx
IN: scratchpad
--- Data stack:
1
IN: scratchpad
--- Data stack:
1
2
IN: scratchpad
--- Data stack:
1
2
3
IN: scratchpad

**デフォルトの名無しさん** · 2019/12/11(水) 16:49:11.52

>>233
約数が3個な数を列挙ってこと？
4の約数は1,2,4の3個だけど
6の約数は1,2,3,6の4個だよ

**デフォルトの名無しさん** · 2019/12/11(水) 17:02:16.44

factorみりゃわかんだろ……
素因数の数だっつーの

**デフォルトの名無しさん** · 2019/12/11(水) 18:01:04.84

素因数の数が3？

**デフォルトの名無しさん** · 2019/12/11(水) 18:12:47.35

こんな具合だろ
seq 100 | ~/bin/factorization.pl| awk -F, "NF==2" | ~/bin/align.pl "=" 4 2,2
6 2,3
9 3,3
10 2,5
14 2,7
......
9995 5,1999
9997 13,769
9998 2,4999

**デフォルトの名無しさん** · 2019/12/11(水) 19:33:15.23

>>224
とりあえず7個を除いて出来た
ここからが長いのかな？

◆QZaw55cn4c · 2019/12/11(水) 19:46:53.70

>>243
そうそう、それに
8 の約数は　1, 2, 4, 8 の 4 個だけれども、これは >>197 のいう半素数ではないんですよね

>>233 は間違っていますね

**デフォルトの名無しさん** · 2019/12/11(水) 20:44:53.42

QZがそういうなら合ってるんじゃないの？

**デフォルトの名無しさん** · 2019/12/11(水) 20:51:31.74

反面教師？

**デフォルトの名無しさん** · 2019/12/11(水) 20:54:12.94

お題
1兆以下の半素数の個数を求めよ

これだとスクリプト系は無理かな？

**デフォルトの名無しさん** · 2019/12/11(水) 21:07:14.36

素数判定で書いた
long long int のC言語でも100億ぐらいを上限でサポートしてる
100億以上の判定は遅くてより優れた判定アルゴリズムが必要だからだ

推定:C系、優れた判定、現代スペックのPCが必要

**デフォルトの名無しさん** · 2019/12/11(水) 21:11:52.35

素数判定する必要は無いよ

◆QZaw55cn4c · 2019/12/11(水) 22:08:07.74

>>252
＞long long int のC言語でも100億ぐらいを上限でサポートしてる

極めて不正確ですよ
C99 後は uint64_t の最大値は 1844'6744'0737'0955'1615=約 1844 京ですね
それに実際に割ってみる方法以外には素数判定があるとでも？そんなものは存在しませんね

**デフォルトの名無しさん** · 2019/12/11(水) 22:15:54.42

何年も前に素数の数を数え時は
シングルスレッドのコードで
5000億までで1時間弱

だから>>251もそのくらいで求まるはず

**デフォルトの名無しさん** · 2019/12/11(水) 22:46:17.34

>>254
詳しそうなので>>251をよろしく！

**デフォルトの名無しさん** · 2019/12/11(水) 23:52:51.14

知り合いから、「子供の宿題なんだけど解けない」と言われた。

「正方形13枚と正三角形を何枚か組み合わせて、凸多面体を作るとき、その頂点数は A 個または B 個である。ただし一つの頂点には同じ枚数の面が集まっているものとし、A < B とする。」

途中経過を適宜表記しつつ解くプログラムを!

**デフォルトの名無しさん** · 2019/12/12(木) 00:28:15.32

A=19
B=38

**デフォルトの名無しさん** · 2019/12/12(木) 00:32:27.86

>>247
3乗して64bitに収まる範囲全て調べたけど
7個から減らなかった

**デフォルトの名無しさん** · 2019/12/12(木) 00:35:14.01

7個になるまで9秒
3乗して64bitに収まる範囲検索2031秒

**デフォルトの名無しさん** · 2019/12/12(木) 00:42:16.14

>>257
プログラミングの宿題じゃないはずなので
解き方は他の板で聞いて

**デフォルトの名無しさん** · 2019/12/12(木) 00:44:42.94

k=30, 33, 42, 52, 74, 75, 84
の時がわからん

**デフォルトの名無しさん** · 2019/12/12(木) 06:06:40.97

>>254
f(x) = 1 if x is prime
f(x) = 0 otherwise
なる f: R -> R が存在しないってなんで言い切れんの？
頭悪すぎ

**デフォルトの名無しさん** · 2019/12/12(木) 07:32:31.88

>>254も>>263も頭悪すぎだ

**デフォルトの名無しさん** · 2019/12/12(木) 07:40:34.14

>>252
素数判定は不要

>>254
uint64_tの上限じゃなくて素数判定可能な範囲
また、実際に割ってみる以外の方法も色々と存在する

>>263
fは存在する
あなたが書いた通りの定義そのままをfとすれば良い
fの存在と素数判定の方法が存在するかどうかは無関係

**デフォルトの名無しさん** · 2019/12/12(木) 07:48:51.75

>>265
お前も馬鹿だな
簡単な関数で f が定義できるならそれはそのまま簡単な計算で素数判定ができることに他ならない
こちらはそんな方法が存在するといっているわけではなく、無配慮にそんなものは存在しないと一刀両断する数学的リテラシーの無さを指摘しているだけなんだが

**デフォルトの名無しさん** · 2019/12/12(木) 08:00:21.80

>>263
｢簡単な関数でfが定義できる｣
非常に重要な言葉を後付けだし
これが最初からあったとしても、
単に簡単に素数判定可能可能を言い換えただけの
何の意味もないレス

リテラシーの無さの指摘の方法が頭が悪い

**デフォルトの名無しさん** · 2019/12/12(木) 08:08:21.73

現代の技術で言えば

特殊な形であれば高速に判定する方法はある
特殊な形でなければ高速に判定する方法は見つかっていない
実際に割ってみる以外の判別方法は存在する

**デフォルトの名無しさん** · 2019/12/12(木) 08:25:56.28

なんでNやZじゃなくてRにしたのか気になる

**デフォルトの名無しさん** · 2019/12/12(木) 12:56:52.27

文脈無視して揚げ足取ってて草
こいつ関数は全射しか存在しないとでも思ってんのか

**デフォルトの名無しさん** · 2019/12/12(木) 14:55:31.47

単純な疑問だよ
なんかこだわりとか深い意味とかあるのかと

｢実数が素である｣とか書かれると、
素数じゃない別の意味を指してると思うわけで

揚げ足取られたと思ったってことは
ミスったと思ったんだよね？

**デフォルトの名無しさん** · 2019/12/12(木) 15:35:29.90

いつもの連投アスペ君じゃん
相手すんなよ

**デフォルトの名無しさん** · 2019/12/12(木) 17:58:27.01

論破されて捨て台詞
の典型でした

**デフォルトの名無しさん** · 2019/12/12(木) 18:19:15.31

5chリテラシーが低いやつもだいぶ減ってはきてるが
このスレはちょいちょい出てくるよなぁ

**デフォルトの名無しさん** · 2019/12/12(木) 18:27:18.35

>>263を書いておいて良くいうよ

**デフォルトの名無しさん** · 2019/12/12(木) 18:39:51.55

論破を禁止します

**デフォルトの名無しさん** · 2019/12/12(木) 19:05:41.21

まーたロンパンジーが湧いてんのかよ

**デフォルトの名無しさん** · 2019/12/12(木) 19:12:35.43

リテラシー君
よほど悔しかったのかな？

◆QZaw55cn4c · 2019/12/12(木) 19:15:24.80

>>256
単に上限値を増やしただけのお題なんて糞だと思います

◆QZaw55cn4c · 2019/12/12(木) 19:17:30.13

>>265
＞また、実際に割ってみる以外の方法も色々と存在する
私の初等的な教科書にはウィルソンの定理の証明が載っていました、>>254 はちょっと勇み足でしたね

>>266
＞無配慮にそんなものは存在しないと一刀両断する数学的リテラシーの無さ
まあ、認めましょう

**デフォルトの名無しさん** · 2019/12/12(木) 19:18:20.83

>>279
そう思うってことは上限値以外の差がわからんて事だね
決定的な違いがある

**デフォルトの名無しさん** · 2019/12/12(木) 20:55:40.30

お題：

(A) 階乗の下z桁がすべて0となるような最小の自然数を、z = 1, 12, 123, 1234, 12345,
123456, 1234567, 12345678, 123456789の各場合について求めよ。存在しないときは
そう表示せよ。

(B) (A)の「下z桁がすべて0となる」という条件を「0である桁がz個となる」に置き換えて解け。

(A)はちょっと考えれば簡単に解ける。(B)はどうか知らない。多桁整数計算ができるJuliaでは
z = 12345までは力任せでも求められたが。

**デフォルトの名無しさん** · 2019/12/12(木) 21:21:19.99

5で割っていく感じだから
整数無視して無限に割ると
おおまかに4で割った感じの数0が並ぶ
126！は4で割って31ぐらい

**デフォルトの名無しさん** · 2019/12/12(木) 22:14:34.83

連投論破野郎が論破されてて草

**デフォルトの名無しさん** · 2019/12/12(木) 23:15:10.76

(A)だけ C++
https://ideone.com/LTC90E

**デフォルトの名無しさん** · 2019/12/12(木) 23:18:50.66

(B)は実際に下限から順番に計算していかないとダメじゃないか？

**デフォルトの名無しさん** · 2019/12/12(木) 23:52:01.98

n ! のゼロの個数は大雑把にはこのくらい
log_10(n !) / 10 + n * 9 / 40

逆算すれば大体の値はわかる
でも最小値だから正確に求めるのには使えない

**デフォルトの名無しさん** · 2019/12/13(金) 18:36:38.11

[お題]
今月(2019/12)の13日は金曜日。
開始年月と終了年月が与えられる。年は西暦、年月区切りは"/"。
開始年月一日から終了年月末日の間で"13日の金曜日"は何回あるか。
現在のグレゴリオ歴ルールで、終了年月末日まで変わらないとする。

制約: 1970 <= 年 <= 300億, 1 <= 月 <= 12, 開始年月 <= 終了年月

1)　2019/1 2019/12 --> 2　　 // 今年は9月と12月が該当
2)　1980/1 2012/12 --> 57
3)　213456/8 453212/4 --> 412380
4)　2000/1 2399/12 -->?
5)　2000/11 2800/9 -->?
6)　123456789/11 7301177081/10 -->?
7)　1970/1 30000000000/12 --> 51599996613

※数十万年程度だと力技だけで回っちゃうので、難易度調整で300億に。
　wikiの"グレゴリオ歴#暦法"　https://ja.wikipedia.org/wiki/%E3%82%B0%E3%83%AC%E3%82%B4%E3%83%AA%E3%82%AA%E6%9A%A6#%E6%9A%A6%E6%B3%95

**デフォルトの名無しさん** · 2019/12/13(金) 20:06:08.48

300億にしても年の曜日は14パターンしか無いのであっという間に解ける気が…
曜日は400年周期だから400年分求めればあとは割り算するだけかな？

**デフォルトの名無しさん** · 2019/12/13(金) 20:15:02.20

C++
https://ideone.com/33NYDY

**デフォルトの名無しさん** · 2019/12/13(金) 20:16:51.71

>>289
今日やらないとダメだよ

**デフォルトの名無しさん** · 2019/12/13(金) 20:32:32.43

>>288
ずるいやり方だが、取り敢えずPowerShellに
　1..12 |? {([DateTime]"2000/$_/13").DayOfWeek -eq "Friday"}
　1..12 |? {([DateTime]"2400/$_/13").DayOfWeek -eq "Friday"}
と打ち込んでみると、
　10
　10
と表示され、2000年と2400年の13日の金曜日は10月で一致する。
これと閏年調整のパターンが400年周期であることを組み合わせると、
13日の金曜日の出現パターンも400年周期になることが分かる。
だから、400年分の表を作ってしまえば、遅いPowerShellでも余裕で解ける。
https://ideone.com/oV7JFk

--- 実行結果 ---
2019/1 2019/12 -> 2
1980/1 2012/12 -> 57
213456/8 453212/4 -> 412380
2000/1 2399/12 -> 688
2000/11 2800/9 -> 1375
123456789/11 7301177081/10 -> 12345678901
1970/1 30000000000/12 -> 51599996613

**デフォルトの名無しさん** · 2019/12/13(金) 20:57:37.29

>>285
正解。例によってRで簡潔に書いた解答例: https://ideone.com/RhEtiC

>>287
やっぱり無理か。ちなみに、Juliaでz = 12345まで力任せで求めるプログラム
https://ideone.com/iqMYLa

--- 実行結果 ---
階乗の0である桁が1個になる最小値は5です。
階乗の0である桁が12個になる最小値は37です。
階乗の0である桁が123個になる最小値は299です。
階乗の0である桁が1234個になる最小値は2414です。
階乗の0である桁が12345個になる最小値は20177です。

**デフォルトの名無しさん** · 2019/12/13(金) 21:27:38.91

>>293

>>287の式に20177を入れると12349.7
かなりいい線行ってる

**デフォルトの名無しさん** · 2019/12/13(金) 22:22:01.90

123456個は存在しない

**デフォルトの名無しさん** · 2019/12/13(金) 23:21:15.30

123457桁になる 30484!
から
下位に0が123457個続く 493840!
までの間に存在しない事を示せば
存在しない証明になりますが
もっと範囲を絞った証明は出来ないですかねえ

確率的には>>287の付近だけ調べれば良いんですが...

**デフォルトの名無しさん** · 2019/12/13(金) 23:47:17.94

全部調べる感じでぶん回してるけど1234567までが長いお(´・ω・`)　(123456は>>295の言う通り出なかった)
1時間ぶん回して900000!まで行ったが…