プログラミングのお題スレ Part9 [無断転載禁止]©2ch.net

**デフォルトの名無しさん** · 2016/12/01(木) 16:58:30.97

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

前スレ
プログラミングのお題スレ Part8©2ch.net
http://echo.2ch.net/test/read.cgi/tech/1444216746/

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

2 名前：デフォルトの名無しさん
　　>>1 使用言語
　　回答本文

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

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

**デフォルトの名無しさん** · 2017/07/23(日) 11:26:56.94

>>375
のほかの実装はでてこないねぇ‥

**デフォルトの名無しさん** · 2017/07/23(日) 12:53:55.68

>>394
使えるコードにするためには、規模がでかくなりすぎるから

**デフォルトの名無しさん** · 2017/07/23(日) 14:15:20.13

C/C++ で最長1000行ぐらいとみて、2日ぐらいあれば、とりあえず動く
土日で仕上がってくるんじゃないかと期待してたんだが

**デフォルトの名無しさん** · 2017/07/23(日) 14:18:11.58

速度が考えられてないコードなんて実用にはならないよ

**デフォルトの名無しさん** · 2017/07/23(日) 14:21:31.79

ていうか、
コードに対する条件とか
サポートする機能とか
条件が無さすぎる

**デフォルトの名無しさん** · 2017/07/23(日) 14:24:57.75

速度‥か‥

どうしてもローテートとかキャリーフラグとかを使いたいから、これはアセンブラの領域になるね
よくみかけるアセンブラ中毒者が今頃爪を研いでいるのだろうか？

**デフォルトの名無しさん** · 2017/07/23(日) 14:25:42.29

>>398
そこは「自分用」だから自由に決めていいんでないかい？

**デフォルトの名無しさん** · 2017/07/23(日) 14:49:34.77

>>394

当日に >>305 で >>390 より10倍以上早いのがでているだろう。
しかも計算量まで書いてある

**デフォルトの名無しさん** · 2017/07/23(日) 14:53:19.94

>>401
お題が違うのでは？

**デフォルトの名無しさん** · 2017/07/23(日) 15:05:28.62

>>402 >>394
確かに違った、すいません。

c++多倍長なら、karatsubaにも対応して300行くらいの以下をパクって使うのも一案

http://sites.google.com/site/indy256/algo_cpp/bigint

**デフォルトの名無しさん** · 2017/07/23(日) 15:09:21.91

>>403
base 10進ならば、表示（operator<<) が楽でいいね、なるほど、それは思いつかなかった

**デフォルトの名無しさん** · 2017/07/23(日) 16:51:44.88

>>399
通常、
処理時間のほとんどが乗算
乗算のほとんどがFFT
アセンブラの出番は当分先

**デフォルトの名無しさん** · 2017/07/23(日) 16:52:43.21

FFTのライブラリをどこからか持ってくるのでもいいけど、
それなら素直に多倍長ライブラリを持ってくれば
ってことになる

**デフォルトの名無しさん** · 2017/07/23(日) 16:54:07.19

今は浮動小数点演算が速いので、
カラツバの出番はあまりない

**デフォルトの名無しさん** · 2017/07/23(日) 18:50:17.74

基数を10のべき乗にするとか(printf()的なものが簡単だから)、乗算はunsigned shortやintとの
乗算に限るとか、除算無しとかいうのは…

プログラムの本体に組み込まれてしまって、再利用可能なライブラリの形で括りだされてる事の
方が少ないかｗ

**デフォルトの名無しさん** · 2017/07/24(月) 00:48:44.80

>>6-285 SageMath
ttp://sagecell.sagemath.org/?q=veftjc

**デフォルトの名無しさん** · 2017/07/24(月) 18:24:00.06

>>408
裁定ラインとしては、乗算は Bigint×Bigint、および除算の実装ですかね、でも足し算の回数での乗算や引き算の回数での除算は嫌ですね

**デフォルトの名無しさん** · 2017/07/24(月) 18:58:12.05

お題：お題スレ3の>>170をファレイ数列を使って解く。
http://peace.2ch.net/test/read.cgi/tech/1390525149/170

**デフォルトの名無しさん** · 2017/07/24(月) 19:10:35.21

>>411 Ruby
def farey_sequence(n)
(1..n-1).map{|i| 1r*i/n}
end

def ans_411(m)
(2..m).map{|i| farey_sequence(i)}.flatten.uniq.sort
end

ans_411 3 #=> [(1/3), (1/2), (2/3)]
ans_411 5 #=> [(1/5), (1/4), (1/3), (2/5), (1/2), (3/5), (2/3), (3/4), (4/5)]

**デフォルトの名無しさん** · 2017/07/24(月) 19:11:05.16

>>296
超高速版が出来ました！
http://ideone.com/FrRkof

一皿9秒が上限であれば、計算オーダーはn
関数自体は何秒でも大丈夫です

コードだけじゃ意味がわからないでしょうけど、とりあえずコードだけ
あまりテストしてないので、バグってたらごめんなさい

**デフォルトの名無しさん** · 2017/07/24(月) 23:00:47.69

オーダーnは凄いな

**デフォルトの名無しさん** · 2017/07/25(火) 05:40:32.49

>>413
うーんわからん。
俺の思考とは別系統かな。
ホントに0秒で解けてるし、素晴らしい。
素直に賞賛。

**デフォルトの名無しさん** · 2017/07/25(火) 11:52:16.57

回転寿しの問題は、部分的な最短経路が全体の最短経路にならないんだよな
だが最短時間はレーン長の2倍程度の再帰回数で出る
そのあと数十億回再帰して総当たりしてもそれより短くならない
最後の皿から逆方向に探索してもおそらく同じ状況

**デフォルトの名無しさん** · 2017/07/25(火) 11:56:47.84

例えば、”122” は最短時間6だが、1周目で2番目の要素”2”をパスしないとそうならない

**411** · 2017/07/26(水) 19:54:35.63

>>412
ファレイ数列の中間数(mediant)を再帰的に生成すると、uniqもsortも要らないのだけど、
mが3や5だと大差無いかｗ

**デフォルトの名無しさん** · 2017/07/26(水) 20:50:49.45

>>411
リンク先が見えません
問題文をもう一回書いてください

**デフォルトの名無しさん** · 2017/07/26(水) 20:52:34.29

と思ったら見れました

ファレイ数列を使って何かを解くわけじゃなくて、
ファレイ数列を求める問題？

**411** · 2017/07/26(水) 23:20:07.89

>>420
元の問題はそういうもの(＝ファレイ数列の両端(0/1と1/1)無し版を求める問題)と
解釈してますです。

**デフォルトの名無しさん** · 2017/07/26(水) 23:26:01.52

#include <list>
#include <iostream>
const int N_MAX = 10;
struct RATIONAL {
int num;
int den;
};
int main() {
std::list < RATIONAL > farey;
RATIONAL zero = {0, 1};
RATIONAL one = {1, 1};
farey.push_back(zero);
farey.push_back(one);
for (int n = 1; n <= N_MAX; n++){
for (std::list < RATIONAL > ::iterator i1 = farey.begin(), i0 = i1++; i1 != farey.end(); i0 = i1, i1++) {
if (i0->den + i1->den <= n) {
RATIONAL m = {i0->num + i1->num, i0->den + i1->den};
farey.insert(i1, m);
}
}
std::cout << n << " : ";
for (std::list < RATIONAL > ::iterator i = farey.begin(); i != farey.end(); i++) {
std::cout << i->num << "/" << i->den << " ";
}
std::cout << "\n";
}
return 0;
}

**デフォルトの名無しさん** · 2017/07/26(水) 23:29:22.49

これから0と1を除けば良いって問題であれば、
表示のループに以下を加えれば
if (i->den != 1)

**デフォルトの名無しさん** · 2017/07/26(水) 23:31:57.11

問題の意味も意図も良くわからん

出題者が「そういうものと解釈しています」とか
出題者が >>418 みたいな回答をバカにする発言とか

なんか非常に感じが悪い

**412** · 2017/07/27(木) 00:12:38.86

そもそも>>412のfarey_sequenceは定義が間違ってたわ
んでもって再帰にすると>>412より遅くなるという
Ruby

class Farey
def self.[](m)
if m == 1
[0r, 1r]
else
succ(m - 1)
end
end

def self.succ(m)
self[m].each_cons(2).inject([0r]){|s, (a, b)|
x = a.denominator + b.denominator
s << 1r*(a.numerator + b.numerator)/x if x == m + 1
s << b
}
end
end

Farey[3] # => [(0/1), (1/3), (1/2), (2/3), (1/1)]
Farey[5] # => [(0/1), (1/5), (1/4), (1/3), (2/5), (1/2), (3/5), (2/3), (3/4), (4/5), (1/1)]

**デフォルトの名無しさん** · 2017/07/27(木) 01:59:46.61

>>411 Java
https://ideone.com/w0q7cN

**デフォルトの名無しさん** · 2017/07/28(金) 18:51:16.80

>>375
のほかの実装はでてこないねぇ‥

**デフォルトの名無しさん** · 2017/07/28(金) 19:19:26.77

昔brainf**kで実装したのあるけどちょっとなぁ

**デフォルトの名無しさん** · 2017/07/28(金) 19:24:55.65

また懐かしい言語を

**デフォルトの名無しさん** · 2017/07/28(金) 19:26:36.86

どうせならチューリングマシンで作ってよ

**デフォルトの名無しさん** · 2017/07/30(日) 10:59:37.37

お題：MathematicaのFareySequence[n,k](引数2つ)に相当するものの実装
ttp://reference.wolfram.com/language/ref/FareySequence.html

**デフォルトの名無しさん** · 2017/07/30(日) 11:53:03.31

>>431
http://ideone.com/m7BnJN
C++。一瞬計算が合わなくてビビったけど、空目だった。
インデックス概念がベーシックなんだな。

**デフォルトの名無しさん** · 2017/07/30(日) 12:00:36.94

っていうか、この関数インデックスに0与えたら何が出力されるんだろう・・・。
早速バグってる気がする。

**デフォルトの名無しさん** · 2017/07/30(日) 12:15:40.85

>>432
バグってた。のでエディトしてFIXした。
所持する数の概念勘違いしてた。

**デフォルトの名無しさん** · 2017/07/30(日) 12:25:59.30

>>422の微妙な変更でいいよね

**デフォルトの名無しさん** · 2017/07/30(日) 12:26:41.34

1個だけ求めるなら、もっといい方法がある？

**デフォルトの名無しさん** · 2017/07/30(日) 12:27:22.69

ていうか、いい加減Fareyはもういいでしょ
他の課題の方が

**デフォルトの名無しさん** · 2017/07/30(日) 12:30:15.15

フィボナッチって素数何だっけ？

**デフォルトの名無しさん** · 2017/07/30(日) 12:44:33.08

1, 1, 2, 3, 5, 8, ...
違うよね

**デフォルトの名無しさん** · 2017/07/30(日) 12:47:15.18

だよねー。>>422ってフィボナッチ使ってない？
あんまり深く考えてないだけど。Orz

**デフォルトの名無しさん** · 2017/07/30(日) 12:47:30.85

じゃあ、任意の二個の数からはじまるフィボナッチ数列で、はじめから連続する素数の数が多い物を探す
って課題で

**デフォルトの名無しさん** · 2017/07/30(日) 12:48:36.68

>>440
フィボナッチではない
wikipediaにのってるレベルの知識で作った

**デフォルトの名無しさん** · 2017/07/30(日) 12:49:05.29

あれ？俺とんちんかんなこと言ってるか？
>>422が数列としてあってるのかよくわからない。Orz
どう考えればいいんだろう。

**デフォルトの名無しさん** · 2017/07/30(日) 12:53:04.00

まぁ、俺のもあってるかどうかはしらんけど。＞＜；

**デフォルトの名無しさん** · 2017/07/30(日) 12:56:43.60

頭悪くてゴメン。爆発しそう。。。

**デフォルトの名無しさん** · 2017/07/30(日) 13:01:05.96

引っ込む。すまんかった。

**デフォルトの名無しさん** · 2017/07/30(日) 14:17:22.65

>>431 Java
http://ideone.com/9AXdRV

**デフォルトの名無しさん** · 2017/07/30(日) 19:16:10.29

>>431 Ruby

def farey(n, k)
return [0r, 1r][k] if n == 1
farey(n - 1, 0..-1).each_cons(2).inject([0r]){|s, (a, b)|
x = a.denominator + b.denominator
s << 1r*(a.numerator + b.numerator)/x if x == n
s << b
}[k]
end

**デフォルトの名無しさん** · 2017/08/03(木) 07:36:01.80

お題：ミンコフスキーの疑問符関数の実装
ttps://en.wikipedia.org/wiki/Minkowski%27s_question_mark_function
ttp://reference.wolfram.com/language/ref/MinkowskiQuestionMark.html

**デフォルトの名無しさん** · 2017/08/03(木) 10:39:36.15

WIKIぺにコード乗ってますが。

**デフォルトの名無しさん** · 2017/08/03(木) 10:48:34.50

>>449のWIKIより。
/* Minkowski's question mark function */
double minkowski(double x) {
long p=x; if ((double)p>x) --p; /* p=floor(x) */
long q=1, r=p+1, s=1, m, n;
double d=1, y=p;
if (x<(double)p||(p<0)^(r<=0)) return x; /* out of range ?(x) =~ x */
for (;;) /* invariants: q*r-p*s==1 && (double)p/q <= x && x < (double)r/s */
{
d/=2; if (y+d==y) break; /* reached max possible precision */
m=p+r; if ((m<0)^(p<0)) break; /* sum overflowed */
n=q+s; if (n<0) break; /* sum overflowed */

if (x<(double)m/n) r=m, s=n;
else y+=d, p=m, q=n;
}
return y+d; /* final round-off */
}

**デフォルトの名無しさん** · 2017/08/05(土) 17:44:11.83

>>375
のほかの実装はでてこないねぇ‥

**デフォルトの名無しさん** · 2017/08/12(土) 18:46:00.57

寿司のオーダーNのやつを理解しようとしたけどまだやってない。
その仕組みと、ほんとに正解してるのかとか。いたら誰が解説して。

**デフォルトの名無しさん** · 2017/08/12(土) 19:04:22.18

>>413です

もちろん合っているつもりのコードです
作者が言っても何の説得力もありませんが

**デフォルトの名無しさん** · 2017/08/12(土) 19:07:04.34

会社に帰ってこない巡回セールスマンだよね
寿司の乗った皿がノード、計算量はO(n!)

**デフォルトの名無しさん** · 2017/08/12(土) 19:10:18.10

それぞれの寿司を食べている期間をレーン上の線分で表します

この線の重なり具合をpileで表しました

効率良く食べられた場合はレーンがpile_max周するまでの間に食べきることが出来ます

170行目の判定がそれで、trueの場合は効率良く食べられない場合です

**デフォルトの名無しさん** · 2017/08/12(土) 19:12:06.32

>>456
もしそれで最適解が得られるなら巡回セールスマンも可能じゃないかな？

**デフォルトの名無しさん** · 2017/08/12(土) 19:17:11.73

巡回セールスマン問題とけたら色々応用範囲アルヨ。
マジでどっかに売り込んでもいいくらい。
天才か。

**デフォルトの名無しさん** · 2017/08/12(土) 19:18:34.85

社会的に言うと交通統制とかもそれじゃないかな？
信号の待ち時間問題。よくしらんけど。

**デフォルトの名無しさん** · 2017/08/12(土) 19:19:17.76

効率良く食べられない方が簡単なのでその場合から

お寿司を以下のグループに分けます
----
各グループのお寿司は、レーンの特定の位置から食べ始めた場合、pile[グループ]周以内で食べ終わることが出来る
このとき、pile_max = Σ pile[グループ]
となる
---
このようなグループの分け方の最小の物が存在します

**デフォルトの名無しさん** · 2017/08/12(土) 19:22:56.16

同じグループのお寿司は連続して食べます
開始時と、各グループのお寿司を食べ終わった後、最初に来るお寿司から食べはじめ、pile[グループ]以内で食べられる食べ方でそのグループを食べ終える
ということを繰り返せば最小の時間で食べ終えることが出来ます

**デフォルトの名無しさん** · 2017/08/12(土) 19:26:29.79

グループ分けした時に1個のグループになった場合は、
効率良く食べられることになります
つまり、pile_max周以下で食べ終えることが出来ます

この時は、コード上にあるダミーのお寿司を追加してから最小時間を求め、ダミーのお寿司を食べてる時間を引けば求められます

**デフォルトの名無しさん** · 2017/08/12(土) 19:28:18.79

うーん、よくわからん
セールスマンの巡回先を一次元にマッピングできれば同じことできそうな
無理か

**デフォルトの名無しさん** · 2017/08/12(土) 19:30:27.01

グループの分け方は少し難しいです

レーンの各整数位置に対して、
お寿司の線の両端にあたる点同士
線の重なりがpile_max未満である区間の点(両端を含む)
を同じグループの点とし、
これらを続けることで最小のグループ分けが出来ます
線の両端の点のグループが、そのお寿司のグループになります

**デフォルトの名無しさん** · 2017/08/12(土) 19:31:42.19

それぞれ、証明は出来ているつもりです

**デフォルトの名無しさん** · 2017/08/12(土) 19:32:49.51

もちろん、一般の巡回問題はこの方法では無理です

**デフォルトの名無しさん** · 2017/08/12(土) 19:37:29.23

全ノードを巡回する最短時間の問題だから、できそうな気がするけどね

**デフォルトの名無しさん** · 2017/08/12(土) 19:39:44.61

372仕様書無しさん2017/08/11(金) 10:31:43.41
フリーランスで検索すると引っかかる零細ITがやっているフリーランスのサイトはだめだ。
高額に見せているけど実際は50万前後
JIET加入した方がいいよ。案件は毎日千件以上末端価格は60万円　平凡な稼働時間の80万円の案件もある。
ユー子も求人をだしてる。名刺も渡せる。ユー子に名刺が渡せるんだぞ。夢のようだ

それらの案件まさぐってHPで転売していたのが零細ITがやるフリーランスサイト

473非決定性名無しさん2017/08/03(木) 15:21:30.71

JIETに加入すれば誰でも3次60万からスタートだ。フリーランスのサイトをやってる
自称エージェントもそこから案件情報を取得しきてる。サイトで60万で釣って40万から55万の
間でやらしている。

446非決定性名無しさん2017/08/02(水) 22:12:48.95

JIETに毎月5千円払えば３次から入場できるだろ？
高額をうたうフリーランスのサイトはだいたい５次から４５万円
JIETで閲覧応募できる末端価格からさらに搾取するのが高額をみせつけるフリーランスサイトでした
高額案件をみせつけるフリーランスサイトも案件の取得はJIETでした

自称エージェントはJIETから流れてくる案件を転売してるだけだった。
JIETに加入すれば誰でも案件に応募することができた。収入が40万50万台にならなくて済む

**デフォルトの名無しさん** · 2017/08/12(土) 19:40:20.63

pile_maxとその位置から下限が得られますが、
>>296 の例では98秒の物以外はすべてその下限になっています
一個その下限になるような例を見つければ答えがわかるのですが、
自力で検索してみればわかると思いますがそのような例はあっさり見つかります

98秒の例は効率良く食べられない場合になります

効率良く食べられる側のなかでも、pileから得られる下限値より大きくなる場合もあります

**デフォルトの名無しさん** · 2017/08/12(土) 19:43:06.99

いずれの場合も、PCを使わなくても手計算で十分可能です

**デフォルトの名無しさん** · 2017/09/15(金) 10:14:33.05

お題:
N次元で１辺のマス目がM個の魔法陣を作る
N>3(任意)、M>=3(任意)の超立方体

**デフォルトの名無しさん** · 2017/09/15(金) 10:20:33.77

参考
http://magcube.la.coocan.jp/magcube/jp/
http://nadamath2012.web.fc2.com/bushi/2004_mas2.pdf

**デフォルトの名無しさん** · 2017/09/17(日) 16:38:38.03

魔方陣は1個作ればいいの？

Mが奇数か4の倍数は簡単
4で割って2余るのは検索するしかないのかな？

**デフォルトの名無しさん** · 2017/09/17(日) 16:55:04.37

バックトラックで組もうかと思ったけど、重そうだったからやめた。
数独より重そう。
それに一列合計をどの数字にするのかちょっとわからなかった。

**デフォルトの名無しさん** · 2017/09/17(日) 23:20:36.13

一列合計は、M*[数字の平均]
になる

つまり
M*(M^N+1)/2

**片山博文MZ** ◆T6xkBnTXz7B0 · 2017/09/18(月) 21:53:54.06

お題: URLから適当なサムネイルを生成するWebプログラム。

**デフォルトの名無しさん** · 2017/09/18(月) 23:06:01.10

お題
0以上90未満の整数nを入力として
タンジェントn°の値が有理数ならば真
そうでなければ偽を返す

**デフォルトの名無しさん** · 2017/09/18(月) 23:29:46.49

bool f(int n){return n==0 || n == 45;}

**デフォルトの名無しさん** · 2017/09/18(月) 23:33:52.62

sed -r -e "s/^(0|45)\$/True/" -e "s/[1-8][0-9]*/False/"

**デフォルトの名無しさん** · 2017/09/19(火) 01:13:30.78

tanの計算しないのはどうかと

**デフォルトの名無しさん** · 2017/09/19(火) 01:57:23.24

計算で有理数かどうか確認？
それは非常に難しいな

by 東大数学科卒

**デフォルトの名無しさん** · 2017/09/19(火) 02:28:36.10

>>480
そう思うなら他者を批判するより行動で示せばいいと思うよ

**デフォルトの名無しさん** · 2017/09/19(火) 03:58:59.82

tan1°が無理数であることの証明すら面倒くせえのに一体どんな回答を求めているんだ

**デフォルトの名無しさん** · 2017/09/19(火) 06:37:25.45

>>483
面倒くさい？
てことは出来るの？
やってみて

**デフォルトの名無しさん** · 2017/09/19(火) 07:41:26.64

>>484
確か京大の過去問にあったでしょう
説明めんどいからは解法は自分で調べて

**デフォルトの名無しさん** · 2017/09/19(火) 08:13:59.26

いや、おれは出来るよ

>>483の実力で出来るのか？と疑問に思っただけ
実力じゃなくてカンニングね

**デフォルトの名無しさん** · 2017/09/19(火) 08:21:42.19

何言ってんだこいつ

**デフォルトの名無しさん** · 2017/09/19(火) 08:49:35.38

問題が悪いな

与えられた有理数rに対し、
tan(πr)が有理数かどうか判別するプログラムを書け

ならテーブルは使えない

**デフォルトの名無しさん** · 2017/09/19(火) 09:06:24.19

また、多倍長精度演算のないC++にはきつい問題を・・・。

**デフォルトの名無しさん** · 2017/09/19(火) 11:09:01.37

>>488
すいません
結果を知っていたらこれでも簡単でした

**デフォルトの名無しさん** · 2017/09/19(火) 11:10:19.22

>>489
多倍長を使っても解決しないでしょ

**デフォルトの名無しさん** · 2017/09/19(火) 12:53:02.28

そもそも出題者はどういう回答を期待してるんだ？
数学の知識無しでは作れないし、数学の知識を使えば>>478になる

**デフォルトの名無しさん** · 2017/09/19(火) 14:38:36.12

tan()の加法定理
　tan（α+β）=（tanα+tanβ）/（1-tanαtanβ）
により

もしtan（α）が有理数なら
tan（nα）　(n = 1,2,3,4・・・)
も全て有理数

このため
整数nにより
tan（n）が無理数なら
nの約数全てによるtan()が無理数

ここで
tan(60)＝√３
が無理数なのは簡単に証明されるため、
tan(1)
も無理数

証明終わり

**デフォルトの名無しさん** · 2017/09/19(火) 14:54:56.80

>>476を解くにはあとtan(18度)が無理数であることを証明しないと