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

**デフォルトの名無しさん** · 2020/03/13(金) 01:42:47.06

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

【出題と回答例】
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/

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

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

**デフォルトの名無しさん** · 2020/06/03(水) 18:36:37.68

趣味のプログラミングと義務のプログラミングは違うくね？

**デフォルトの名無しさん** · 2020/06/03(水) 18:37:26.73

>>801
3ヶ月はやりすぎやろーｗって思って覗いたら力作だった
今後も頑張って！

◆QZaw55cn4c · 2020/06/03(水) 19:17:52.62

>>801
！！
js 習いたくなりました！

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/03(水) 19:29:40.86

戦艦もののシミュレーションゲームとか作って、ゲーム会社に売り込んでみたら？

**デフォルトの名無しさん** · 2020/06/03(水) 19:29:44.33

>>808
お前は先に解析やれよ

**デフォルトの名無しさん** · 2020/06/03(水) 19:29:56.09

>>801
おお cool

**デフォルトの名無しさん** · 2020/06/03(水) 19:30:52.69

久々に良いものを見た。
そういえばしばらく仕事以外のプログラミングから遠ざかっていたな...

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/03(水) 19:54:45.25

5chに天才ハッカー登場！

**デフォルトの名無しさん** · 2020/06/03(水) 22:36:49.52

お題を出される方は当然自身の答えを持ち合わせてから出すんですよね？
答え合わせはどうだったのですか？

**デフォルトの名無しさん** · 2020/06/03(水) 23:19:23.57

多分俺日本で一番技術力高いと思うわ
日本一になった今思うことは問題が簡単過ぎてつまらん。ライバルがいなくてつまらん
俺はいったい何をすればいいのか

**デフォルトの名無しさん** · 2020/06/03(水) 23:21:04.09

Adobeのようにソフトを市場にうちだすんだ
あなたはお金を得られて嬉しい
ユーザーは高い技術がうみだした良い製品を使えて嬉しい

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/03(水) 23:27:43.96

>>815
人工知能チップと人工知能ドライバーモデルを開発して人工知能OSを設計する。

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/03(水) 23:29:22.65

>>815
逆コンパイラを作る。

**デフォルトの名無しさん** · 2020/06/03(水) 23:32:03.80

お題：
説明変数に測定誤差が生じるとき最小二乗法は一致推定量とはならない
正しい尤度方程式を導出し最尤推定を行えば一致推定になるが、驚くべきことに測定誤差があるとき最尤推定は確率収束が最も早い推定にはならない
これをネイマン・スコット問題という
最も収束が早い推定にはどのようなものがあるかを述べ、実装を行え

**デフォルトの名無しさん** · 2020/06/03(水) 23:34:25.50

>>815
>>801を3次元化してSAOガンゲイルオンライン作れそう？

**デフォルトの名無しさん** · 2020/06/04(木) 11:03:46.04

お題：
文脈依存文法を実装せよ

**デフォルトの名無しさん** · 2020/06/04(木) 19:17:36.11

>>821
コテハン付け忘れてんぞハゲ

**デフォルトの名無しさん** · 2020/06/04(木) 21:00:37.65

>>822
ヒント：線形有界オートマトン

**デフォルトの名無しさん** · 2020/06/04(木) 21:50:03.46

世界的なゲームプログラマーと言えば、
スーパー鳥スロットの作者、jagarikin だろw

jagarikin
https://twitter.com/jagarikin/status/1189024376915521536
https://twitter.com/5chan_nel (5ch newer account)

**デフォルトの名無しさん** · 2020/06/04(木) 22:56:02.12

>>795
お題：文字列を読み取って同じ構造のフォルダを作ってください
入力:
フグ亜目
├─マンボウ科
│　├─ヤリマンボウ属
│　│　└─ヤリマンボウ
│　├─マンボウ属
│　│　├─ウシマンボウ
│　│　├─マンボウ
│　│　└─カクレマンボウ
│　└─クサビフグ属
│　　　└─クサビフグ
└─ハリセンボン科
　　└─ハリセンボン属
　　　　└─ハリセンボン

**デフォルトの名無しさん** · 2020/06/07(日) 00:53:20.20

>>795
お題：ソートされたパスの一覧を構造化して出力してください
入力：
フグ亜目
フグ亜目\ハリセンボン科
フグ亜目\ハリセンボン科\ハリセンボン属
フグ亜目\ハリセンボン科\ハリセンボン属\ハリセンボン
フグ亜目\マンボウ科
フグ亜目\マンボウ科\クサビフグ属
フグ亜目\マンボウ科\クサビフグ属\クサビフグ
フグ亜目\マンボウ科\マンボウ属
フグ亜目\マンボウ科\マンボウ属\ウシマンボウ
フグ亜目\マンボウ科\マンボウ属\カクレマンボウ
フグ亜目\マンボウ科\マンボウ属\マンボウ
フグ亜目\マンボウ科\ヤリマンボウ属
フグ亜目\マンボウ科\ヤリマンボウ属\ヤリマンボウ

出力：
フグ亜目
　　ハリセンボン科
　　　　ハリセンボン属
　　　　　　ハリセンボン
　　マンボウ科
　　　　クサビフグ属
　　　　　　クサビフグ
　　　　マンボウ属
　　　　　　ウシマンボウ
　　　　　　カクレマンボウ
　　　　　　マンボウ
　　　　ヤリマンボウ属
　　　　　　ヤリマンボウ

**デフォルトの名無しさん** · 2020/06/07(日) 01:34:43.22

>>826 Python
https://ideone.com/B07fvO

**デフォルトの名無しさん** · 2020/06/07(日) 01:39:43.13

>>826
align.pl
https://ideone.com/eVvf7e

sh && perl
https://i.imgur.com/v3CQyhy.png

**540** · 2020/06/07(日) 03:37:39.04

>>826 Perl5

for (<DATA>) {
　chomp;
　@a = split/\\/;
　sub {
　　$k = shift @a;
　　$_[0]->{$k} = {} unless exists $_[0]->{$k};
　　__SUB__->($_[0]->{$k}) if @a;
　}->(\%h);
}
sub ($s, $r) {
　for (keys %$r) {
　　print "$s$_\n";
　　__SUB__->($s.'　', $r->{$_})
　}
}->('', \%h);
__DATA__
フグ亜目
フグ亜目\ハリセンボン科
フグ亜目\ハリセンボン科\ハリセンボン属
フグ亜目\ハリセンボン科\ハリセンボン属\ハリセンボン
フグ亜目\マンボウ科
フグ亜目\マンボウ科\クサビフグ属
フグ亜目\マンボウ科\クサビフグ属\クサビフグ
フグ亜目\マンボウ科\マンボウ属
フグ亜目\マンボウ科\マンボウ属\ウシマンボウ
フグ亜目\マンボウ科\マンボウ属\カクレマンボウ
フグ亜目\マンボウ科\マンボウ属\マンボウ
フグ亜目\マンボウ科\ヤリマンボウ属
フグ亜目\マンボウ科\ヤリマンボウ属\ヤリマンボウ

**デフォルトの名無しさん** · 2020/06/07(日) 04:14:38.16

>>829 一行目にこれが要る
use feature qw{current_sub signatures};

**デフォルトの名無しさん** · 2020/06/07(日) 04:16:35.18

>>829
ｆor (keys %$r) {
↓
for (sort keys %$r) {

**540** · 2020/06/07(日) 11:06:41.99

>>830 ハイハイ、修正修正

use feature qw{current_sub signatures};
for (<DATA>) {
　chomp;
　@a = split/\\/;
　sub ($hr) {
　　$k = shift @a;
　　$hr->{$k} //= {};
　　__SUB__->($hr->{$k}) if @a;
　}->(\%h);
}
sub ($s, $r) {
　for (sort keys %$r) {
　　print "$s$_\n";
　　__SUB__->($s.'　', $r->{$_});
　}
}->('', \%h);
__DATA__
フグ亜目
フグ亜目\ハリセンボン科
フグ亜目\ハリセンボン科\ハリセンボン属
フグ亜目\ハリセンボン科\ハリセンボン属\ハリセンボン
フグ亜目\マンボウ科
フグ亜目\マンボウ科\クサビフグ属
フグ亜目\マンボウ科\クサビフグ属\クサビフグ
フグ亜目\マンボウ科\マンボウ属
フグ亜目\マンボウ科\マンボウ属\ウシマンボウ
フグ亜目\マンボウ科\マンボウ属\カクレマンボウ
フグ亜目\マンボウ科\マンボウ属\マンボウ
フグ亜目\マンボウ科\ヤリマンボウ属
フグ亜目\マンボウ科\ヤリマンボウ属\ヤリマンボウ

**デフォルトの名無しさん** · 2020/06/08(月) 01:00:33.65

>>826
Kotlin
https://paiza.io/projects/l95d6I0l3E22Ja0pr50_wg

**デフォルトの名無しさん** · 2020/06/11(木) 14:12:11.60

お題: 年月日の日付が YYYY-MM-DD 形式の文字列(2020-06-11 など)で与えられるので、その日付から月の第何週か、及び第何曜日(第2月曜日、など)かを求める。

月の第何週かを求める場合、週の初めの曜日が日曜日と月曜日で違う答えになるが、これは両方を求めて分かるように二つ出力しても良いし、片方だけでも良い。

ヒント: Kotlin 7 スレの43。ただしこのヒントはほとんど答えなので自分で考えて解きたい人は見ない方が良い。
https://mevius.5ch.net/test/read.cgi/tech/1588748438/43

**デフォルトの名無しさん** · 2020/06/11(木) 15:10:12.05

お題: 表面積が1平米である凸八面体のうち、体積が最大となる立体について、その体積の値を15桁以上の精度で計算せよ。

**デフォルトの名無しさん** · 2020/06/11(木) 18:02:24.18

なぜ凸をつけた？

**デフォルトの名無しさん** · 2020/06/11(木) 18:15:04.87

1÷√√34992

**デフォルトの名無しさん** · 2020/06/11(木) 18:25:55.50

正八面体は最大ではないよ

**デフォルトの名無しさん** · 2020/06/11(木) 18:30:55.39

まぢか

**デフォルトの名無しさん** · 2020/06/11(木) 18:53:27.65

6角柱とたまたま同じになった

**デフォルトの名無しさん** · 2020/06/11(木) 18:59:28.26

あと考えられるのは5角形4枚と4角形4枚か？
面倒で計算したくない

**デフォルトの名無しさん** · 2020/06/11(木) 19:14:23.40

１週間誰も解けなかったらやってみる

**デフォルトの名無しさん** · 2020/06/11(木) 20:45:21.07

冷やっこ食ってて思ったんだけど
一辺が41cmの立方体の
表面積が 41×41×6＝10086cm^2なので
角っこの頂点２ヶ所を
ちょんちょんって切り落とした立体なら八面ではあるって事かな？
三角錐二個の部分の表面積計算めんどくさそう

**デフォルトの名無しさん** · 2020/06/11(木) 22:59:47.78

1÷√√34992＝0.0731152229418051367121788278776110586200…
これより大きいものをとりあえず見つける方針で

**デフォルトの名無しさん** · 2020/06/12(金) 01:51:54.36

0.0743448680932299869829269776944062867445462710353362468

見つかった

**デフォルトの名無しさん** · 2020/06/12(金) 01:53:02.92

形状を求めよ
の方が難しい

**>>843** · 2020/06/12(金) 05:02:56.70

>>843
コナミ音ゲーのビートストリームの箱猫にゃん
みたいのは難しいので中学生が電卓で出来るレベルで
立方体から直方体切り出しの六角柱で考える
https://i.imgur.com/Mjl42Cr.jpg

表面積約1平米の凸六面体の最大体積は？

まず一辺wが41cmの立方体の表面積w×w×6=10086[cm^2]
10086から86を引き算して10000に近づけます
do-whileループでxを求める
do{
a=10086-(x*x*2);
x+=0.00001;
}while(a>10000);
直方体切取幅x=6.557449999920918[cm]
切取直方体上下部分だけの表面積x×x×2=86.0000387051257[cm^2]

凸六面体（六角柱）の表面積=
9999.999961294874[cm^2]≠約1平米

元の立方体の体積41^3=68921[cm^3]

切取った直方体の体積x×x×w=1763.0061705599765[cm^3]

凸六面体（六角柱）の体積は
41^3-(x×x×w)=67157.99382944002[cm^3]
≒0.06715799382944002[m^3]
15桁以上
1割るルートルート34992ってなんだろう？

**デフォルトの名無しさん** · 2020/06/12(金) 05:12:49.16

>>847
訂正、凸六面体じゃなくて凸八面体ね

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

なんでわざわざ凸じゃないのを選ぶ？

**デフォルトの名無しさん** · 2020/06/12(金) 08:46:39.31

最大でもないし

**デフォルトの名無しさん** · 2020/06/12(金) 08:48:13.07

最大じゃないのは私も人の事を言えないので...

>>845が合ってるか間違ってるか
だけ教えていただけると

**デフォルトの名無しさん** · 2020/06/12(金) 08:49:21.18

8面体が解決したら次は9面体

8面も最大の証明は難しいけど

**デフォルトの名無しさん** · 2020/06/12(金) 11:36:58.81

>>851
値は15桁まで合ってると思いますよ。

プログラムのお題なんでコードを載せていただけると幸いですが。
何らか組んだんですよね？

**デフォルトの名無しさん** · 2020/06/12(金) 12:49:26.30

しまった
係数 (1/3) がdoubleのままだったんでdouble精度になってしまった

**デフォルトの名無しさん** · 2020/06/12(金) 12:50:37.41

家に帰ったらdouble版をアップします
ひたすら乱数探索

**デフォルトの名無しさん** · 2020/06/12(金) 15:14:32.23

>>849
えー凸にしたら10面になっちゃう
https://imgur.com/SthLxF6.jpg

**デフォルトの名無しさん** · 2020/06/12(金) 17:08:36.60

>>856
それも凸じゃないだろ

**デフォルトの名無しさん** · 2020/06/12(金) 17:09:01.45

あ、すまん、そういうことかw

**デフォルトの名無しさん** · 2020/06/12(金) 17:33:15.44

´+ω+`

**デフォルトの名無しさん** · 2020/06/12(金) 18:47:56.69

>>826
Ruby で

text = <<'TEXT' # \ を解釈しない
フグ亜目
フグ亜目\ハリセンボン科
フグ亜目\ハリセンボン科\ハリセンボン属
フグ亜目\ハリセンボン科\ハリセンボン属\ハリセンボン
フグ亜目\マンボウ科
フグ亜目\マンボウ科\クサビフグ属
フグ亜目\マンボウ科\クサビフグ属\クサビフグ
フグ亜目\マンボウ科\マンボウ属
フグ亜目\マンボウ科\マンボウ属\ウシマンボウ
フグ亜目\マンボウ科\マンボウ属\カクレマンボウ
フグ亜目\マンボウ科\マンボウ属\マンボウ
フグ亜目\マンボウ科\ヤリマンボウ属
フグ亜目\マンボウ科\ヤリマンボウ属\ヤリマンボウ
TEXT

text.each_line.inject( [ ] ) do | prev, line | # 1行ずつ処理する
line.chomp! # 末尾の改行を削除する
words = line.split( "\\" )

words.zip( prev ).each_with_index do |ary, idx|
if ary[ 0 ] != ary[ 1 ] # 1つ前の要素と、異なっていたら
puts "　　" * idx + ary[ 0 ] # インデント
break
end
end

words
end

**デフォルトの名無しさん** · 2020/06/12(金) 19:40:29.35

>>835
https://ideone.com/Jl6zc2

**デフォルトの名無しさん** · 2020/06/12(金) 19:41:47.16

0.07434486809322999110989619368561104639139465256280075733904874421666917664911993460556242342418454035329397732967289108376252183488881295875664070934129575634827025621624683049348057877455258195173727352356013738382644505372784969303844994902036116602401810098334065219605813513311533511174260033990459633153548001338432276618035960031832154258304694684989102438831680515365572871146875031327576991273892059152446734535128273692394822907010528329236904909366291990338992899295703583566824654618008406013296199027300917662863617783075767465771504776626571228664315942541587700314169499095063254202509152441676340208925338471982724370712866011599553438930619611265418988185782338154673314730488687674137034907531619005947724301353643581914084747714077649289547685930833278053943700887036591347108420548965969583488831450816524757283812848039192907796866305854763801905888387806167102792468923

**デフォルトの名無しさん** · 2020/06/12(金) 20:27:51.93

ルートとか使うと何の値なの

**デフォルトの名無しさん** · 2020/06/12(金) 20:28:24.27

>>861のコードはわかるけど、>>862の桁数はどうやったんだろう

**デフォルトの名無しさん** · 2020/06/12(金) 20:32:58.85

´°Д°｀いったいどんな形なんだ

**デフォルトの名無しさん** · 2020/06/12(金) 20:54:30.86

>>864
scalarを自作8192bit floatライブラリにして計算しました
全部はれなかったのではれたぶんだけ
探索方法は全く同じです
頂点座標は体積最大値の半分くらいの桁数で求まるようです

今回は速度が問題にならなかったので力業のままで

8面体は自力で面の構成を考えたけど
9面以上はコンピューターに考えさせないとキビシイかな？

**デフォルトの名無しさん** · 2020/06/12(金) 21:00:40.76

>>861
お見事です

ま、頂点の座標があったら図を描くのも難しくはないですけどね
http://i.imgur.com/JxLdedg.gif

**デフォルトの名無しさん** · 2020/06/12(金) 21:12:57.61

>>867
はヘェ～スゲーな～
どんな勉強したらこんなの解けるんだろ
理屈の解説は>>861のコードを読み解くしか無いですかね？
頂点を求める考え方や手順の解説お願いします。

**デフォルトの名無しさん** · 2020/06/12(金) 21:20:14.56

6角柱がハズレだとわかった後
6分で思い付いた形
すごいでしょ？
さすが数オリ経験者の東大数学科卒

**デフォルトの名無しさん** · 2020/06/12(金) 22:01:32.78

>>869
>840-841は6分かも知れないけど
出題から3時間後>>837で答えてるのは
答えを知ってたのでは？
何処で習うんだか、というか有名な立体でプラトンとか名前付いてたりするんかな

**デフォルトの名無しさん** · 2020/06/12(金) 22:38:07.50

>>866
＞8面体は自力で面の構成を考えたけど
＞9面以上はコンピューターに考えさせないとキビシイかな？

そうねえ。
面の構成まで計算できたら汎用性があっていいし、何より「本当にコレで最大なのか？」ってとこで悩まなくてもいいですね

**デフォルトの名無しさん** · 2020/06/13(土) 02:39:55.37

>>869
すごいって言われたい症候群
http://mevius.5ch.net/test/read.cgi/tech/1573948822/304,686-687,844,891,892
コンプレックスの裏返しで自演
心療内科行け

**デフォルトの名無しさん** · 2020/06/13(土) 09:06:18.02

私にもファンが出来たか
5か月前の書き込みを覚えていてくれるとは

では次5か月後にあらわれます

**デフォルトの名無しさん** · 2020/06/13(土) 11:29:01.03

このスレにキチガイがいるな

**540** · 2020/06/13(土) 11:52:25.11

>>834 Perl5、こういうﾌﾟﾛｸﾞﾗﾑでは日にち曜日計算を自分で考える前に、いつも使っているﾓｼﾞｭｰﾙの利用をすぐ思い浮かんでしまう

for (qw{02 03 07 10 14 17 21 24 28 30 31}) {
　$date = "2020-05-$_";
　$t = Time::Piece->strptime($date, '%Y-%m-%d');
　$n = 1 + int($t->mday / 7);
　$dw = $n . $t->day(qw{日月火水木金土});
　$t1 = $t - ONE_DAY * ($t->mday - 1);
　$wk = 1 + int(($t->mday + $t1->_wday - 1) / 7);
　print "$date: 第$wk週目の第$dw曜日です。\n";
}

実行結果
~ $ perl 17_834_date_week_day.pl
2020-05-02: 第1週目の第1土曜日です。
2020-05-03: 第2週目の第1日曜日です。
2020-05-07: 第2週目の第2木曜日です。
2020-05-10: 第3週目の第2日曜日です。
2020-05-14: 第3週目の第3木曜日です。
2020-05-17: 第4週目の第3日曜日です。
2020-05-21: 第4週目の第4木曜日です。
2020-05-24: 第5週目の第4日曜日です。
2020-05-28: 第5週目の第5木曜日です。
2020-05-30: 第5週目の第5土曜日です。
2020-05-31: 第6週目の第5日曜日です。

**540** · 2020/06/13(土) 14:11:48.60

>>875
先頭に
use Time::Piece;
use Time::Seconds;
が抜けてたてorz

**デフォルトの名無しさん** · 2020/06/13(土) 20:11:04.26

>>834
Java
https://paiza.io/projects/8qXLWkOSbm0Wi9qDDgyTZA

**デフォルトの名無しさん** · 2020/06/14(日) 16:08:51.18

お題
フィボナッチ数列のN番目の項を求めよ。
ただし非常に大きなNにも対応するために計算量をO(logN)に抑える事。

**デフォルトの名無しさん** · 2020/06/14(日) 16:39:40.62

>>878
最速のフィボナッチ数計算を考える | 雑記帳
https://blog.miz-ar.info/2019/01/fast-fibonacci/

行列でも使えばいいの？

**デフォルトの名無しさん** · 2020/06/14(日) 16:56:58.68

>>878 perl
use bigint;

$ap = $am = $bp = $b0 = 1;
$a0 = $bm = 0;

for($m = $N = $ARGV[0]; $m ;$m /= 2) {
($ap, $a0, $am) = ($ap * $bp + $a0 * $b0, $ap * $b0 + $a0 * $bm, $a0 * $b0 + $am * $bm) if $m % 2;
($bp, $b0, $bm) = ($bp * $bp + $b0 * $b0, $bp * $b0 + $b0 * $bm, $b0 * $b0 + $bm * $bm)
}
print "F_$N = $a0\n";

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/14(日) 17:02:48.95

お題:

https://google.com/robots.txt

をダウンロードするプログラム。

**デフォルトの名無しさん** · 2020/06/14(日) 17:10:27.96

>>878 Perl
# 簡略化
use bigint;

$am = $b0 = 1;
$a0 = $bm = 0;

for($m = $N = $ARGV[0]; $m ;$m /= 2) {
$a0b0 = $a0 * $b0;
$b0b0 = $b0 * $b0;
($a0, $am) = ($a0b0 + $am * $b0 + $a0 * $bm, $a0b0 + $am * $bm) if $m % 2;
($b0, $bm) = ($b0b0 + 2 * $b0 * $bm, $b0b0 + $bm * $bm)
}
print "F_$N = $a0\n";

**デフォルトの名無しさん** · 2020/06/14(日) 17:17:33.08

>>881
Java
https://paiza.io/projects/BGptgg_BylJmzHj6K6Savw

**蟻人間** ◆T6xkBnTXz7B0 · 2020/06/14(日) 17:41:46.79

>>881 C++/Win32
#include <windows.h>
#include <urlmon.h>
#include <stdio.h>

int main(void)
{
CoInitialize(NULL);
HRESULT hr = URLDownloadToFileA(NULL,
"http://google.com/robots.txt";,
"robots.txt",
0,
NULL);
if (SUCCEEDED(hr)) puts("OK");
else puts("No good");
CoUninitialize();
return 0;
}

**デフォルトの名無しさん** · 2020/06/14(日) 17:45:17.97

なおNが大きくなると扱う数の桁数がNに比例して増える
よって>>882の乗算の回数はO(logN)であるが、
計算量はO(logN)に抑えることはできない
おそらくO(N^2)とか、乗算を効率化してもO(N*(logN)^2)あたりではないか

**download.html** · 2020/06/14(日) 17:50:45.87

>>881 HTML
<a href="https://google.com/robots.txt"; download="https://google.com/robots.txt";>https://google.com/robots.txtのダウンロード</a>

**デフォルトの名無しさん** · 2020/06/15(月) 03:05:29.95

お題: (835の類題)表面積が1平米である九面体のうち、体積が最大となる立体について、その体積の値を15桁以上の精度で計算せよ。

※ 立体の形状をあらかじめ限定しないアルゴリズムがあればなお良い。

**デフォルトの名無しさん** · 2020/06/15(月) 03:19:03.34

>>884 標準コマンド
Windows 10 だったら多くのOSと同様に標準で curl が使えるから
curl -LO https://google.com/robots.txt

**デフォルトの名無しさん** · 2020/06/15(月) 03:46:31.62

>>881

#!/bin/sh

wget https://google.com/robots.txt

**デフォルトの名無しさん** · 2020/06/15(月) 03:55:40.44

>>881
Kotlin script

java.net.URL("https://google.com/robots.txt";).openConnection().inputStream.bufferedReader().lineSequence().forEach(::println)

script じゃなくてもほぼ同じ。main() の中に入れるだけ。

**デフォルトの名無しさん** · 2020/06/15(月) 06:06:01.39

>>881 Ruby
require 'open-uri'
$><<open('https://google.com/robots.txt').read

**540** · 2020/06/15(月) 21:52:43.85

>>881 Perl5

use LWP::Simple;
use LWP::Protocol::https;
print mirror('https://google.com/robots.txt', './robots.txt'), "\n";

実行結果
~ $ perl 17_881_https_mirror_robots_txt.pl
200
~ $ ls -o robots.txt
-rw-r--r-- 1 user 7035 Jun 9 03:00 robots.txt

**デフォルトの名無しさん** · 2020/06/19(金) 07:19:50.29

>>659
https://ideone.com/dJhsj8
10日以上回してようやく229の次が出た
先頭から30個は以下。A098044 の記載とも一致
7, 13, 19, 37, 43, 79, 163, 223, 229, 608981812891, 608981812951, 608981812993, 608981813507, 608981813621, 608981813819, 608981813837, 608981813861, 608981813929, 608981813941, 608981814019, 608981814143, 608981814247, 608981814823, 608981814829, 608981815027, 608981815051, 608981815069, 608981816791, 608981816803, 608981816893

**デフォルトの名無しさん** · 2020/06/19(金) 07:48:33.21

>>893
VPSは１０コアぐらい使いましたか？
それにしても興味深い

**デフォルトの名無しさん** · 2020/06/19(金) 19:45:31.37

>>894
いや、コード見てもらったらわかるようにシングルスレッドです。
工夫して並列化したらコアの数だけ高速化できる余地はあるので
誰かやってみてほしいですよ。

計算を続けると6千億台に1万個以上解が見つかります。

**デフォルトの名無しさん** · 2020/06/19(金) 20:32:41.02

連続した数調べるならエラトステネスの篩したほうが速いんじゃないの

**デフォルトの名無しさん** · 2020/06/19(金) 20:58:48.38

だね
シングルスレッドで1時間くらいで>>893

**デフォルトの名無しさん** · 2020/06/19(金) 21:30:26.89

篩を使った方法は>>645(修正>>695)
こんなん１時間でいける？

**デフォルトの名無しさん** · 2020/06/19(金) 21:45:20.37

>>896
えっともしかしてそれは
6000億ちょっとの要素数の配列を扱うと言ってますか？

**デフォルトの名無しさん** · 2020/06/19(金) 21:47:33.79

>>899
言ってません

**デフォルトの名無しさん** · 2020/06/19(金) 22:02:05.91

6000億ちょっとの要素数の配列を扱います

**デフォルトの名無しさん** · 2020/06/19(金) 23:21:23.40

エラとステネスの篩はビット単位で判定させれば8分の1の要素数で済む
もっと効率化すれば32分の1くらいになったはず
まぁそれでも億単位の要素を使うわけだけど

◆QZaw55cn4c · 2020/06/19(金) 23:38:14.59

私なら篩を一定の大きさに制限して窓として扱い、
試し割り法と併用しますね
https://ideone.com/6Ww9nq
https://mevius.5ch.net/test/read.cgi/tech/1580471646/245

**デフォルトの名無しさん** · 2020/06/20(土) 01:08:51.20

>>893
そのまま解を検索してみたけど、610968213803 から先は解が見つからなくなった
次はいくつになるかをお題にしてもいいかもね

**デフォルトの名無しさん** · 2020/06/20(土) 01:10:45.21

ほとんどの6000億幅が取れない人は、"区間ふるい"で.

「1億の6千回」か「10億の600回」でやるのが現実的(平方分割的感性が必要？)
実際に1億でやってみたら、２時間ほどかかった。

ideonは時間の関係で1億配列にして、
最初と(ログから)最後の1億をやっている。
https://ideone.com/Ikn4qW

**デフォルトの名無しさん** · 2020/06/20(土) 01:11:48.05

>>896
いや
篩も試してみたけどやっぱりBPSW法が速かったよ
それと単純に割り算の繰り返しで試すより面白そうだったってのもある