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

**デフォルトの名無しさん** · 2020/07/14(火) 13:53:46.47

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

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

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

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

**デフォルトの名無しさん** · 2020/10/11(日) 17:59:24.67

>>629
Perl
https://paiza.io/projects/wfQXB4arB2zbu4DNYEiU5w

eval が使えて楽なので、つい。

**デフォルトの名無しさん** · 2020/10/11(日) 18:12:54.16

>>602
Java
https://paiza.io/projects/2bTy71EhJJcVemamRx8IPA

**デフォルトの名無しさん** · 2020/10/11(日) 19:16:59.76

>>635
Kotlin
https://paiza.io/projects/O1eUwROYKRQepDq3N-VyUA

paiza.io はHTMLで出力させてブラウザに表示できるのでカラーバーも出してみました。

**デフォルトの名無しさん** · 2020/10/11(日) 19:36:20.38

これは良いこと聞いた

**デフォルトの名無しさん** · 2020/10/11(日) 19:48:34.94

>>640 Common Lisp (SBCL, 処理系依存)
https://ideone.com/nImAgj

**デフォルトの名無しさん** · 2020/10/11(日) 20:05:13.72

>>663
かわいい

**デフォルトの名無しさん** · 2020/10/11(日) 20:30:57.80

>>640 brainfuck

>+++++++[<+++++++>-]<+.--.++.--.+.-.+.-.

◆QZaw55cn4c · 2020/10/11(日) 23:30:32.44

>>605
C(C++) https://mevius.5ch.net/test/read.cgi/tech/1434079972/70

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/11(日) 23:57:26.40

お題: 整数0～10の一様分布の乱数を百回発生させて、その統計結果を基に各整数を横軸、回数を縦軸とする棒グラフ(ヒストグラム)を描画し(棒の最大幅32字)、ほぼ一様であることを確かめなさい。

**デフォルトの名無しさん** · 2020/10/11(日) 23:59:14.77

paiza.io は、text, html, json で出力できるのか

html も、部分的なタグだけを書いても動くのか

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/12(月) 00:01:01.03

>>667 訂正。

(誤)各整数を横軸、回数を縦軸とする棒グラフ
(正)各整数を縦軸、回数を横軸とする棒グラフ

**デフォルトの名無しさん** · 2020/10/12(月) 01:06:07.97

haskell

import Data.List
import System.Random

trial records = do
x <- randomRIO (0, 1) :: IO Double
return $ (truncate $ x * 11) : records

experience = iterate (>>= trial) (return [0..10])
showrecord (x,y) = let
clsstr = id
$ reverse
$ take 2
$ ( ++ (repeat '0'))
$ reverse
$ show x
barstr = replicate y '*'
in concat [ clsstr," ", barstr,"\n"]

main = do
records <- experience !! 100
putStr $ id
$ concat
$ map showrecord
$ map (\x -> (head x, length $ tail x))
$ group
$ sort
$ records

**デフォルトの名無しさん** · 2020/10/12(月) 01:14:47.03

>>640
Kotlin
https://paiza.io/projects/32aoQfPBUWZ0tnk2E5fg2g

ちょっと変わった感じにしようと思ってこうした。

**デフォルトの名無しさん** · 2020/10/12(月) 08:18:41.49

>>667
二項分布になるんじゃないかなあ

**デフォルトの名無しさん** · 2020/10/14(水) 15:32:12.46

>>649 オーディオ系はシステム依存度が高いから余りサンプルが出てこないね。　WindowsやMacOSなら簡単そうだが
1kHz の音源ファイルを使うのは違反かな？

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/14(水) 17:27:12.50

>>649 C/C++ with Win32
#include <windows.h>
int main(void)
{
Beep(1000, 1000);
}

**デフォルトの名無しさん** · 2020/10/14(水) 17:57:41.59

答え出てきても試しようがない

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/14(水) 18:19:22.34

>>649 Python with Windows
import winsound
winsound.Beep(1000, 1000);

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/14(水) 18:24:27.45

>>649 Python with Mac
import os
os.system('play -n synth 1 sin 1000');

**デフォルトの名無しさん** · 2020/10/14(水) 18:36:34.21

シェル経由でSoX playコマンド呼んでるだけやんけ
Pythonもmacも関係ねぇわ

**デフォルトの名無しさん** · 2020/10/14(水) 19:17:28.35

OSにつく前置詞はwithじゃなくてonだろハゲ

**デフォルトの名無しさん** · 2020/10/14(水) 19:30:29.09

試すと音が出て煩そうだから試す気になれない。

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/14(水) 19:50:02.82

テレビの自主規制音と同じ音だよ。
https://ja.wikipedia.org/wiki/%E8%87%AA%E4%B8%BB%E8%A6%8F%E5%88%B6%E9%9F%B3

**デフォルトの名無しさん** · 2020/10/15(木) 00:45:55.69

>>678 iOS じゃ、play コマンドがなさそうなんだよな。

**デフォルトの名無しさん** · 2020/10/15(木) 01:18:37.79

結局このスレで盛り上がるのは環境によらず言語によらず、でも多少は骨のあるようなお題なんだよな
へぇその言語ならそんなにキレイに書けるのかとか普段自分が使ってない言語ではどうやって実現できるのか眺めたりするのが面白い
音鳴らすなんてお題として面白いわけがない

**デフォルトの名無しさん** · 2020/10/15(木) 01:40:28.46

ヒント：出題者

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

蟻人間さんへのお題：
文字列を受け取って、RFC 5322にしたがってSubjectヘッダを出力してください(MUST、SHOULDを満たすこと)
文字コードはUTF-8、エンコーディングはBエンコードを使用してください

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/15(木) 13:28:00.67

>>685 Python
text = 'Re: Re: This is a subject text. This is a subject text. This is a subject text.'
if text[0:4] == 'Re: ':
＿＿while text[0:4] == 'Re: ':
＿＿＿＿text = text[4:]
＿＿text = 'Re: ' + text
import base64
text = 'Subject: ' + base64.b64encode(text.encode('utf-8')).decode('ascii');
ary = []
while len(text) > 78:
＿＿ary.append(text[0:78])
＿＿text = text[78:]
ary.append(text)
import platform
if platform.system() == 'Windows':
＿＿for item in ary:
＿＿＿＿print(item)
else:
＿＿for item in ary:
＿＿＿＿print(item, end="\r\n")

**蟻人間** ◆T6xkBnTXz7B0 · 2020/10/15(木) 13:49:36.74

>>686 訂正。

(誤) while len(text) > 78:
(正) while len(text) >= 78:

**デフォルトの名無しさん** · 2020/10/15(木) 14:30:46.93

>>686
いろいろ違うので・・・メールのソース眺めたが良いです

**デフォルトの名無しさん** · 2020/10/15(木) 19:14:52.90

頭の悪いお題を連投してるやつはやっぱり頭悪いんだな

**デフォルトの名無しさん** · 2020/10/15(木) 19:22:03.88

そんなことないと思うけどな
あまり他人の悪口言わない方が良いよ

**デフォルトの名無しさん** · 2020/10/15(木) 19:28:33.00

くっそみたいなお題延々と出し続けてウザいじゃん

**デフォルトの名無しさん** · 2020/10/15(木) 19:29:59.48

>>629なんかはいいと思うんだけどな
eval使えないHaskellerにはいい練習問題になってる

**デフォルトの名無しさん** · 2020/10/15(木) 19:38:14.90

ウェブアプリつくれだのGUIアプリつくれだの酷い事散々言いまくってるじゃん
◆QZaw55cn4cの方がまだマシ

**デフォルトの名無しさん** · 2020/10/15(木) 19:40:03.02

まぁ音出せとかカラーバー出せとかはクソお題だな

**デフォルトの名無しさん** · 2020/10/16(金) 19:48:26.76

お題
自然数nを 2進数にした時、 1が連続する個数の最大値を求める。
例
20201016 -> 5

**デフォルトの名無しさん** · 2020/10/16(金) 20:25:28.88

haskell

import Data.List

maxConsOne = maximum . map length . filter ((==1).head) . group . unfoldr (\n ->if n == 0 then Nothing else Just (mod n 2, div n 2))

main = print $ maxConsOne $ 20201016
----
5

**デフォルトの名無しさん** · 2020/10/16(金) 21:03:20.03

>>695 ruby
f = ->(n) {n.to_s(2).chars.chunk(&:itself).map {|k, a| k == '1' ? a.size : 0}.max}

**デフォルトの名無しさん** · 2020/10/16(金) 21:35:06.98

>>695 JavaScript
const f = (n = 0) => Math.max(...(n.toString(2).match(/1+/g) ?? ['']).map(m => m.length))

console.log(f(20201016)) //=> 5

**253** · 2020/10/16(金) 23:46:14.30

>>695 Perl5

use List::Util max;
print max map{length} sprintf('%b', 20201016) =~ /(1+)/g;

実行結果
$ perl 18_695_bin1.pl
5

**デフォルトの名無しさん** · 2020/10/17(土) 00:12:43.05

>>695 Ruby

p 20201016.to_s(2).split("0").max.size

**デフォルトの名無しさん** · 2020/10/17(土) 00:24:53.85

haskell
無名関数というインチキ使ってしまったので再挑戦バージョン

import Data.List

main = print $ maximum . map sum . group . map (flip mod 2) . takeWhile (/=0) . iterate (flip div 2) $ 20201016

◆QZaw55cn4c · 2020/10/17(土) 00:46:45.25

>>695
https://ideone.com/hQRZsi

**デフォルトの名無しさん** · 2020/10/17(土) 01:58:38.41

>>695 Ruby
n = 20201016
nl = 0
while n > 0
nl = ((nc = (n.odd?)? nc + 1 : 0) > nl)? nc : nl
n >>= 1
end
puts nl

**デフォルトの名無しさん** · 2020/10/17(土) 01:59:32.64

>>702
125の時1になった

**デフォルトの名無しさん** · 2020/10/17(土) 02:41:53.51

これだからQZ

◆QZaw55cn4c · 2020/10/17(土) 02:53:09.41

>>695
C++: https://ideone.com/7H5HBx

>>702 にはロジックに欠陥がありました

>>704
ご指摘ありがとうございます。

**703** · 2020/10/17(土) 09:22:46.57

>>704
自分の >>703 も
- nl = 0
+ nc = nl = 0

**デフォルトの名無しさん** · 2020/10/17(土) 12:19:37.13

>>695 C++
https://wandbox.org/permlink/udZVzwy6pQKbcPa5

**デフォルトの名無しさん** · 2020/10/17(土) 13:46:27.15

>>695 C
https://ideone.com/DwvfL2

**デフォルトの名無しさん** · 2020/10/17(土) 18:57:16.23

>>695 Ruby
n = 20201016
i = 0
while n != 0
n &= n << 1
i += 1
end
puts i

**デフォルトの名無しさん** · 2020/10/17(土) 19:11:03.54

>>710
うまいやり方だな。それには気付かなかった。

**デフォルトの名無しさん** · 2020/10/17(土) 19:29:10.40

>>695
Kotlin script

2進数文字列に変換してから 0 を区切り文字にして分割してその中の一番長い文字列の長さを出すというやり方。

fun max1(n: Int) = n.toString(2).split('0').map { it.length }.max()

println(max1(20201016))

これで 5 が出る。

同じやり方を perl5 で(特にモジュール使わずに)やるとこんな感じ。

perl -e '$n=20201016; $m = (sort { $b <=> $a } map { length($_) } (split /0+/, sprintf("%b", $n)))[0]; print "$m\n"'

同じく5が出る。

**デフォルトの名無しさん** · 2020/10/17(土) 20:05:12.33

>>695 JavaScript

let num = 20201016
let ans = num.toString(2).split`0`.sort().pop().length

console.log(ans) //=> 5

**デフォルトの名無しさん** · 2020/10/17(土) 20:50:55.84

>>695 octave
https://ideone.com/dO1WZp
f = @(n) max(cellfun(@numel, strsplit(dec2bin(n), '0')));
f(20201016)

**デフォルトの名無しさん** · 2020/10/17(土) 20:54:06.81

>>695
Perl
https://paiza.io/projects/VO9_CFHPOVbO5iI8uDc__A

**デフォルトの名無しさん** · 2020/10/17(土) 23:22:21.08

いつもRubyの人からなにかしら学ぶことがある

**デフォルトの名無しさん** · 2020/10/18(日) 00:42:32.49

ハトにエサを与えないで下さい

**デフォルトの名無しさん** · 2020/10/18(日) 01:06:12.01

QZは自分の頭が悪いのをわざわざ周知してもらうためにクソコードを晒すのか
ドMだな

**デフォルトの名無しさん** · 2020/10/18(日) 03:42:01.18

>>710
なるほどな、感心した
n>>1でもいいわけだな

**デフォルトの名無しさん** · 2020/10/18(日) 04:21:53.35

>>700 のほうが速いと思うが

**デフォルトの名無しさん** · 2020/10/18(日) 04:36:21.54

rubyよく分からないけど>>710の方が100倍くらいは速そうに思える

**デフォルトの名無しさん** · 2020/10/18(日) 05:53:44.51

うん。速そうだね。数値計算しかしてないし。

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

20201016みたいに桁数が少ないならそうだけど桁数が増えるとルーブが答え×桁数回だからアルゴリズムそのものは遅い
桁数が少ないうちは内部関数しか使ってない利点で早いとは思うけど

**デフォルトの名無しさん** · 2020/10/18(日) 12:07:22.14

>>723

>>700,710 をベンチマークで比較してみた
https://ideone.com/t0ZYLg
https://ideone.com/DUEWmC
20201016 では 5~6倍
巨大な数だと 7000倍以上

◆QZaw55cn4c · 2020/10/18(日) 15:53:00.51

>>693
残念ながら、私の出題もイマイチ人気がありませんでした‥‥

**デフォルトの名無しさん** · 2020/10/18(日) 21:20:25.30

>>685
Java
https://paiza.io/projects/frlfP3US1doeAO_HZBnLWQ

**デフォルトの名無しさん** · 2020/10/18(日) 21:43:52.67

>>669
Java
https://paiza.io/projects/0xfQrdW_qpJoUbOO6ux6wg

**デフォルトの名無しさん** · 2020/10/18(日) 22:06:12.69

>>667
Java
https://paiza.io/projects/x3j0kbjqsSoVGNJX5qun2g

**デフォルトの名無しさん** · 2020/10/18(日) 22:14:01.64

>>695
Java
https://paiza.io/projects/Eo4-fSXmdnoqjyHOMsi92Q

**デフォルトの名無しさん** · 2020/10/18(日) 22:45:56.62

>>724
x = n
while x != 0
x &= x << 1
i += 1
end
ruby知らないから適当に書いたけど
こんな感じにしないと二回目以降nが0でやってるっぽい
こうすると22倍くらいの差だった

**デフォルトの名無しさん** · 2020/10/18(日) 22:54:22.98

連続してるビット数が多いときはこの方が速くなったけど
すごいごちゃごちゃになった

i = 0
s = 1
x = n
while x != 0
t = x & (x << s)
if t == 0 then
if s == 1 then
i += 1
break
end
s = 1
else
x = t
i += s
s = (s << 1)
end
end

**デフォルトの名無しさん** · 2020/10/18(日) 22:56:30.98

>>724
GJ

**デフォルトの名無しさん** · 2020/10/19(月) 06:30:07.92

>>695
echo "obase=2; 20201016" | bc | sed 's/./&\n/g' | uniq -c | perl -lane 'if ($M < $F[0]) {$M = $F[0]}; END{print $M}'
5

**724** · 2020/10/19(月) 15:26:08.35

>>730
ご指摘の通りでした
https://ideone.com/DUEWmC
を修正しました

>>710 は 1の連続が多いと弱いですね

**デフォルトの名無しさん** · 2020/10/19(月) 17:41:48.20

お題
2^(1/3)を百桁表示せよ

出力
1.2599210498948731647672106072782283505702514647015079800819751121552996765139594837293965624362550941543102560356156652593990240

参考　
a(n+3) = 3*a(n+2) + 3*a(n+1) + an
a(1) = 3, a(2)=3, a(3)=15
で決まる数列について
min (a(n)/a(n+1), a(n+1)/a(n+2), a(n+3)/a(n+3) )
. < 2^(1/3)-1
. max (a(n)/a(n+1), a(n+1)/a(n+2), a(n+2)/a(n+3) )
で lim a(n)/a(n+1) = 2^(1/3)-1となる

**デフォルトの名無しさん** · 2020/10/19(月) 21:36:44.38

>>733
echo "obase=2; 20201016" | bc | sed 's/0\+/\n/g' | wc -L
でよくね？

**デフォルトの名無しさん** · 2020/10/19(月) 22:00:17.27

>>735 Bash

$ echo "scale=100;e(l(2) / 3)" | bc -l
1.25992104989487316476721060727822835057025146470150
79800819751121552996765139594837293965624362550940

**デフォルトの名無しさん** · 2020/10/19(月) 22:30:56.87

>>735 dc
dcには平方根の演算子はあるが立方根はないので
c^(1/3) = sqrt(c/sqrt(c/sqrt(c/...)))
と変形して適当に反復

100k[2z333>av/]dsaxp

1.2599210498948731647672106072782283505702514647015079800819751121552\
996765139594837293965624362550941

**デフォルトの名無しさん** · 2020/10/19(月) 23:54:32.40

>>735
参考を誰も参考にしてくれないのにワロタ
参考を元に作ったコード
Haskell

https://ideone.com/ywJPjx

ちなみにrubyのpi.rbが同じテクニック使ってます
そっちは定数係数でない三項間関係の漸化式パージョン

**デフォルトの名無しさん** · 2020/10/20(火) 02:36:31.61

>>736
たしかに

**デフォルトの名無しさん** · 2020/10/20(火) 04:08:13.14

おっせえじゃん

**デフォルトの名無しさん** · 2020/10/20(火) 21:07:35.93

お題：ぼくがかんがえたさいきょうのちょうわへいきん
2変数調和平均 f(x, y) = 2 / (1 / x + 1 / y)
但し xまたはyがゼロの時、f(x, y) == 0
の計算を考える。
Cのdoubleに相当する精度で計算する際に、理論値(≒多倍長精度で計算したものを
double相当に丸めたもの)と一致しない割合(エラー率)をモンテカルロ法で求めよ。
x, yは、正規分布N(0, 1)の絶対値に従うものとする。

発展的なお題：
前掲の調和平均の式を変形し、エラー率の変化を調べよ。例えば、
f(x, y) = 2 * x * y / (x + y)
のエラー率は、前掲の式と比較してどうだろうか？
また、2 / (1 / x + 1 / y)よりエラー率が低い式は在るだろうか？

発展的なお題その2：
x, yがベンフォードの法則に従う場合はどうなるか。
正規分布N(0, 1)の絶対値の場合と比べて、エラー率の順位は入れ替わるだろうか？

**デフォルトの名無しさん** · 2020/10/22(木) 01:17:43.95

...but nobody answered

**253** · 2020/10/23(金) 09:17:13.74

正規乱数をモンテカルロ法で算出して
doubleと多倍長浮動小数点でそれぞれ計差した2変数調和平均に
差がある率を求めれば良いの？
正規乱数なら書かなくともﾗｲﾌﾞﾗﾘにあるんだよね。
あるいは誤り率の計算そのものにモンテカルロ法を使うようなやり方なのかな

**デフォルトの名無しさん** · 2020/10/23(金) 14:48:32.25

次のプリーズ

**デフォルトの名無しさん** · 2020/10/23(金) 20:02:36.14

お題
乱数を使ってシミュレーションしてください

ジョーカーを除くトランプ52枚から1枚抜き取って箱に入れました
その後、残った51枚から3枚引いて確認したら3枚ともダイヤでした
では最初に抜き取った箱の中が1枚がダイヤである確率はいくらか？

**デフォルトの名無しさん** · 2020/10/23(金) 20:27:24.17

>>746 Ruby

DECK = [*0..51]

def hoge
a, *b = DECK.sample(4) until b&.all?{|v| v % 4 == 0}
a % 4 == 0
end

def fuga(n)
1.0 * n.times.count{hoge} / n
end

p fuga(1e6) # => 0.20414
p 10 / 59 # => 0.20408163265306123

**デフォルトの名無しさん** · 2020/10/23(金) 20:46:02.34

途中で送信してしまった

99%信頼区間は 0.20414±0.0010382 => [0.2031018, 0.2051782]
理論値は 10/49 = 0.20408163265306123

**デフォルトの名無しさん** · 2020/10/23(金) 20:49:11.91

確認前後での確率変動を出せ、っていうお題じゃねえのコレ

**デフォルトの名無しさん** · 2020/10/24(土) 01:58:03.71

お題
ジョーカー除く５２枚のカードがある
スートを表すS,H,D,Cと、A,K,Q,J,X,9,8,7,6,5,4,3,2との２文字からなる文字列でカードの柄を表すこととする。
当然ながら５２枚のうちどの２枚も同じ柄でない

この一組のカードを使ってポーカーゲームをする。
すなわち、初手のカード５枚を配られたプレイヤーは、そのうち０枚以上５枚以下のカードを場に捨てて、同じ枚数のカードを残りの山札４７枚から引くことができる。
このときどのカードが選ばれるかは無作為であり、プレイヤーは選択できない。

さて、プレイヤーは、初手のカード５枚に対して、最終的に手に残る５枚のカードをポーカーのルールで評価して得られる配点の期待値が最も高くなるように捨て札を決めたい。
入力としてカード５枚が与えられたとき、そのような捨て札の組み合わせを求めるプログラムを記述せよ。

ここで、ポーカーの役と配点の対応はプログラムのパラメータとして組み込んで良い。

**デフォルトの名無しさん** · 2020/10/24(土) 21:09:47.85

>>750

// 任意の役の配点を0とする
printf(0);

**デフォルトの名無しさん** · 2020/10/24(土) 23:43:41.92

>>746
Java
https://paiza.io/projects/rAHsohR9Q8vhh4kLZnfKAw

**デフォルトの名無しさん** · 2020/10/25(日) 09:27:41.37

>>750 Perl5
https://ideone.com/stLkXf

役のレートは某ゲームのものを使用

ノーペアのとき
全部変えるより1枚残したほうが有利という結果が出てちょっと驚いてる

**デフォルトの名無しさん** · 2020/10/25(日) 10:49:49.59

3が1枚来てて、3が3枚残ってるなら、
・手元の3を捨てると、3が来ても何も起こらない
・手元に3を残すと、3が来るとワンペアになる

つまり、残りのカードを引けるかどうか、になる
>>746の亜種

**デフォルトの名無しさん** · 2020/10/25(日) 11:44:16.63

746はシミュレーションも何もただの引っ掛け問題でしょ
モンティホール問題の方が引っ掛け問としては出来が良い

**デフォルトの名無しさん** · 2020/10/25(日) 12:58:23.56

引っかけも何も条件付確率の初歩の初歩の問題だろ

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

52枚トランプの束で２枚目から４枚目までがダイアであるような場合に限定してその内先頭もダイアであるような場合は全体の何バーセントかを問う問題

**742** · 2020/10/26(月) 01:17:42.67

>>744
正規乱数自体は既存のライブラリで(Pythonならnumpy.random.randnとか)。
例えば、最初のお題(2/(1/x+1/y)で計算するやつ)は、34%位のエラー率になります。

**253** · 2020/10/27(火) 00:51:04.61

そうすると俺にはﾓﾝﾃｶﾙﾛ法をわざわざ使うところが
どこかよく分からないんだよね…