12月14日

这几天继续研究概率dp，先做的是Collecting Bugs，这道题一道入门题，题意是一个软件有s个子系统，会产生n种bug，一天发现一个bug,这个bug属于一个子系统，属于一个分类    每个bug属于某个子系统的概率是1/s,属于某种分类的概率是1/n    问发现n种bug,每个子系统都发现bug的天数的期望。用dp[i][j]表示已经找到了i中bug，j个系统的bug，达到目标天数的期望值，从dp[n][s]=0开始推，最推出来出来dp[0][0]，一共有4种情况，加起来就是 dp[i][j]=(i*(s-j)*dp[i][j+1]+(n-i)*j*dp[i+1][j]+(n-i)*(s-j)*dp[i+1][j+1]+n*s)/(n*s-i*j)。然后做One Person Game，题意是有三个骰子，分别有k1,k2,k3个面。每次掷骰子，如果三个面分别为a,b,c则分数置0，否则加上三个骰子的分数之和。当分数大于n时结束。求游戏的期望步数。初始分为0,。这个题的公式真心不好推，需要列公式然后解方程组，用dp[i]表示分数达到i时的期望步数，当i>=n时，dp[i]=0，要出来dp[0]。dp[i]=和号(dp[i+k]*Pk)+E[0]*P0+1。其中Pk表示三个骰子和为k的概率，其中P0就是回到0个概率。 然后解方程组求dp[0]。设dp[i]=A[i]*dp[0]+B[i]，将 dp[i+k]=A[i+k]*dp[0]+B[i+k]代入上述方程右边得到：dp[i]=和号(pk*A[i+k]*dp[0]+pk*B[i+k])+dp[0]*p0+1=(和号(pk*A[i+k])+p0)dp[0]+和号(pk*B[i+k])+1；有A[i]=(∑(pk*A[i+k])+p0)，B[i]=和号(pk*B[i+k])+1，当i>n时dp[i]=A[i]*dp[0]+B[i]=0，故A[i>n]=B[i>n]=0。 所以可以先递推求得A[0]和B[0]。那么 dp[0]=B[0]/(1-A[0])。然后做Aeroplane chess这道题，题意是有n个格子，掷色子的掷出的数目就是你一次到移动格数。其中有m个飞行通道可以让你直接从第xi格飞到第yi格。问你走到终点的期望是多少。用dp[i]表示当前在第i格走到终点的期望。dp[i]=1+ (dp[i+1]+dp[i+2]+dp[i+3]+dp[i+4]+dp[i+5]+dp[i+6])/6。当i是飞行通道的起点的时候特判一下令dp[xi]=dp[yi]。最后做的是Activation这道题，题意是有n个人排队等着在官网上激活游戏。Tomato排在第m个。队首位置有四种情况1，注册失败，不影响队列顺序 ，概率为p1。2，连接失败，队首的人排到队尾，概率为p2。3，注册成功，队首离开队列，概率为p3。4，服务器崩溃，激活停止，概率为p4。求Tomato的位置在K以内，而且服务器崩溃的概率。用dp[i][j]表示队列总共有i个人，站在第j个位置时，所求事件发生的概率。这样就有转移方程，dp[i][1]=dp[i][1]*p1+dp[i][i]*p2+p4，dp[i][j]=dp[i][j]*p1+dp[i][j-1]*p2+dp[i-1][j-1]*p3+p4(j<=k)，dp[i][j]=dp[i][j]*p1+dp[i][j-1]*p2+dp[i-1][j-1]*p3(j>k)。整理一下令  p21=p2/(1-p2)  p31=p3/(1-p1)  p41=p4/(1-p1)，得dp[i][1]=dp[i][i]*p21+p41，dp[i][j]=dp[i][j-1]*p21+dp[i-1][j-1]*p31+p41(j<=k)，dp[i][j]=dp[i][j-1]*p21+dp[i-1][j-1]*p31(j>k)。dp[i-1][j-1]在递推的时候都可以解决，我们令后面的部分为c[j]，dp[i][j]就可以由dp[i][j-1]推出，但是dp[i][1]却是是由dp[i][i]推出来的，这里要迭代一下，先解出dp[i][i]，进而可以得到dp[i][1]，剩下的递推出来了 dp[i][i]=dp[i][i]*p21^i+c[2]^p21^(i-1)……c[i]+p41*p21。最后得到dp[n][m]就是答案。还需要特判一下崩溃的概率是不是小于等于0，如果小于等于便是输出0,。