プログラミングのお題スレ Part17
レス数が950を超えています。1000を超えると書き込みができなくなります。
最大じゃないのは私も人の事を言えないので...
>>845が合ってるか間違ってるか
だけ教えていただけると 8面体が解決したら次は9面体
8面も最大の証明は難しいけど >>851
値は15桁まで合ってると思いますよ。
プログラムのお題なんでコードを載せていただけると幸いですが。
何らか組んだんですよね? しまった
係数 (1/3) がdoubleのままだったんでdouble精度になってしまった 家に帰ったらdouble版をアップします
ひたすら乱数探索 >>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 0.07434486809322999110989619368561104639139465256280075733904874421666917664911993460556242342418454035329397732967289108376252183488881295875664070934129575634827025621624683049348057877455258195173727352356013738382644505372784969303844994902036116602401810098334065219605813513311533511174260033990459633153548001338432276618035960031832154258304694684989102438831680515365572871146875031327576991273892059152446734535128273692394822907010528329236904909366291990338992899295703583566824654618008406013296199027300917662863617783075767465771504776626571228664315942541587700314169499095063254202509152441676340208925338471982724370712866011599553438930619611265418988185782338154673314730488687674137034907531619005947724301353643581914084747714077649289547685930833278053943700887036591347108420548965969583488831450816524757283812848039192907796866305854763801905888387806167102792468923 >>861のコードはわかるけど、>>862の桁数はどうやったんだろう >>864
scalarを自作8192bit floatライブラリにして計算しました
全部はれなかったのではれたぶんだけ
探索方法は全く同じです
頂点座標は体積最大値の半分くらいの桁数で求まるようです
今回は速度が問題にならなかったので力業のままで
8面体は自力で面の構成を考えたけど
9面以上はコンピューターに考えさせないとキビシイかな? >>861
お見事です
ま、頂点の座標があったら図を描くのも難しくはないですけどね
http://i.imgur.com/JxLdedg.gif >>867
はヘェ〜スゲーな〜
どんな勉強したらこんなの解けるんだろ
理屈の解説は>>861のコードを読み解くしか無いですかね?
頂点を求める考え方や手順の解説お願いします。 6角柱がハズレだとわかった後
6分で思い付いた形
すごいでしょ?
さすが数オリ経験者の東大数学科卒 >>869
>840-841は6分かも知れないけど
出題から3時間後>>837で答えてるのは
答えを知ってたのでは?
何処で習うんだか、というか有名な立体でプラトンとか名前付いてたりするんかな >>866
>8面体は自力で面の構成を考えたけど
>9面以上はコンピューターに考えさせないとキビシイかな?
そうねえ。
面の構成まで計算できたら汎用性があっていいし、何より「本当にコレで最大なのか?」ってとこで悩まなくてもいいですね 私にもファンが出来たか
5か月前の書き込みを覚えていてくれるとは
では次5か月後にあらわれます >>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日曜日です。 >>875
先頭に
use Time::Piece;
use Time::Seconds;
が抜けてたてorz お題
フィボナッチ数列のN番目の項を求めよ。
ただし非常に大きなNにも対応するために計算量をO(logN)に抑える事。 >>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"; >>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"; >>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;
} なおNが大きくなると扱う数の桁数がNに比例して増える
よって>>882の乗算の回数はO(logN)であるが、
計算量はO(logN)に抑えることはできない
おそらくO(N^2)とか、乗算を効率化してもO(N*(logN)^2)あたりではないか お題: (835の類題)表面積が1平米である九面体のうち、体積が最大となる立体について、その体積の値を15桁以上の精度で計算せよ。
※ 立体の形状をあらかじめ限定しないアルゴリズムがあればなお良い。 >>884 標準コマンド
Windows 10 だったら多くのOSと同様に標準で curl が使えるから
curl -LO https://google.com/robots.txt >>881
Kotlin script
java.net.URL("https://google.com/robots.txt").openConnection().inputStream.bufferedReader().lineSequence().forEach(::println)
script じゃなくてもほぼ同じ。main() の中に入れるだけ。 >>881 Ruby
require 'open-uri'
$><<open('https://google.com/robots.txt').read >>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 >>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 >>893
VPSは10コアぐらい使いましたか?
それにしても興味深い >>894
いや、コード見てもらったらわかるようにシングルスレッドです。
工夫して並列化したらコアの数だけ高速化できる余地はあるので
誰かやってみてほしいですよ。
計算を続けると6千億台に1万個以上解が見つかります。 連続した数調べるならエラトステネスの篩したほうが速いんじゃないの 篩を使った方法は>>645(修正>>695)
こんなん1時間でいける? >>896
えっともしかしてそれは
6000億ちょっとの要素数の配列を扱うと言ってますか? エラとステネスの篩はビット単位で判定させれば8分の1の要素数で済む
もっと効率化すれば32分の1くらいになったはず
まぁそれでも億単位の要素を使うわけだけど >>893
そのまま解を検索してみたけど、610968213803 から先は解が見つからなくなった
次はいくつになるかをお題にしてもいいかもね ほとんどの6000億幅が取れない人は、"区間ふるい"で.
「1億の6千回」か「10億の600回」でやるのが現実的(平方分割的感性が必要?)
実際に1億でやってみたら、2時間ほどかかった。
ideonは時間の関係で1億配列にして、
最初と(ログから)最後の1億をやっている。
https://ideone.com/Ikn4qW >>896
いや
篩も試してみたけどやっぱりBPSW法が速かったよ
それと単純に割り算の繰り返しで試すより面白そうだったってのもある お題: 整数を記入可能な空欄の横3×縦3マスが格子状に並んでいる。左からh+1番目で上からk+1番目のマスに書かれている整数を「C(h, k)」と表すことにする。
ただし、空欄のときやマスの範囲外の場合は「C(h, k)」はゼロとする。
C(h, k)の位置に整数mを鉛筆で記入する手続きを「P(h, k, m)」と表すことにする。
手続きにおいてはC(h, k)は入れ子を許し、演算子+と-は通常の足し算と引き算を意味するものとする。
(1) 次の手続きのリストを順に実行し、結果を表示せよ。
P(0, 0, 1), P(1, 0, C(0, 0)+2), P(2, 0, C(1, 0)-1),
P(C(1, 0)-1, C(0, 0)+C(2, 0)-2, 5), P(2, 2, C(2, 1)-2).
空欄はゼロで表すものとする。 >>906
前言撤回。>>905のように篩でやったほうが格段に速い。そりゃそうだよね。 お題: CSVをパースせよ
例
入力:
1,abc,"def
ghi",4
2,abc,"def
ghi",4
出力:
[[1, "abc", "def\nghi", 4], [2, "abc", "def\nghi", 4]] >>909 Ruby
require 'csv'
$><<CSV.parse($<.read) >>887
60面体まで出来ました
60面体にもなると
5角形と6角形だけからなるものだけでも
極大点が10個以上にもなって
本当に最大が探せてるのかいまいち確信は持てないんですが...
収束は速いので形の概要が分かれば高精度化は簡単です
14面体以上は全て
5角形が12面、残りが6角形の多面体が最大です お題:多面体の面の構成を、以下のルールに従う面列で表現する。
頂点のリストで表現された面の集合(例:正四面体なら {[1,2,3],[1,3,4],[1,4,2],[4,3,2]})で多面体が与えられたとき、
それに対応する多面体の面列表現を出力せよ。
1. 3〜9角形の面の表現を t,q,p,x,h,o,n とする。10角形以上の面がある場合はエラーとしてよい。
2. 適当な面を選び、その面の表現をまず面列の先頭に置く。
3. 選ばれていない面が残っている場合は、それまでの面列の後にハイフンを置き、続けて、既に選ばれている面の外周に接する面を順に列挙する。
選ばれていない面がなくなるまでこのプロセスを繰り返す。
4. ハイフンの前後の面列の、それぞれ先頭の文字に対応する面は、互いに接しているものを選ぶ。
5. ハイフンで区切られた面列に、pppp のように繰り返し部分がある場合は、p4 のように繰り返し回数を数字で表す。
繰り返しの対象が単一の面でない場合は対象の前後に括弧を置く。この括弧は入れ子にできる。(t(xp)2)2 = txpxptxpxp
(この括弧の中にハイフンを含めることはできない)
6. 可能な表現が複数ある場合は、表現に含まれるアルファベットの個数が最も少ない表現を採用する。
例) 三角錐: t-t3, 六角柱: x-q6-x, 正十二面体: p-p5-p5-p
切頂二十面体(サッカーボール形): p-x5-(xp)5-(px)5-x5-p >>912
そのルールでは任意の多面体を表現出来ないと思うんですが だね
トポロジー的に球面限定だし
その限定でも
4. が複数可能性があるから
>>912の表現では一意に決まらない >>909
これってエクセープサれた改行文字で渡されるのが標準的じゃない?
規格としてあるのかわからないけど 最初から区切り文字をエスケープする方法まで決まっていたらxmlなんて出てこなかったかもな >>911
面の数が可変のものを作ってしまったか・・・
20面体、32面体、42面体、72面体あたりが面白い形状だと思うよ 14面体以上は全て5角形6角形のみで
あまり面白くない
対称性が良いか悪いかだけ
現在82面体
手抜きで一部計算オーダーが非常に悪い所があって進まなくなってきた
(面数の5乗オーダー)
100面まで到達するにはコードを変えないとダメかも 対称性が良ければ良いほど最大解となるかというとそうではない
八面体の解をみてもそれはわかるし
四十二面体の対称性の高いのは切稜十二面体だと思うが、解はそれとは違う
面の数の5乗というのは総当たりをしているのかな
高速化するには焼きなまし法等を使うといいのかと思うがどうやるんだろう 1面少ない体積最大の立体から1面足した形を探索
なので最大が突然変異すれば最大を見逃すかもしれない
ランダムよりは効率が良さそう
構造計算のオーダーが4乗
ここはおそらく2乗くらいに減らせると思う
極大の収束は良くある微分を使うやつ
収束時も毎回構造計算してるけど
(ほぼ)同じ構造であるのでオーダーは1乗に減らせる 形自体はもう面白くないので
改善はしないかも
むしろ体積が小さい極大の方がいろんな形になって面白い >>921
数値はこれと同じ?
9 0.076898933926867766
10 0.078734752898039751
11 0.080055026399577983
12 0.08168837182418255218049
13 0.082432267303420834
14 0.083349245941114841
16 0.084742718358283536
17 0.085264872589057683
20 0.086626966830007951
32 0.089493100466131958
33 0.089603827451613424
42 0.090574499972086386 お題: 3つ以上の言語で動く"Hello World"(標準出力 改行有無不問) >>925
#! /usr/bin/perl
print "Hello World";
#!/usr/bin/python
print "Hello World"
<?php
print "Hello World";
?> >>925
$ cc hw.c && ./a.out
Hello, World!
$ bash hw.c
Hello, World!
$ perl hw.c
Hello, World!
$ cat hw.c
#if 0 /*
/ 2>\/dev\/null
print() {
echo Hello, World!
}
2>\/dev\/null /;
#*/
#endif
#include <stdio.h>
#define print main() { printf ("Hello, World!\n"); } char *p =
print "Hello, World!\n"; >>925 面白いね。俺の回答はつまらんけど
$ cat hello_world
print("Hello World!")'
$ruby hello_world
Hello World!
$perl hello_world
Hello World!
$python hello_world
Hello World! お題: CSV を、JSON に変換せよ
あ,い
1,2
[
["あ","い"],
["1","2"]
]
Convert CSV to JSON
https://www.convertcsv.com/csv-to-json.htm あ,い
1,2
2, 3
↓
[
{ "あ": 1, "い": 2 },
{ "あ": 2, "い": 3 }
]
とかじゃないんだな
❯ node
> const csv2json = text => text.trim().split('\n').map(line => line.split(','))
undefined
> text = `
あ,い
1,2`
> csv2json(text)
[ [ 'あ', 'い' ], [ '1', '2' ] ] >>930 Ruby
require 'csv'
require 'json'
puts CSV.parse($<.read).to_json # => [["あ","い"],["1","2"]] >>930
sh というか awk というか
( echo あ,い ; echo 1,2 ) | awk -F, 'BEGIN{print "["} { printf "[\"%s\"", $1; for (i = 2; i <= NF; i++) printf ",\"%s\"", $i; print "]"} END{print "]"}' >>930
>>934はレコードごとのカンマが抜けてた。ちょっと修正。
( echo あ,い ; echo 1,2 ) | awk -F, 'BEGIN{print "["} {if(NR>1)print","; printf "[\"%s\"", $1; for (i = 2; i <= NF; i++) printf ",\"%s\"", $i; printf "]"} END{print "\n]"}' >>909
>>930
何気にPythonの仕様が気になったので。
import csv
import io
import json
with io.StringIO("""
1,abc,"def
ghi",4
2,abc,"def
ghi",4
""") as f:
__print([r for r in csv.reader(f) if r])
with io.StringIO("""
あ,い
1,2
""") as f:
__print(json.dumps([r for r in csv.reader(f) if r], ensure_ascii=False)) お題:"チューリングマシン"を使って一行CSVパーサーを作る
チューリングマシンも自作すること
【CSVの仕様】
・カラム値は「a」、「,」、「"」、の組み合わせ、または「」(値無し)しか取らない
・カラム値に「,」が含まれる場合は「"」で囲む必要がある、含まれない場合はどっちでもいい
・カラム値に「"」が含まれる場合は「\"」とエスケープする
入力例) "aa,a",aa,,"aa\"aa"
チューリングマシンからの出力例、Pythonの場合) [None, 'a', 'a', ',', 'a', None, 'X', 'a', 'a', 'X', 'X', None, 'a', 'a', None, '"', 'a', 'a', None]
(Noneは無視、'X'は区切りポイントを表している)
【開発規模】Pythonで15行程度 チューリングマシンなのにPythonとは?
意味不明だ 正方形の面積を求める問題に対して
この世に正方形なんて存在しないから求められません
とか言っちゃう系? 比喩がとんちんかんだぞ
今回は「チューリングマシンも自作すること」って明言されてますやん 自作チューリングマシンを Pythonの改行15回程度で作れって?神だわ 動作を規定すればよくて
実際に動く機械を作る必要はない と思ったが
Pythonでエミュレーターを作れって事か? 別に無限にメモリがある必要は無いけど
実際に動くどんなプログラムにも限界がある >>946
メモリが有限だったらチューリングマシンとはいいましぇん>< しかもCSVパーサーって
つまりはPython15行でExcelを再現しろと? >>938の生成言語
CSV → COL | COL ',' CSV
COL → '"' WCOM '"' | NCOM
WCOM → ε | 'a' WCOM | ',' WCOM | '\' '"' WCOM
NCOM → ε | 'a' NCOM | '\' '"' NCOM
状態セル10個前後でいけそう >>938 Perl 状態セル4個
$ cat input
"aa,a",aa,,"aa\"aa"
$ perl odai938 < input
[None, 'a', 'a', ',', 'a', None, 'X', 'a', 'a', 'X', 'X', None, 'a', 'a', None, '"', 'a', 'a', None]
$ cat odai938
%maton = (
'N' => { 'a' => ['N', q('a')] , ',' => ['N', q('X')] , '"' => ['Q', q(None)] , '\\' => ['E', q(None)] },
'Q' => { 'a' => ['Q', q('a')] , ',' => ['Q', q(',')] , '"' => ['N', q(None)] , '\\' => ['F', q(None)] },
'E' => { 'a' => ['N', q('a')] , ',' => ['N', q(',')] , '"' => ['N', q('"') ] , '\\' => ['N', q(None)] },
'F' => { 'a' => ['Q', q('a')] , ',' => ['Q', q(',')] , '"' => ['Q', q('"') ] , '\\' => ['Q', q(None)] },
);
@res = ();
$state = 'N';
for $c ( split //, <> ) {
($next, $out) = @{$maton{$state}{$c}};
last unless defined $next;
push @res, $out;
$state = $next;
}
printf "[%s]\n", join ', ', @res; レス数が950を超えています。1000を超えると書き込みができなくなります。