【题目】

把N个骰子扔在地上，所有骰子朝上一面的点数之和为S。输入N，打印出S的所有可能的值出现的概率。

【分析】

典型的动态规划题目。

设n个骰子的和为s出现的次数记为f(n,s)，其中n=[1-N],s=[n-6n]。

n=1, s=[1-6], f(n,s)=1;

n=[2-N], s=[n-6n], f(n,s)= f(n-1,s-1)+ f(n-1,s-2)+ f(n-1,s-3)+ f(n-1,s-4)+ f(n-1,s-5)+ f(n-1,s-6) = sum(f(n-1,s-t)) ，其中t=[1-6] and t<s。

最终点数之和为S出现的次数为f(N,S)，概率为f(N,S)/6N。

【代码】

 C++ Code 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
  
// 43_PrintSumProbabilityOfDices.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <iostream>
#include <cmath>
using namespace std;

/*
DP
f(n,s): the count of sum s of n dices, where n=[1-N],s=[n-6n]
n=1, s=[1-6], f(n,s)=1;
n=[2-N], s=[n-6n], f(n,s)= f(n-1,s-1)+ f(n-1,s-2)+ f(n-1,s-3)+ f(n-1,s-4)+ f(n-1,s-5)+ f(n-1,s-6) = sum(f(n-1,s-t)), where t=[1-6]&&t<s。
the final result is: f(N,S)/6^N
*/
 // N>14 f will overflow

int f[MAXN][MAXS];

void PrintSumProbabilityOfDices(int N)
{
     || N > MAXN)
        return;
    // init f array
; i < MAXN; ++i)
        ; j < MAXS; ++j)
            f[i][j] = ;

int n, s, t;
    // base case n=1
;
    for (s = n; s <= FACES * n; ++s)
        f[n][s] = ;

// n>=2
; n <= N; n++)
    {
        for (s = n; s <= FACES * n; ++s)
        {
            ; t <= FACES && t < s; ++t)
            {
                f[n][s] += f[n - ][s - t];
            }
        }
    }

// print the probability or counts
    n = N;
    double total = pow((double)FACES, n);
    for (s = n; s <= FACES * n; ++s)
    {
        cout << s << "--->" << f[n][s] << endl;
        //cout<<f[n][s]/total<<endl;
    }
}

void test_base(int n)
{
    cout << n << " dices." << endl;
    PrintSumProbabilityOfDices(n);
    cout << "----------------------\n";
}

void test_main()
{
    ; n <= MAXN; n++)
        test_base(n);
}

int _tmain(int argc, _TCHAR *argv[])
{
    test_main();
    ;
}

【运行结果】

win7+vs2010运行结果：

只贴出了n=1到8的结果。

1 dices.
1--->1
2--->1
3--->1
4--->1
5--->1
6--->1
----------------------
2 dices.
2--->1
3--->2
4--->3
5--->4
6--->5
7--->6
8--->5
9--->4
10--->3
11--->2
12--->1
----------------------
3 dices.
3--->1
4--->3
5--->6
6--->10
7--->15
8--->21
9--->25
10--->27
11--->27
12--->25
13--->21
14--->15
15--->10
16--->6
17--->3
18--->1
----------------------
4 dices.
4--->1
5--->4
6--->10
7--->20
8--->35
9--->56
10--->80
11--->104
12--->125
13--->140
14--->146
15--->140
16--->125
17--->104
18--->80
19--->56
20--->35
21--->20
22--->10
23--->4
24--->1
----------------------
5 dices.
5--->1
6--->5
7--->15
8--->35
9--->70
10--->126
11--->205
12--->305
13--->420
14--->540
15--->651
16--->735
17--->780
18--->780
19--->735
20--->651
21--->540
22--->420
23--->305
24--->205
25--->126
26--->70
27--->35
28--->15
29--->5
30--->1
----------------------
6 dices.
6--->1
7--->6
8--->21
9--->56
10--->126
11--->252
12--->456
13--->756
14--->1161
15--->1666
16--->2247
17--->2856
18--->3431
19--->3906
20--->4221
21--->4332
22--->4221
23--->3906
24--->3431
25--->2856
26--->2247
27--->1666
28--->1161
29--->756
30--->456
31--->252
32--->126
33--->56
34--->21
35--->6
36--->1
----------------------
7 dices.
7--->1
8--->7
9--->28
10--->84
11--->210
12--->462
13--->917
14--->1667
15--->2807
16--->4417
17--->6538
18--->9142
19--->12117
20--->15267
21--->18327
22--->20993
23--->22967
24--->24017
25--->24017
26--->22967
27--->20993
28--->18327
29--->15267
30--->12117
31--->9142
32--->6538
33--->4417
34--->2807
35--->1667
36--->917
37--->462
38--->210
39--->84
40--->28
41--->7
42--->1
----------------------
8 dices.
8--->1
9--->8
10--->36
11--->120
12--->330
13--->792
14--->1708
15--->3368
16--->6147
17--->10480
18--->16808
19--->25488
20--->36688
21--->50288
22--->65808
23--->82384
24--->98813
25--->113688
26--->125588
27--->133288
28--->135954
29--->133288
30--->125588
31--->113688
32--->98813
33--->82384
34--->65808
35--->50288
36--->36688
37--->25488
38--->16808
39--->10480
40--->6147
41--->3368
42--->1708
43--->792
44--->330
45--->120
46--->36
47--->8
48--->1
----------------------

注意，当N>14时，f数组中的若干值会溢出。

【参考】

http://zhedahht.blog.163.com/blog/static/254111742009101524946359/

http://blog.csdn.net/wuzhekai1985/article/details/6640449

43. 动态规划求解n个骰子的点数和出现概率(或次数)[Print sum S probability of N dices]的更多相关文章

  1. 【剑指offer】面试题43：n个骰子的点数

    第一种思路是,每一个骰子的点数从最小到最大,如果为1-6,那么全部的骰子从最小1開始,我们如果一种从左向右的排列,右边的最低,索引从最低開始,推断和的情况. def setTo1(dices, sta ...

  2. 《剑指offer》 面试题43 n个骰子的点数 （java）

    引言:写这篇文章的初衷只是想做个笔记,因为这道题代码量有点大,有点抽象,而书上并没有详细的注释.为了加深印象和便于下次复习,做个记录. 原题:把n个骰子扔到地上,所有骰子朝上一面的点数之后为s. 输入 ...

  3. 43：n个骰子的点数

    算法分析: 算法1.基于递归求色子点数,时间效率不高 现在我们考虑如何统计每一个点数出现的次数.要向求出n个骰子的点数和,可以先把n个骰子分为两堆:第一堆只有一个,另一个有n-1个.单独的那一个有可能 ...

  4. （剑指Offer）面试题43：n个骰子的点数

    题目: 把n个骰子仍在地上,所有骰子朝上一面的点数之和为s.输入n,打印出s的所有可能的值出现的概率. 思路: s可能出现的值的范围为:n--6*n 1.全排列 回溯法枚举n个骰子(6面)的全排列,然 ...

  5. 剑指 Offer 60. n个骰子的点数 + 动态规划 + 空间优化

    剑指 Offer 60. n个骰子的点数 Offer_60 题目详情 题解分析 package com.walegarrett.offer; /** * @Author WaleGarrett * @ ...

  6. 剑指Offer面试题43（Java版）：n个骰子的点数

    题目:把n个骰子仍在地上.全部骰子朝上一面的点数之和为s,输入n,打印出s的全部可能的值出现的概率. 解法一:基于递归求骰子的点数,时间效率不够高 如今我们考虑怎样统计每个点数出现的次数. 要向求出n ...

  7. 《剑指offer(第二版)》面试题60——n个骰子的点数

    一.题目描述 把n个骰子仍在地上,所有的骰子朝上的一面的点数之和为s,输入n,打印出s所有可能的值出现的概率. 二.题解 <剑指offer>上给出的两种方法,尤其是代码,晦涩难懂且没有注释 ...

  8. 【编程题目】n 个骰子的点数

    67.俩个闲玩娱乐(运算).2.n 个骰子的点数.把 n 个骰子扔在地上,所有骰子朝上一面的点数之和为 S.输入 n,打印出 S 的所有可能的值出现的概率. 思路:用递归把每个骰子的可能情况变量,记录 ...

  9. 【面试题043】n个骰子的点数

    [面试题043]n个骰子的点数 题目:     把n个骰子扔在地上,所有骰子朝上一面的点数之和为s, 输入n,打印出s的所有可能的值出现的概率.   n个骰子的总点数,最小为n,最大为6n,根据排列组 ...

随机推荐

  1. Maven 初学（一）基本概念

    Pom中有三个主要元素 Groupid,artifactid,version goupid 是一个组织唯一的标识  例如 com.ibm.www artifactid  是一个工程呢ID        ...

  2. CentOS7 安装 Mono

    官网参考:http://www.mono-project.com/docs/getting-started/install/linux/#centos-7-fedora-19-and-later-an ...

  3. AMD&CMD

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  4. BZOJ-1015 StarWar星球大战 并查集+离线处理

    1015: [JSOI2008]星球大战starwar Time Limit: 3 Sec Memory Limit: 162 MB Submit: 4105 Solved: 1826 [Submit ...

  5. 【poj2342】 Anniversary party

    http://poj.org/problem?id=2342 (题目链接) 题意 没有上司的舞会... Solution 树形dp入门题. dp[i][1]表示第i个节点的子树当节点i去时的最大值,d ...

  6. hihoCoder 1195 高斯消元.一

    传送门 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Ho:喂不得了啦,那边便利店的薯片半价了! 小Hi:啥?! 小Ho:那边的便利店在打折促销啊. 小Hi:走走走, ...

  7. 使用servletAPI三种方式简单示例

    一.直接实现Action接口或集成ActionSupport类(推荐) public class HelloAction implements Action { @Override public St ...

  8. 位图索引:原理（BitMap index）

    http://www.cnblogs.com/LBSer/p/3322630.html 位图(BitMap)索引 前段时间听同事分享,偶尔讲起Oracle数据库的位图索引,顿时大感兴趣.说来惭愧,在这 ...

  9. 二、Ubuntu14.04下安装Hadoop2.4.0 （伪分布模式）

    在Ubuntu14.04下安装Hadoop2.4.0 (单机模式)基础上配置 一.配置core-site.xml /usr/local/hadoop/etc/hadoop/core-site.xml ...

  10. MYSQL远程登录权限设置(转)

    Mysql默认关闭远程登录权限,如下操作允许用户在任意地点登录: 1. 进入mysql,GRANT ALL PRIVILEGES ON *.* TO 'root'@'%' IDENTIFIED BY ...