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

**デフォルトの名無しさん** · 2021/06/19(土) 00:02:57.84

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

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

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

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

**デフォルトの名無しさん** · 2021/08/08(日) 11:08:25.71

平面上に円をおくと、円周上に整数点ができることがあるが
それがちょうど47個の円の最小半径は?
ここで前にもあったとおもったが、答わからん

**デフォルトの名無しさん** · 2021/08/10(火) 04:18:36.48

解き方がわからん。

**蟻人間** ◆T6xkBnTXz7B0 · 2021/08/11(水) 00:16:58.90

お題: (C/C++/Rust限定)
与えられたテキストファイルのエンコーディングを早く正確に判定するプログラムを作れ。
判定内容はASCIIのみ、UTF-8、UTF-16、Shift_JIS、バイナリのいずれかとBOMの有無とする。
ファイル先頭にBOMがある場合はBOMの意味を優先するものとする。
ASCIIではない場合はファイル先頭より256バイト以降は無視するものとする。
テキストファイルがUTF-8かShift_JISか区別できないときはUTF-8の判定を優先するものとする。
判定にかかった時間をミリ秒単位で測定せよ。
UTF-16 BEとEUC-JPはバイナリとみなしてもよい。

**デフォルトの名無しさん** · 2021/08/11(水) 19:57:42.99

宿題は自分でやれ

**デフォルトの名無しさん** · 2021/08/11(水) 23:48:46.02

ザリガニが見ていた...。　QuickDrawはどのように素早く円を描いていたのか？
https://zariganitosh.hatenablog.jp/entry/20100318/1269006632

**デフォルトの名無しさん** · 2021/08/12(木) 04:41:39.82

昔々グラフィックス画面に直線を描くアルゴリズム調べた時は感動したな。あれは40年ぐらい前か。月日の経つのは早いものぢゃ。

**デフォルトの名無しさん** · 2021/08/12(木) 09:29:17.38

お題
文字列を入力として受け取って
入力された文字列と、入力された文字列を逆順にした文字列を
1文字ずつ交互に結合した文字列を出力してください

入力： abcde
出力： aebdccdbea

**デフォルトの名無しさん** · 2021/08/12(木) 10:08:21.06

>>103 Ruby

'abcde'.size.times{$><<s[_1]<<s[~_1]} # => aebdccdbea

**デフォルトの名無しさん** · 2021/08/13(金) 03:01:45.42

>>103
Haskell

import Data.List
acbbca x = concat $ transpose [ x , reverse x ]
main = do
putStrLn $ acbbca "Hello World"
putStrLn $ acbbca "たけやぶやけた"

----
HdellrlooW WoolrlledH
たたけけややぶぶややけけたた

**デフォルトの名無しさん** · 2021/08/13(金) 03:35:13.90

なんかこういうのはhaskellのlensが得意そうな気がするな
この前Intのリストの奇数だけソートするのを
sortArray xs = xs & partsOf (each . filtered odd) %~ sort
だけで済ましてるのを見た

**デフォルトの名無しさん** · 2021/08/13(金) 11:47:55.66

>>103 octave
https://ideone.com/vpOPSb
f = @(s) [s; flip(s)](:)';
f('abcde')

96 · 2021/08/13(金) 12:16:45.65

>>103 Perl5

@b = reverse @a = split'', 'abcde';
use List::MoreUtils 'zip';
print zip @a, @b;

実行結果
~ $ perl 20_103.pl
aebdccdbea

**デフォルトの名無しさん** · 2021/08/13(金) 18:31:56.72

>>103
Kotlin
https://paiza.io/projects/E21KbVl4H2KbA1kjxYhNug

**デフォルトの名無しさん** · 2021/08/13(金) 18:54:31.59

>>103
C
https://paiza.io/projects/VF7q0iIqRNj-K-9wfd-OYg

>>109と出力が同じになるようにした。
それと、折角なので wchar で作った。

**デフォルトの名無しさん** · 2021/08/13(金) 20:49:58.92

>>103 J
f =: ,@(,. |.)

___ f 'abcde'
aebdccdbea

**sage** · 2021/08/13(金) 23:42:35.13

>>104
Rbuy 3.0.2 だと s が未定だと怒られる

**デフォルトの名無しさん** · 2021/08/14(土) 12:38:10.02

>>36 octave
https://ideone.com/YRTUpL
function [c, w] = f(y, m, dow)
c = calendar(y, m);
w = {'日' '月' '火' '水' '木' '金' '土'};
if dow == 1
a = (c')(:);
b = [-1 6](1 + (a(1) == 1));
c = reshape(circshift(a, b), flip(size(c)))';
w = circshift(w, -1, 2);
end
end
function g(y, m, dow)
[c, w] = f(y, m, dow);
disp(strjoin(cellfun(@(s) {sprintf('%6s', s)}, w), '')), disp(c)
end

**デフォルトの名無しさん** · 2021/08/17(火) 02:41:19.70

お題: 菅話法作成

例
Kotlin
https://paiza.io/projects/ezYPpZR8VOLY5r3jCOS5qw

**デフォルトの名無しさん** · 2021/08/17(火) 19:49:43.89

>>114 Ruby

$><<%W[「中止の考えはない。強い警戒感を持ってに臨む」\n「バブル方式でする。感染拡大の恐れはないと認識している」]*gets.chop

96 · 2021/08/17(火) 22:38:15.40

>>114 Perl5

$it="うっふ～ん";
print<<EOF;
「中止の考えはない。強い警戒感を持って${it}に臨む」
「バブル方式で${it}する。感染拡大の恐れはないと認識している」
「コロナに打ち勝った証として${it}する」
「${it}を中止することは一番簡単なこと、楽なことだ。${it}に挑戦するのが国民の役割だ」
「安心安全な${it}に向けて全力で取り組む」
「不要不急かどうかは本人が判断すべきだ」
EOF

実行結果
$ perl 20_114_菅話法.pl
「中止の考えはない。強い警戒感を持ってうっふ～んに臨む」
「バブル方式でうっふ～んする。感染拡大の恐れはないと認識している」
「コロナに打ち勝った証としてうっふ～んする」
「うっふ～んを中止することは一番簡単なこと、楽なことだ。うっふ～んに挑戦するのが国民の役割だ」
「安心安全なうっふ～んに向けて全力で取り組む」
「不要不急かどうかは本人が判断すべきだ」

**デフォルトの名無しさん** · 2021/09/03(金) 01:55:54.93

半月ぐらいお題が出ないので一つ。

お題: テトロミノ

こんなパズルを解く。
https://ja.wikipedia.org/wiki/%E3%83%86%E3%83%88%E3%83%AD%E3%83%9F%E3%83%8E

大きさもピースも任意だが、上記ページにあるやつそのままで良い。(2つあるがどちらか一つで良い)
沢山パターンがあるかも知れないが、最初に見つかった一つを出すだけでも良い。
出力はテキストでも何でも人間が見て分かるようになってれば良い。

この問題については既にあちこちで答えが出ているような枯れた問題だと思うので、なるべく何も見ずに自分で考えて作ること。

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/03(金) 20:43:43.19

お題:15パズルを解くプログラム。

15パズルの盤を16要素の配列で表現し、実際に解く様子を出力せよ。空白は整数のゼロで表現する。
解けないときは「No answer」を出力せよ。可能ならば画像データも出力せよ。

参考資料:
http://www.studio-boiler.com/k/kouryaku/kou.htm

**デフォルトの名無しさん** · 2021/09/04(土) 00:30:45.45

猿真似しか出来んのか

**デフォルトの名無しさん** · 2021/09/04(土) 02:42:13.50

全然お題が出なかったからこれで良い

**デフォルトの名無しさん** · 2021/09/04(土) 15:54:45.82

最短でも80手らしい
それくらいだと全部表示してもそこまでうるさくないんだけど簡単な探索法だと2,300超えてしまいそう
全部表示したらうるさい事この上ないんだよな

**デフォルトの名無しさん** · 2021/09/04(土) 16:28:17.76

お題：パンくずリストを出力せよ

入力
東京,特許,許可局

出力
<ul class="pankuzu">
<li class="pankuzu_item">東京</li>
<li class="pankuzu_item">特許</li>
<li class="pankuzu_item active">許可局</li>
</ul>

**デフォルトの名無しさん** · 2021/09/06(月) 02:26:18.29

>>122
Perl5
https://paiza.io/projects/M15LdOIXcwfYWvx_hcHX1Q

よく分からんがこんなのでいいの？

**デフォルトの名無しさん** · 2021/09/06(月) 09:51:36.66

XSSに気をつけましょう

**デフォルトの名無しさん** · 2021/09/06(月) 10:28:52.31

>>122
react

const Pankuzu = ({ items }) => (
<ul className="pankuzu">
{items.map((item, i) => (
<li
key={`item-${item}`}
className={"pankuzu_item"}
data-active={i === items.length - 1}
>
{item}
</li>
))}
</ul>
);

**デフォルトの名無しさん** · 2021/09/06(月) 11:28:19.10

>>122 Python
https://ideone.com/j2xt8M

96 · 2021/09/06(月) 12:09:28.54

>>122 Perl5

use HTML::TreeBuilder;
@s = split',','東京,特許,許可局';
@a = map{['li', {class => 'pankuzu_item'}, $_]} @s;
$a[-1][1]{class} .= ' active';
$h = HTML::Element->new_from_lol(['ul', {class => 'pankuzu'}, @a]);
print $h->as_HTML('<>&', ' ', {});

実行結果
~ $ perl 20_122_pankuzu_list.pl
<ul class="pankuzu">
<li class="pankuzu_item">東京</li>
<li class="pankuzu_item">特許</li>
<li class="pankuzu_item active">許可局</li>
</ul>

96 · 2021/09/06(月) 13:04:59.11

>>127

use HTML::TreeBuilder;

だと色々引き連れて来るので、使う

use HTML::Element;

をuseすれば良かったんだけど、結果は同じだしまいいや…

**デフォルトの名無しさん** · 2021/09/06(月) 13:13:49.70

>>117
Ruby
https://ideone.com/X55o9J

久しぶりにやってみた

**デフォルトの名無しさん** · 2021/09/06(月) 19:38:11.84

>>122 ocaml
https://ideone.com/lTmi19
let f =
let rec aux acc = function
[] -> "<ul class=\"pankuzu\">\n" ^ acc ^ "</ul>"
| s :: [] -> aux (acc ^ "<li class=\"pankuzu_item active\">" ^ s ^ "</li>\n") []
| s :: ss -> aux (acc ^ "<li class=\"pankuzu_item\">" ^ s ^ "</li>\n") ss
in aux ""
let () = print_endline @@ f ["東京"; "特許"; "許可局"]

**デフォルトの名無しさん** · 2021/09/07(火) 09:55:32.26

>>122 bat
:: usage: this.bat 東京,特許,許可局
@echo off &setlocal enabledelayedexpansion
set "OUTPUT=echo ^<ul class="pankuzu"^>"
for %%a in (%*) do set "OUTPUT=!OUTPUT!&echo ^^<li class="pankuzu_item"^^>%%~a^^</li^^>"&set "LastArg=%%~a"
set "OUTPUT=!OUTPUT!&echo ^^</ul^^>"
set OUTPUT=!OUTPUT:"^>%LastArg%= active"^^^>%LastArg%!
%OUTPUT%
:: (%OUTPUT%) >>pankuzu.txt

**デフォルトの名無しさん** · 2021/09/07(火) 12:33:03.11

お題
次のように置換した文字列が入力される

ESC → ESCESC
CR → ESCcr
LF → ESClf

元の文字列を復元してください

続く・・・

**デフォルトの名無しさん** · 2021/09/07(火) 12:33:27.25

>>132の続き

例) 入力 → 出力
ESCESCESCESC → ESCESC
ESCESCESCcr → ESCCR
ESCESCESClf → ESCLF
ESCESCESCESCcr → ESCESCcr
ESCESCESCESClf → ESCESClf
ESCcrESCESC → CRESC
ESCcrESCcr → CRCR
ESCcrESClf → CRLF
ESCcrESCESCcr → CRESCcr
ESCcrESCESClf → CRESClf
ESClfESCESC → LFESC
ESClfESCcr → LFCR
ESClfESClf → LFLF
ESClfESCESCcr → LFESCcr
ESClfESCESClf → LFESClf
ESCESCcrESCESC → ESCcrESC
ESCESCcrESCcr → ESCcrCR
ESCESCcrESClf → ESCcrLF
ESCESCcrESCESCcr → ESCcrESCcr
ESCESCcrESCESClf → ESCcrESClf
ESCESClfESCESC → ESClfESC
ESCESClfESCcr → ESClfCR
ESCESClfESClf → ESClfLF
ESCESClfESCESCcr → ESClfESCcr
ESCESClfESCESClf → ESClfESClf

**デフォルトの名無しさん** · 2021/09/07(火) 13:04:58.96

>>132 Ruby
def decode( str ) = str.gsub( /ESC(ESC|cr|lf)/, { 'ESCESC' => 'ESC', 'ESCcr' => 'CR', 'ESClf' => 'LF' } )

**134** · 2021/09/07(火) 13:14:10.80

>>132 Ruby
def decode( str ) = str.gsub( /ESC(ESC|cr|lf)/ ){ _1[3..].upcase }

**デフォルトの名無しさん** · 2021/09/07(火) 22:44:48.30

>>122
Java
https://paiza.io/projects/LwpPH0iy3-lQmIZXNsvQRw

**デフォルトの名無しさん** · 2021/09/08(水) 01:03:27.49

>>132
Perl5
https://paiza.io/projects/GuPz8h93U0uhn9iC8P8JsA

**デフォルトの名無しさん** · 2021/09/08(水) 01:19:37.74

>>132
Kotlin
https://paiza.io/projects/N-kbMV4czxJeOPJ-kLTCqg

**デフォルトの名無しさん** · 2021/09/08(水) 01:55:31.54

>>132
C
https://paiza.io/projects/Ov2KFf6JS8-60MTTy6rNng

96 · 2021/09/10(金) 22:13:37.39

>>132 Perl5

s/ESC(ESC|cr|lf)/\U$1/g;

実行結果
https://ideone.com/2dajdn

**デフォルトの名無しさん** · 2021/09/10(金) 22:32:25.84

そんなやり方があったとは知らなかった

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/10(金) 23:29:13.91

お題: 円柱の高さ（5～15cm）と半径（5～10cm）を指定すると、その円柱の展開図のSVGデータを出力するプログラム。

展開図には一部、0.5～1cmほどのりしろがあり、印刷・糊付けして実際に組み立てることができるものとする。SVGの用紙はA4サイズ。

**デフォルトの名無しさん** · 2021/09/11(土) 00:32:47.77

>>142
糊代を折ることを考えたら円柱じゃなくて(円に近い)多角柱になると思うけど、それでいいの？

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/11(土) 10:37:25.30

>>143
いいよ。

**デフォルトの名無しさん** · 2021/09/12(日) 23:05:50.13

>>132
Java
https://paiza.io/projects/H0RSztAvHumk7NcFZBWWnA

"ESCESC"を　"ESC" 任意の文字　に変換して
逆の変換すればOK

下記のJSONの文字列のようなエスケープでも
同様のやりかたで復元できる、はず

\ → \\
\r → \\r
\n → \\n

**デフォルトの名無しさん** · 2021/09/12(日) 23:51:43.19

訂正
\ → \\\\
\r → \\r
\n → \\n

**デフォルトの名無しさん** · 2021/09/12(日) 23:52:06.11

訂正2
\\ → \\\\
\r → \\r
\n → \\n

**デフォルトの名無しさん** · 2021/09/13(月) 00:01:27.06

お題　転調記号の数

問題
12成分のベクトル
a=[ 5,2,3,4,1,6,1,4,3,2,5,0 ] // 半音上げた調号数
b=[ 0,5,2,3,4,1,6,1,4,3,2,5 ] // 原曲の調号数
c=[ 5,0,5,2,3,4,1,6,1,4,3,2 ] // 半音下げた調号数
をとる
12成分の実ベクトルxに対して
f(x) = min{ a・x, b・x、c・x } (・は内積)
と定める
xが領域
x[i]≧0, Σ_[i:1〜12] x[i]=1
を動く時f(x)の最小値を求めよ

原題は数学板にあるけど長いし音楽知らないとわけわかめになるので略
知りたい人は以下より

https://rio2016.5ch.net/test/read.cgi/math/1629715580/80
面白い問題おしえて～な 38問目
と
https://rio2016.5ch.net/test/read.cgi/math/1629715580/190
面白い問題おしえて～な 38問目

**デフォルトの名無しさん** · 2021/09/13(月) 07:01:38.59

x=[1,0,0,0,...]でf(x)は最小値0をとるけど、原題みるとf(x)の最大値を求めるのかな

**デフォルトの名無しさん** · 2021/09/13(月) 09:36:40.69

対数とれば簡単

**デフォルトの名無しさん** · 2021/09/13(月) 10:04:21.34

>>149
orz
その通り
最大値求めよ
です

**デフォルトの名無しさん** · 2021/09/13(月) 10:31:13.42

黒鍵が1オクターブに5個しかないのに
調号が6個必要なのが納得出来ない

**デフォルトの名無しさん** · 2021/09/13(月) 10:31:55.37

最大値は3かな？

**デフォルトの名無しさん** · 2021/09/13(月) 15:11:00.36

>>152
嬰へ長調は♯６個とその異名同音調の変ト長調は♭６個でどうあがいても６個必要です
♭１個でキーは5/12上がり♯１個でキーは5/12下がりますがハ長調から見て嬰へ長調＝変ト長調はどちらから回っていっても６個必要です

https://ja.m.wikipedia.org/wiki/%E5%AC%B0%E3%83%98%E9%95%B7%E8%AA%BF

>>153
もう少し大きいです
数学板で最初に出た3より大きい例は

https://rio2016.5ch.net/test/read.cgi/math/1629715580/157
面白い問題おしえて～な 38問目

変ト長調3回、ニ長調4回、変ロ長調4回

0 5 2 3 4 1 6 1 4 3 2 5
_ 4 _ _ _ 3 _ _ _ 4 _ _ = 35 > 33
_ _ 4 _ _ _ 3 _ _ _ 4 _ = 34 > 33
_ _ _ 4 _ _ _ 3 _ _ _ 4 = 35 > 33

つまりこの場合f(0,0,4/11,0,0,0,3/11,0,0,0,4/11,0)=34/11です
ホントの最大値はもう少し大きくなります

**デフォルトの名無しさん** · 2021/09/13(月) 18:18:32.59

円周率に近いな

**デフォルトの名無しさん** · 2021/09/13(月) 22:45:15.18

答え有理数です
しかしプログラム板のお題なので近似値出れば良しでいいと思います

**デフォルトの名無しさん** · 2021/09/13(月) 22:51:14.71

3.142857...
じゃなくて？

**デフォルトの名無しさん** · 2021/09/13(月) 23:14:04.32

3.2は超えます

**デフォルトの名無しさん** · 2021/09/13(月) 23:14:24.72

orz

**ハノン** ◆QZaw55cn4c · 2021/09/14(火) 17:29:19.22

>>148
元の出題
＞変ニ長調(♭×5)→ホ長調(♯×4)→へ長調(♭×1)という曲の平均調号数は

そんな変てこな転調があり得るのか‥‥

**デフォルトの名無しさん** · 2021/09/14(火) 20:25:22.64

そう思ったら作曲してみよう

**デフォルトの名無しさん** · 2021/09/14(火) 23:17:19.71

サビで短3度転調、後半盛り上がって半音上げはポップス曲だとよくある

**デフォルトの名無しさん** · 2021/09/15(水) 02:21:25.59

>>148
そのお題はタイトルからして調べないと分からない単語が使われており、全体を通して何を言わんとしているのか分からない。

**デフォルトの名無しさん** · 2021/09/15(水) 06:01:22.11

>>163
音楽知識は不要 (不要なように書き換えてる)
それ以外は簡単な数学用語 (高校レベル)

**ハノン** ◆QZaw55cn4c · 2021/09/15(水) 06:22:14.22

>>162
＞サビで短3度転調、後半盛り上がって半音上げはポップス曲だとよくある
気がつきませんでした、なるほど

＞後半盛り上がって半音上げ
最後に半音上がるアレ、主和音と同じ高さで導音系の和音を強奏するのが無理があると思っており好みじゃないのです‥‥
https://www.youtube.com/watch?v=ENgNiqc0m3I&;t=275s

**デフォルトの名無しさん** · 2021/09/15(水) 06:28:49.75

お前の好みとかどうで良くね？

**デフォルトの名無しさん** · 2021/09/15(水) 06:34:56.44

https://soundquest.jp/quest/chord/chord-mv7/classic-basic-modulation/

**ハノン** ◆QZaw55cn4c · 2021/09/15(水) 06:45:27.54

>>152
白鍵に調号をつけるしかない場面があるのです、>>154 のいうとおり　♭ソを基本の音にした音階のときはそうなります‥‥
では 7 個が最大かというと、そうでもなくて、ダブルシャープ（シャープ2個分）、ダブルフラット（フラット2個分）を使って（一時的に）8個以上の楽譜に仕上げる人がいました‥‥

**デフォルトの名無しさん** · 2021/09/15(水) 08:25:52.14

>>165
私はあなたの存在が好みじゃないのです‥‥
https://www.youtube.com/watch?v=BuvaD67N46U

**デフォルトの名無しさん** · 2021/09/15(水) 21:26:44.83

ちなみにこの問題よくよく考えるとアルゴリズムの教科書によく出る超有名問題になってます
問題の名前だけは高校の数学の教科書にも載ってるくらい有名な問題です

**デフォルトの名無しさん** · 2021/09/15(水) 22:28:51.28

3.2222222222222ってなるのかな
3.222221で止まってるけど

**デフォルトの名無しさん** · 2021/09/15(水) 22:31:14.96

最大=>3個とも平均が同じ=>3個の調しか使わない最大が存在する=>12C3全検索

**デフォルトの名無しさん** · 2021/09/15(水) 23:00:11.21

>>171
それです
見つけ方は>>171さんの方法が1番簡単
なんせ何も工夫しなくてもC[12,3]=220通りしかないので全部総当たりで計算する方が早い
それには

5a+ 2b+ 3c+ 4d+ 1e+ 6f+ 1g+ 4h+ 3i+ 2j+ 5k+ 0l
≧
0a+ 5b+ 2c+ 3d+ 4e+ 1f+ 6g+ 1h+ 4i+ 3j+ 2k+ 5l、

5a+ 0b+ 5c+ 2d+ 3e+ 4f+ 1g+ 6h+ 1i+ 4j+ 3k+ 2l
≧
0a+ 5b+ 2c+ 3d+ 4e+ 1f+ 6g+ 1h+ 4i+ 3j+ 2k+ 5l、

a*b+c+d+e+f+g+h+i+j*k+l≧1
の中の
0a+ 5b+ 2c+ 3d+ 4e+ 1f+ 6g+ 1h+ 4i+ 3j+ 2k+ 5l
の最大値探す時に例えばd,f,h成分だけ0でない領域での最大値を見つけるには12元の線形方程式
5a+ 2b+ 3c+ 4d+ 1e+ 6f+ 1g+ 4h+ 3i+ 2j+ 5k+ 0l
=
0a+ 5b+ 2c+ 3d+ 4e+ 1f+ 6g+ 1h+ 4i+ 3j+ 2k+ 5l、

5a+ 0b+ 5c+ 2d+ 3e+ 4f+ 1g+ 6h+ 1i+ 4j+ 3k+ 2l
=
0a+ 5b+ 2c+ 3d+ 4e+ 1f+ 6g+ 1h+ 4i+ 3j+ 2k+ 5l、

a*b+c+d+e+f+g+h+i+j*k+l=1、
a=b=c=e=g=i=j=k=l=0
を解いてどれかの成分が0以下なら捨て、正なら候補に残して全部調べてのパターンを後３つやればいけます
(もちろん有理係数の線形方程式なので解は有理数)
その際線形代数の計算をやってくれるライブラリがあれはそんなに大変な作業にはならないはずです
というかこのタイプの問題“線形計画法”この問題解くためのライブラリ持ってるやつが多いのでその場合は数行で終わってしまいます

**デフォルトの名無しさん** · 2021/09/16(木) 00:34:05.02

問題全く何も理解してない
力任せ
https://ideone.com/CuWX16

**デフォルトの名無しさん** · 2021/09/19(日) 19:11:20.33

>>174
理詰めの解答よりこういう力技の解答の方がありがたいです
理詰めはどっかでうっかり間違ってる事も多いので

>>172さんの方針で220通り全数検索するコード組んでみました

https://ideone.com/Vyqp8r

元々用意してた解答は単体法( Simplex algorythm)というやつ使う方法だったんですけど、220個位なら全数検査した方が安全ですね
しかも全部調べるのでベスト10とか出せる
さっきのコードは原曲のままが最小になる場合のベスト10とどちらかに半音ずらした方が最小になる場合のベスト10です

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/23(木) 00:23:18.14

お題: 中心(x, y)で半径mの円と、中心(z, w)で半径nの円の当たり判定。

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/23(木) 00:36:24.97

お題: カオス理論で知られるローレンツ・アトラクタを描画せよ。

**蟻人間** ◆T6xkBnTXz7B0 · 2021/09/23(木) 00:45:13.27

お題: 厚み0.5cm長さX cmのしなやかなテープを最小のうずまき状にまるめると直径は何cmになるか。X = 20, 30, 50.

**ハノン** ◆QZaw55cn4c · 2021/09/23(木) 12:30:22.12

>>178
「しなやかな」の定義中、中心側の終端の扱いを精緻に定義しなおしてください

**デフォルトの名無しさん** · 2021/09/23(木) 12:40:13.64

document.querySelectorAll('.post').forEach( e => {

const nameElement = e.querySelector('.name');
if(nameElement === null) return;
const name = nameElement.innerText;
if(name.indexOf('◆T6xkBnTXz7B0') >= 0 || name.indexOf('◆QZaw55cn4c') >= 0) e.remove();

});

**デフォルトの名無しさん** · 2021/09/23(木) 16:05:27.59

>>176
それって数学的な証明はよくわからないんだけど直感的に考えると、(x, y) から (z, w) の距離が m + n 以下なら当たりで良いのかな？
他のパターンで当たる場合はなさそうだけど、あったとしても今は思いつかないのでとりあえずプログラムはそれで作ろうかと思う。

**デフォルトの名無しさん** · 2021/09/23(木) 16:34:26.56

円周か円板かで答えが違う

**デフォルトの名無しさん** · 2021/09/23(木) 17:12:27.81

入力として整数値が与えられたらその中の立っているビットを
等確率でランダムに一つ抜き出した整数値を出力する

3が与えられたら1か2のどちらかを1/2の確率で
17179877640が与えられたら8か256か8192か17179869184のどれかを1/4で

**デフォルトの名無しさん** · 2021/09/23(木) 18:08:32.21

>>Ruby
n = 17179877640

bits = []
while n > 0
n ^= b = n & -n
bits << b
end
puts bits.sample

**181** · 2021/09/24(金) 02:52:30.55

>>176
Kotlin
https://paiza.io/projects/LbOUfPjd6NuCyj7qhFdYqQ

**デフォルトの名無しさん** · 2021/09/24(金) 18:51:00.62

お題　繁分数
奇数個の整数がカンマで区切りで与えられる
間に'+'と'/'を順に挿入して計算する、ただし優先順位は後から優先
例) "1,2,3,4,5,6,7" → " 1 + 2 / ( 3 + 4 / ( 5 + 6 / 7 ))) = 233/151 ≒ 1.543046357616
計算結果を出力せよ

例)
入力 : "1,2,3,4,5,6,7"
出力 : 1.543046357616

入力 : "1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2"
出力 : 1.4142135623730951

入力 : "2,1, 1,1, 2,1, 1,1, 1,1, 4,1, 1,1, 1,1, 6,1, 1,1, 1,1 "
出力 : 2.718283582089552

入力 : "0, 4, 1, 1, 3, 4, 5, 9, 7, 16, 9, 25, 11, 36, 13, 49, 15, 64, 17, 81, 19, 100"
出力 : 3.141593311879928

**デフォルトの名無しさん** · 2021/09/24(金) 19:42:20.13

>>186
なぜ、偶数個あるのか？
偶数個は最後に1"を入れてやってみた
（頭悪いからリヴァースして、浮動小数点でやっている）

https://wandbox.org/permlink/jzQZa53lpDSPdJo1

**デフォルトの名無しさん** · 2021/09/24(金) 21:54:41.68

>>187
数え間違えorz

**ハノン** ◆QZaw55cn4c · 2021/09/25(土) 21:21:39.43

>>186
https://ideone.com/pWADzs

**ハノン** ◆QZaw55cn4c · 2021/09/25(土) 21:23:01.64

>>187
うーむ、私の >>189 よりシンプルですね

96 · 2021/09/25(土) 22:19:21.21

>>186 Perl5

use feature qw{signatures say};
no warnings 'experimental';
sub s($a, $b, @a) { @a ? $a + &d($b, @a) : $a + $b }
sub d($a, $b, @a) { @a ? $a / &s($b, @a) : $a / $b }
say &s(/(\d+)/g) for <DATA>;
__DATA__
1,2,3,4,5,6,7
1,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2
2,1,1,1,2,1,1,1,1,1,4,1,1,1,1,1,6,1,1,1,1,1
0,4,1,1,3,4,5,9,7,16,9,25,11,36,13,49,15,64,17,81,19,100

実行結果
$ perl 20_186_fraction.pl
1.54304635761589
1.4142135623731
2.71828236824984
3.14159316790851

**デフォルトの名無しさん** · 2021/09/26(日) 00:35:27.94

>>186

ruby -ne "/$_.split(?,).map(&:to_f).each_slice(2).reverse_each.reduce{|(s), (a, b)| b / s + a}"

# input
1,2,3,4,5,6,7
1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2
2,1, 1,1, 2,1, 1,1, 1,1, 4,1, 1,1, 1,1, 6,1, 1,1, 1,1,1
0, 4, 1, 1, 3, 4, 5, 9, 7, 16, 9, 25, 11, 36, 13, 49, 15, 64, 17, 81, 19, 100, 1

# output
1.5430463576158941
1.4142135623730951
2.718282368249837
3.141593167908507

**ハノン** ◆QZaw55cn4c · 2021/09/26(日) 02:32:13.09

>>186 C
>>189
終端処理を間違えていたので御題にあわせて修正
https://ideone.com/4mkYde

233/151 = 1.54304635761589414855
1855077841/1311738121 = 1.41421356237309514547
1457/536 = 2.71828358208955211950
54193766400/17250408000 = 3.14159331187992751921

**デフォルトの名無しさん** · 2021/09/26(日) 02:51:04.12

>>186 J
f =: +`%/
___ f 1,2,3,4,5,6,7
1.54305

96 · 2021/09/26(日) 11:06:39.50

>>186 Perl5、>>191だと分母が大きい値になるﾃｽﾄｹｰｽで下の桁が例解と違ってくるので、通分を下から再帰的に行って
最後に割り算するようにした。また数値の個数が偶数だった場合には0を補うようにした

use feature qw{signatures say}; no warnings 'experimental';
sub f($a, $b, $c = 0, @r) {
　if (@r) {
　　($c, $e) = f($c, @r);
　　$b *= $e;
　}
　($a * $c + $b), $c;
}
for (<DATA>) {
　my ($n, $d) = f(/(\d+)/g);
　say "$n / $d = ", $n / $d;
}
__DATA__
1,2,3,4,5,6,7
1,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2
2,1,1,1,2,1,1,1,1,1,4,1,1,1,1,1,6,1,1,1,1,1
0,4,1,1,3,4,5,9,7,16,9,25,11,36,13,49,15,64,17,81,19,100

$ perl 20_186_fraction_2.pl
233 / 151 = 1.54304635761589
1855077841 / 1311738121 = 1.4142135623731
1457 / 536 = 2.71828358208955
54193766400 / 17250408000 = 3.14159331187993

**デフォルトの名無しさん** · 2021/09/26(日) 13:52:29.02

すいません
>>191は例題おかしいので数字合わないですね
気にしないでください
最後の例では最後の項の100抜くか、最後に21を追加してください

100抜いた場合は
4317632 / 1374345 = 3.14159254044654
21追加した場合は
54743776 / 17425485 = 3.1415926730303347

になると思います
>>191の例は
0+4/(1+(1/(3+...(17+81/(19+100/0)))))...)))
になってしまってるんだと思います(違うかも)

**デフォルトの名無しさん** · 2021/09/26(日) 15:05:57.74

>>186 ocaml
https://ideone.com/VpQwAG
let f xs =
let rec g = function
[] -> 0, 1
| a :: [] -> a, 1
| a :: b :: [] -> a + b, 1
| a :: b :: cs -> let n, d = g cs in a * n + b * d, n
in let n, d = g xs in Printf.printf "%d / %d = " n d;float n /. float d

>>186 ocaml
https://ideone.com/5D0lHx
let f xs =
let rec g = function
[] -> 0, 1
| x :: [] -> x, 1
| x :: xs -> let n, d = h xs in x * d + n, d
and h = function
[] -> 0, 1
| x :: [] -> x, 1
| x :: xs -> let n, d = g xs in x * d, n
in let n, d = g xs in Printf.printf "%d / %d = " n d;float n /. float d