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

**デフォルトの名無しさん** · 2022/11/13(日) 19:00:36.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/

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

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

**蟻人間** ◆T6xkBnTXz7B0 · 2023/05/14(日) 14:43:17.58

>>663
無視です。

**デフォルトの名無しさん** · 2023/05/14(日) 14:44:40.86

>>658
PowerShell 3以上

function bceil([UInt]$n)
{
　　if (!$n) {return 0}
　　$m, $n = [UInt]1, ($n - 1)
　　while ($n) {$m, $n = ($m -shl 1), ($n -shr 1)}
　　$m
}

function GCD([UInt]$a, [UInt]$b)
{
　　while ($b) {$a, $b = $b, ($a % $b)}
　　$a
}

2..20 + 40, 60 |% {
　　if (($a = bceil $_) -eq $_) {return "x/$_ = x*(1/$a)"}
　　$d= GCD $a $_
　　$b, $c = (($a - $_) / $d), ($a / $d * $_)
　　"x/$_ = x*(1/$a) + x*($b/$c)"
}

**デフォルトの名無しさん** · 2023/05/14(日) 14:44:57.54

[実行結果]
x/2 = x*(1/2)
x/3 = x*(1/4) + x*(1/12)
x/4 = x*(1/4)
x/5 = x*(1/8) + x*(3/40)
x/6 = x*(1/8) + x*(1/24)
x/7 = x*(1/8) + x*(1/56)
x/8 = x*(1/8)
x/9 = x*(1/16) + x*(7/144)
x/10 = x*(1/16) + x*(3/80)
x/11 = x*(1/16) + x*(5/176)
x/12 = x*(1/16) + x*(1/48)
x/13 = x*(1/16) + x*(3/208)
x/14 = x*(1/16) + x*(1/112)
x/15 = x*(1/16) + x*(1/240)
x/16 = x*(1/16)
x/17 = x*(1/32) + x*(15/544)
x/18 = x*(1/32) + x*(7/288)
x/19 = x*(1/32) + x*(13/608)
x/20 = x*(1/32) + x*(3/160)
x/40 = x*(1/64) + x*(3/320)
x/60 = x*(1/64) + x*(1/960)

**デフォルトの名無しさん** · 2023/05/15(月) 20:46:01.76

>>662
R
ideone.com/W0gYyh

合っているか分からないが、対称な盤面を除去すればWikipedia英語版に載っている
138通りに一致するから、合っていそう。

**デフォルトの名無しさん** · 2023/05/15(月) 20:54:26.02

>>662

(3^9)/4
= 4920.75

じゃね？

**デフォルトの名無しさん** · 2023/05/15(月) 21:03:02.31

>>668

(((3^9)-3)/4)+3
= 4923.

だったか

**デフォルトの名無しさん** · 2023/05/15(月) 21:10:44.54

>>667の情報を見て思ったが、日本では先攻が〇、英語版では先攻が×。
その時点でパターンが変わる。

**658** · 2023/05/15(月) 22:33:44.26

>>660
除数が2の整べき乗よりちょっと大きい場合は、引き算にした方が精度が出たりします。
tps://sagecell.sagemath.org/?z=eJx9j8FqwkAQQO-B_MOClyyahlmN8eLHSGlQ6DZL2A9wEw-KJ-uhlF7FioUIvbT0YD9mCMlnOIlGb8LusPN23uxsGEeSSSVHeswmUkWxpqxDO7QtqR5U_PQ4BDFgLUYR0zWmKabz4vuY_23RfJa7eZG9oXnHZGlbeggemQ4EvBI8CDB5vePZloonL9rRVTmaNZosn27QfJX_KzSzYrcsFwc6UGe4tO4K3g6fo5F2wK-B3-vx5jWq6Yo2XRAkRzTj9Ll7cepcBIKUmwN9FzyCttVixcdvfjiW-03-k2GSNP-jta1nvY4MHUbB1dzT_ApFBcUZngDsTonA&lang=python&interacts=eJyLjgUAARUAuQ==

**658** · 2023/05/15(月) 22:37:15.73

すぐ消えるかも知れない方のリンク(中身は同じ)
tps://sagecell.sagemath.org/?q=tkqhls

**658** · 2023/05/16(火) 04:50:29.58

補足：>>671の「引き算」は、>>658での「b)残り(除数の逆数-a)」を負数にする事に
相当します。

**デフォルトの名無しさん** · 2023/05/16(火) 09:51:49.00

>>673
こんな結果が出ましたが、この検証の仕方は合っていますか？
shortlink
https://sagecell.sagemath.org/?q=lftfab
permalink
https://sagecell.sagemath.org/?z=eJyFUs1OwkAQvpPwDg1c2AA2u_3BC89iiJFIYqUh-wBuMQTFixpjjCYcIIBiwHiRGMSHGQu9-QrOtlvUJmLSTme_mW929_tabdQdzXGdCt_Xao5bb3BcFfCtplOOu-U29nbLlG1rWQ0jNC-h2YRme_W88F8HIIbBqL2aXIO4Aa-TTrmN2iHPZf7tzJB0Kp3iZarjPjlaInK8TkvgXWzgxvO5bAdxCWLiH_VBjIP3cxDHq1EnOJligpOpGm0wkq8e1Cs8R60QsEyTxLthj8HyWEAQOSw-jk2KihOuWYkh5ZtD7SLVEUSK8YOimi2b5KPUNCxGkjTLxqqsSAGy2up25k8XwUPff5mA58Ua4TMI77q-slHQMBQ50TlZg0yCLAFSCdIYXHviVHb3a4d7kaRl20yIHO88VLq2W8urJ4zBsIemRCda3j36b1gdx76cfcxOMYLogegiC0Tnc979w8Sx3-2AJ4LFHDwcMlz274LRXP4J3DZ3aDnUUVfCl4jC2W881DnKN9iCPONvXuRQtAg92mSSElVOlLrKb0JveUhVSloh76VK9HcpA-IeRAvVypAv4NdijA==&lang=python&interacts=eJyLjgUAARUAuQ==

**デフォルトの名無しさん** · 2023/05/16(火) 09:59:14.32

元のお題に「分ける」とあるのてマイナスなしなんじゃないの？

**デフォルトの名無しさん** · 2023/05/16(火) 15:16:43.16

回分かどうか判定するプログラムを書きなさい。回分とは逆から読んでも同じ文である。
例　step on no pets
使用言語　C

**デフォルトの名無しさん** · 2023/05/16(火) 18:13:35.98

回文ね

**デフォルトの名無しさん** · 2023/05/16(火) 18:38:20.59

>>676 補足
入力した文が回分かどうか判定するプログラムを書きなさい。

**デフォルトの名無しさん** · 2023/05/16(火) 20:23:20.23

>>676,678
これ後から絵文字やらクラスタやらどんどん補足していくパターン

**デフォルトの名無しさん** · 2023/05/17(水) 02:03:28.90

>>678
これなんの補足をしたつもりだったんだ…？

**デフォルトの名無しさん** · 2023/05/17(水) 10:28:21.37

>>676

Haskell

isPalindrome x = x == ( reverse x )

main = do
print $ isPalindrome "step on no pets"

**667** · 2023/05/17(水) 23:14:17.75

>>662
最終盤面が何通りか理論的に求めてみた。○を先手とし、最終盤面での○と×の個数をそれぞれa, bとすると、(a, b) = (3, 2), (3, 3), (4, 3), (4, 4), (5, 4) の5種類に場合分けされる。

(1) (a, b) = (3, 2) のとき
・○は縦・横・斜めの8本のどれかに並べて置くので8通り
・×は残り6マスのどれかに置くので ₆C₂ = 15通り
・以上より、盤面数は 8 × 15 = 120通り

(2) (a, b) = (3, 3) のとき
・○は縦・横・斜めの8本のどれかに並べて置くので8通り
・×は残り6マスのどれかに置くので ₆C₃ = 20通りだが、○が縦/横に並ぶ場合は、×が縦/横の残り2本のどちらかに並ぶ2通りを除外する
・以上より、盤面数は 8 × 20 － 6 × 2 = 148通り

(3) (a, b) = (4, 3) のとき
・○は3個を縦・横・斜めの8本のどれかに並べ、1個を残り6マスのどれかに置くので、8 × 6 = 48通り
・×は残り5マスのどれかに置くので ₅C₃ = 10通りだが、○が縦/横に並ぶ場合は、×が縦/横の残り1本に並ぶ1通りを除外する
・以上より、盤面数は 48 × 10 － 6 × 6 × 1 = 444通り

(4) (a, b) = (4, 4) のとき
・×は3個を縦・横・斜めの8本のどれかに並べ、1個を残り6マスのどれかに置くので、8 × 6 = 48通り
・○は残り5マスのどれかに置くので ₅C₄ = 5通りだが、○が縦/横に並ぶ場合は、×3個が縦/横の残り1本に並び1個が残り2マスのどちらかに並ぶ2通りを除外する
・以上より、盤面数は 48 × 5 － 6 × 6 × 2 = 168通り

(5) (a, b) = (5, 4) のとき
・×は9マスのどれかに置くので ₉C₄ = 126通りだが、3個を縦・横・斜めの8本のどれかに並べ1個を残り6マスのどれかに置く 8 × 6 = 48通りを除外する
・○は残り5マスに置くので1通り
・以上より、盤面数は (126 － 48) × 1 = 78通り

**667** · 2023/05/17(水) 23:14:48.75

(1)～(5)を合計すれば、盤面数は 120 + 148 + 444 + 168 + 78 = 958通り

よって>>667のプログラムでの計算値が正しいことが検証された。>>667のプログラムの行列Sを調べると、(1)～(5)の各場合の盤面数も上記の値と一致した。上のように場合分けして盤面を組み立てていく手順をプログラムで記述するのは複雑で面倒なので、>>667のようにゲームを実際に行ってみるプログラムの方が簡潔に書けて良い。

**デフォルトの名無しさん** · 2023/05/18(木) 14:03:59.67

>> 676 の回答例を載せときます
#include <stdio.h>
#include <string.h>
void kaibun_judge(char *str);
int main(void)
{
char ss[80];

printf("文字列を入力してください ");
gets(ss);
kaibun_judge(ss);

}
void kaibun_judge(char *str)
{
char *p, *q, dmy[strlen(str)];
p = str;
q = dmy;

for (int i = strlen(str)-1;i >= 0;i--) {
dmy[i] = *p++;
}
if (strcmp(str, dmy) == 0) printf("回文です\n");
else printf("回分ではありません\n");

}

**667** · 2023/05/18(木) 21:14:06.06

>>682の(4)の訂正
(誤)○が縦/横に並ぶ場合は、×3個が
(正)×が縦/横に並ぶ場合は、○3個が

**デフォルトの名無しさん** · 2023/05/20(土) 00:00:46.48

>>684
よくあるバグ入でかつ分かりやすい速度向上の余地を持つ良い例だ

**デフォルトの名無しさん** · 2023/05/20(土) 05:48:17.89

>>684
そのプログラムは char 一つで1文字を表す文字コード以外ではまともに動かない。
80文字以上の入力でバッファーオーバーランする。
dmy[] 配列の最後の次のメモリに \0 以外のゴミが入っていたら strcmp() が想定した通りに動かない。

**デフォルトの名無しさん** · 2023/05/20(土) 05:50:36.82

C使い慣れてないなら他の文字列が楽々使えてメモリ管理についてほとんど考えなくても良い言語使って作った方が良いのではないか？

**デフォルトの名無しさん** · 2023/05/20(土) 06:06:53.76

>>662
Kotlin

馬鹿正直に○×と交互に打ちながら途中で縦横斜め何れかで3つ揃ったらそこで終了という感じで全パターンやらせたら
255168 パターンになったのだが、これで合っているのか謎だ。対象になるパターンを抜くとすると4で割って 63792 になる。
それと出力が多すぎて恐らくどのサイトでも実行されると結果がまともに出ない。
paiza.io にソース置いておくが実行するとパターン 1651 の途中で Timeout して中途半端に終わっている。

https://paiza.io/projects/2eV3pNvLM7kdKO0gYfHc-w

**デフォルトの名無しさん** · 2023/05/20(土) 06:14:37.20

よく見たら >>676 は使用言語がCだけになってたんだな。
じゃあC使うのは仕方がないが、Cだと色々と気を付けなければならないことが多いので要注意だ。

**デフォルトの名無しさん** · 2023/05/20(土) 17:36:52.34

>>676
C

Locale 考慮して wchar_t 使って書いた。日本語読ませても処理できる。
但しやってることは文字列の順序入れ替えて比較しているだけなので>>684と考え方は同じ。
(書きながら一々メモリ確保して逆順に並べなくても出来ると気付いたがそれは後で公開する)。
入力元が Unicode になっていた場合は濁点半濁点の合成文字みたいなのもあって MacOS でやった時とその他 OS で違ってしまうかも知れないが、その辺は考慮されていない。

https://paiza.io/projects/YXLsX3FqAsVF6DVDxye7bA

**デフォルトの名無しさん** · 2023/05/20(土) 17:50:35.91

>>676
C

>>691の改良版。リバースした文字列をわざわざ作らず半分を前後から比較するだけ。

https://paiza.io/projects/VPN0R2nAEGYugRnNLiVBfQ

**692** · 2023/05/20(土) 18:08:07.96

>>692は改行のみが入力された時の反応がおかしかたのでちょっとだけこっそり修正した。
同様に>>691もおかしかったがそっちは面倒なので放置。

**デフォルトの名無しさん** · 2023/05/21(日) 09:10:07.45

>>684 gets()のかわりにfgets()を使うと改行を含んでしまうので\0で埋めてやらないといけない。

**デフォルトの名無しさん** · 2023/05/21(日) 11:15:13.97

文字列がらみはもはやライブラリ使わないと無理やろ
ユニコードのコーディングルールはもう複雑すぎて素人が自作で作ってもまずうまくいかない
この手のルールは数学的にこうなってるはずだ、こうなってないとおかしいとかいうロジックが何一つ通用しない、単なる取り決めるなのでただひたすら覚えるしかないけど、もはや色んな言語全部取り込むために魑魅魍魎の世界

**デフォルトの名無しさん** · 2023/05/22(月) 00:54:31.79

>>676
U+10000以上の文字、異体字セレクタ、結合文字に対応したプログラムはPowerShellでは
簡単に書ける。ただし、IDS（漢字構成記述文字列）には対応していない。
tps://tio.run/##XZJPS8MwGMbv/RQxVGjRFnvxMPA4hjdhw8sYoxvZVumS2WQqzoKtf5kHD8KGG4i3zcNQbx7ULxM63E3wA8ykm1sxl6bv8@T3Pm/TBjlEHq0h151OK01cZg7BgDLPQweaSnWlpQCxVAa2AITxPp9xScl2nWNbes0s8xxc3cYVUkilMojl0BFLu6iOMEvjZh15NiOeRIGTVdCKQWoRrAmkP0crvqL8a21tJpobe8TBQKV5wzJNQ6Wmi3CV1QrxOeIhu1wDwg6ECTYQo0BgMBEk1IDrAE5GPR4Mv9qvP4/nAQ9ueNj7fnr74MFg9pCFhToUbqiDWeMN07Ti1PGbXBrUxp2X6P1Wl2Rt3B3Eez2vFgsLUzyjUwHailoUqNlEIr8PkEvRX8HalKXFIRhdnI07V9GoG113U0KCS2lyepmUGEzkifoPQhPD8OCZByEP2zz45EGfh3cyYkK@lzEpMMpoXzCWcefX6svPqe5kd8XvIC4iZ5dcNJ3@Ag

これはOutputの下の方に表示されているようにLinuxでの実行結果だが、Windowsと同様に
動作しOSに依存しない。CではOSに依存しないで書くのは難しそう。

(旧式)はUTF-16のワードを逆順に並べただけなので、U+10000以上のU+29E3D、異体字
セレクタ使用のU+82B1 U+E0101、結合文字のU+304B U+309Aが含まれる2番目の例では
正しく判定できていないが、(新式)はできている。異体字は対応しているフォント（例えば
Windowsの標準フォントでは游ゴシック）でないと正しく表示されない。

**デフォルトの名無しさん** · 2023/05/22(月) 17:00:26.90

>>696 ワイド文字対応だけでもすごいのにそこまでやるとは

**デフォルトの名無しさん** · 2023/05/23(火) 00:39:57.91

>>696をIDS（漢字構成記述文字列）対応にした。ビャンビャン麵のビャンという漢字の
IDSを正しく表示できるのは、源ノ等幅、源ノ明朝、源ノゴシックのバージョン1.xxxと
Visual Studio Code, Edgeなどの組み合わせに限られる。正しく表示されるとこうなる。
://i.imgur.com/TEDafUA.png
バージョン2.xxxxではビャンという漢字に単独のコードU+30EDEが割り当てられたので、
IDSでは表示できなくなった。

function strrev($s)
{
　　$m = 2, 2, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2
　　$p = $q = ""
　　$n = 0
　　[Globalization.StringInfo]::GetTextElementEnumerator($s) |% {
　　　　if (0x2FF0..0x2FFB -contains $_[0]) {
　　　　　　if (!$n) {$n = 1}
　　　　　　$n += $m[$_[0] - 0x2FF0]
　　　　}
　　　　$p += $_
　　　　if (!$n -or !--$n) {$p, $q = "", ($p + $q)}
　　}
　　$p + $q
}

foreach ($s in "⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麵で麵⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心", "⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麵で麺⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心") {
　　$t = strrev $s
　　"元文字列: $s"
　　"逆文字列: $t"
　　("回文ではありません。", "回文です。")[$s -ceq $t]
　　""
}

**デフォルトの名無しさん** · 2023/05/23(火) 00:40:48.32

[実行結果]
元文字列: ⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麵で麵⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心
逆文字列: ⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麵で麵⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心
回文です。

元文字列: ⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麵で麺⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心
逆文字列: ⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心麺で麵⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心⿺辶⿳穴⿰月⿰⿲⿱幺長⿱言馬⿱幺長刂心
回文ではありません。

**デフォルトの名無しさん** · 2023/05/23(火) 06:35:45.73

文字列はカオスだわ
https://jp.emeditor.com/forums/topic/mac%E3%81%A7%E5%85%A5%E5%8A%9B%E3%81%95%E3%82%8C%E3%81%9F%E3%83%8F%E3%83%B3%E3%82%B0%E3%83%AB%E6%96%87%E5%AD%97%E3%81%AE%E8%A1%A8%E7%A4%BA%E3%81%AB%E5%AF%BE%E5%BF%9C%E3%81%97%E3%81%A6%E6%AC%B2/

>>698
これは対応が難しい？

https://tio.run/##pVPBattAEL3rKzZCBZtYwnYhB0MvhST0VkjoxZigmHWsIq8UadOGpgbLcpzWJZik0NAE7B4KTgihLbnkUvdb3EWl6cWGfoA7u4ol1fTWZdGO5r15M5od2dZz7LhVbJrTaWWHlKlhEeRSx8HPUoqblvYkBEuh6AGSZWEXV01rUzeNFzrnamvUMcjWI1KxSoXCKqbreJcum7iGCV0mOzXs6NRyuBR6eQ/tCSFlAy2CZP1OWqpL0lzqfDJ3DULyGb7vi53/5w65NpffjotVCNjZ/69bKPBlVFAqu5tfWclqmjgfIrVsEaobxIUPK2ZL6QR7FrGgEHCLYnL1v1DwLULJtaKIRSoKxUsRKabDx3HqhjSnjFTLQQuqGuawM7MGZFCKh8BrOhSpz3rEfbzr823PLSX6rj61DIIUt6jmNE1VXM3EZItWS@K6LAfr5SokcBGQZBtTF4EMsUAJ25Bavr06Zd75z86X3x9aHvPesObp5OLmK/MG4cEdEXoObB7E/NfM70RP7hoPW@NhbzzcnzM4NukPjib97rfvjcNJ/@giNgeROTjmJicH3Y/M7zL/mjU/M/8SjVrHo3Zv1L78ddOIbASIwK@BCxFy8jbF8N4NKDQm8svBvv/j3UFwdRK8OikAIsfQbaOdhGgMpeTgrAcYdIR5n5jXZM0O84bMO2PNt6LgGH4vp4vQarWMt0EjHg6YcrgN5fHaE/iJ4R7X9U0TT6d/AA

**デフォルトの名無しさん** · 2023/05/23(火) 21:04:58.28

>>700
こっちのPowerShellでやれば全部OK
https://github.com/PowerShell/PowerShell
Name Value
---- -----
PSVersion 7.3.4
PSEdition Core

.NET 7凄すぎ

**デフォルトの名無しさん** · 2023/05/27(土) 01:36:20.91

お題：現在の日時の文字列が与えられるのでその文字列の先頭から秒数分を切り出して出力せよ

**デフォルトの名無しさん** · 2023/05/27(土) 02:11:47.73

>>702
それって秒の部分だけが残れば良いってこと? たとえば 2023/05/27 02:03:45 なら最後の 45 だけになれば良いの?
もしそうなら sed 使ってこうすればできるな。

sed 's/^.*://'

YYYY/MM/DD hh:mm:ss 形式ではなく秒が最後に出て来ない形式だったとしても最後のコロンの直後の数字が秒で、
その秒の数字文字列の後に数字でない文字が必ず来るのだとすれば (例えば Sat May 27 02:05:47 JST 2023 の
ような形式だとすれば)こうすれば良い。

sed 's/^.*://;s/[^0-9].*$//'

**デフォルトの名無しさん** · 2023/05/27(土) 02:42:42.27

サンプルケース書いてくれ

**デフォルトの名無しさん** · 2023/05/28(日) 09:45:16.03

お題　ナイトツアー

5×5の狭いチェスボードのどこかにナイトのコマを置きそこからナイトの可能な動きでコマを動かします
ただし同じマスは2回通れません
このとき全マス25マスを全て通る手順の数を数え上げて下さい
しかし実際全部数えるのは流石に無駄なので対称性を利用して
♯{A1→B3}×8
+ ♯{A3→B5}×8 + ♯{A3→C4}×8
+ ♯{B2→C4}×8
+ ♯{C3→D5}×8
はよいとします

例
♯{A1→B3} = 152

--ナイトの動き--
桂馬の動きと同じ、水平方向の移動量+垂直方向の移動量 = 3となる移動(もちろん盤の中)

**デフォルトの名無しさん** · 2023/05/28(日) 11:22:07.87

>>705

訂正
> --ナイトの動き--
水平方向の移動量+垂直方向の移動量 = 3となる移動で水平、垂直な移動はなし(もちろん盤の中)

**デフォルトの名無しさん** · 2023/05/28(日) 23:12:57.01

>>705
対称性を考慮するのは面倒なので全部数え上げる。Rで書いたら5秒以内に到底終わらなかったので、
ideone.com/y5Dwy7
C++に移植した。盤面表示はideoneでは途中で切れてしまうので非表示にした（関数moveの
コメント部分をコメントでなくせば表示される）。
ideone.com/6Muw9F

対称性を考慮すれば、手順数は1728 / 8 = 216通りになる。

**蟻人間** ◆T6xkBnTXz7B0 · 2023/05/29(月) 01:40:51.78

お題: 水の流れる複数のタンクにおいて水の流れをシミュレーションしたいので、以下に定義する入力文字列でタンクと
水の流れなどを記述する。入力文字列は次のような形式の行の並びとする。

「(タンクの名前), (タンクの底面積[m^2]), (タンクの底からの排水口までの高さ[m]), (次のタンクの名前)」

各行は各タンクの情報とその次のタンクを表す。入力文字列の最初の行は最初に水が入るタンクに対応する。最初のタンクには外部から水が0.5[m^3/min]の速度で入る。タンクの形状は円錐か四角錐とし、底面は平面状とする。
各タンクには排水口が付いているかも知れない。排水口があるのならそれはタンクの底より高く、タンクの側面に付いている。
排水口は次のタンクにつながっている。次のタンクはその前のタンクより低い位置にある。各タンクにおいて排水口に水が届いていないときは水は排水されない。
タンクとタンクの間の長さと水の落差にかかる時間は無視して構わない。次のタンクがN/Aならそのタンクに排水口と排水先はない。

シミュレーション開始後、1分ごとに各タンクに何[m^3]の水が入っているかを出力せよ。タンクCが満タンなら終了せよ。

(入力)
A, 3, 0.5, B
B, 2, 0.7, C
C, 1, 0.4, N/A

**デフォルトの名無しさん** · 2023/05/29(月) 01:46:57.43

>>705
Kotlin

馬鹿正直に再帰で全パターンやっているだけ。
途中経過は出力しない。既に置かれたマスは1bitで記録されるようにして 5 * 5 = 25 bit で Int (32bit) に収まるのでそれを利用している。
恐らく Boolean の配列を使って書くとそこまで最適化されなさそうなのでこうした。

https://paiza.io/projects/O4Ejxs76h64GSeXYCVs2mg

結果は同じ 1728 になったので多分これで合っていると思う。

**デフォルトの名無しさん** · 2023/05/29(月) 23:38:05.56

>>705
>>707にちょっと付け加えたら対称性を考慮して数え上げることができた。
://ideone.com/NSyJcK (R)
://ideone.com/S1VIUA (C++)

Rでは5秒以内に終わらなかった。

**蟻人間** ◆T6xkBnTXz7B0 · 2023/05/30(火) 20:06:49.36

お題: x-y平面上に摩擦のないすべり台がある。重力加速度を9.8[m/s^2]とし、重力はy軸とは逆向きに作用している。
この滑り台は動かないように固定されている。横位置x>=0[m]についてf(x)がすべり台の面の高さ[m]を表す。
体重35kgの子どもがx=0の場所から滑り台を滑り始めた。次の場合の5秒後の子どもの位置ベクトルと速度ベクトルを求めよ。空気抵抗は無視する。

(i) f(x)=10-sqrt(x)。
(ii) f(x)=10-sqrt(sqrt(x))。

※ xの平方根をsqrt(x)とする。

**デフォルトの名無しさん** · 2023/06/01(木) 09:25:39.43

お題

アルファベットが入力されます
> [A, A, B, B, B, C, C, D, D, D, E]

出力上限が入力されます
> 7

出力上限以内のアルファベットを出力してください
> [A, A, B, B, B, C, C]

ただし、出力上限で切った場合に同じアルファベットが分割される場合は出力上限を超えて
同じアルファベットが続く限り出力してください

アルファベット: [A, A, B, B, B]
出力上限: 3
出力: [A, A, B, B, B]

**デフォルトの名無しさん** · 2023/06/01(木) 21:49:43.70

>>712
C#でループを使わずに出力
www.ideone.com/Xk3ACt

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

お題：矩形に整形された全角文字のテキストが与えられる。左上から時計回りに
渦巻き状に読んだ文字列を、元の矩形と同文字数で桁数がより多く最も近い矩形に
左上から時計回りに渦巻き状に並べて出力せよ。

[例]
プログラ　　→　　プログラミン
題スレミ　　　　　レス題おのグ
おのグン

[問題]
メロスは激怒した。必ず、かの邪智暴虐の王を除かなけ
ては、人一倍に敏感であった。きょう未明メロスは村れ
し、内気な妹と二人暮しだ。この妹は、村の或る律をば
対の嫁の衣裳やら祝宴の御馳走やらを買いに、は気出な
に六花ロスには竹馬の友があった。セリヌンテるな発ら
悪十、メてみるつもりなのだ。久しく逢わなィば一しぬ
邪。え。ねスは、まちの様子を怪しく思っかウる牧、と
もいゆた訪ロのは当りまえだが、けれどたっス市人野決
ど無れいらメい全体が、やけに寂しいも。たでにをを意
れもそ歩かに暗市不安になって来た。、ひのあや、越し
け房、られちの、ん衆をつかまえ。のなっだるっ近えた
。女はぶこうちくだい若たっ逢で路んんそか。て々山。
た。スら、るま無んだ、もスロメなきだりら今来、越メ
来いロぶをい、はでりかばいせの夜、かし、はた花えロ
て無メを友ててち落も日に既うも。るいて訪此の婿、ス
しも。路のい歩。るあでみし楽がのく行てねのだと十に
暮母る大そ。るいてしを工石、で市のスクラシ。し里は
で、あの都らかれそ、め集い買を々品のそ、ず先ては政
んもでのなか近間も式婚結。たいてっなに事るえ迎な治
遊父はにスロメ。た来てっやに市のスクラシの此たれが
と羊、き吹を笛。るあで人牧の村、はスロメ。ぬらかわ

**デフォルトの名無しさん** · 2023/06/02(金) 00:02:20.01

>>714
この部分の意味がわからない。どうしろというのか？

> 元の矩形と同文字数で桁数がより多く最も近い矩形

**デフォルトの名無しさん** · 2023/06/02(金) 00:10:19.23

>>715
[例]で言えば、3行×4桁=12文字が矩形に並ぶのは1, 2, 3, 4, 6, 12桁のときだから、
4桁より多く4桁に最も近いのは6桁となる。

**デフォルトの名無しさん** · 2023/06/03(土) 18:04:48.61

>>712
Kotlin
https://paiza.io/projects/o3Y5W0EXe0CxGFZTAXhtAQ

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/04(日) 14:05:22.78

お題: 郵便切手の問題。

切り離されていない、一列に並んだn枚の切手を考える。
一枚の切手の上に全ての切手を折り込む時、
左端の切手を表向きで一番上に折り込む方法は何通りか計算せよ。

**638** · 2023/06/04(日) 19:16:56.02

>>714 Perl5
https://www.ideone.com/XpTafZ

**638** · 2023/06/04(日) 20:52:10.72

>>712 Perl5

use feature qw{signatures say};
sub f($n, $c, @s) {
　if (1 < $n) {
　　$c, f($n - 1, @s);
　} elsif (1 == $n) {
　　$c, f(0, $c, @s);
　} else {
　　$c, f(0, $c, splice @s, 1) if $c eq $s[0];
　}
}

say f 7, qw[A A B B B C C D D D E];
say f 3, qw[A A B B B];

※見易くするためインデントは全角スペースに置換してあります。

実行結果

$ perl 21_712.pl
AABBBCC
AABBB

**638** · 2023/06/04(日) 23:24:30.12

>>712 Perl5、もっと簡潔なやり方あるんだった…orz

use feature qw{signatures say};
sub f($n, $s) {
　my $r = '^.{' . ($n - 1) . '}(.)\1*';
　$s =~ /$r/;
　$&;
}
say f 7, 'AABBBCCDDDE';
say f 3, 'AABBB';

実行結果
$ perl 21_712_rex.pl
AABBBCC
AABBB

**デフォルトの名無しさん** · 2023/06/05(月) 17:20:26.45

お題　逆順に

文字列が与えられます
使われている文字のリストを作って並べその“反対の”文字を対応させます
例えば登場文字列が
HELLO WORLD
なら出てくる文字をコード順に並べると
" DEHLORW"
であるので
E → O, L →H, ' '→W
となります(空白、句読点等一切区別せず一文字と考えます)
この対応で入力された文字列を変換してください

例　
HELLO WORLD →
LOHHEW EDHR

⑨⑫⑦②⑤⑭⑰⑳㉔①⑩⑬㉒⑲④⑮⑪⑧㉓⑥③⑯㉑⑱ →
⑯⑬⑱㉓⑳⑪⑧⑤①㉔⑮⑫③⑥㉑⑩⑭⑰②⑲㉒⑨④⑦

**デフォルトの名無しさん** · 2023/06/05(月) 19:56:23.90

>>722
Java
https://paiza.io/projects/ThIWBUU33EC_Tv88KPog8A

**638** · 2023/06/05(月) 20:41:22.90

>>720 もっとすんなり書けるんだった… orz

　　$c, f(0, $c, splice @s, 1) if $c eq $s[0];
　　↓
　　$c, f(0, @s) if $c eq $s[0];

関数全体のコードとしてはこうなる：

sub f($n, $c, @s) {
　if (1 < $n) {
　　$c, f($n - 1, @s);
　} elsif (1 == $n) {
　　$c, f(0, $c, @s);
　} else {
　　$c, f(0, @s) if $c eq $s[0];
　}
}

すっきりした形になりました

**638** · 2023/06/05(月) 21:31:53.30

>>722 Perl5

use feature qw{signatures say};
use utf8; binmode STDOUT => 'utf8';
use Encode 'encode';

sub f($s) {
　@s = split '', $s;
　%h = map{unpack('H*', encode('utf-8', $_)) => $_} @s;
　@i = reverse @c = sort keys %h;
　%o = map{$h{$c[$_]} => $h{$i[$_]}} 0..$#c;
　say map{$o{$_}} @s;
}

f 'HELLO WORLD';
f '⑨⑫⑦②⑤⑭⑰⑳?①⑩⑬?⑲④⑮⑪⑧?⑥③⑯?⑱';

※見易くするためインデントは全角スペースに置換してあります。

実行結果:
$ perl 21_722.pl
LOHHEW EDHR
⑯⑬⑱?⑳⑪⑧⑤①?⑮⑫③⑥?⑩⑭⑰②⑲?⑨④⑦

**638** · 2023/06/05(月) 21:35:18.29

>>722 Perl5 ｽﾏｿ、ﾋﾞｭｰﾜｰから書き込んだら丸数字の一部が文字化けしたのでﾌﾞﾗｳｻﾞから再書き込み

use feature qw{signatures say};
use utf8; binmode STDOUT => 'utf8';
use Encode 'encode';

sub f($s) {
　@s = split '', $s;
　%h = map{unpack('H*', encode('utf-8', $_)) => $_} @s;
　@i = reverse @c = sort keys %h;
　%o = map{$h{$c[$_]} => $h{$i[$_]}} 0..$#c;
　say map{$o{$_}} @s;
}

f 'HELLO WORLD';
f '⑨⑫⑦②⑤⑭⑰⑳㉔①⑩⑬㉒⑲④⑮⑪⑧㉓⑥③⑯㉑⑱';

実行結果
$ perl 21_722.pl
LOHHEW EDHR
⑯⑬⑱㉓⑳⑪⑧⑤①㉔⑮⑫③⑥㉑⑩⑭⑰②⑲㉒⑨④⑦

**デフォルトの名無しさん** · 2023/06/05(月) 22:09:27.65

>>722
PowerShell 4以上 (U+10000以上の文字、結合文字、異体字セレクタに対応)

foreach ($s in "HELLO WORLD", "⑨⑫⑦②⑤⑭⑰⑳㉔①⑩⑬㉒⑲④⑮⑪⑧㉓⑥③⑯㉑⑱", "魚と花󠄁か゚𩸽") {
　　$c = [Globalization.StringInfo]::GetTextElementEnumerator($s) | sort | unique
　　$h = @{}
　　0..($c.length - 1) |% {$h[$c[$_]] = $c[-$_ - 1]}
　　-join $h[[Globalization.StringInfo]::GetTextElementEnumerator($s)]
}

[実行結果]
LOHHEW EDHR
⑯⑬⑱㉓⑳⑪⑧⑤①㉔⑮⑫③⑥㉑⑩⑭⑰②⑲㉒⑨④⑦
と魚花󠄁𩸽か゚

**638** · 2023/06/06(火) 01:24:29.34

>>726 しくった
utf-8文字を可変長16進数にunpackして辞書的順番でsortするので、例題では意図した通りの結果を得られたが、
半角と全角などutf-8でバイト数の異なる文字が混在する文字列の場合には、順番が違ってくるだろうな…

**638** · 2023/06/06(火) 02:40:16.03

>>726 Perlで文字の数値コードを返すord()関数は半角文字の場合にASCIIコードを返すが、
全角Unicode文字ではコードポイントを返す、ということなので、素直にそれを使い、
文字列としてっではなく数値としてsortし

　%h = map{unpack('H*', encode('utf-8', $_)) => $_} @s;
　@i = reverse @c = sort keys %h;
　　↓を
　%h = map{ord($_) => $_} @s;
　@i = reverse @c = sort{$a <=> $b} keys %h;

と修正するのがよさげ。連投ｽﾏｿ

**デフォルトの名無しさん** · 2023/06/06(火) 20:31:50.40

>>714 octave
https://ideone.com/MXrAYt
m2c = @(m, r, c) mat2cell(m, repelem(r, size(m, 1) / r), repelem(c, size(m, 2) / c));
s2c = @(s) m2c(s, 1, 3);
function v = e2v(e, v)
if nargin < 2
v = e2v(e, []);
elseif ~isempty(e)
v = e2v(rot90(e(2:end, :)), [v e(1, :)]);
end
end
function ind = eddyind(varargin)
if nargin == 2
ind = e2v(reshape(1:times(varargin{1:2}), varargin{1:2}));
elseif nargin == 1
ind = eddyind(size(varargin{1}, 1), size(varargin{1}, 2));
end
end
s = 省略
c = s2c(s);
n = numel(c);
m = find(rem(n, 1:n) == 0 & size(c, 2) < 1:n)(1);
d = cell(n / m, m);
d(eddyind(d)) = c(eddyind(c));
disp(cell2mat(d));

**デフォルトの名無しさん** · 2023/06/06(火) 21:42:34.58

>>722 ruby
https://ideone.com/yqb6sd
f = ->(s) {s.chars.uniq.sort.join.tap{|u| break s.tr(u, u.reverse)}}
p f.('HELLO WORLD')
p f.('⑨⑫⑦②⑤⑭⑰⑳㉔①⑩⑬㉒⑲④⑮⑪⑧㉓⑥③⑯㉑⑱')
#p f.('HELLO WORLD-') # trのパターンになっちゃう文字があるとダメ

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

>>722
>>727では大・小・全角・半角文字が区別されないのを修正した。PowerShell 2でも正しく動作するようになった。

foreach ($s in "HELLO WORLD", "⑨⑫⑦②⑤⑭⑰⑳㉔①⑩⑬㉒⑲④⑮⑪⑧㉓⑥③⑯㉑⑱", "魚と花󠄁か゚𩸽", "HeLLO ＷoRＬD") {
　　$c = [Globalization.StringInfo]::GetTextElementEnumerator($s) | sort -unique {[int[]][char[]]$_}
　　$h = New-Object Collections.Hashtable
　　0..($c.length - 1) |% {$h[$c[$_]] = $c[-$_ - 1]}
　　-join $h[[Globalization.StringInfo]::GetTextElementEnumerator($s)]
}

[実行結果]
LOHHEW EDHR
⑯⑬⑱㉓⑳⑪⑧⑤①㉔⑮⑫③⑥㉑⑩⑭⑰②⑲㉒⑨④⑦
と魚花󠄁𩸽か゚
oLeeRＷ HODＬ

**デフォルトの名無しさん** · 2023/06/07(水) 23:16:05.89

お題
A,B,Cがランダムにいくつか入力されます

隣り合う文字が異なるときそれらは同じグループです
A,B → [A,B]

隣り合う文字が同じ時それらは異なるグループです
A,A → [A],[A]

A,B,A,A,A,B → [A,B,A],[A],[A,B]

もっとも長いグループを出力してください

入力
B,A,B,C,B,B,B,A,C,C,A,B,A,A,C,A,B,C,C,C,B,C,A,B,C,A,A,A,C,B

出力
[C,B,C,A,B,C,A]

**638** · 2023/06/08(木) 01:33:27.84

>>733　Perl5

use feature qw{:5.16 signatures say};
@a = ([]);
sub ($c, @s) {
　push @{$a[-1]}, $c;
　push @a, [] if $c eq $s[0];
　__SUB__->(@s) if @s;
}->(qw{B A B C B B B A C C A B A A C A B C C C B C A B C A A A C B});
$c = (sort{@$a <=> @$b} @a)[-1];
$" = ','; say "[@$c]";

※見易くするためインデントを全角スペースに置換してあります。

実行結果
$ perl 21_733.pl
[C,B,C,A,B,C,A]

**デフォルトの名無しさん** · 2023/06/08(木) 10:35:46.63

>>712
Kotlin

抜き出す数の指定が 0 以下の場合は空のリストになるようにした。

https://paiza.io/projects/oE3GVsi4b0sNqSXQ78_7gg

**デフォルトの名無しさん** · 2023/06/08(木) 21:24:06.92

>>733
PowerShell

$s = "B,A,B,C,B,B,B,A,C,C,A,B,A,A,C,A,B,C,C,C,B,C,A,B,C,A,A,A,C,B"
@($s -replace "(.),(?=\1)", "`$1/" -split "/" | group length | sort -descending {[int]$_.name})[0] |% {$_.group}

[実行結果]
C,B,C,A,B,C,A

**デフォルトの名無しさん** · 2023/06/08(木) 22:05:38.22

>>733 ocaml
https://ideone.com/jfFx3Y
open List
let (<<) f g x = f (g x)
let flip f x y = f y x
let f =
let aux x = function
[] :: ys -> [x] :: ys
| (z :: _) as y :: ys -> if x = z
then [x] :: y :: ys
else (x :: y) :: ys
in rev << hd << rev << map snd << sort compare << map (fun x -> (length x, x)) << fold_left (flip aux) [[]]
let g = print_newline << iter print_char << f
let () = g ['B';'A';'B';'C';'B';'B';'B';'A';'C';'C';'A';'B';'A';'A';'C';'A';'B';'C';'C';'C';'B';'C';'A';'B';'C';'A';'A';'A';'C';'B']

>>733 ruby
f = ->(a) {a.each_with_object([[]]) {|x, acc| acc[-1][-1] != x ? acc[-1] << x : acc << [x]}.sort_by(&:size)[-1]}
p f.(%w(B A B C B B B A C C A B A A C A B C C C B C A B C A A A C B))

**デフォルトの名無しさん** · 2023/06/08(木) 22:14:54.10

お題：1234567以下の素数 a, b, c で a + b + c = 1234567 かつ a ≤ b ≤ c を満たす組み合わせは何通りあるか計算せよ。

**デフォルトの名無しさん** · 2023/06/09(金) 04:35:40.14

>>738
https://ideone.com/4YftKh

力業っぽくなったな、3で割って1余る整数計しか解けない

**デフォルトの名無しさん** · 2023/06/09(金) 08:47:59.66

実質力技しかないやろな
割ったあまりとかで多少減らせる程度にしかならない

**デフォルトの名無しさん** · 2023/06/09(金) 08:55:35.77

入力nに対してn以下の素数はn オーダーある
最初に素数list upするのにn オーダー、そこから足してnを見つけるのはどうやってもn²オーダーよりは早くできそうにない

**デフォルトの名無しさん** · 2023/06/09(金) 20:08:26.01

>>738 c
https://ideone.com/Dbl6mc
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
void sieve_of_eratosthenes(int n, int **isprime, int **primes_beg, int **primes_end) {
int i, m, j, nbytes = (1 + n) * sizeof **isprime, *map = malloc(nbytes), *beg = malloc(nbytes), *end;
for (i = 2; i <= n; i++) map[i] = 1;
for (i = 2, m = sqrt(n); i <= m; i++)
if (map[i]) for (j = i * i; j <= n; j += i) map[j] = 0;
for (i = 2, end = beg; i <= n; i++) if (map[i]) *end++ = i;
*isprime = map, *primes_beg = beg, *primes_end = end;
}
int main() {
int n = 1234567, *isprime, *primes_beg, *primes_end;
sieve_of_eratosthenes(n, &isprime, &primes_beg, &primes_end);

int *a, *b, c, count = 0;
for (a = primes_beg; a < primes_end && *a <= n / 3; a++)
for (b = a; b < primes_end && *b <= (c = n - *a - *b); b++)
if (isprime[c]) count++;
printf("count = %d\n", count);
return 0;
}

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/09(金) 20:55:26.51

お題: 与えられた4頂点が長方形または正方形を形成するかどうか判定せよ。長方形または正方形のときtrueを返すものとする。頂点の順序は無視する。

例)
{ (0, 0), (1, 0), (1, 1), (0, 1) } → true.
{ (0, -1), (0, 1), (1, 0), (-1, 0) } → true.
{ (0, -2), (1, 0), (0, 2), (-1, 0) } → false.

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/09(金) 22:16:01.89

お題: コマンドロボット。

ロボットはx-y平面の原点(0, 0)に居る。
コマンドIDはロボットによる何かのアクションを表す0以上4以下の整数で、次のように定義する。

コマンド0は何もしない。
コマンド1は位置が(-1, 0)だけ動く。
コマンド2は位置が(+1, 0)だけ動く。
コマンド3は位置が(0, -1)だけ動く。
コマンド4は位置が(0, +1)だけ動く。

入力として、コマンドIDのリストが与えられる。コマンドを順番に実行したときの最終位置を出力せよ。

例)
(2, 3, 0, 3, 1) → (0, -2).
(1, 1, 4, 0, 0, 2) → (-1, 1).

**デフォルトの名無しさん** · 2023/06/09(金) 22:46:11.50

>>743
R
www.ideone.com/hNR3VK

**デフォルトの名無しさん** · 2023/06/09(金) 22:55:53.04

>>739 >>742
正解。

出題者によるプログラム
www.ideone.com/0R9KaI

**デフォルトの名無しさん** · 2023/06/09(金) 23:13:25.78

>>744
R
www.ideone.com/eH07Wk

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/09(金) 23:20:35.60

お題: コマンドロボット2。

ロボットはx-y平面の原点(0, 0)に居る。
コマンドIDはロボットによる何かのアクションを表す0以上4以下の整数で、次のように定義する。

コマンド0は何もしない。
コマンド1は位置が(-1, 0)だけ動く。
コマンド2は位置が(+1, 0)だけ動く。
コマンド3は位置が(0, -1)だけ動く。
コマンド4は位置が(0, +1)だけ動く。

入力として、最終位置が与えられる。原点から最終位置まで移動する最短のコマンドIDリストを生成せよ。
最終位置のx成分とy成分がどちらも非ゼロならx成分から先に移動せよ。

例)
(2, -3) → (2, 2, 3, 3, 3).
(-1, 3) → (1, 4, 4, 4).

**デフォルトの名無しさん** · 2023/06/09(金) 23:59:17.57

>>748
R

FormatArray <- function(a) paste0("(", paste0(a, collapse = ", "), ")")
P <- list(c(2, -3), c(-1, 3), c(0, 0))
for (p in P) cat(FormatArray(p), " → ", FormatArray(rep((sign(p) + c(3, 7)) / 2, abs(p))), "\n", sep = "")

[実行結果]
(2, -3) → (2, 2, 3, 3, 3)
(-1, 3) → (1, 4, 4, 4)
(0, 0) → ()

**デフォルトの名無しさん** · 2023/06/10(土) 01:30:29.55

>>722
Kotlin

○付き数字はブラウザやフォントの設定によっては表示が文字ごとに大きかったり小さかったりするかも知れない。

https://paiza.io/projects/I5rPp0OBjWg6RdMlWu3t5g

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/10(土) 11:10:34.93

お題: x=0～πの範囲(step: π/8)でsin(x)を数値微分(dx=0.001)すると、その値sin'(x)がほぼcos(x)に一致することを示せ。
xの各値に対する計算結果(x, sin(x), sin'(x), cos(x))を一行に並べて表示せよ。
各数値は小数点以下3桁まで表示するものとする。

**デフォルトの名無しさん** · 2023/06/10(土) 11:11:16.11

>>748 ruby

def moves(x, y)
[].fill(x < 0 ? 1 : 2, ...x.abs) + [].fill(y < 0 ? 3 : 4, ...y.abs)
end
=> :moves
moves(2, -3)
=> [2, 2, 3, 3, 3]
moves(-1, 3)
=> [1, 4, 4, 4]
moves(0, 0)
=> []

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/10(土) 12:03:39.21

お題: 物理シミュレーション。

運動方程式によると、位置の微分が速度になり、速度の微分が加速度となっている。
逆に言えば加速度の積分と初速から現在の速度が求まり、速度の積分と初期位置から現在の位置が求まる。

数直線上を移動可能な点Pがある。Pの初期位置は原点0[cm]とする。Pの初速はゼロ[cm/s]とする。

入力として加速度[cm/s^2]と加速時間[s]のペアのリストが与えられる。時刻t=0[s]から入力に従って点Pが加速する。
dt=0.1として区分求積法により加速度と速度を数値積分し、1秒ごとのPの位置と速度を5秒後まで出力せよ。各数値は小数点以下2桁まで表示せよ。

例) ((1, 1), (-1, 1), (2, 0.5), (-4, 0.25), (3, 1.25), (-8, 1))

**デフォルトの名無しさん** · 2023/06/10(土) 15:57:09.17

>>738 c
・>>742から若干の省メモリ化
　判別用配列の要素をintからcharへ
　素数の個数を数えてから一覧用領域を確保
・組み合わせループの条件を>>746さんを参考に改善
https://ideone.com/CixK7I

>>738 c
・上記からさらに省メモリ化
　判別用配列はint配列だが判別はビットごと
https://ideone.com/JDcFjV

**754** · 2023/06/10(土) 16:17:23.88

>>754
誤： for (i = 2; i <= n; i++) if (at(map, i)) nprimes++;
正： for (i = 2, nprimes = 0; i <= n; i++) if (at(map, i)) nprimes++;

**蟻人間** ◆T6xkBnTXz7B0 · 2023/06/10(土) 20:13:01.15

>>751
C++ 回答例
https://ideone.com/QdLtnJ

**デフォルトの名無しさん** · 2023/06/10(土) 23:02:02.89

>>744
Java
https://paiza.io/projects/Y3CoVhKtAgw03ZnnwF8GEA

**デフォルトの名無しさん** · 2023/06/10(土) 23:37:30.81

>>749
Java
https://paiza.io/projects/79057xQfec28jB_BAD40oQ?language=java

**デフォルトの名無しさん** · 2023/06/11(日) 00:07:26.45

お題：数学的に植物を描画せよ

**デフォルトの名無しさん** · 2023/06/11(日) 00:37:53.77

>>753
www.ideone.com/jv00La

Rで書いたらdouble型の精度の問題で期待通りに動作しなかったので、decimal型を使えるC#で書いた。
例のリストで解析解と一致する数値解が得られるように、速度は台形則、位置は台形則と中点則の平均
（結局は刻み幅が半分での台形則に一致するが）で計算した。参考のためstdinに解析解を書いておいた。

**638** · 2023/06/11(日) 02:58:42.46

力ずく問題は苦手なインタプリタ―だが、速度に配慮した書き方をすればこの規模だと数分程度で解ける…

use feature qw{:5.16 signatures say};
$m = 1234567;
%h = map{$_ => ++$k} 2..$m;
for $i (2..int(sqrt $m)) {
　if (exists $h{$i}) {
　　for ($i..(int($m / $i)+1)) { delete $h{$i * $_} }
　}
}
@p = sort{$a <=> $b} keys %h;
%o = map{$p[$_] => $_} 0..$#p;
for $i (0..$#p-2) {
　$p1 = $p[$i];
　last if $m <= 3 * $p1;
　for $j ($i+1..$#p-1) {
　　$p2 = $p[$j]; $s = $p1 + $p2;
　　$r = $m - $s; last if $r <= $p2;
　　$n++ if exists $o{$r};
　}
}
print "n = $n\n";

実行結果 (CPUは Core-i7 8995u、Mem: 32GB)
$ time perl 21_738_prime+1234567.pl
n = 151055501

real 5m48.035s
user 5m44.468s
sys 0m3.250s

出題者の解と合わん…＞＜

**638** · 2023/06/11(日) 03:01:19.54

>>761 アンカーしくった＞＜
>>738 Perl5

**デフォルトの名無しさん** · 2023/06/11(日) 20:59:06.26

>>761
インタプリタでもRだともっと速い。以下のプログラムを https: //tio.run で実行すると
9.96秒で完了し、解は>>746の151060650通りと合う。>>761は制限時間 (1分) 以内に
完了しない。Rは自前のループをなるべく書かずにベクトル演算にすれば割と速い。

n <- 1234567L
isprime <- c(FALSE, rep(TRUE, n - 1))
for (i in 2:sqrt(n)) if (isprime[i]) isprime[seq(i + i, n, i)] <- FALSE
p <- which(isprime)
q <- cumsum(isprime)

k <- 0L
for (i in p[1:q[n %/% 3L]]) {
　　m <- n - i
　　k <- k + sum(isprime[m - p[q[i]:q[m %/% 2L]]])
}
cat(k, "通り\n", sep = "")