ACM训练日记—1月4日

        新的一年的第一篇半周总结。昨天刚结束了常微分考试，考的我拔凉拔凉的，常微分老师对我们这一级真的是关照，相比15级的常微分期末题，难度提升好多。。QAQ，还是好好准备物理，英语吧。

       虽然这半周考试有点紧张，但还是有有点收获的，其实看了一些博客，关于概率的主要知识点主要就是求期望和概率dp，其他还有结合组合数学的。整理一下水的几道相关题吧。

      LightOJ 1038  Race to 1 Again

      题意：给出一个数n,我们可以选择 1－n中可以被n整除的数字，然后用n出得到一个新的数字n1；然后在找所有n１的因子，用n１除，直到得到１；问除的次数的期望值；

     一道比较明显的概率dp，n一定与所有因子期望有关，设n=(p1^a1)...(pn^kn)。dp[n]=((dp[1]+1)/k(dp[p1]+1)/k)+....+(dp[pk]+1)/k+(dp[n]+1)/k，（k是因子个数，包括本身与1），整理下dp[n]=(dp[1]+dp[p1]+...dp[pk]+k)/k。然后循环预处理就可以了，注意时间限制，要做一个sqrt()的处理。

      LightOJ 1274 Beating the Dataset

      题意：给出n表示有n个单词（单词只有两个yes与no），给出s表示有s个字母（yes又3个字母，no有2个字母），下面用1表示yes，用0表示no，问给出n和s，求所有01串将末尾去掉，头上补上1形成的新串与原来的串有多少个不一样，求所有组合的平均值。

     这个题竟然有两种解法，可以分别用概率dp和组合数学做。 

     1，组合数学方法

来自这位大佬的证明http://blog.csdn.net/doris1104/article/details/52950700

设1有x个，0有y个，那么解二元方程组，x=s-2*n，y=3*n-s。

那么最高位的时候只有为0的才会不一样: C(x,x+y-1)（补充：第一个位置放0，剩下x+y-1个放x个1种情况）
其他位置只要跟前面一个位置不一样就可以了，那么就是说当前和前面一个位置一个0一个1 ， 共C(x-1 , x+y-2)，因为1，0可以交换所以要乘2（补充：首端放1，剩下x+y-1个位置相当于用x+y-2个位置放n-1个1，然后用一个0进行插空，以此形成变化，所以是(x+y-1)*C(x-1,x+y-2)）。

平均值就是(2*(x+y-1)*C(x-1,x+y-2)+C(x,x+y-1)) / C(x,x+y)

(2*x*C(x,x+y-1)+C(x,x+y-1)) / C(x,x+y)
=(2*x+1)*C(x,x+y-1) /(C(x,x+y-1)+C(x-1,x+y-1))
=(2*x+1)*C(x,x+y-1) /((1+x/y)*C(x,x+y-1))
=(2*x*y+y)/(x+y)这就是最终公式了。

      2，概率dp   来自这位大佬的讲解  http://blog.csdn.net/qq_31759205/article/details/54646913

题解：期望dp，逆推。dp[i][j][k]表示当前状态为第i行，已经输出了j个YES，这一行要输出YES(0)/NO(1)
因此可以得出，①如果第i要输出YES，那么i+1行如果输出NO，则i的期望值为i+1的期望值加上+1再乘以输出NO的概率；i+1行如果输出YES，则i的期望值为i+1的期望值乘以输出YES的概率。②如果第i行要输出NO，做法和①相同
递推关系式为
dp[i][j][0]=dp[i+1][j+1][0]*p1+(dp[i+1][j][1]+1)*p2
dp[i][j][1]=(dp[i+1][j+1][0]+1*p1)*p1+dp[i+1][j][1]*p2

       LightOJ - 1079 Just another Robbery

       题意：抢银行，现在给出一个概率p,和银行的个数n，接下去给出每个银行可以抢到的钱,还有抢劫这个银行被抓的概率，问在被抓的概率小于等于p的情况下,最多抢到多少钱。

       这其实是一道做过的题， 其实就是01背包，dp[j] 代表抢劫到 j 元钱,逃跑的概率;那么dp[j] = max( ( dp[j - m[i]] ) * p[i])。

p[i]=1-p[i];
memset(dp,0,sizeof(dp));
dp[0]=1;
for(i=0;i<n;i++)
{
      for(j=sum;j>=m[i];j--)
      {
          dp[j]=max(dp[j],dp[j-m[i]]*p[i]);
      }
}

先整理这些吧，新年快乐！