剑指offer-74、n个骰⼦的点数

16 阅读4分钟

题目描述

把 n 个骰⼦扔在地上,所有骰⼦朝上⼀⾯的点数之和为 s 。输⼊ n ,打印出 s 的所有可能的值出现的概率。

你需要⽤⼀个浮点数数组返回答案,其中第 i 个元素代表这 n 个骰⼦所能掷出的点数集合中第 i ⼩的那个的概率。

示例1:

输⼊: 1 输出: [0.16667,0.16667,0.16667,0.16667,0.16667,0.16667]

示例2

输⼊: 2 输出:[0.02778,0.05556,0.08333,0.11111,0.13889,0.16667,0.13889,0.11111,0.08333,0.05556,0.02778]

思路及解答

暴力递归

枚举所有骰子组合。递归计算每个骰子的点数,统计所有可能的和

public class Solution {
    public double[] dicesProbability(int n) {
        // 骰子点数范围:n到6n,共5n+1种可能
        int[] counts = new int[5 * n + 1];
        
        // 递归统计所有可能的和
        backtrack(n, 0, counts);
        
        // 计算概率
        double total = Math.pow(6, n);
        double[] res = new double[counts.length];
        for (int i = 0; i < counts.length; i++) {
            res[i] = counts[i] / total;
        }
        return res;
    }
    
    private void backtrack(int remain, int sum, int[] counts) {
        if (remain == 0) {
            counts[sum - remain]++; // 统计和出现的次数
            return;
        }
        
        // 当前骰子可以是1到6点
        for (int i = 1; i <= 6; i++) {
            backtrack(remain - 1, sum + i, counts);
        }
    }
}

如果使⽤暴⼒法,每⼀个骰⼦扔到 1 - 6 的概率都是 1/6,如果有 n 个 骰⼦,先不看重复的情况,⼀共有 6n6^n 种情况,点数的范围是 n ~ 6n ,也就是 5n+1 种。

  • 时间复杂度:O(6ⁿ),每个骰子有6种可能,n个骰子组合数为6ⁿ
  • 空间复杂度:O(n),递归调用栈的深度

以上的计算复杂度实在太⾼,我们不能接受。

动态规划(推荐)

其实,这道题可以⽤动态规划来处理, 1 个骰⼦的情况是已知的,⽽ 2 个骰⼦的情况呢? 2 个骰⼦的情况,可以使⽤ 1 个骰⼦的情况推出, 3 个骰⼦的情况,可以使⽤ 2 个骰⼦的结果推出...

dp[i][j]表示i个骰子和为j的出现次数

执行过程示例(n=2):

初始化dp[1][1]=1, dp[1][2]=1, ..., dp[1][6]=1
计算dp[2][2] = dp[1][1] = 1
dp[2][3] = dp[1][2] + dp[1][1] = 2
...
dp[2][12] = dp[1][11] (不存在) + ... + dp[1][6] = 1
public class Solution {
    public double[] dicesProbability(int n) {
        // dp[i][j]:i个骰子和为j的出现次数
        int[][] dp = new int[n + 1][6 * n + 1];
        
        // 初始化:1个骰子的情况
        for (int j = 1; j <= 6; j++) {
            dp[1][j] = 1;
        }
        
        // 填充状态转移表
        for (int i = 2; i <= n; i++) {            // 骰子数量从2到n
            for (int j = i; j <= 6 * i; j++) {     // 和的范围:i到6i
                for (int k = 1; k <= 6 && k < j; k++) { // 最后一个骰子的点数
                    dp[i][j] += dp[i - 1][j - k];
                }
            }
        }
        
        // 计算概率
        double total = Math.pow(6, n);
        double[] res = new double[5 * n + 1];
        for (int j = n; j <= 6 * n; j++) {
            res[j - n] = dp[n][j] / total;
        }
        return res;
    }
}
  • 时间复杂度:O(n²),外层循环n次,内层循环最多6n次
  • 空间复杂度:O(n²),需要二维数组存储中间结果

空间优化动态规划

通过观察发现当前状态只依赖前一个状态,可以优化空间。通过滚动数组减少空间使用

public class Solution {
    public double[] dicesProbability(int n) {
        // 使用两个一维数组交替更新
        int[] prev = new int[6 * n + 1];
        int[] curr = new int[6 * n + 1];
        
        // 初始化第一个骰子
        for (int j = 1; j <= 6; j++) {
            prev[j] = 1;
        }
        
        // 动态规划填充
        for (int i = 2; i <= n; i++) {
            Arrays.fill(curr, 0); // 清空当前数组
            for (int j = i; j <= 6 * i; j++) {
                for (int k = 1; k <= 6 && k < j; k++) {
                    curr[j] += prev[j - k];
                }
            }
            // 交换数组,准备下一轮
            int[] temp = prev;
            prev = curr;
            curr = temp;
        }
        
        // 计算概率
        double total = Math.pow(6, n);
        double[] res = new double[5 * n + 1];
        for (int j = n; j <= 6 * n; j++) {
            res[j - n] = prev[j] / total;
        }
        return res;
    }
}
  • 时间复杂度:O(n²),外层循环n次,内层循环最多6n次
  • 空间复杂度:O(n),只保留前一轮的结果,空间从O(n²)降到O(n)

数学公式法(多项式展开)

和为s的概率对应于(x+x²+...+x⁶)ⁿ展开式中xˢ的系数

public class Solution {
    public double[] dicesProbability(int n) {
        // 初始化多项式系数:1个骰子时各项系数为1
        int[] coeff = new int[6 * n + 1];
        for (int j = 1; j <= 6; j++) {
            coeff[j] = 1;
        }
        
        // 多项式乘法:计算(1x+1x²+...+1x⁶)^n
        for (int i = 2; i <= n; i++) {
            int[] newCoeff = new int[6 * i + 1];
            // 多项式乘法计算
            for (int j = 1; j <= 6; j++) {
                for (int k = i - 1; k <= 6 * (i - 1); k++) {
                    newCoeff[k + j] += coeff[k];
                }
            }
            coeff = newCoeff;
        }
        
        // 计算概率
        double total = Math.pow(6, n);
        double[] res = new double[5 * n + 1];
        for (int j = n; j <= 6 * n; j++) {
            res[j - n] = coeff[j] / total;
        }
        return res;
    }
}
  • 时间复杂度:O(n²),多项式乘法计算
  • 空间复杂度:O(n),存储多项式系数