面白い問題おしえて～な 28問目

面白い問題おしえて～ ..

57:１３２人目の素数さん
18/11/03 17:47:45.86 OVkXWZOI.net
>>54
なるほど！そうすればできるね。
,var n
R.<x> = CC['x']
g = (x/6)/(1-(x^2+x^3+x^4+x^5+x^6)/6)/(1-x) #x^nの係数は総和がn以下になる確率
pfd = g.partial_fraction_decomposition()
p = sum(c.numerator()/c.denominator().subs(x=0)*(-1/c.denominator().subs(x=0))^n for c in pfd[1]) #pは総和がn以下になる確率
sage: p
-0.924602908258674*e^(-0.0450727249412852*n)
- (0.0141317451913899 + 0.0178673357104935*I)*e^(-(0.329151925064039 - 1.18446726809051*I)*n)
- (0.0141317451913899 - 0.0178673357104935*I)*e^(-(0.329151925064039 + 1.18446726809051*I)*n)
- 0.0172163389039878*e^(-(0.361410021965218 - 3.14159265358979*I)*n)
- (0.0149586312272785 + 0.00913725707805078*I)*e^(-(0.363486436096737 - 2.18509474953750*I)*n)
- (0.0149586312272785 - 0.00913725707805078*I)*e^(-(0.363486436096737 + 2.18509474953750*I)*n)
+ 1.00000000000000
sage: 1-p.subs(n=49)
0.101575341438976 + 5.79026428787119e-24*I

次ページ

続きを表示

1を表示