プログラミングのお題スレです。
前スレ
プログラミングのお題スレ 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/
宿題は宿題スレがあるのでそちらへ。
プログラミングのお題スレ Part9 [無断転載禁止]©2ch.net
■ このスレッドは過去ログ倉庫に格納されています
2016/12/01(木) 16:58:30.97ID:gTkHDluD
371デフォルトの名無しさん
2017/07/14(金) 14:31:43.15ID:DwybRUfK 競プロみたいな相手方の環境使う物だと標準と準標準の差はでかい
自分の環境なら導入すればいいだけだが
自分の環境なら導入すればいいだけだが
>>370
個体発生は系統発生を繰り返す
個体発生は系統発生を繰り返す
373デフォルトの名無しさん
2017/07/15(土) 12:42:06.98ID:odVkuNfb >>361
厳密解を出しているのなら、チャレンジ
(わかって近似値解狙いなら気にしないで)
"14432" と "887654329"
両方とも既出の"貪欲つぶし"(?)数列
"14432"は 20秒 (ゼロインデック順で02341)
"887654329"は 80秒(同123456708)でいける。
厳密解を出しているのなら、チャレンジ
(わかって近似値解狙いなら気にしないで)
"14432" と "887654329"
両方とも既出の"貪欲つぶし"(?)数列
"14432"は 20秒 (ゼロインデック順で02341)
"887654329"は 80秒(同123456708)でいける。
374デフォルトの名無しさん
2017/07/15(土) 14:59:21.20ID:OEoVgGO0 >>373
http://ideone.com/cBzPSj
C++。それ解くとほかの問題が解けなくなる。
厳密解のつもりだったが、ちょっと自分の領分超えてるなぁ。
うまくいかないものだ。
真実が奥の方にあると貪欲法は弱いな。Orz
http://ideone.com/cBzPSj
C++。それ解くとほかの問題が解けなくなる。
厳密解のつもりだったが、ちょっと自分の領分超えてるなぁ。
うまくいかないものだ。
真実が奥の方にあると貪欲法は弱いな。Orz
375デフォルトの名無しさん
2017/07/16(日) 16:33:07.01ID:8ZBD9z9c お題:
自分用多倍長整数演算関数
…って思ったけど、処理系の標準ではないとか、仕事でGNU MP使っては駄目とかの
制約で、簡易的なもの(乗算くらいまでとか)を書いた事ある人は少なくないと見た。
自分用多倍長整数演算関数
…って思ったけど、処理系の標準ではないとか、仕事でGNU MP使っては駄目とかの
制約で、簡易的なもの(乗算くらいまでとか)を書いた事ある人は少なくないと見た。
376デフォルトの名無しさん
2017/07/16(日) 18:30:49.57ID:8+Akms5T 多倍長整数演算がサポートされている言語を使う
終わり
終わり
378デフォルトの名無しさん
2017/07/16(日) 20:34:05.63ID:yctBkD01 掛け算の実装がキモだろう。
ここがボトルネックになるはず。
ここができると円周率とか、ルート計算も高速化できるはず。
ここがボトルネックになるはず。
ここができると円周率とか、ルート計算も高速化できるはず。
>>378
うん,FFTを使うそうだが‥いまいちよくわからない
うん,FFTを使うそうだが‥いまいちよくわからない
380375
2017/07/16(日) 20:41:46.81ID:8ZBD9z9c381デフォルトの名無しさん
2017/07/17(月) 22:48:25.95ID:5edeqhg+ >>296
手計算で計算出来るレベルにまで計算量を減らせた
もちろん数学的な裏付け付きで
ある条件を見たせば一瞬で求まる
"123456789123456789" > 98秒
残念ながら、これだけはその条件を満たしてない
手計算で計算出来るレベルにまで計算量を減らせた
もちろん数学的な裏付け付きで
ある条件を見たせば一瞬で求まる
"123456789123456789" > 98秒
残念ながら、これだけはその条件を満たしてない
382デフォルトの名無しさん
2017/07/18(火) 06:37:17.84ID:nFCFlf58 >>381
22とか2323もその条件を満たしてない感じ?
22とか2323もその条件を満たしてない感じ?
383デフォルトの名無しさん
2017/07/18(火) 07:37:05.09ID:Ew0RSScO 22 は微妙
2323 は大丈夫
2323 は大丈夫
384デフォルトの名無しさん
2017/07/18(火) 07:41:44.22ID:Ew0RSScO まだコードになってないんで、
コードになったらアップします
寿司を食べる時間 < レーンの回転周期
という前提をつけちゃおうと思ったけど、
つけない方が良さそうですね
寿司を食べる時間がレーンの回転周期の整数倍の寿司は
ちょっと特別な処理が必要
コードになったらアップします
寿司を食べる時間 < レーンの回転周期
という前提をつけちゃおうと思ったけど、
つけない方が良さそうですね
寿司を食べる時間がレーンの回転周期の整数倍の寿司は
ちょっと特別な処理が必要
385デフォルトの名無しさん
2017/07/18(火) 08:01:32.43ID:Ew0RSScO 整数倍の寿司が無いもので
条件に当てはまらない最小は
2222
かな
条件に当てはまらない最小は
2222
かな
386デフォルトの名無しさん
2017/07/18(火) 08:49:51.93ID:YLlwVFMJ388デフォルトの名無しさん
2017/07/19(水) 18:12:29.37ID:Np9hKHT2 >>296
>>373
http://ideone.com/B9vl8l
C++。結局、i7-6700のmem2G使って7分で解けた。
どうしようもない位遅いな。
でも一応題意には添えたと思う。
もう見たくない・・・。Orz
高速化するにはインラインアセンブリ使うか、スレッド分割できるようなアルゴリズムかんがえるか。
よくわからんけど、数学で頑張ってる人に期待だ。
>>373
http://ideone.com/B9vl8l
C++。結局、i7-6700のmem2G使って7分で解けた。
どうしようもない位遅いな。
でも一応題意には添えたと思う。
もう見たくない・・・。Orz
高速化するにはインラインアセンブリ使うか、スレッド分割できるようなアルゴリズムかんがえるか。
よくわからんけど、数学で頑張ってる人に期待だ。
389386
2017/07/20(木) 01:57:16.81ID:Q7XnESC/ 100のべき乗に変更
ttp://sagecell.sagemath.org/?q=mciykc
ttp://sagecell.sagemath.org/?q=mciykc
390デフォルトの名無しさん
2017/07/21(金) 15:21:18.13ID:7e+pM3K/ >>296
http://ideone.com/mXPglY
C++。試しに再起化してみたら処理速度倍になった。
自分の環境では3分ちょいで解ける。
相変わらずメモリ馬鹿食いするけど。
もう俺には無理。
俺の中では終了でーす。Orz
http://ideone.com/mXPglY
C++。試しに再起化してみたら処理速度倍になった。
自分の環境では3分ちょいで解ける。
相変わらずメモリ馬鹿食いするけど。
もう俺には無理。
俺の中では終了でーす。Orz
391デフォルトの名無しさん
2017/07/22(土) 08:54:36.28ID:OQXA8cUK392デフォルトの名無しさん
2017/07/22(土) 08:55:28.19ID:OQXA8cUK このスレが無くならないうちに解明します
393デフォルトの名無しさん
2017/07/22(土) 10:43:30.03ID:apsnR2Z0394デフォルトの名無しさん
2017/07/23(日) 11:26:56.94ID:ipiEUPYV >>375
のほかの実装はでてこないねぇ‥
のほかの実装はでてこないねぇ‥
395デフォルトの名無しさん
2017/07/23(日) 12:53:55.68ID:7fREas1L >>394
使えるコードにするためには、規模がでかくなりすぎるから
使えるコードにするためには、規模がでかくなりすぎるから
396デフォルトの名無しさん
2017/07/23(日) 14:15:20.13ID:ipiEUPYV C/C++ で最長1000行ぐらいとみて、2日ぐらいあれば、とりあえず動く
土日で仕上がってくるんじゃないかと期待してたんだが
土日で仕上がってくるんじゃないかと期待してたんだが
397デフォルトの名無しさん
2017/07/23(日) 14:18:11.58ID:7fREas1L 速度が考えられてないコードなんて実用にはならないよ
398デフォルトの名無しさん
2017/07/23(日) 14:21:31.79ID:7fREas1L ていうか、
コードに対する条件とか
サポートする機能とか
条件が無さすぎる
コードに対する条件とか
サポートする機能とか
条件が無さすぎる
399デフォルトの名無しさん
2017/07/23(日) 14:24:57.75ID:ipiEUPYV 速度‥か‥
どうしてもローテートとかキャリーフラグとかを使いたいから、これはアセンブラの領域になるね
よくみかけるアセンブラ中毒者が今頃爪を研いでいるのだろうか?
どうしてもローテートとかキャリーフラグとかを使いたいから、これはアセンブラの領域になるね
よくみかけるアセンブラ中毒者が今頃爪を研いでいるのだろうか?
400デフォルトの名無しさん
2017/07/23(日) 14:25:42.29ID:ipiEUPYV >>398
そこは「自分用」だから自由に決めていいんでないかい?
そこは「自分用」だから自由に決めていいんでないかい?
401デフォルトの名無しさん
2017/07/23(日) 14:49:34.77ID:TcY6qE9r402デフォルトの名無しさん
2017/07/23(日) 14:53:19.94ID:ipiEUPYV >>401
お題が違うのでは?
お題が違うのでは?
403デフォルトの名無しさん
2017/07/23(日) 15:05:28.62ID:TcY6qE9r >>402 >>394
確かに違った、すいません。
c++多倍長なら、karatsubaにも対応して300行くらいの以下をパクって使うのも一案
http://sites.google.com/site/indy256/algo_cpp/bigint
確かに違った、すいません。
c++多倍長なら、karatsubaにも対応して300行くらいの以下をパクって使うのも一案
http://sites.google.com/site/indy256/algo_cpp/bigint
404デフォルトの名無しさん
2017/07/23(日) 15:09:21.91ID:ipiEUPYV >>403
base 10進ならば、表示(operator<<) が楽でいいね、なるほど、それは思いつかなかった
base 10進ならば、表示(operator<<) が楽でいいね、なるほど、それは思いつかなかった
405デフォルトの名無しさん
2017/07/23(日) 16:51:44.88ID:7fREas1L406デフォルトの名無しさん
2017/07/23(日) 16:52:43.21ID:7fREas1L FFTのライブラリをどこからか持ってくるのでもいいけど、
それなら素直に多倍長ライブラリを持ってくれば
ってことになる
それなら素直に多倍長ライブラリを持ってくれば
ってことになる
407デフォルトの名無しさん
2017/07/23(日) 16:54:07.19ID:7fREas1L 今は浮動小数点演算が速いので、
カラツバの出番はあまりない
カラツバの出番はあまりない
408デフォルトの名無しさん
2017/07/23(日) 18:50:17.74ID:5CSy1R8t 基数を10のべき乗にするとか(printf()的なものが簡単だから)、乗算はunsigned shortやintとの
乗算に限るとか、除算無しとかいうのは…
プログラムの本体に組み込まれてしまって、再利用可能なライブラリの形で括りだされてる事の
方が少ないかw
乗算に限るとか、除算無しとかいうのは…
プログラムの本体に組み込まれてしまって、再利用可能なライブラリの形で括りだされてる事の
方が少ないかw
409デフォルトの名無しさん
2017/07/24(月) 00:48:44.80ID:XgJeE+LA >>6-285 SageMath
ttp://sagecell.sagemath.org/?q=veftjc
ttp://sagecell.sagemath.org/?q=veftjc
410デフォルトの名無しさん
2017/07/24(月) 18:24:00.06ID:UuUAOyUA >>408
裁定ラインとしては、乗算は Bigint×Bigint、および除算の実装ですかね、でも足し算の回数での乗算や引き算の回数での除算は嫌ですね
裁定ラインとしては、乗算は Bigint×Bigint、および除算の実装ですかね、でも足し算の回数での乗算や引き算の回数での除算は嫌ですね
411デフォルトの名無しさん
2017/07/24(月) 18:58:12.05ID:5ve8i6tz412デフォルトの名無しさん
2017/07/24(月) 19:10:35.21ID:nJVItCRy >>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)]
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)]
413デフォルトの名無しさん
2017/07/24(月) 19:11:05.16ID:7nQ6Z7f9 >>296
超高速版が出来ました!
http://ideone.com/FrRkof
一皿9秒が上限であれば、計算オーダーはn
関数自体は何秒でも大丈夫です
コードだけじゃ意味がわからないでしょうけど、とりあえずコードだけ
あまりテストしてないので、バグってたらごめんなさい
超高速版が出来ました!
http://ideone.com/FrRkof
一皿9秒が上限であれば、計算オーダーはn
関数自体は何秒でも大丈夫です
コードだけじゃ意味がわからないでしょうけど、とりあえずコードだけ
あまりテストしてないので、バグってたらごめんなさい
414デフォルトの名無しさん
2017/07/24(月) 23:00:47.69ID:fjGi9Yh0 オーダーnは凄いな
415デフォルトの名無しさん
2017/07/25(火) 05:40:32.49ID:ubbfnjuS416デフォルトの名無しさん
2017/07/25(火) 11:52:16.57ID:bLUUDw7G 回転寿しの問題は、部分的な最短経路が全体の最短経路にならないんだよな
だが最短時間はレーン長の2倍程度の再帰回数で出る
そのあと数十億回再帰して総当たりしてもそれより短くならない
最後の皿から逆方向に探索してもおそらく同じ状況
だが最短時間はレーン長の2倍程度の再帰回数で出る
そのあと数十億回再帰して総当たりしてもそれより短くならない
最後の皿から逆方向に探索してもおそらく同じ状況
417デフォルトの名無しさん
2017/07/25(火) 11:56:47.84ID:bLUUDw7G 例えば、”122” は最短時間6だが、1周目で2番目の要素”2”をパスしないとそうならない
418411
2017/07/26(水) 19:54:35.63ID:6H34MdHA419デフォルトの名無しさん
2017/07/26(水) 20:50:49.45ID:s8dUUqTb420デフォルトの名無しさん
2017/07/26(水) 20:52:34.29ID:s8dUUqTb と思ったら見れました
ファレイ数列を使って何かを解くわけじゃなくて、
ファレイ数列を求める問題?
ファレイ数列を使って何かを解くわけじゃなくて、
ファレイ数列を求める問題?
421411
2017/07/26(水) 23:20:07.89ID:6H34MdHA422デフォルトの名無しさん
2017/07/26(水) 23:26:01.52ID:lPM9zwS7 #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;
}
#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;
}
423デフォルトの名無しさん
2017/07/26(水) 23:29:22.49ID:lPM9zwS7 これから0と1を除けば良いって問題であれば、
表示のループに以下を加えれば
if (i->den != 1)
表示のループに以下を加えれば
if (i->den != 1)
424デフォルトの名無しさん
2017/07/26(水) 23:31:57.11ID:lPM9zwS7425412
2017/07/27(木) 00:12:38.86ID:qteH6K3e そもそも>>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)]
んでもって再帰にすると>>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)]
426デフォルトの名無しさん
2017/07/27(木) 01:59:46.61ID:GuEy9AL1427デフォルトの名無しさん
2017/07/28(金) 18:51:16.80ID:XBSdfIgC >>375
のほかの実装はでてこないねぇ‥
のほかの実装はでてこないねぇ‥
428デフォルトの名無しさん
2017/07/28(金) 19:19:26.77ID:mqZJq6H+ 昔brainf**kで実装したのあるけどちょっとなぁ
429デフォルトの名無しさん
2017/07/28(金) 19:24:55.65ID:WViVOgsq また懐かしい言語を
430デフォルトの名無しさん
2017/07/28(金) 19:26:36.86ID:WViVOgsq どうせならチューリングマシンで作ってよ
431デフォルトの名無しさん
2017/07/30(日) 10:59:37.37ID:A7gIx2b1 お題:MathematicaのFareySequence[n,k](引数2つ)に相当するものの実装
ttp://reference.wolfram.com/language/ref/FareySequence.html
ttp://reference.wolfram.com/language/ref/FareySequence.html
432デフォルトの名無しさん
2017/07/30(日) 11:53:03.31ID:EQKnHSgY433デフォルトの名無しさん
2017/07/30(日) 12:00:36.94ID:EQKnHSgY っていうか、この関数インデックスに0与えたら何が出力されるんだろう・・・。
早速バグってる気がする。
早速バグってる気がする。
434デフォルトの名無しさん
2017/07/30(日) 12:15:40.85ID:EQKnHSgY435デフォルトの名無しさん
2017/07/30(日) 12:25:59.30ID:B3p9Yl5S >>422の微妙な変更でいいよね
436デフォルトの名無しさん
2017/07/30(日) 12:26:41.34ID:B3p9Yl5S 1個だけ求めるなら、もっといい方法がある?
437デフォルトの名無しさん
2017/07/30(日) 12:27:22.69ID:B3p9Yl5S ていうか、いい加減Fareyはもういいでしょ
他の課題の方が
他の課題の方が
438デフォルトの名無しさん
2017/07/30(日) 12:30:15.15ID:EQKnHSgY フィボナッチって素数何だっけ?
439デフォルトの名無しさん
2017/07/30(日) 12:44:33.08ID:B3p9Yl5S 1, 1, 2, 3, 5, 8, ...
違うよね
違うよね
440デフォルトの名無しさん
2017/07/30(日) 12:47:15.18ID:EQKnHSgY だよねー。>>422ってフィボナッチ使ってない?
あんまり深く考えてないだけど。Orz
あんまり深く考えてないだけど。Orz
441デフォルトの名無しさん
2017/07/30(日) 12:47:30.85ID:B3p9Yl5S じゃあ、任意の二個の数からはじまるフィボナッチ数列で、はじめから連続する素数の数が多い物を探す
って課題で
って課題で
442デフォルトの名無しさん
2017/07/30(日) 12:48:36.68ID:B3p9Yl5S443デフォルトの名無しさん
2017/07/30(日) 12:49:05.29ID:EQKnHSgY444デフォルトの名無しさん
2017/07/30(日) 12:53:04.00ID:EQKnHSgY まぁ、俺のもあってるかどうかはしらんけど。><;
445デフォルトの名無しさん
2017/07/30(日) 12:56:43.60ID:EQKnHSgY 頭悪くてゴメン。爆発しそう。。。
446デフォルトの名無しさん
2017/07/30(日) 13:01:05.96ID:EQKnHSgY 引っ込む。すまんかった。
447デフォルトの名無しさん
2017/07/30(日) 14:17:22.65ID:t+CfDp82448デフォルトの名無しさん
2017/07/30(日) 19:16:10.29ID:LizATlBz >>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
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
449デフォルトの名無しさん
2017/08/03(木) 07:36:01.80ID:cLWzUq7C お題:ミンコフスキーの疑問符関数の実装
ttps://en.wikipedia.org/wiki/Minkowski%27s_question_mark_function
ttp://reference.wolfram.com/language/ref/MinkowskiQuestionMark.html
ttps://en.wikipedia.org/wiki/Minkowski%27s_question_mark_function
ttp://reference.wolfram.com/language/ref/MinkowskiQuestionMark.html
450デフォルトの名無しさん
2017/08/03(木) 10:39:36.15ID:ONmyLPuf WIKIぺにコード乗ってますが。
451デフォルトの名無しさん
2017/08/03(木) 10:48:34.50ID:ONmyLPuf >>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 */
}
/* 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 */
}
452デフォルトの名無しさん
2017/08/05(土) 17:44:11.83ID:40G0sflG >>375
のほかの実装はでてこないねぇ‥
のほかの実装はでてこないねぇ‥
453デフォルトの名無しさん
2017/08/12(土) 18:46:00.57ID:953va2dM 寿司のオーダーNのやつを理解しようとしたけどまだやってない。
その仕組みと、ほんとに正解してるのかとか。いたら誰が解説して。
その仕組みと、ほんとに正解してるのかとか。いたら誰が解説して。
454デフォルトの名無しさん
2017/08/12(土) 19:04:22.18ID:Bi4KH0eW455デフォルトの名無しさん
2017/08/12(土) 19:07:04.34ID:4r/z/Qd5 会社に帰ってこない巡回セールスマンだよね
寿司の乗った皿がノード、計算量はO(n!)
寿司の乗った皿がノード、計算量はO(n!)
456デフォルトの名無しさん
2017/08/12(土) 19:10:18.10ID:Bi4KH0eW それぞれの寿司を食べている期間をレーン上の線分で表します
この線の重なり具合をpileで表しました
効率良く食べられた場合はレーンがpile_max周するまでの間に食べきることが出来ます
170行目の判定がそれで、trueの場合は効率良く食べられない場合です
この線の重なり具合をpileで表しました
効率良く食べられた場合はレーンがpile_max周するまでの間に食べきることが出来ます
170行目の判定がそれで、trueの場合は効率良く食べられない場合です
457デフォルトの名無しさん
2017/08/12(土) 19:12:06.32ID:4r/z/Qd5 >>456
もしそれで最適解が得られるなら巡回セールスマンも可能じゃないかな?
もしそれで最適解が得られるなら巡回セールスマンも可能じゃないかな?
458デフォルトの名無しさん
2017/08/12(土) 19:17:11.73ID:6XNTCj+p 巡回セールスマン問題とけたら色々応用範囲アルヨ。
マジでどっかに売り込んでもいいくらい。
天才か。
マジでどっかに売り込んでもいいくらい。
天才か。
459デフォルトの名無しさん
2017/08/12(土) 19:18:34.85ID:6XNTCj+p 社会的に言うと交通統制とかもそれじゃないかな?
信号の待ち時間問題。よくしらんけど。
信号の待ち時間問題。よくしらんけど。
460デフォルトの名無しさん
2017/08/12(土) 19:19:17.76ID:Bi4KH0eW 効率良く食べられない方が簡単なのでその場合から
お寿司を以下のグループに分けます
----
各グループのお寿司は、レーンの特定の位置から食べ始めた場合、pile[グループ]周以内で食べ終わることが出来る
このとき、pile_max = Σ pile[グループ]
となる
---
このようなグループの分け方の最小の物が存在します
お寿司を以下のグループに分けます
----
各グループのお寿司は、レーンの特定の位置から食べ始めた場合、pile[グループ]周以内で食べ終わることが出来る
このとき、pile_max = Σ pile[グループ]
となる
---
このようなグループの分け方の最小の物が存在します
461デフォルトの名無しさん
2017/08/12(土) 19:22:56.16ID:Bi4KH0eW 同じグループのお寿司は連続して食べます
開始時と、各グループのお寿司を食べ終わった後、最初に来るお寿司から食べはじめ、pile[グループ]以内で食べられる食べ方でそのグループを食べ終える
ということを繰り返せば最小の時間で食べ終えることが出来ます
開始時と、各グループのお寿司を食べ終わった後、最初に来るお寿司から食べはじめ、pile[グループ]以内で食べられる食べ方でそのグループを食べ終える
ということを繰り返せば最小の時間で食べ終えることが出来ます
462デフォルトの名無しさん
2017/08/12(土) 19:26:29.79ID:Bi4KH0eW グループ分けした時に1個のグループになった場合は、
効率良く食べられることになります
つまり、pile_max周以下で食べ終えることが出来ます
この時は、コード上にあるダミーのお寿司を追加してから最小時間を求め、ダミーのお寿司を食べてる時間を引けば求められます
効率良く食べられることになります
つまり、pile_max周以下で食べ終えることが出来ます
この時は、コード上にあるダミーのお寿司を追加してから最小時間を求め、ダミーのお寿司を食べてる時間を引けば求められます
463デフォルトの名無しさん
2017/08/12(土) 19:28:18.79ID:4r/z/Qd5 うーん、よくわからん
セールスマンの巡回先を一次元にマッピングできれば同じことできそうな
無理か
セールスマンの巡回先を一次元にマッピングできれば同じことできそうな
無理か
464デフォルトの名無しさん
2017/08/12(土) 19:30:27.01ID:Bi4KH0eW グループの分け方は少し難しいです
レーンの各整数位置に対して、
お寿司の線の両端にあたる点同士
線の重なりがpile_max未満である区間の点(両端を含む)
を同じグループの点とし、
これらを続けることで最小のグループ分けが出来ます
線の両端の点のグループが、そのお寿司のグループになります
レーンの各整数位置に対して、
お寿司の線の両端にあたる点同士
線の重なりがpile_max未満である区間の点(両端を含む)
を同じグループの点とし、
これらを続けることで最小のグループ分けが出来ます
線の両端の点のグループが、そのお寿司のグループになります
465デフォルトの名無しさん
2017/08/12(土) 19:31:42.19ID:Bi4KH0eW それぞれ、証明は出来ているつもりです
466デフォルトの名無しさん
2017/08/12(土) 19:32:49.51ID:Bi4KH0eW もちろん、一般の巡回問題はこの方法では無理です
467デフォルトの名無しさん
2017/08/12(土) 19:37:29.23ID:4r/z/Qd5 全ノードを巡回する最短時間の問題だから、できそうな気がするけどね
468デフォルトの名無しさん
2017/08/12(土) 19:39:44.61ID:2Yw2XYfL 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千円払えば3次から入場できるだろ?
高額をうたうフリーランスのサイトはだいたい5次から45万円
JIETで閲覧応募できる末端価格からさらに搾取するのが高額をみせつけるフリーランスサイトでした
高額案件をみせつけるフリーランスサイトも案件の取得はJIETでした
自称エージェントはJIETから流れてくる案件を転売してるだけだった。
JIETに加入すれば誰でも案件に応募することができた。収入が40万50万台にならなくて済む
フリーランスで検索すると引っかかる零細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千円払えば3次から入場できるだろ?
高額をうたうフリーランスのサイトはだいたい5次から45万円
JIETで閲覧応募できる末端価格からさらに搾取するのが高額をみせつけるフリーランスサイトでした
高額案件をみせつけるフリーランスサイトも案件の取得はJIETでした
自称エージェントはJIETから流れてくる案件を転売してるだけだった。
JIETに加入すれば誰でも案件に応募することができた。収入が40万50万台にならなくて済む
469デフォルトの名無しさん
2017/08/12(土) 19:40:20.63ID:Bi4KH0eW pile_maxとその位置から下限が得られますが、
>>296 の例では98秒の物以外はすべてその下限になっています
一個その下限になるような例を見つければ答えがわかるのですが、
自力で検索してみればわかると思いますがそのような例はあっさり見つかります
98秒の例は効率良く食べられない場合になります
効率良く食べられる側のなかでも、pileから得られる下限値より大きくなる場合もあります
>>296 の例では98秒の物以外はすべてその下限になっています
一個その下限になるような例を見つければ答えがわかるのですが、
自力で検索してみればわかると思いますがそのような例はあっさり見つかります
98秒の例は効率良く食べられない場合になります
効率良く食べられる側のなかでも、pileから得られる下限値より大きくなる場合もあります
470デフォルトの名無しさん
2017/08/12(土) 19:43:06.99ID:Bi4KH0eW いずれの場合も、PCを使わなくても手計算で十分可能です
■ このスレッドは過去ログ倉庫に格納されています
ニュース
- 【LIVE】国分太一 騒動後初の公の場 司法記者クラブで会見 [ひかり★]
- 生クリームだけの真っ白なクリスマスケーキ 大手メーカーが販売、その理由は…フルーツなしで価格は半額以下に [おっさん友の会★]
- 【文春】元TOKIO・国分太一(51)「女性スタッフ2名への“わいせつ事案”」日テレ事情聴取の全貌が分かった! ★2 [Ailuropoda melanoleuca★]
- 【山上裁判】安倍氏が狙わた理由 旧統一教会の関係者が「安倍氏は『われわれの味方』」と宣伝していた [1ゲットロボ★]
- 【文春】スクープ撮! 超人気ゴルファー・都玲華(21) “30歳上”石井忍コーチ(51)と路上レッスン禁断愛 [冬月記者★]
- ラサール石井氏「こういう空気が戦前の空気」と警鐘「すでに来てしまっている」「大政翼賛会も始まっています」 [少考さん★]
- 【高市悲報】「餃子の王将」社長射殺、初公判で被告側が無罪主張。遺族が法廷でブチギレ [834922174]
- 【速報】国分太一会見 [115996789]
- ひるおび「国連の常任理事国は中華民国。中華人民共和国は中華民国を正式に引き継いでない」高市 [931948549]
- ミヤネ屋、国分太一会見
- なあ、中出しセックスしても、そこで使用されたマンコの内部は洗われずに放置されるケースも往々にしてあるわけだろ?
- 【悲報】自民党、通行人に叩かれるwwwwwwwwwwwwwww [834922174]
