zoj3329 One Person Game
来源:互联网 发布:80端口提供什么服务 编辑:程序博客网 时间:2024/06/07 08:38
There is a very simple and interesting one-person game. You have 3 dice, namelyDie1, Die2 and Die3. Die1 hasK1 faces. Die2 has K2 faces.Die3 has K3 faces. All the dice are fair dice, so the probability of rolling each value, 1 toK1, K2, K3 is exactly 1 /K1, 1 / K2 and 1 / K3. You have a counter, and the game is played as follow:
- Set the counter to 0 at first.
- Roll the 3 dice simultaneously. If the up-facing number of Die1 isa, the up-facing number of Die2 is b and the up-facing number ofDie3 is c, set the counter to 0. Otherwise, add the counter by the total value of the 3 up-facing numbers.
- If the counter's number is still not greater than n, go to step 2. Otherwise the game is ended.
Calculate the expectation of the number of times that you cast dice before the end of the game.
Input
There are multiple test cases. The first line of input is an integer T (0 <T <= 300) indicating the number of test cases. Then T test cases follow. Each test case is a line contains 7 non-negative integersn, K1, K2, K3,a, b, c (0 <= n <= 500, 1 < K1,K2, K3 <= 6, 1 <= a <= K1, 1 <=b <= K2, 1 <= c <= K3).
Output
For each test case, output the answer in a single line. A relative error of 1e-8 will be accepted.
Sample Input
20 2 2 2 1 1 10 6 6 6 1 1 1
Sample Output
1.1428571428571431.004651162790698
我们可以推出公式dp[i]=sum{pk*dp[i+k]}+p0*dp[0]+1,dp[i]表示当前分数是i是时的到游戏结束的期望,那么dp[0],就是要求的,其实,我们可以发现一个归律,求期望都是从后住前推的,比如,这里,是由i+k推到i的,求概率的时候是刚好相反的,但我们发现求出的这个式子,是个环形的,所以要变形,因为,每一个都是和dp[0]相关的,我们可以设,dp[i]=a[i]*dp[0]+b[i],那么dp[0]不就是,b[i]/(1-a[i])了么,我们,把这个式子代入上式就可以得到,dp[i]=sum{pk*a[i+k]}+p0,b[i]=sum{pk*b[i+k]}+1,这样,就可以马上求结果来了!
#include <iostream>#include <stdio.h>#include <string.h>using namespace std;#define MAXN 550double pa[MAXN],pb[MAXN],dp[MAXN],p[MAXN];int main(){ int tcase,i,j,k,n,a,b,c,k1,k2,k3; double p0; scanf("%d",&tcase); while(tcase--) { scanf("%d%d%d%d%d%d%d",&n,&k1,&k2,&k3,&a,&b,&c); p0=1.0/k1/k2/k3; memset(pa,0,sizeof(pa)); memset(pb,0,sizeof(pb)); memset(dp,0,sizeof(dp)); memset(p,0,sizeof(p)); for(i=1;i<=k1;i++) for(j=1;j<=k2;j++) for(k=1;k<=k3;k++) if(i!=a||j!=b||k!=c)//不同时相等 p[i+j+k]+=p0; int temp=k1+k2+k3; for(i=n;i>=0;i--) { pa[i]=p0; pb[i]=1; for(k=1;k<=temp;k++) { pa[i]+=pa[i+k]*p[k]; pb[i]+=pb[i+k]*p[k]; } } printf("%.15f\n",pb[0]/(1.0-pa[0])); } return 0;}
- zoj3329 One Person Game
- ZOJ3329 One Person Game
- zoj3329---One Person Game
- ZOJ3329-One Person Game
- One Person Game(ZOJ3329)
- 【概率dp】zoj3329 One Person Game
- zoj3329 One Person Game(期望dp)
- ZOJ3329 One Person Game 概率dp
- [ZOJ3329]One Person Game(概率dp)
- zoj3329——One Person Game(概率DP)
- ZOJ3329-One Person Game(概率DP求数学期望)
- zoj3329 One Person Game 概率dp求期望
- zoj3329--One Person Game(概率dp第六弹:形成环的dp,带入系数,高斯消元)
- ZOJ One Person Game
- ZOJ_3329 One Person Game
- One Person Game ZOJ
- One Person Game ZOJ
- zoj 3329 One Person Game
- JavaScript面向对象继承方法
- Qt QSortFilterProxyModel示例代码, 使用方法
- 移动支付成诈骗新宠 汇付宝提醒手机购物需谨慎
- 文档 用 微软雅黑
- 【xinfanqie】快速选择有线网络 提高网络速度
- zoj3329 One Person Game
- 三、对象到底如何创建
- Hadoop 新 MapReduce 框架 Yarn 详解
- 25个增强iOS应用程序性能的提示和技巧(初级篇)
- Android的monkey测试
- ARC属性的定义
- 更换svn账户
- SQL SERVER 函数组合实现oracle的LPAD函数功能
- HDU_2032 杨辉三角