高校数学の質問スレ Part437

[過去ﾛｸﾞ] 高校数学の質問スレ Part437 (1002ﾚｽ)
上下前次1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞
抽出解除必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞あぼーん

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索歴削→次ｽﾚ栞削→次ｽﾚ過去ﾛｸﾞﾒﾆｭｰ

177: １３２人目の素数さん [sage] 2024/07/21(日) 21:34:30.71 ID:3skoHSGy

>>166（Wolfram言語での想定解）
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= calc[n_,m_] :=(
If[m<1||m>4,Return[0]];
f[x_] :=(
{p1,p2,p3,p4}={1/10,2/10,3/10,4/10};
{b1,b2,b3,b4}=x;
Binomial[n,b1] p1^b1 Binomial[n-b1,b2] p2^b2 Binomial[n-b1-b2,b3] p3^b3 p4^b4
);
x=Select[Partition[Flatten@Table[{x1,x2,x3,n-x1-x2-x3},{x1,Range[0,n]},{x2,Range[0,n]},{x3,Range[0,n]}],4],#[[4]]>=0&];
Total[f /@ Select[x,Total@Boole@Table[#[[i]]>0,{i,Range[4]}]==m&]]
)
In[2]:= Table[calc[2,i],{i,1,4}]
3   7
Out[2]= {--, --, 0, 0}
10  10
In[3]:= Table[calc[3,i],{i,1,4}]
1   3  3
Out[3]= {--, -, --, 0}
10  5  10
In[4]:= Table[calc[4,i],{i,1,4}]
177   2111  303  36
Out[4]= {----, ----, ---, ---}
5000  5000  625  625
In[5]:= Table[calc[5,i],{i,1,4}]
13   141  561   18
Out[5]= {----, ---, ----, ---}
1000  500  1000  125
In[6]:= Table[calc[6,i],{i,1,4}]
489    18571   2859  297
Out[6]= {------, ------, ----, ----}
100000  100000  5000  1250
In[7]:= Table[calc[7,i],{i,1,4}]
187    6111   54831   819
Out[7]= {------, -----, ------, ----}
100000  50000  100000  2500
In[8]:= Table[calc[8,i],{i,1,4}]
36177    4040351   1590309  319977
Out[8]= {--------, --------, -------, ------}
50000000  50000000  3125000  781250
In[9]:= Table[calc[9,i],{i,1,4}]
14117    1344309   23183889  753543
Out[9]= {--------, --------, --------, -------}
50000000  25000000  50000000  1562500
In[10]:= Table[calc[10,i],{i,1,4}]
22173     7205087   4179633   1364751
Out[10]= {---------, ---------, --------, -------}
200000000  200000000  10000000  2500000

小数表示
{0.3, 0.7, 0., 0.}

{0.1, 0.6, 0.3, 0.}

{0.0354, 0.4222, 0.4848, 0.0576}

{0.013, 0.282, 0.561, 0.144}

{0.00489, 0.18571, 0.5718, 0.2376}

{0.00187, 0.12222, 0.54831, 0.3276}

{0.00072354, 0.080807, 0.508899, 0.409571}

{0.00028234, 0.0537724, 0.463678, 0.482268}

{0.000110865, 0.0360254, 0.417963, 0.5459}

0.5を超えるのは10人以上。

http://rio2016.5ch.net/test/read.cgi/math/1721071007/177

182: １３２人目の素数さん [sage] 2024/07/21(日) 23:23:25.56 ID:3skoHSGy

>>175
可読性を度外視して1行に纏めた。自分でも解読するのに難渋しそう。

*  血液型頻度比はA:O:B:AB＝4:3:2:1,   n 人に m 種類の血液型が存在する確率  *)
solve[n_,m_] := Total[Binomial[n,#[[1]]] (1/10)^#[[1]] Binomial[n-#[[1]],#[[2]]] (2/10)^#[[2]] Binomial[n-#[[1]]-#[[2]],#[[3]]] (3/10)^#[[3]] (4/10)^#[[4]]& /@ (Select[Select[Partition[Flatten@Table[{x1,x2,x3,n-x1-x2-x3},{x1,Range[0,n]},{x2,Range[0,n]},{x3,Range[0,n]}],4],#[[4]]>=0&],Total@Boole@Table[#[[i]]>0,{i,Range[4]}] == m &])]

実行結果

In[1]:= solve[n_,m_] := Total[Binomial[n,#[[1]]] (1/10)^#[[1]] Binomial[n-#[[1]],#[[2]]] (2/10)^#[[2]] Binomial[n-#[[1]]-#[[2]],#[[3]]] (3/10)^#[[3]] (4/10)^#[[4]]& /@ (Select[Select[Partition[Flatten@Table[{x1,x2,x3,n-x1-x2-x3},{x1,Range[0,n]},{x2,Range[0,n]},{x3,Range[0,n]}],4],#[[4]]>=0&],Total@Boole@Table[#[[i]]>0,{i,Range[4]}] == m &])]

In[2]:= solve[40,3]

4660364912515525380760230058213856169
Out[2]= ---------------------------------------
312500000000000000000000000000000000000

尿瓶チンパフェチの嫁なしPhimoseくん（別名、自演認定厨）はWolfram言語はわからないみたいだな。
初学の俺と達人のコードが区別できないみたいだし。
嫁なしPhimoseくんが自演認定厨であることは度々、識者より指摘されている。
俺の想定解よりエレガントな解が投稿されていることを理解できないようだ。

http://rio2016.5ch.net/test/read.cgi/math/1721071007/182

上下前次1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ

ぬこの手ぬこTOP 1.712s*