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

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

[前50を表示]
150:デフォルトの名無しさん
19/12/07 00:44:47.39 7zlttp1T.net
>>101 Perl5
no warnings 'experimental::signatures';
use feature qw(say signatures);
sub Y($f) { sub {$f->(Y($f))->(@_)} } # Y combinator w/ recursive
$s = sub ($n) {$n and $n + $s->($n - 1)}; # calc series recursive
say Y(sub ($f) { $s })->(100);

151:デフォルトの名無しさん
19/12/07 01:33:07.92 7zlttp1T.net
>>101 Perl5
use List::Util 'sum';
print sum 1..100;

152:デフォルトの名無しさん
19/12/07 03:39:00.62 GrS1V5od.net
お題
1以上、20以下の整数がランダムに10個与えられます。
与えられた整数の集合から1つ以上の整数を抜き出し、残った整数の集合からまた1つ以上の整数を抜き出します。
この時2つの抜き出した整数の集合のそれぞれの和が同値であり、かつ最小のものを出力して下さい。

例
与えられた整数の集合 : 2 3 4 5 6 7 8 9 10 11
期待される出力 : 5
この場合、与えられた整数の集合から
「2,3」「5」の整数を抜き出しています。

153:デフォルトの名無しさん
19/12/07 03:47:40.06 GrS1V5od.net
元ネタ
ｽﾚﾘﾝｸ(news4vip板:18番)

154:デフォルトの名無しさん
19/12/07 04:00:27.84 G4M88S/c.net
・「同値」の意味が謎
・解がない時の処理が謎(e.g. input: 9 10 11 12 13 14 15 16 17 18
・同値の使い方が間違っているのでそもそも「集合」の元が重複しないものとしてよいものか悩ましい

155:デフォルトの名無しさん
19/12/07 04:23:11 GrS1V5od.net
>>152
同値は単に同じ(等しい)値という意味で使っています。

解が存在しない場合は「ないよー」と出力して下さい。

与えられる整数は重複する可能性があります。

156:デフォルトの名無しさん
19/12/07 05:16:02 HQTo5ewj.net
ならば結局こういうことでよいのかね.
もとのスレを見ても出題者本人の主張が不明瞭なうえに二転三転していて気持ち悪いが.

与えられた数列を {a_n} に対して,
{a_n} の異なる項からなる任意の部分列の内それぞれの和が等しくなるものを {b_n}, {c_n} として
Σb_n (= Σc_n) が最小となる {b_n}, {c_n} を求めよ.

そして今回は b, c の項数をそれぞれ 2, 1 に限るものとすると.

157:デフォルトの名無しさん
19/12/07 05:38:25.75 GrS1V5od.net
もう少し例を載せるべきでした。
すいません。
例えば
 >>152
さんの解が存在しないとしているものですが、
9 10 11 12 13 14 15 16 17 18
を与えられた場合の出力は
21となります。(10,11と9,12)

入力が
1 1 1 1 1 1 1 1 1 1
の場合の出力は
1となります。(1と1)

158:デフォルトの名無しさん
19/12/07 06:50:35.03 HQTo5ewj.net
>>155
なるほど、概ね理解した

159:デフォルトの名無しさん
19/12/07 09:57:16.56 WrheNqRo.net
>>150
取り敢えずRで力任せ。これでも瞬時に終わるので工夫の必要なし。
URLﾘﾝｸ(ideone.com)

160:デフォルトの名無しさん
19/12/07 15:33:28.08 WrheNqRo.net
実は >>157は「残った整数の集合から」の条件を忘れていて、元の集合から抜き出すと
勘違いして書いてしまったプログラム。
が、改めて条件を考えてみると、既に抜き出された数と同じものを選んでしまうのは、
「片方の部分集合の要素が2個以上で、もう片方の部分集合の要素が3個以上の場合」(A)
に限られる。例えば、2+9=11と2+3+6=11。この場合、2回選んでしまった2を取り除いた
部分集合は、和9=9と3+6=9が11より小さく、かつ(A)の場合に該当しないので2回選んで
しまった数は存在しない。
だから結局、>>157のプログラムのままで正解が得られることになる。

161:デフォルトの名無しさん
19/12/07 16:28:15.65 tj55yZgB.net
へなへななお題へなへなな回答

162:デフォルトの名無しさん
19/12/07 20:45:08.83 HU7sPj+p.net
>>150
この問題から関連して考えてたんだけど
5を
1 1 1 1 1
1 1 1 2
1 1 3
1 2 2
1 4
2 3
5
このように分割する方法は7通り、10だと42通りあった
数が与えられたとき何通りの分け方があるかどうやって計算するのか考えてたけど
そんな簡単な問題ではないようだった
とりあえず数が与えられたら上の例のように全ての分け方を列挙する
というのを問題としてみる

163:デフォルトの名無しさん
19/12/07 22:36:30.68 WrheNqRo.net
>>160
Rでお決まりの再帰呼び出し。
URLﾘﾝｸ(ideone.com)

164:デフォルトの名無しさん
19/12/07 23:22:48 eT8T+vHJ.net
分割数でググれば考え方いっぱい出てくるよ

165:デフォルトの名無しさん
19/12/08 13:36:47.70 jvur7pXC.net
>>160
C++
URLﾘﾝｸ(ideone.com)

166:デフォルトの名無しさん
19/12/08 13:57:09.43 FOSx0Jk/.net
>>154
最後の文はどこから出てきたの？

167:デフォルトの名無しさん
19/12/08 13:58:38.12 xElyalHo.net
>>164
元のスレ

168:デフォルトの名無しさん
19/12/08 14:21:04.90 FOSx0Jk/.net
しらんがな

169:デフォルトの名無しさん
19/12/08 14:22:15.55 FOSx0Jk/.net
その制限が無い方がアルゴリズムとしておもしろい
制限があるとつまらなすぎる

170:デフォルトの名無しさん
19/12/08 17:19:46.37 xElyalHo.net
知らんがなって
なら聞くなアホ

171:デフォルトの名無しさん
19/12/08 17:53:26.22 FOSx0Jk/.net
すまん
不満は >>150に言ったつもり
そんな重要な事を省略すんなって

172:デフォルトの名無しさん
19/12/08 19:17:22.87 DgjgjjxW.net
別に項数の制限はつけていませんよ
「項数の制限が無い」という事を省略したのに怒っているのならすいません。
ですが制限があるのなら普通に問題文に加えますし、別に書くほどの事では無いかな～と

173:デフォルトの名無しさん
19/12/08 19:37:55.14 KCeBLlvA.net
>>150
URLﾘﾝｸ(ideone.com)
C++。総当たりです。スカイレークのi7で12秒くらいかかります。
久しぶりにまじめに総当たりを書いた気がしました。

174:デフォルトの名無しさん
19/12/08 19:44:34.35 KCeBLlvA.net
>>150
URLﾘﾝｸ(ideone.com)
オマケで、答えが見える版を置いておきます。C++。

175:デフォルトの名無しさん
19/12/08 20:22:49.16 KCeBLlvA.net
ちなみにオーダーは大体O(N!)位です。（笑

176:デフォルトの名無しさん
19/12/08 20:35:11.71 KCeBLlvA.net
ギャグですけど、並列化は比較的簡単なのでそれで時間短縮はできます。
底の値をシェアードにすると早く終わります。Nになってると思うんだけど。

177:デフォルトの名無しさん
19/12/08 20:48:51.01 KCeBLlvA.net
一回を関数に切り出して実行した場合、
一回のイテレーションが大体100回のループに収まるはずなのでザクザクおわります。
多分。

178:デフォルトの名無しさん
19/12/08 20:49:16.37 KCeBLlvA.net
一回を関数に切り出して実行した場合、
一回のイテレーションが大体100回のループに収まるはずなのでザクザクおわります。
多分。

179:デフォルトの名無しさん
19/12/08 20:59:32.19 KCeBLlvA.net
ぐあ、重複した・・・。

180:デフォルトの名無しさん
19/12/08 21:17:04.79 FKbRmDMb.net
>>150
これは問題の設定がつまらないな。1～20の中から10個を選んで元の集合を作るから、
結果に1個か2個の和しかほとんど現れず、集合の最初の方をパッと見ただけで
暗算でも分かってしまう。1～5000の中から10個を選ぶ設定にすると、
結果がなしだったり、3個の和と

181:4個の和だったり、2個の和と6個の和だったり、変化に富んで面白くなる。https://ideone.com/kmihyT 例えば、リンク先にある入力: [63, 70, 269, 949, 1337, 2670, 3538, 3764, 4183, 4320] 出力: Σ[3764, 4183] = Σ[63, 70, 269, 1337, 2670, 3538] = 7947 なんてパッと見では思いつかないから、コンピュータに探させる意義がある。

182:デフォルトの名無しさん
19/12/08 21:27:36.39 KCeBLlvA.net
>>178
異様に早いなーと思ったら、言語にコンビネーションあるんかいな。
裏山シー。

183:デフォルトの名無しさん
19/12/08 22:01:27.05 h14g0YSH.net
サンプルだから人間が簡単に検証できるようにしてるんでしょ
普通それぐらいはわかりそうなもんだけど、>>165みたいに項数だと思う奴とか >>178みたいにイチャモンつける奴とか世の中広いわw

184:デフォルトの名無しさん
19/12/08 22:35:33.47 FKbRmDMb.net
>>180
そんなことは分かっているよ。だから、お題通りの1～20の場合も >>157でちゃんと回答した。
その上で、もっと面白い場合の追加を提案してみただけ。

185:デフォルトの名無しさん
19/12/09 00:06:03 QbXWD96q.net
>>150
N!より速い方法ある？

186:デフォルトの名無しさん
19/12/09 00:53:44 rq2SBWAq.net
>>182
動的計画法？

187:デフォルトの名無しさん
19/12/09 01:15:42 2eMu76Ef.net
URLﾘﾝｸ(ideone.com)
全ての和を計算して並べ替えるだけ
多分最も愚直な方法

188:デフォルトの名無しさん
19/12/09 01:57:50 2eMu76Ef.net
bit演算で面倒なことやってたけどpairっての使えば良かったのか

189:デフォルトの名無しさん
19/12/09 02:04:04.39 vzskLW//.net
>>150　URLﾘﾝｸ(ideone.com)
By PyPy、　ノーマルpythonでは力業の(N=20)が8秒くらいかな、
力業が　2^N * N
最小値だけなら、N*数列の合計 = N^2 * 数の平均(/2) ででる(みたい？)
(自信ががないDP解)

190:デフォルトの名無しさん
19/12/09 02:51:44 ElWitvQQ.net
>>180
日本語が読め無い馬鹿発見

191:デフォルトの名無しさん
19/12/09 04:54:54.63 wE9bCkNR.net
>>181
わかってたら
> これは問題の設定がつまらないな。
なんていう物言いにはならんだろ
 >>187
夜中まで必死だな…
何に必死なのかよくわからんけどw

192:デフォルトの名無しさん
19/12/09 05:00:22.33 ElWitvQQ.net
必死なのはお前だろ
お前一人だけ日本語すらまともに読めてない馬鹿だって気づけよ

193:デフォルトの名無しさん
19/12/09 06:12:33.29 wE9bCkNR.net
うわっ、アホが無駄に絡んできたよw

194:デフォルトの名無しさん
19/12/09 06:42:21 QCNDYaVq.net
明け方からどんだけ必死なんだよ
以降、劣等感の塊のID:wE9bCkNRくんが全レスしてくれるってよ！

195:デフォルトの名無しさん
19/12/09 06:46:34 PLlkWb6P.net
こいつ少し上の方でレスバしてたアホやろ？
さんざん馬鹿にされて悔しい思いしたから早朝にちょろっと顔出してるんやろ

196:デフォルトの名無しさん
19/12/09 07:23:51 RwnUxfkW.net
単芝ガイジ君、情けなさ過ぎて草

197:デフォルトの名無しさん
19/12/09 12:30:02.62 G+LM1RHL.net
>>192
自己紹介乙ww

198:デフォルトの名無しさん
19/12/09 15:19:21 gONUrOAf.net
URLﾘﾝｸ(ideone.com)
C++面白いな

199:デフォルトの名無しさん
19/12/09 15:48:30 gONUrOAf.net
(sum[i].second & sum[i + 1].second) == 0
この比較はいらないのかな
これが重なってるならより小さい重なってない組合せが必ず存在するか

200:デフォルトの名無しさん
19/12/09 21:07:53.73 l5WymCFL.net
お題：2つの素数（2つは同じでもよい）の積で表される数は半素数と呼ばれる。
1万以下の半素数をすべて表示せよ。

201:デフォルトの名無しさん
19/12/09 21:20:18.06 QbXWD96q.net
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};

202:デフォルトの名無しさん
19/12/09 21:22:25.12 QbXWD96q.net
91じゃなかった

203:
19/12/09 21:29:40.49 ZryjKmS+.net
>>197
ｽﾚﾘﾝｸ(tech板)
2625個、でしょうか？

204:
19/12/09 21:35:00.48 ZryjKmS+.net
やりなおします
 >>197
ｽﾚﾘﾝｸ(tech板:62番)
2625個

205:デフォルトの名無しさん
19/12/09 22:42:10.29 QbXWD96q.net
ふるいで5000以下の素数を求めて
2重ループで列挙
が速いかな

206:デフォルトの名無しさん
19/12/09 22:48:18.91 l5WymCFL.net
>>201
正解。Rでエラトステネスの篩を使って求めるプログラムを一応貼っておく。
URLﾘﾝｸ(ideone.com)

207:デフォルトの名無しさん
19/12/10 06:55:28 cIwr+d9F.net
>>197 J

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

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

208:デフォルトの名無しさん
19/12/10 07:06:32 qBy9puuu.net
問題の一区分である素数判定、並びに範囲内の素数列挙するコード
URLﾘﾝｸ(ideone.com)
~/bin/is_prime.exe 2 10000
とやれば1万までの素数が列挙され、
~/bin/is_prime.exe 2017
と1つ引数与えればそれだけ判定
引数無いとURLの用にOFする限界付近まで全部

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

209:デフォルトの名無しさん
19/12/10 08:49:47 92MPgAr5.net
5000までの素数で十分だって言ってるのに

210:デフォルトの名無しさん
19/12/10 09:32:39 WOcT9SPT.net
>>197
お題：このお題の回答を論理式で表すとどうなるでしょうか。

211:デフォルトの名無しさん
19/12/10 09:53:12.04 gKYhlG5V.net
>>207
それはプログラミングのお題なのか？

212:デフォルトの名無しさん
19/12/10 13:07:12.61 bINIS1ks.net
また数(ry

213:デフォルトの名無しさん
19/12/10 13:10:37.69 FDwwVytW.net
出題者はいろんな言語の表記方法を知りたいだけか？
数学やアルゴリズム的には全然おもしろくないのばかり

214:デフォルトの名無しさん
19/12/10 15:09:39 zIz8I18p.net
>>197
URLﾘﾝｸ(ideone.com)
C++。>>205 の素数判定パクりました。楽すぎ。
と、思ったらこれ、俺の回答間違ってる。

215:デフォルトの名無しさん
19/12/10 15:15:40 zIz8I18p.net
>>197
URLﾘﾝｸ(ideone.com)

C++。こうかいな。

216:デフォルトの名無しさん
19/12/10 16:37:43 hI+yeapE.net
>>210
お前が面白い問題出せばいいじゃん
たぶん誰も解かないだろうけどwww

217:デフォルトの名無しさん
19/12/10 17:22:00.22 Ajx0JUvY.net
過去スレからお題引っ張りたいんだけど、有料会員じゃないからむりぽ

218:デフォルトの名無しさん
19/12/10 17:45:48.68 qBy9puuu.net
2chscとかいうのが無料サルベージに向いていると聞いたことがある

219:デフォルトの名無しさん
19/12/10 17:47:09.91 ClyY78bX.net
>>214
普通のブラウザで見ても出ないんだっけ？

220:デフォルトの名無しさん
19/12/10 18:57:11.15 W3sLZ8lM.net
>>213
問題を出して人に解かせるのはあまり好きじゃない
解く方が好き

221:デフォルトの名無しさん
19/12/10 19:01:43.67 W3sLZ8lM.net
過去の良問があればおしえろください

222:デフォルトの名無しさん
19/12/10 19:30:15.36 Ajx0JUvY.net
>>216
見れた( ﾟДﾟ)

223:デフォルトの名無しさん
19/12/10 20:07:19.37 ClyY78bX.net
>>219
そうか。普通のブラウザだとエロ広告が激しく付くからそれで過去スレ見せる料金なんとかしてるのかもね。
>>217
そんなあなたにとっておきのお題をひとつ。
お題: 面白いお題を作れ。

224:デフォルトの名無しさん
19/12/10 20:25:27.92 6QYDHDQi.net
じゃあ四角形を全部違う大きさの円で埋める

225:デフォルトの名無しさん
19/12/10 21:05:07.05 0RQ6ozIG.net
>>207
答えは高々有限個の整数でしかないんだから論理式にはならない

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

C#, Julia, PowerShellでも類似の書き方ができる（>>206の通り素数は5000までにした）。
C#　URLﾘﾝｸ(ideone.com)
Julia　URLﾘﾝｸ(ideone.com)
PowerShell　URLﾘﾝｸ(ideone.com)
R　URLﾘﾝｸ(ideone.com)

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は遅いと言われているが、
スクリプト言語の中ではかなり速い。

227:デフォルトの名無しさん
19/12/10 22:54:55 ZImsJVBi.net
>>217
じゃあ、これ解いてみる?

整数x, y, z, kに関する方程式x^3 + y^3 + z^3 = kの解を、k = 1から100までの場合について求めよ。
URLﾘﾝｸ(engineer.fabcross.jp)

228:デフォルトの名無しさん
19/12/10 23:25:26 qBy9puuu.net
k = 64, z = 4の時
任意の整数を+-反転した組が x,yの解であり、その個数は無限

229:デフォルトの名無しさん
19/12/10 23:27:52 W3sLZ8lM.net
なぜ
k=1, z=1
じゃない？

230:デフォルトの名無しさん
19/12/10 23:47:24 qBy9puuu.net
あとから追加されそうな条件の
仮に全部0以上の整数とした時に
5*5*5>125,100>4*4*4のメモ代り

231:デフォルトの名無しさん
19/12/10 23:59:19 RjwvfByt.net
k=1から100のどれかに対してじゃなくて、
k=1から100のそれぞれすべてに対して求めるんじゃろ…

232:デフォルトの名無しさん
19/12/11 00:11:02 10jfhd7e.net
URLﾘﾝｸ(ideone.com)

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

233:デフォルトの名無しさん
19/12/11 00:13:04 10jfhd7e.net
C#の10000以下より速い！

234:デフォルトの名無しさん
19/12/11 00:30:48 10jfhd7e.net
>>224
k=1の時からいきなり難しいなあ

235:デフォルトの名無しさん
19/12/11 00:40:32 10jfhd7e.net
>>224
ん？
解を全て求めるのではなく
各kに対して1個解を求めればいいの？

236:デフォルトの名無しさん
19/12/11 09:11:06 aadkbL3F.net
>>197 seq, factor, awk

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

237:デフォルトの名無しさん
19/12/11 09:24:12 ztpKOEip.net
>>233
awkのとここれどういう意味？わたし女騎士だけど教えて！

238:デフォルトの名無しさん
19/12/11 10:31:05 dG8VWZ74.net
>>234
女騎士？

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

239:デフォルトの名無しさん
19/12/11 12:05:48.97 dG8VWZ74.net
>>197
Kotlin
URLﾘﾝｸ(paiza.io)

240:デフォルトの名無しさん
19/12/11 13:38:08.03 ivhCTlPt.net
>>233
素因数分解しちゃえばいいのか
サイコー

241:デフォルトの名無しさん
19/12/11 13:42:56.98 QbvBtpFM.net
>>233
やってみたら
--- Data stack:
って出力が10000行並ぶだけなんだけど…
なんかオプションいる？

242:デフォルトの名無しさん
19/12/11 14:03:27.78 jagg9gKF.net
普通にできたけど
何のシェル使ってるの？

243:デフォルトの名無しさん
19/12/11 14:35:31 QbvBtpFM.net
bash。macで。

244:デフォルトの名無しさん
19/12/11 14:38:38 jagg9gKF.net
俺もbashもだけど
seq 100くらいなら動くの？

245:デフォルトの名無しさん
19/12/11 15:36:13.32 QbvBtpFM.net
$ 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

246:デフォルトの名無しさん
19/12/11 16:49:11 1E1+DBtw.net
>>233
約数が3個な数を列挙ってこと？
4の約数は1,2,4の3個だけど
6の約数は1,2,3,6の4個だよ

247:デフォルトの名無しさん
19/12/11 17:02:16.44 hUZCfnLs.net
factorみりゃわかんだろ……
素因数の数だっつーの

248:デフォルトの名無しさん
19/12/11 18:01:04.84 10jfhd7e.net
素因数の数が3？

249:デフォルトの名無しさん
19/12/11 18:12:47.35 ivhCTlPt.net
こんな具合だろ
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

250:デフォルトの名無しさん
19/12/11 19:33:15.23 ReYSFEXH.net
>>224
とりあえず7個を除いて出来た
ここからが長いのかな？

251:
19/12/11 19:46:53.70 6E3wj7zP.net
>>243
そうそう、それに
8 の約数は　1, 2, 4, 8 の 4 個だけれども、これは >>197 のいう半素数ではないんですよね
 >>233 は間違っていますね

252:デフォルトの名無しさん
19/12/11 20:44:53.42 WjX/lCwK.net
QZがそういうなら合ってるんじゃないの？

253:デフォルトの名無しさん
19/12/11 20:51:31.74 7CYZ1E2N.net
反面教師？

254:デフォルトの名無しさん
19/12/11 20:54:12.94 7CYZ1E2N.net
お題
1兆以下の半素数の個数を求めよ
これだとスクリプト系は無理かな？

255:デフォルトの名無しさん
19/12/11 21:07:14.36 ivhCTlPt.net
素数判定で書いた
long long int のC言語でも100億ぐらいを上限でサポートしてる
100億以上の判定は遅くてより優れた判定アルゴリズムが必要だからだ
推定:C系、優れた判定、現代スペックのPCが必要

256:デフォルトの名無しさん
19/12/11 21:11:52.35 7CYZ1E2N.net
素数判定する必要は無いよ

257:
19/12/11 22:08:07.74 6E3wj7zP.net
>>252
＞long long int のC言語でも100億ぐらいを上限でサポートしてる
極めて不正確ですよ
C99 後は uint64_t の最大値は 1844'6744'0737'0955'1615=約 1844 京ですね
それに実際に割ってみる方法以外には素数判定があるとでも？そんなものは存在しませんね

258:デフォルトの名無しさん
19/12/11 22:15:54.42 10jfhd7e.net
何年も前に素数の数を数え時は
シングルスレッドのコードで
5000億までで1時間弱
だから >>251もそのくらいで求まるはず

259:デフォルトの名無しさん
19/12/11 22:46:17.34 10jfhd7e.net
>>254
詳しそうなので >>251をよろしく！

260:デフォルトの名無しさん
19/12/11 23:52:51.14 VagVWTF0.net
知り合いから、「子供の宿題なんだけど解けない」と言われた。
「正方形13枚と正三角形を何枚か組み

261:合わせて、凸多面体を作るとき、その頂点数は A 個または B 個である。ただし一つの頂点には同じ枚数の面が集まっているものとし、A < B とする。」途中経過を適宜表記しつつ解くプログラムを!

262:デフォルトの名無しさん
19/12/12 00:28:15.32 UvCIGzud.net
A=19
B=38

263:デフォルトの名無しさん
19/12/12 00:32:27.86 UvCIGzud.net
>>247
3乗して64bitに収まる範囲全て調べたけど
7個から減らなかった

264:デフォルトの名無しさん
19/12/12 00:35:14.01 UvCIGzud.net
7個になるまで9秒
3乗して64bitに収まる範囲検索2031秒

265:デフォルトの名無しさん
19/12/12 00:42:16.14 UvCIGzud.net
>>257
プログラミングの宿題じゃないはずなので
解き方は他の板で聞いて

266:デフォルトの名無しさん
19/12/12 00:44:42.94 UvCIGzud.net
k=30, 33, 42, 52, 74, 75, 84
の時がわからん

267:デフォルトの名無しさん
19/12/12 06:06:40.97 wmbGd8tm.net
>>254
f(x) = 1 if x is prime
f(x) = 0 otherwise
なる f: R -> R が存在しないってなんで言い切れんの？
頭悪すぎ

268:デフォルトの名無しさん
19/12/12 07:32:31.88 UvCIGzud.net
>>254も >>263も頭悪すぎだ

269:デフォルトの名無しさん
19/12/12 07:40:34.14 UvCIGzud.net
>>252
素数判定は不要
 >>254
uint64_tの上限じゃなくて素数判定可能な範囲
また、実際に割ってみる以外の方法も色々と存在する
 >>263
fは存在する
あなたが書いた通りの定義そのままをfとすれば良い
fの存在と素数判定の方法が存在するかどうかは無関係

270:デフォルトの名無しさん
19/12/12 07:48:51.75 wmbGd8tm.net
>>265
お前も馬鹿だな
簡単な関数で f が定義できるならそれはそのまま簡単な計算で素数判定ができることに他ならない
こちらはそんな方法が存在するといっているわけではなく、無配慮にそんなものは存在しないと一刀両断する数学的リテラシーの無さを指摘しているだけなんだが

271:デフォルトの名無しさん
19/12/12 08:00:21.80 UvCIGzud.net
>>263
｢簡単な関数でfが定義できる｣
非常に重要な言葉を後付けだし
これが最初からあったとしても、
単に簡単に素数判定可能可能を言い換えただけの
何の意味もないレス
リテラシーの無さの指摘の方法が頭が悪い

272:デフォルトの名無しさん
19/12/12 08:08:21.73 UvCIGzud.net
現代の技術で言えば
特殊な形であれば高速に判定する方法はある
特殊な形でなければ高速に判定する方法は見つかっていない
実際に割ってみる以外の判別方法は存在する

273:デフォルトの名無しさん
19/12/12 08:25:56.28 UvCIGzud.net
なんでNやZじゃなくてRにしたのか気になる

274:デフォルトの名無しさん
19/12/12 12:56:52.27 wmbGd8tm.net
文脈無視して揚げ足取ってて草
こいつ関数は全射しか存在しないとでも思ってんのか

275:デフォルトの名無しさん
19/12/12 14:55:31.47 2PtCb74v.net
単純な疑問だよ
なんかこだわりとか深い意味とかあるのかと
｢実数が素である｣とか書かれると、
素数じゃない別の意味を指してると思うわけで
揚げ足取られたと思ったってことは
ミスったと思ったんだよね？

276:デフォルトの名無しさん
19/12/12 15:35:29.90 zRSUYply.net
いつもの連投アスペ君じゃん
相手すんなよ

277:デフォルトの名無しさん
19/12/12 17:58:27.01 2C/x7O5N.net
論破されて捨て台詞
の典型でした

278:デフォルトの名無しさん
19/12/12 18:19:15.31 RcQND7/k.net
5chリテラシーが低いやつもだいぶ減ってはきてるが
このスレはちょいちょい出てくるよなぁ

279:デフォルトの名無しさん
19/12/12 18:27:18.35 cffvEFbc.net
>>263を書いておいて良くいうよ

280:デフォルトの名無しさん
19/12/12 18:39:51.55 JPTMcftT.net
論破を禁止します

281:デフォルトの名無しさん
19/12/12 19:05:41.21 f2lx0DlC.net
まーたロンパンジーが湧いてんのかよ

282:デフォルトの名無しさん
19/12/12 19:12:35.43 cffvEFbc.net
リテラシー君
よほど悔しかったのかな？

283:
19/12/12 19:15:24.80 bwSkQtN1.net
>>256
単に上限値を増やしただけのお題なんて糞だと思います

284:
19/12/12 19:17:30.13 bwSkQtN1.net
>>265
＞また、実際に割ってみる以外の方法も色々と存在する
私の初等的な教科書にはウィルソンの定理の証明が載っていました、>>254 はちょっと勇み足でしたね
 >>266
＞無配慮にそんなものは存在しないと一刀両断する数学的リテラシーの無さ
まあ、認めましょう

285:デフォルトの名無しさん
19/12/12 19:18:20.83 cffvEFbc.net
>>279
そう思うってことは上限値以外の差がわからんて事だね
決定的な違いがある

286:デフォルトの名無しさん
19/12/12 20:55:40.30 uIjgN9ep.net
お題：
(A) 階乗の下z桁がすべて0となるような最小の自然数を、z = 1, 12, 123, 1234, 12345,
123456, 1234567, 12345678, 123456789の各場合について求めよ。存在しないときは
そう表示せよ。
(B) (A)の「下z桁がすべて0となる」という条件を「0である桁がz個となる」に置き換えて解け。

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

287:デフォルトの名無しさん
19/12/12 21:21:19.99 qa03h3dy.net
5で割っていく感じだから
整数無視して無限に割ると
おおまかに4で割った感じの数0が並ぶ
126！は4で割って31ぐらい

288:デフォルトの名無しさん
19/12/12 22:14:34.83 WrPVO4K2.net
連投論破野郎が論破されてて草

289:デフォルトの名無しさん
19/12/12 23:15:10.76 UvCIGzud.net
(A)だけ C++
URLﾘﾝｸ(ideone.com)

290:デフォルトの名無しさん
19/12/12 23:18:50.66 UvCIGzud.net
(B)は実際に下限から順番に計算していかないとダメじゃないか？

291:デフォルトの名無しさん
19/12/12 23:52:01.98 UvCIGzud.net
n ! のゼロの個数は大雑把にはこのくらい
log_10(n !) / 10 + n * 9 / 40
逆算すれば大体の値はわかる
でも最小値だから正確に求めるのには使えない

292:デフォルトの名無しさん
19/12/13 18:36:38.11 0Lb/VhZx.net
[お題]
今月(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の"グレゴリオ歴#暦法"　URLﾘﾝｸ(ja.wikipedia.org)

293:デフォルトの名無しさん
19/12/13 20:06:08.48 KAf60mjk.net
300億にしても年の曜日は14パターンしか無いのであっという間に解ける気が…
曜日は400年周期だから400年分求めればあとは割り算するだけかな？

294:デフォルトの名無しさん
19/12/13 20:15:02.20 QuzLYYNu.net
C++
URLﾘﾝｸ(ideone.com)

295:デフォルトの名無しさん
19/12/13 20:16:51.71 QuzLYYNu.net
>>289
今日や

296:らないとダメだよ

297:デフォルトの名無しさん
19/12/13 20:32:32.43 7JZ8Vcg3.net
>>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でも余裕で解ける。
URLﾘﾝｸ(ideone.com)
--- 実行結果 ---
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

298:デフォルトの名無しさん
19/12/13 20:57:37.29 7JZ8Vcg3.net
>>285
正解。例によってRで簡潔に書いた解答例: URLﾘﾝｸ(ideone.com)
>>287
やっぱり無理か。ちなみに、Juliaでz = 12345まで力任せで求めるプログラム
URLﾘﾝｸ(ideone.com)
--- 実行結果 ---
階乗の0である桁が1個になる最小値は5です。
階乗の0である桁が12個になる最小値は37です。
階乗の0である桁が123個になる最小値は299です。
階乗の0である桁が1234個になる最小値は2414です。
階乗の0である桁が12345個になる最小値は20177です。

299:デフォルトの名無しさん
19/12/13 21:27:38.91 QuzLYYNu.net
>>293
>>287の式に20177を入れると12349.7
かなりいい線行ってる

300:デフォルトの名無しさん
19/12/13 22:22:01.90 QuzLYYNu.net
123456個は存在しない

301:デフォルトの名無しさん
19/12/13 23:21:15.30 QuzLYYNu.net
123457桁になる 30484!
から
下位に0が123457個続く 493840!
までの間に存在しない事を示せば
存在しない証明になりますが
もっと範囲を絞った証明は出来ないですかねえ
確率的には >>287の付近だけ調べれば良いんですが...

302:デフォルトの名無しさん
19/12/13 23:47:17.94 r78PztnN.net
全部調べる感じでぶん回してるけど1234567までが長いお(´・ω・`)　(123456は >>295の言う通り出なかった)
1時間ぶん回して900000!まで行ったが…

303:デフォルトの名無しさん
19/12/13 23:49:02.36 r2S5lo/C.net
数学でワープできる分手計算のほうが早い

304:デフォルトの名無しさん
19/12/13 23:49:38.96 QuzLYYNu.net
kwsk

305:デフォルトの名無しさん
19/12/13 23:50:22.63 r2S5lo/C.net
>283

306:デフォルトの名無しさん
19/12/13 23:52:09.33 QuzLYYNu.net
(A)は解決済み
(B)の話

307:デフォルトの名無しさん
19/12/13 23:53:42.34 QuzLYYNu.net
数学で絞れるのは >>296くらいかと

308:デフォルトの名無しさん
19/12/13 23:56:32.85 QuzLYYNu.net
下限の n! を計算するのは速い
そこから順番に1個ずつnを増やしていって
0を数えていかないとならない
見つからなければ上限までずっと

309:デフォルトの名無しさん
19/12/13 23:59:51.46 QuzLYYNu.net
ちなみに私、東大の数学科卒なんで

310:デフォルトの名無しさん
19/12/14 00:07:59.08 ruM/1CH5.net
ワープ
(* 4 1234567) 4938268
(/ 4938268 5.0) 987653.6
(/ 987653 5.0) 197530.6
(/ 197530 5.0) 39506.0
(/ 39506 5.0) 7901.2
(/ 7901 5.0) 1580.2
(/ 1580 5.0) 316.0
(/ 316 5.0) 63.2
(/ 63 5.0) 12.6
(/ 12 5.0) 2.4
(+ 987653 197530 39506 7901 1580 316 63 12 2) 1234563

311:デフォルトの名無しさん
19/12/14 00:14:15.42 03SKvnKb.net
>>305
だからそれは >>282の(A)の話じゃろ…

312:デフォルトの名無しさん
19/12/14 00:15:48.22 ruM/1CH5.net
Bの話か

313:デフォルトの名無しさん
19/12/14 02:58:36 03SKvnKb.net
階乗の0である桁が1234567個になる最小値は1542987です。
お、出た…4時間か

314:デフォルトの名無しさん
19/12/14 03:55:22 DhSlagHw.net
お題: 目覚まし時計の現在時刻とアラーム時刻が入力されるので、待機時間を出力しなさい

例:
現在時刻 = 10:00
アラーム時刻 = 16:00
待機時間 = 06:00

現在時刻 = 22:00
アラーム時刻 = 04:00
待機時間 = 06:00

315:デフォルトの名無しさん
19/12/14 09:12:53.23 ruM/1CH5.net
>>309
これだな
URLﾘﾝｸ(paiza.io)
~/bin/sleep.pl

316:10:00 2909 00:48:29 のように使ってる

317:デフォルトの名無しさん
19/12/14 09:59:37.97 Lp9CT80F.net
>>304
進振り何点だった？　連投アスペ君

318:デフォルトの名無しさん
19/12/14 22:12:52.90 8CJijD11.net
>>309
この手の処理はPowerShellが得意だな。
function WaitingTime($now, $alarm)
{
　　$wait = [DateTime]$alarm - [DateTime]$now
　　if ($wait -lt 0) {$wait += [TimeSpan]"12:00" + [TimeSpan]"12:00"}
　　"現在時刻 = $now"
　　"アラーム時刻 = $alarm"
　　"待機時刻 = $wait"
　　""
}
WaitingTime 10:00 16:00
WaitingTime 22:00 4:00
WaitingTime 12:34:56 1:23:45
--- 実行結果 ---
現在時刻 = 10:00
アラーム時刻 = 16:00
待機時刻 = 06:00:00
現在時刻 = 22:00
アラーム時刻 = 4:00
待機時刻 = 06:00:00
現在時刻 = 12:34:56
アラーム時刻 = 1:23:45
待機時刻 = 12:48:49

319:デフォルトの名無しさん
19/12/15 00:16:05 K3IRPfAN.net
>>312
キャスト出来るんだな
楽そう

320:デフォルトの名無しさん
19/12/15 06:09:40 K3IRPfAN.net
お題: あなたは研究所で助手をしています
ある日、博士からこんな手紙が届きました

　助手よ。この数値の名前を教えて欲しい
　16776960

あなたは博士に返事を書こうと思います
その返事の内容をお答えください

321:デフォルトの名無しさん
19/12/15 07:44:00 VJoyKt1C.net
>>314 直感
Yellow

322:デフォルトの名無しさん
19/12/15 10:09:42 /SOvb5b/.net
黄色だな

323:デフォルトの名無しさん
19/12/15 11:56:51.84 wptEtnNg.net
お題
与えられた自然数の0を全て右端に寄せる
例
2019 -> 2190
102030 -> 123000
123 -> 123

324:デフォルトの名無しさん
19/12/15 12:02:23.68 fpSJINfx.net
>>309
Ruby で、
require 'active_support/time'
def calculate_time_span( now, alarm )
current = Time.now
now_time = Time.parse( now, current ) #=> 2019-12-15 22:00:00 +0900
alarm_time = Time.parse( alarm, current ) #=> 2019-12-15 04:00:00 +0900
diff_time = alarm_time - now_time
# アラームが翌日の場合
diff_time = alarm_time.tomorrow - now_time if diff_time < 0
Time.at( diff_time ).utc.strftime( '%H:%M' )
end
inputs = %w(10:00 16:00 22:00 04:00) # (現在時刻, アラーム時刻)の組
inputs.each_slice( 2 ) do |now, alarm| # 2個ずつ処理する
puts calculate_time_span( now, alarm )
end

次ページ