HDU 4191

原创 2013年12月05日 18:52:28

题意:一场网球比赛,告诉你你能赢得一分的概率,输出赢得一局,一轮,和一场的比赛的概率,规则和正式比赛一模一样。

正常情况:dp[i][j] = p*dp[i+1][j] + (1 - p) * dp[i][j+1];

对于可能一直打的情况: dp[i][i] = p * p / (1 - p * (1 - p) * 2);

三次dp 都利用上次的计算结果

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#define DB double
using namespace std;

const int MAXN = 20;
DB dp[MAXN][MAXN];
void init(){
    for(int i = 0;i <= 10; i++)
        for(int j = 0; j <= 10; j++)
            dp[i][j] = 0.0;
}
int main()
{
 //   freopen("test.txt","r",stdin);
    DB p;
    while(scanf("%lf",&p) != EOF){
        if(p == -1.0) break;
        init();
        for(int i = 0; i < 3; i++) dp[4][i] = 1.0;
        dp[3][3] = p * p / (1 - p * (1 - p) * 2);
        for(int i = 3; i >= 0; i--)
            for(int j = 3; j >= 0;j--){
                if(i == 3 && j == 3) continue;
                dp[i][j] = p * dp[i+1][j] + (1 - p) * dp[i][j + 1];
            }
        DB gp = dp[0][0];
        init();
        for(int i = 0; i < 6; i++) dp[7][i] = 1.0;
        dp[6][6] = p * p / (1 - p * (1 - p) * 2);
        for(int i = 6; i >= 0; i--)
            for(int j = 6; j >= 0;j--){
                if(i == 6 && j == 6) continue;
                dp[i][j] = p * dp[i+1][j] + (1 - p) * dp[i][j + 1];
            }
        DB p7 = dp[0][0];
    //    cout<<"p7:"<<p7<<endl;
        init();
        dp[6][6] = p7;
        dp[6][5] = gp + (1.0 - gp) * dp[6][6];
        dp[5][6] = gp * dp[6][6];
        for(int i = 0; i <= 4;i++) dp[6][i] = 1.0;
        for(int i = 5; i >= 0; i--)
            for(int j = 5;j >= 0;j--){
         //       DB p1,p2;
                dp[i][j] = gp * dp[i+1][j] + (1.0 - gp) * dp[i][j+1];
            //    if(i + 1 >= j + 2 && i + 1 >= 6) p1 = 1.0;
           //     else if(i + 3 <= j && j <= 6) p1 = 0.0;
          //      else p1 = dp[i+1][j];

           //     if(i  >= j + 3 && i >= 6) p2 = 1.0;
           //     else if(i + 2 <= j + 1 && j + 1 >= 6) p2 = 0.0;
           //     else p2 = dp[i][j+1];
           //     dp[i][j] = gp * p1 + (1.0 - gp) * p2;
            }
        DB sp = dp[0][0];
        init();
        for(int i = 0; i < 2; i++) dp[2][i] = 1.0;
        for(int i = 1; i >= 0; i--)
            for(int j = 1; j >= 0; j--){
                dp[i][j] = dp[i+1][j] * sp + dp[i][j+1] * (1.0 - sp);
            }
        DB mp = dp[0][0];
        printf("%.11f %.11f %.11f\n",gp,sp,mp);
    }
    return 0;
}


相关文章推荐

spoj 4191. Sky Code(容斥)

http://www.spoj.com/problems/MSKYCODE/ 给出一个集合,含有n个元素,每次任意从中取出4个使得他们的gcd是1,问有多少种取法。 可以先考虑...

boost thread warning C4191和algorithm warning C4242

今天抽出时间来整理一下项目中的警告信息,发现有几个警告很神奇,经过一番查证,终于找到了问题所在。第一个警告:boost/thread/win32/thread_primitives.hpp(315):...
  • natahew
  • natahew
  • 2017年06月16日 09:47
  • 304

hdu1002大数相加

  • 2014年12月16日 12:44
  • 877B
  • 下载

acm课件 HDU 算法大全

  • 2013年08月08日 20:57
  • 6.87MB
  • 下载

HDU:1257 最少拦截系统(贪心+动态规划DP||LIS?)

最少拦截系统 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total S...
  • zugofn
  • zugofn
  • 2016年08月11日 09:49
  • 703

HDU2501解题报告

  • 2014年07月10日 10:44
  • 12KB
  • 下载

hdu 3341(ac自动机+状态压缩)

  • 2013年10月29日 17:03
  • 23KB
  • 下载

hdu 4777 树状数组+合数分解

Rabbit Kingdom Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...

hdu+acm课件

  • 2014年06月26日 20:33
  • 8.66MB
  • 下载

HDU1019(2028)解题报告

  • 2014年07月11日 10:05
  • 555B
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 4191
举报原因:
原因补充:

(最多只允许输入30个字)