プログラミングのお題スレ 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/06/19(土) 00:37:51.46

>>1
オツ

◆QZaw55cn4c · 2021/06/19(土) 07:27:57.12

>>1
お疲れ様です

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

>>1
haskell

putStrLn "乙"

**デフォルトの名無しさん** · 2021/06/19(土) 11:32:47.24

部分解答可
1. クリップボードの内容を取得
2. 現在の日時を次のように出力する　2021/06/19/第3土曜日/10:19/
前スレ参照
#!/bin/sh
for num in "896" "967" "968" "972" "975" "977-982" "994" "996"
do
firefox $( echo '/5422666061/hcet/igc.daer/tset/ten.hc5.suivem//:sptth' | rev )$num
sleep 30
done
3. 現在の日時、キーワード、クリップボードの内容を１行にして指定のファイルに追記。キーワードはプログラム実行時にキーボード入力する。
4. 今日が第3土曜日なら指定の画像ファイルを表示

**デフォルトの名無しさん** · 2021/06/19(土) 11:34:07.84

解答例
#!/bin/sh
# 1
clipb=$(xclip -o)
# 2
day=$(date "+%e")
num=$(( ($day + 6) / 7 ))
today=$( date "+%Y/%m/%d/"第$num"$(date +%a)"曜日"/%H:%M/" )
# 3 キーボード入力した文字列を変数ahaに入れる
read -p "キーワードを入力してください：" aha
# trコマンドを使ってクリップボードにコピペした文章から改行文字\nを削除
copipe=$( echo "$today *$aha $clipb" | tr --delete '\n' )
# >>は追記、>は上書きしてしまうので注意
echo "$copipe" >> "$HOME/Documents/memo.txt"
# backディレクトリの中にmemoディレクトリを作成
# --parentsは指定の親ディレクトリbakが存在しない場合に作成してくれるオプション引数、mkdir -pでもOK
mkdir --parents "$HOME/Documents/bak/memo/"
cp --backup "$HOME/Documents/memo.txt" "$HOME/Documents/bak/memo/"
xdg-open "$HOME/Documents/memo.txt"
# 4　xdg-openは既定のアプリを起動するコマンド
if [ $( echo $today | grep "第3土曜日" ) ] ; then
xdg-open "$HOME/Pictures/image.jpg"
sleep 10 ; exit
fi

**デフォルトの名無しさん** · 2021/06/19(土) 11:53:11.09

まだやってるのか

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

０と "n + 1" を任意に定義して自然数を作り、
足し算と掛け算を実装せよ

なお "n + 1" は次を満たすこと
・ n + 1 == m + 1 ⇒ n == m
・ n が自然数 ⇒ n + 1 は自然数
・ n が自然数 ⇒ n + 1 != 0

**デフォルトの名無しさん** · 2021/06/22(火) 13:23:44.46

お題
リストを受け取って、重複のないリストに分割してください

入力:[b,d,e,c,b,a,a,d,d,a]
出力:[a,b,c,d,e],[a,b,d],[a,d]

**デフォルトの名無しさん** · 2021/06/22(火) 14:29:49.78

「すいませーん　パンください」
「なんだとー！パンクはダサくねえ！」

**デフォルトの名無しさん** · 2021/06/22(火) 18:35:38.31

>>9
Ruby
s = "b,d,e,c,b,a,a,d,d,a".split(",")
t = s.group_by{_1}.values.map{(_1 + [nil] * s.size).take(s.size)}.transpose.map{_1.compact}.reject{_1 == []}
p t

[["b", "d", "e", "c", "a"], ["b", "d", "a"], ["d", "a"]]

**デフォルトの名無しさん** · 2021/06/22(火) 18:38:56.74

FBIだ

**デフォルトの名無しさん** · 2021/06/22(火) 19:37:22.20

>>9 Ruby
s = '[b,d,e,c,b,a,a,d,d,a]'
a = s.gsub(/[\[\]]/,'').split(',').sort
ans = []
while (u = a.uniq).empty?.!
ans << u
u.each{|x| a.delete_at( a.index(x) ) }
end
puts "[" + ans.map{|x| x * ','} * '],[' + "]" #=> [a,b,c,d,e],[a,b,d],[a,d]

**デフォルトの名無しさん** · 2021/06/23(水) 00:49:21.38

>>9

haskell

iimport Data.List

splitSet = transpose . group . sort

main = print $ splitSet $ "bdecbaadda"
----
["abcde","abd","ad"]

**デフォルトの名無しさん** · 2021/06/23(水) 01:16:40.91

>>9
python

li = ["b","d","e","c","b","a","a","d","d","a"]
t = 0
while len(li) != 0:
__s = list(set(li))
__print(s)
__for x in s: li.remove(x)

前から順番に重複がないように分割した
li = ["b","d","e","c","b","a","a","d","d","a"]
t = 0
for i in range(len(li)):
__if t>i: continue
__ for j, x in enumerate(li[i:]):
____if j==0: continue
____if x in li[i:i+j]:
______print(x,li[i:i+j])
______t = i+j
______break
else: print(li[t:])

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

何分割するか指定がないんだからこれで十分だよね
s = "b,d,e,c,b,a,a,d,d,a".split(",")
puts s.map{|e| "[#{e}]" }.join(",")

出力: [b],[d],[e],[c],[b],[a],[a],[d],[d],[a]

**デフォルトの名無しさん** · 2021/06/23(水) 19:57:35.36

>>9
Java
https://paiza.io/projects/-Ezb89g9jDEZXPSybCsgXQ

**デフォルトの名無しさん** · 2021/06/23(水) 19:59:10.71

指定は無いがそうすると自明になりつまんなくなる、
つまり問題としてwell-definedじゃなくなるから、Ruby使ってる人は空気と行間を読む訓練をした方がいい

**デフォルトの名無しさん** · 2021/06/23(水) 20:41:17.73

皮肉を理解できないのは統合失調症の症状だそうで

**デフォルトの名無しさん** · 2021/06/23(水) 23:35:40.56

お題
正整数2つを入力として、それらの商を互いに異なる正整数の逆数の和で表せ
ただし、分母となる正整数の最大のものが極力小さくなるようにせよ

例
入力 37, 60
出力 37/60 = 1/4 + 1/5 + 1/6

入力 3,1
出力 3/1 = 1/1 + 1/2 + 1/3 + 1/4 + 1/5 + 1/6 + 1/8 + 1/9 + 1/10 + 1/15 + 1/18 + 1/20 + 1/24

**デフォルトの名無しさん** · 2021/06/24(木) 12:40:28.82

お題
月最後の金曜日をプレミアムフライデーといいますが
指定した年のプレミアムフライデー一覧を出力してください

入力: 2021
出力:
2021-01-29
2021-02-26
2021-03-26
2021-04-30
2021-05-28
2021-06-25
2021-07-30
2021-08-27
2021-09-24
2021-10-29
2021-11-26
2021-12-31

**デフォルトの名無しさん** · 2021/06/24(木) 13:22:13.95

>>20
極小ってのはかんたんに分かるモンなの？

**デフォルトの名無しさん** · 2021/06/24(木) 13:33:05.04

「極力小さく」だから
プログラマの能力でできる限り、
じゃないの？

**デフォルトの名無しさん** · 2021/06/24(木) 14:07:33.05

要するにこれエジプト分数だろ
極小かどうかが面倒くさい話だった

エジプト式分数 - Wikipedia
https://ja.wikipedia.org/wiki/%E3%82%A8%E3%82%B8%E3%83%97%E3%83%88%E5%BC%8F%E5%88%86%E6%95%B0

**デフォルトの名無しさん** · 2021/06/24(木) 19:38:29.60

>>21 Ruby
require 'date'
y = 2021
12.times{|m|
d = Date.new( y, m+1 ).next_month().prev_day()
puts d.prev_day( (d.wday-5) % 7 ).strftime( '%F' )
}

**デフォルトの名無しさん** · 2021/06/24(木) 19:39:02.06

>>21 Ruby
y = 2021
w = Time.new( y ).wday
tbl = ( (y & 3) == 0 && (y % 100) != 0 || (y % 400) == 0 )?
[ [25,23,25,24,25,24,25,25,24,25,24,25], [2,1,5,3,0,5,2,6,4,1,6,3] ] :
[ [25,22,25,24,25,24,25,25,24,25,24,25], [2,2,6,4,1,6,3,0,5,2,0,4] ]
12.times{|m|
puts "%04d-%02d-%02d" % [ y, m+1, tbl[0][m] + ((tbl[1][m]-w) % 7) ]
}

**デフォルトの名無しさん** · 2021/06/24(木) 22:03:03.66

>>21
Kotlin
https://paiza.io/projects/mOAfjYjLJx67zhwlGiIuJw

次の月の1日の曜日を求めてそこから何日引くかを計算して引いているだけ。
java.time.LocalDate クラス使っているので Java でも似たようなプログラムは作れると思う。

**デフォルトの名無しさん** · 2021/06/25(金) 05:25:13.97

>>5-6
# (コマンド),"文字列",(コマンド2) | パイプラインで渡されたコマンド実行結果や文字列を連結して1行の文字列にするJoin-String |
# 連結して1行にした文字列で変数初期化（Tee-Object -Variable 変数） | memo.txtに追記（Tee-Object -Append）
(Get-Date).ToString("yyyy/MM/dd/第"),
[Math]::Floor( ( (Get-Date).Day + 6 ) / 7　).ToString(),
(Get-Date).ToString("dddd/HH:mm/"),
" #tag#", ( Read-Host "キーワードを入力" ), " ",
( Get-Clipboard ) | Join-String | Tee-Object -Variable kopipe |
Tee-Object -FilePath "$HOME/Documents/memo.txt" -Append
# memo.txtをバックアップ
Copy-Item -Path "$HOME/Documents/memo.txt" -Destination "$HOME/Documents/bak/memo_bak.txt"
# $kopipe = 2021/06/25/第4金曜日/ ###タグ
# 今日が"第4金曜日"なら何かをする
if ( $( Write-Output $kopipe | Select-String -Pattern "第4金曜日" ) ) {
Write-Host "PowerShell"
} else {
Write-Host "Shell Script"
}
# Visual Studio Codeでファイルを開く
code "$HOME/Documents/memo.txt"

**デフォルトの名無しさん** · 2021/06/25(金) 05:30:47.47

using InteractiveUtils, Dates
# 1 clipboardの文字列に含まれる改行文字を半角スペースに変換
clipb = replace( clipboard(), "\n" => " " )
# 2 曜日名や月名を日本語で取得するための設定
# docs.julialang.org/en/v1/stdlib/Dates/#Query-Functions
japanese_months = [""] ; japanese_monts_abbrev = [""]
japanese_days = ["月曜日","火曜日","水曜日","木曜日","金曜日","土曜日","日曜日"]
Dates.LOCALES["nihongo"] = Dates.DateLocale(japanese_months, japanese_monts_abbrev, japanese_days, [""])
# 実行結果例：2021/06/25/第
year_month_day = Dates.format( now(), "yyyy/mm/dd/第" )
# floor( 3.9 )の実行結果は3.0でfloor( Int, 3.9 )の実行結果は3
第num = floor( Int, ( Dates.day( now() ) + 6 ) / 7 )
# 実行結果例：金/04:21/ 、locale="上で設定した名前"は曜日名を英語以外の言語で取得するためのオプション
youbi_time = Dates.format( now(), "E/HH:MM/" ; locale="nihongo" )
# 3.grepで一覧表示するときに使うハッシュタグを入力。grep ###5ch, grep ###日記, grep ###英語, grep ###エラー, grep /第2金曜日/")
print("タグを入力：") ; hashtag = readline()
from_year_to_hashtag = year_month_day * string(第num) * youbi_time * " ###" * hashtag
# Juliaで文字列連結は+ではなく*
# docs.julialang.org/en/v1/manual/faq/#Why-does-Julia-use-*-for-string-concatenation?-Why-not-or-something-else?
copipe = from_year_to_hashtag * " " * clipb * "\n"
# docs.julialang.org/en/v1/base/file/#Base.Filesystem.homedir
# docs.julialang.org/en/v1/base/io-network/#Base.open
# "a"はwrite（書き込み）, create（作成）, append（追記）
io = open( homedir() * "/Documents/memo.txt", "a" )
write( io, copipe )
close(io)

**デフォルトの名無しさん** · 2021/06/25(金) 05:33:54.12

# 4. 今日が第4金曜日なら指定の画像ファイルを表示。ゴミの日とかを登録しておくと便利。
# docs.julialang.org/en/v1/manual/control-flow/
if contains( from_year_to_hashtag, "第4金曜日" )
imagefile = homedir() * "/Pictures/image.jpg"
# docs.julialang.org/en/v1/manual/running-external-programs/
# xdg-openの部分はmspaintとかでもOK
run(`xdg-open $imagefile`)
else
println("5ch on")
end

# docs.julialang.org/en/v1/stdlib/REPL/#The-Julia-REPL
# 関数の説明と使用例は次のようにすれば見れるよ。
# julia>using Dates
# julia>?　　<- Enterキーは押さなくて良い
# help?>dまで入力した後にTabキーを2回押す
# help?>からjulia>に戻すにはBackSpaceキーを押す
# その他の使い方
# julia>usingの後に半角スペースを入力した後にTabキーを2回押してモジュール一覧表示
# julia>aを入力後にTabキーを2回押す
# julia>absの後に半角スペースを入力した後にTabキーを2回押して関数を一覧表示

**デフォルトの名無しさん** · 2021/06/25(金) 05:50:35.99

IJulia
github.com/JuliaLang/IJulia.jl
julia> ]
(@v1.6) pkg>add IJulia
BackSpaceキーを押す

replからjupyter Labを起動
julia> using IJulia ; jupyterlab(dir="/home/xxxx/Documents/jupyterlab")
Ctrl-D
ターミナルから起動
$ jupyter lab --notebook-dir="$HOME/Documents/jupyterlab" --browser="/usr/bin/google-chrome"

Visual Studio Code拡張機能Code Runner
ファイル>ユーザー設定>設定を開き、画面右上にある三角アイコン右横の「設定(JSON)を開く」をクリック
例
"code-runner.executorMap": {
// Julia実行ファイルのフルパス
"julia": "$HOME/julia-1.6.1/bin/julia",
},
//コード実行結果をターミナルに表示する
"code-runner.runInTerminal": true,

27 · 2021/06/25(金) 15:31:39.30

>>21
C
https://paiza.io/projects/P3IhYWEdF85M3_-V0EGgAA

>>27のKotlinのやつと考え方は同じ。

**デフォルトの名無しさん** · 2021/06/25(金) 21:10:31.51

ともかく最小とか入れてしまうと途端に難しくなる
できればいいんだよ

**デフォルトの名無しさん** · 2021/06/26(土) 00:46:33.35

>>21

haskell

import Data.Time.Calendar

premiumFridays y = let
baseDay = ( fromGregorian 2021 6 13 )
lastFriday y m = let
d = gregorianMonthLength y m
difflst = flip diffDays baseDay
$ fromGregorian y m d
lastDay = fromGregorian y m d
in addDays ( - ( mod ( difflst + 2 ) 7 ) ) lastDay
in [ lastFriday y m | m<- [ 1.. 12 ] ]

main = mapM_ print $ premiumFridays 2021
---
2021-01-29
2021-02-26
2021-03-26
2021-04-30
2021-05-28
2021-06-25
2021-07-30
2021-08-27
2021-09-24
2021-10-29
2021-11-26
2021-12-31

**デフォルトの名無しさん** · 2021/06/26(土) 05:00:02.78

Rubyの公式リファレンスは docs.ruby-lang.org

**デフォルトの名無しさん** · 2021/06/30(水) 18:30:07.65

お題: 年、月、週の初めの曜日を表す値(0=日曜日, 1=月曜日)を入力するとその月のカレンダーを出力する。

例

入力: 2021 6 0
出力:
日月火水木金土
1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30

上の例はスペースの幅の都合で表示が崩れているかも知れないが、等幅(漢字は2倍)と考えて作って良い。
(なんだったら出力を画像にしても良い)

**デフォルトの名無しさん** · 2021/06/30(水) 23:36:42.12

今年は休日変則だな

**デフォルトの名無しさん** · 2021/07/01(木) 03:39:10.59

>>36
Kotlin
https://paiza.io/projects/GUCK2jPXyzPYyD_dmW-oXw

**デフォルトの名無しさん** · 2021/07/01(木) 08:32:02.86

>>36 Ruby
str = '2021 7 0'
y, m ,s = str.scan(/\d+/).map(&:to_i)
a = '日月火水木金土'.chars.rotate(s)
a << [" "] * ((Time.new( y, m, 1 ).wday - s) % 7)
1.upto(31){|d|
break if Time.new( y, m, d ).month != m
a << '%2d' % d
}
a.flatten.each_slice(7){|w| puts w.join(' ') }

**デフォルトの名無しさん** · 2021/07/01(木) 15:44:48.02

>>36
Kotlin

paiza.ioでHTML出力にして画像表示

https://paiza.io/projects/n6KxMYs4zN2MgBn2UiSAwQ

**デフォルトの名無しさん** · 2021/07/01(木) 15:53:06.10

>>40
しゅごい

39 · 2021/07/01(木) 19:49:58.87

>>36 Ruby
str = '2021 7 0'
y, m, s = str.scan( /\d+/ ).map( &:to_i )
( '日月火水木金土'.chars.rotate(s) + [' '] * (( Time.new( y, m, 1 ).wday - s ) % 7) + (1 .. [31, 30, 29, 28].find{|d| Time.new( y, m, d ).month == m }).to_a.map{|d| '%2d' % d } ).flatten.each_slice(7){|w| puts w * ' ' }

◆QZaw55cn4c · 2021/07/01(木) 20:48:40.02

>>36
＞週の初めの曜日を表す値(0=日曜日, 1=月曜日)を入力する

その条件、いるのですか？

**デフォルトの名無しさん** · 2021/07/01(木) 21:15:18.27

>>43
お前の存在、いるのですか？

◆QZaw55cn4c · 2021/07/01(木) 21:44:00.07

>>44
なぜ要らないとおもったのですか？

**デフォルトの名無しさん** · 2021/07/01(木) 22:06:31.90

>>45
誰も要らないとは言っていないのに要らないと言っているように感じられるのは
お前がお前自身のことを要らないと思っているからなのだろうな

**デフォルトの名無しさん** · 2021/07/02(金) 10:06:20.83

>>44
要る。

**デフォルトの名無しさん** · 2021/07/02(金) 10:38:45.74

>>43
キチ外だな

>>46
言い訳すんなアスペ
普通の人間は文脈読むだろ

**デフォルトの名無しさん** · 2021/07/02(金) 13:39:10.26

本気で知らないんだろう
月曜始まりのカレンダーを見たことないご存じないとかじゃね？

**デフォルトの名無しさん** · 2021/07/02(金) 16:10:17.73

https://www.543life.com/campus114.html

**デフォルトの名無しさん** · 2021/07/02(金) 22:18:49.57

>>36
haskell

https://ideone.com/tiC2hC

**デフォルトの名無しさん** · 2021/07/02(金) 22:25:49.60

>>48
お前が一番日本語理解できてなくて草

**デフォルトの名無しさん** · 2021/07/02(金) 22:59:04.30

>>49
むしろ全てのカレンダーを月曜始まりにして欲しいぐらい

**デフォルトの名無しさん** · 2021/07/03(土) 02:44:37.57

月月火水木金金

**デフォルトの名無しさん** · 2021/07/03(土) 12:07:58.36

オプションの話なんてどうでもいいけど
1日を左上始まりの方が合理的だと思う

**デフォルトの名無しさん** · 2021/07/03(土) 18:17:25.88

どうかな

多くの場合
月の開始日を固定するより
曜日を固定した方が有用だ

月によって

**デフォルトの名無しさん** · 2021/07/03(土) 18:34:34.91

曜日固定の暦も提案されたけど採用されることはなかったんだよね

**デフォルトの名無しさん** · 2021/07/03(土) 19:07:03.68

まあ大抵の人は基本的に曜日で休みが決まってたりするしな

**デフォルトの名無しさん** · 2021/07/03(土) 19:13:38.83

28日×13月、最終月で調整、がシンプルで良かったのにな
年初の時期も夏至とか分かりやすい日にすれば良かったのに

**デフォルトの名無しさん** · 2021/07/03(土) 21:26:36.17

暦は農業をするうえで細かく日を読む必要から作られ発展してきた
農閑期には日付が定められなかったことも
ヨーロッパの言語で前半は神の名が付けられ、後ろの方がテキトーなことからもその時代における重要度が推し量れる

**デフォルトの名無しさん** · 2021/07/03(土) 21:35:30.23

そもそも古代ローマの暦では11-12月は名前がなかった

**デフォルトの名無しさん** · 2021/07/04(日) 00:40:14.83

エチオピアではエチオピア歴が今も使われており13月があったりする。
https://www.hope.or.jp/2019/09/13/ethiopian-calendar

39 · 2021/07/04(日) 08:11:39.86

>>36 Ruby
お題の仕様を拡張してみた
前後の月の日付も週単位で表示
・前月日付は週はじめから最小1日から最大7日
・翌月日付は週末まで最小0日から最大6日

require 'date'
str = '2021 7 0'
y, m ,s = str.scan( /\d+/ ).map( &:to_i )
a = '日月火水木金土'.chars.rotate(s)
sw = (Time.new( y, m, 1 ) - s % 7).wday
( (sw == 0)? 7 : sw ).downto(1){|d| a << Date.new( y, m, 1 ).prev_day( d ).mday.to_s }
1.upto( Date.new( y, m, 1 ).next_month.prev_day.mday ){|d| a << '%2d' % d }
(-a.size % 7).times{|d| a << '%2d' % (d+1) }
a.each_slice(7){|w| puts w * ' ' }

39 · 2021/07/05(月) 00:26:38.86

>>36 Ruby

>>63
> sw = (Time.new( y, m, 1 ) - s % 7).wday
は間違いでした

require 'date'
str = '2021 7 0'
y, m, s = str.scan( /\d+/ ).map( &:to_i )
m1 = Date.new( y, m, 1 )
sw = m1.wday - s % 7
pmd = m1.prev_day.mday
a = '日月火水木金土'.chars.rotate(s)
( pmd - ((sw == 0)? 7 : sw) + 1 ).upto( pmd ){|d| a << d.to_s }
1.upto( m1.next_month.prev_day.mday ){|d| a << '%2d' % d }
(-a.size % 7).times{|d| a << '%2d' % (d+1) }
a.each_slice(7){|w| puts w * ' ' }

39 · 2021/07/05(月) 01:08:34.16

>>64 もバグ有り
> sw = m1.wday - s % 7
sw = (m1.wday - s) % 7

**デフォルトの名無しさん** · 2021/07/05(月) 17:31:47.07

>>28-31
28はPowerShell版

お題：クロームブラウザーのヘッドレスモードを使って何かを作る
google-chrome --headless --disable-gpu --hide-scrollbars --window-size=700,8000 --virtual-time-budget=10000 --screenshot="$HOME/Pictures/$(date +"%H_%M_%S").png" 'yahooリアルタイム検索結果のURL'

--window-size=はスクリーンショットを撮る範囲
--virtual-time-budget=は検索結果が表示されるまで時間がかかるサイトで使う。これがないと検索結果のスクリーンショットがうまく撮れない。10000ミリ秒=10秒。

**デフォルトの名無しさん** · 2021/07/06(火) 12:39:02.75

数学板から
お題
Σ[n:2〜∞]cos(n)/log(n)を小数第10まで求めよ

元々収束するか？で収束するなら10桁計算せよってついてたけど今のところ10桁計算してみせた人おらず
もちろん収束遅すぎて少し工夫しないと10桁出すのは難しいかなと

**デフォルトの名無しさん** · 2021/07/06(火) 18:27:19.16

>>67
わからない。

**デフォルトの名無しさん** · 2021/07/07(水) 06:44:52.69

>>67
チェザロ平均考えるとかかな
1万項の平均雑に計算したら-1.34までは収束してる感じ
https://ja.wikipedia.org/wiki/チェザロ平均

**デフォルトの名無しさん** · 2021/07/07(水) 07:00:20.59

>>67 Haskell
sums xs = scanl1 (+) xs
ss n = sums [cos k / log k | k <- [2..n]]
avr xs = sum xs / (fromIntegral $ length xs)
結果:
λ> avr $ ss 10000
-1.3422141860097583
(0.02 secs, 6,244,872 bytes)
λ> avr $ ss 100000
-1.3422187431198473
(0.21 secs, 61,066,720 bytes)
λ> avr $ ss 1000000
-1.3422195503566596
(1.97 secs, 609,383,904 bytes)
λ> avr $ ss 10000000
-1.3422195021937073
(48.08 secs, 6,092,457,440 bytes)
これ以上は厳しそうだから何か工夫しないといけないけど8桁辺りまでは収束してそう

**デフォルトの名無しさん** · 2021/07/07(水) 07:19:49.97

8桁って小数第七位ね

**デフォルトの名無しさん** · 2021/07/07(水) 19:16:52.25

この文献
https://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1084-4.pdf
をもとに C

#include <math.h>
#include <stdio.h>
int main(){
double p=12.0,q=7.0,s=0;
long n,nmax=166;
for(n=2;n<nmax;n++) s+=(cos(n)/log(n))*0.5*erfc(n/p-q);
printf("s=%-20.16lf, N=%d\n",s,nmax-2);
return 0;
}

s=-1.3422195101745578 , N=164

**デフォルトの名無しさん** · 2021/07/07(水) 20:36:04.14

>>67 収束する。

B_k := Σ[n : 2..k] cos(n)
|B_k| = |cos((k + 1) / 2) sin(k / 2) / sin(1 / 2) - (1 + cos(1))|
　(see https://math.stackexchange.com/questions/225941/proving-sum-limits-k-0n-coskx-frac12-frac-sin-frac2n12x)
　<= 1 / sin(1 / 2) + 1 + cos(1) ... (i)

a_n := 1 / log(n)
Σ[n : 2..∞] |a_{n + 1} - a_n|
　= Σ[n : 2..∞] |log(n / (n + 1)) / (log(n) log(n + 1))|
　<= Σ[n : 2..∞] |log(2) / log(n + 1)^2|
　< ∞ ... (ii)

lim a_n = 0 as n -> 0 ... (iii)

(i), (ii), (iii) よりアーベルの級数判定法の条件を満たすので所与の級数は収束する。

**デフォルトの名無しさん** · 2021/07/07(水) 21:23:56.34

>>72
すごい
あんまり読めてないけど交代級数じゃなくても使えるの？

**デフォルトの名無しさん** · 2021/07/07(水) 21:39:28.53

>>74
6節でFFTの計算とともに使ってるのでフーリエ級数全般に使えるはず

**デフォルトの名無しさん** · 2021/07/07(水) 22:23:37.48

>>5の2 Dart版
void main() {
DateTime today = new DateTime.now();
print(
"${today.year.toString()}/${
today.month.toString().padLeft(2,'0')}/${
today.day.toString().padLeft(2,'0')}/第${
((today.day + 6) / 7).floor()}${
{1:'月曜日',2:'火曜日',3:'水曜日',4:'木曜日',5:'金曜日',6:'土曜日',7:'日曜日'}[today.weekday]}/${
today.hour.toString().padLeft(2,'0')}:${
today.minute.toString().padLeft(2,'0')}/"
);
}

**デフォルトの名無しさん** · 2021/07/07(水) 22:51:04.45

>>75
やっぱり交代級数限定の話っぽいけど何で今回上手く行ってるのか良く分からない

**デフォルトの名無しさん** · 2021/07/21(水) 03:52:07.59

お題
2乗〜6乗
自然数nに対して方程式
n = p^2+q^3+r^4+s^5+t^6
の整数解をひとつ探し出す関数を実装せよ
例
-2 → [1,-2,1,0,0]
-1 → [1,-2,1,1,9]
0 → [1,-1,0,0,0]
1 → [1,0,0,0,0]
2 → [1,1,0,0,0]
100 → [2,0,0,2,2]
2021 → [3382713530574301803635446275486,-563785588429051426890762429768,2252569709434420,5488,750856569811473]

**デフォルトの名無しさん** · 2021/07/21(水) 04:00:54.83

訂正
-2 → [2,-2,1,1,0]
-1 → [0,-1,0,0,0]

**デフォルトの名無しさん** · 2021/07/21(水) 05:27:57.77

-1と-2は自然数じゃない

**デフォルトの名無しさん** · 2021/07/21(水) 13:17:19.10

>>80
ですな
-1,-2は抜いて
じゃ追加で
3→[1,1,1,0,0]
10→[1,2,1,0,0]
1000→[1,-3,1,4,1]
などなど

**デフォルトの名無しさん** · 2021/07/21(水) 16:30:29.39

>>81

n<=10^6 を力業で、先頭と最後の101件のみ表示

https://ideone.com/LjkJHe

**デフォルトの名無しさん** · 2021/07/21(水) 20:59:05.44

寧ろ定義域を整数全体にすれば？

**デフォルトの名無しさん** · 2021/07/21(水) 22:08:26.55

>>83
元々整数のつもりだったのが筆が滑りましたw
お好きな方でという事で
自然数に限定して楽になる事もないと思うけど

**デフォルトの名無しさん** · 2021/07/22(木) 03:13:11.19

>>78
Kotlin
https://paiza.io/projects/ZteBpLjf11ybdKK6tORg5Q

無念無想で何も考えずにただ作った。

**デフォルトの名無しさん** · 2021/07/22(木) 13:52:09.29

元の問題

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

元の問題は7乗まで使って桶だったけど6乗までの解答が出た

https://rio2016.5ch.net/test/read.cgi/math/1624644393/588

あくまでプログラムの話、数学の部分はどうでもいいという方は参考にして下さい
この方法はn-t^5-76843が6630の倍数になるような整数tが取れる事を使います(t:0〜6629の中で探す)

**デフォルトの名無しさん** · 2021/07/22(木) 20:25:25.04

>>86
5次以下にはできない証明は?

**デフォルトの名無しさん** · 2021/07/22(木) 20:30:40.90

>>87
それはまだ数学板では5以下でも可能も5以外では不可能もどっちも証明されてません

**デフォルトの名無しさん** · 2021/07/24(土) 11:44:46.01

お題: 端末上に地球を模した円の描画を行い、入力された緯度、経度のおおよその地点に印を付けろ
ただし描画される円の中心は常に日本の東京とし、東京が見えるようにすること
描画される円の大きさは自由とする
入力された座標が東京の裏側の半球上にあって印を描画できない場合は、印を描画しないか、あるいは円（地球）を透過させて印を描画する

地球の半径=6,371 km
東京の経度=139°41′30″
東京の緯度=35°41′22″

＊出題者は自分で実装してません

**デフォルトの名無しさん** · 2021/07/24(土) 11:48:58.50

>>89
画面の上を地球上のどの方位にするとか、平行投影なのか透視投影なのか画角どうするかとか決めないといけないんでないの？

**デフォルトの名無しさん** · 2021/07/24(土) 11:53:24.04

>>89
>>90
画面の上端を北とします
投影方法と画角はよくわからないので自由とします＾＾；

**デフォルトの名無しさん** · 2021/07/24(土) 12:18:19.24

地球は球体じゃないから無理じゃない？

**デフォルトの名無しさん** · 2021/07/25(日) 16:25:19.73

地形データってどっかから取って来れるの？そこも自作しろと？それとも地形描画なしのツルツルの地球で良い？

**デフォルトの名無しさん** · 2021/07/26(月) 04:06:01.70

いいよ

**デフォルトの名無しさん** · 2021/07/27(火) 08:42:26.61

>>5の2　JavaScript版
let today = new Date()
console.log(
today.getFullYear().toString() + '/' +
today.getMonth().toString().padStart(2,'0') + '/' +
today.getDate().toString().padStart(2,'0') + '/' +
'第' + ( Math.floor( ( today.getDate() + 6 ) / 7 ) ).toString() +
['日曜日', '月曜日', '火曜日', '水曜日', '木曜日', '金曜日', '土曜日'][today.getDay()] + '/' +
today.getHours().toString().padStart(2,'0') + ':' +
today.getMinutes().toString().padStart(2,'0') + '/'
)

96 · 2021/07/28(水) 13:12:24.23

>>20 Perl5、DPで解けないかなと考えてたけど、すぐには思いつかずnaiveな解法

use List::Util 'max';
my $m;
sub f {
　my ($a, $b, $c) = @_;
　return () if $m and $m <= $c;
　my $aa = $a * $c - $b;
　return ($m = $c) if 0 == $aa;
　my @s;
　if (0 < $aa) {
　　my $bc = $b * $c;
　　@s = f($aa, $bc, max(int($bc / $aa), $c + 1));
　　@s = ($c, @s) if @s;
　}
　my @t = f($a, $b, max(int($b / $a), $c + 1));
　if (@s and @t) {
　　$s[-1] < $t[-1] ? @s: @t;
　} else {
　　@s ? @s: @t;
　}
}
$m = 0; my @a = f(37, 60, 1);
printf "37 / 60 = %s\n", join ' + ', map{"1/$_"} grep{$_} @a;
#$m = 0; my @a = f(3, 1, 1);

実行結果
~ $ perl 20_20_re_expan.pl
37 / 60 = 1/4 + 1/5 + 1/6

なお、この解法だと 3/1 = は時間がかかりすぎて解けない

**デフォルトの名無しさん** · 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

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

let g xs = Printf.printf "%.16f\n" @@ f xs
let () =
g [];g [1];g [1;2];g [1;2;3];g [1;2;3;4];g [1;2;3;4;5;6];
g [1;2;3;4;5;6;7];
g [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];
g [2;1;1;1;2;1;1;1;1;1;4;1;1;1;1;1;6;1;1;1;1;1];
g [0;4;1;1;3;4;5;9;7;16;9;25;11;36;13;49;15;64;17;81;19;100];
g [0;4;1;1;3;4;5;9;7;16;9;25;11;36;13;49;15;64;17;81;19];
g [0;4;1;1;3;4;5;9;7;16;9;25;11;36;13;49;15;64;17;81;19;100;21]
↓
0 / 1 = 0.0000000000000000
1 / 1 = 1.0000000000000000
3 / 1 = 3.0000000000000000
5 / 3 = 1.6666666666666667
9 / 7 = 1.2857142857142858
59 / 37 = 1.5945945945945945
233 / 151 = 1.5430463576158941
1855077841 / 1311738121 = 1.4142135623730951
4178 / 1537 = 2.7182823682498372
66628546560 / 21208521600 = 3.1415931679085070
12434780160 / 3958113600 = 3.1415925404465401
315324149760 / 100370793600 = 3.1415926730303347

**ハノン** ◆QZaw55cn4c · 2021/09/26(日) 19:57:37.78

>>197
＞[] -> 0, 1
＞| x :: [] -> x, 1
＞| x :: xs -> let n, d = h xs in x * d + n, d

>>193 を書いた私としてはなんとなくわかってしまうのが‥‥悲しいですね

**ハノン** ◆QZaw55cn4c · 2021/09/26(日) 19:58:46.43

つまり >>193 は言語として冗長だ、という点が悲しい、という意味です、誤解を招かないように補充します