| For All Solvers
OMC114

OMC114(E)

ユーザー解説 by tria

 うまい方法が思い付かなくても計算力次第で何とかなることもあるよということで,計算力を使って解く方法を紹介します.


 まず,(nk)\dbinom{n}{k} で,0kn0\leq k\leq n のときは通常の二項係数,k<0k\lt0k>nk\gt n のときは 00 を表すとします.nn が負の場合はここでは扱いません.
 多項式 f(x)f(x)xnx^n の係数を [xn]f(x)[x^n]f(x) で表します(多項式でなくてもいいのですがここでは多項式としておきます).例えば,f(x)=5x214x+3f(x)=5x^2-14x+3 のときは, [x0]f(x)=3, [x1]f(x)=14, [x3]f(x)=0[x^0]f(x)=3,\ [x^1]f(x)=-14,\ [x^3]f(x)=0 などとなります.
 この記法を用いると,多項式 f(x),g(x)f(x),g(x) に対し, k=0n[xk]f(x)[xnk]g(x)=[xn](f(x)g(x))\displaystyle\sum_{k=0}^n[x^k]f(x)[x^{n-k}]g(x)=[x^n] (f(x)g(x)) となります(証明は割愛するので各自確かめてください).
 また,[xk](1+x)p=(pk)[x^k] (1+x)^p=\dbinom{p}{k} であることを使うと
k=0n(pk)(qnk)=k=0n[xk](1+x)p[xnk](1+x)q=[xn+k](1+x)p+q=(p+qn)\begin{aligned} \sum_{k=0}^n\binom{p}{k}\binom{q}{n-k} &=\sum_{k=0}^n[x^k] (1+x)^p[x^{n-k}] (1+x)^q\\ &=[x^{n+k}] (1+x)^{p+q} =\binom{p+q}{n} \end{aligned} となります(k>pk\gt p でも [xk](1+x)p=(pk)[x^k] (1+x)^p=\dbinom{p}{k} は成り立っていることに注意!).
 これは Vandermondeの恒等式 と呼ばれており,あとで使います.


 問題の解説に入ります.解説の都合上,棒を並び替えるという状況を,各マスに整数を書き込むという状況に置き換えます.スコアの期待値は,各マスでの LlRrLlRr の期待値の合計になります.左から k+1k+1 番目のマスに m+1m+1 が書き込まれているとして,他の 9999 マスへの整数の書き込み方 99!99! 通りについての LlRrLlRr の期待値を E(k,m)E(k,m) とします.
 11 から 100100 までに m+1m+1 より小さい整数は mm 個あり,大きい整数は 99m99-m 個あります.また, m+1m+1 の左には kk マスあり,右には 99k99-k マスあります.
 m+1m+1 より小さい mm 個の整数のうち ll 個が m+1m+1 より左,mlm-l 個が右にあるとします.このとき, LlRrLlRr の値は l(kl)(ml)(99km+l)l(k-l)(m-l)(99-k-m+l) となり,このような書き込み方は (ml)(99mkl)k!(99k)!\dbinom{m}{l}\dbinom{99- m}{k-l}k!(99-k)! 通りあります(l>kl\gt k のときは LlRrLlRr の値が負になりますが,二項係数の部分が 00 になるので帳尻は合っています).

 これにより, E(k,m)=199!l=0m(ml)(99mkl)k!(99k)! l(kl)(ml)(99km+l)=k!(99k)!99!l=0mm!l!(ml)!(99m)!(kl)!(99km+l)!l(kl)(ml)(99km+l)=m(m1)(99m)(98m)k!(99k)!99!l=1m1(m2)!(l1)!(ml1)!(97m)!(kl1)!(98km+l)!=m(m1)(99m)(98m)k!(99k)!99!l=1m1(m2l1)(97mkl1)=m(m1)(99m)(98m)k!(99k)!99!(95k2)=m(m1)(99m)(98m)k(k1)(99k)(98k)99989796\begin{aligned} \displaystyle E(k,m) &=\frac{1}{99!}\sum_{l=0}^m\dbinom{m}{l}\dbinom{99- m}{k-l}k!(99-k)!\ l(k-l)(m-l)(99-k-m+l) \\ &=\frac{k!(99-k)!}{99!}\sum_{l=0}^m\frac{m!}{l!(m-l)!}\frac{(99-m)!}{(k-l)!(99-k-m+l)!}l(k-l)(m-l)(99-k-m+l) \\ &=m(m-1)(99-m)(98-m)\frac{k!(99-k)!}{99!}\sum_{l=1}^{m-1}\frac{(m-2)!}{(l-1)!(m-l-1)!}\frac{(97-m)!}{(k-l-1)!(98-k-m+l)!} \\ &=m(m-1)(99-m)(98-m)\frac{k!(99-k)!}{99!}\sum_{l=1}^{m-1}\binom{m-2}{l-1}\binom{97-m}{k-l-1} \\ &=m(m-1)(99-m)(98-m)\frac{k!(99-k)!}{99!}\binom{95}{k-2} \\ &=\frac{m(m-1)(99-m)(98-m)k(k-1)(99-k)(98-k)}{99\cdot98\cdot97\cdot96} \end{aligned} と計算できます.
 求める期待値 EEE=k=0991100m=099E(k,m)=110099989796(k=099k(k1)(99k)(98k))2\begin{aligned} E &=\sum_{k=0}^{99}\frac{1}{100}\sum_{m=0}^{99}E(k,m)\\ &=\frac{1}{100\cdot99\cdot98\cdot97\cdot96}\left(\sum_{k=0}^{99}k(k-1)(99-k)(98-k)\right)^2 \end{aligned} となります.

 次に,Sn=k=0nk(k1)(nk)(nk1)S_n=\displaystyle\sum_{k=0}^nk(k-1)(n-k)(n-k-1) を求めましょう.SnS_n を計算してみると, S0=S1=S2=S3=0,S4=4,S5=24,S6=84,S7=224,S8=504,S_0=S_1=S_2=S_3=0,\quad S_4=4,\quad S_5=24,\quad S_6=84,\quad S_7=224,\quad S_8=504,\ldots となるので, Sn=(n+1)n(n1)(n2)(n3)30S_n=\dfrac{(n+1)n(n-1)(n-2)(n-3)}{30} と予想できます(予想のコツは階差ではなく比を見ることです).
 ここで,SnS_n は高々 55 次の多項式であることがわかっているので,この予想は正しいです.
 また,次のように求めることもできます.x<1|x|\lt1 のとき 11x=n=0xn\displaystyle\frac{1}{1-x}=\sum_{n=0}^\infty x^n であることを思い出しましょう.すると,[xn]11x=1\displaystyle[x^n]\frac{1}{1-x}=1 と書けます([xn]f(x)[x^n]f(x) の説明で多項式でなくてもよいと言っていましたが,これのことです.収束半径内であれば項の順序を変えてもよいので,畳み込みの式も成り立っています).微分することで [xn]2(1x)3=(n+1)(n+2),[xn]120(1x)6=(n+1)(n+2)(n+3)(n+4)(n+5)\displaystyle[x^n]\frac{2}{(1-x)^3}=(n+1)(n+2),\quad [x^n]\frac{120}{(1-x)^6}=(n+1)(n+2)(n+3)(n+4)(n+5) を得ます(級数と微分は入れ替えてよいときと悪いときがありますが,今回はよいときです).
 これを使うと, Sn=k=0n[xk2]2(1x)3[xnk2]2(1x)3=[xn4]4(1x)6=(n+1)n(n1)(n2)(n3)30\begin{aligned} S_n &=\sum_{k=0}^n[x^{k-2}]\frac{2}{(1-x)^3}[x^{n-k-2}]\frac{2}{(1-x)^3}\\ &=[x^{n-4}]\frac{4}{(1-x)^6}=\dfrac{(n+1)n(n-1)(n-2)(n-3)}{30} \end{aligned} と計算できます.

 最後に,Sn=(n+1)n(n1)(n2)(n3)30S_n=\dfrac{(n+1)n(n-1)(n-2)(n-3)}{30} を代入すると, E=S99210099989796=10099989796900=10038336\begin{aligned} E&=\frac{S_{99}^2}{100\cdot99\cdot98\cdot97\cdot96} &=\frac{100\cdot99\cdot98\cdot97\cdot96}{900}=\mathbf{10038336} \end{aligned} となります.


 解説なので長くなりましたが,これぐらいの計算なら慣れれば10分ほどで終わります.