【题目】
描述描述:
骰子是一个六面分别刻有一到六点的立方体,每次投掷骰子,理论上得到 1 1 1 到 6 6 6 的概率都是 1 / 6 1/6 1/6
有骰子一颗,连续投掷 n n n 次,问点数总和大于等于 x x x 的概率是多少
输入格式:
输入一行 2 2 2 个整数,分别表示 n , x n,x n,x,其中 1 1 1 ≤ n n n ≤ 24 , 0 24,0 24,0 ≤ x x x < 150 150 150
输出格式:
1 1 1 行,一个分数,要求用最简的形式精确的表达投掷 n n n 次骰子,总数大于等于 x x x 的概率。如果是 0 / 1 0/1 0/1,则输出 0 0 0,如果 1 / 1 1/1 1/1 输出 1 1 1
样例数据:
输入
3 9
输出
20/27
【分析】
题解:概率DP入门题
用 f i , j f_{i,j} fi,j 表示第 i i i 次掷骰子,总点数为 j j j 的方案总数
那么很显然,最后的总方案数就是 ∑ i = x 6 ∗ n f n , i \sum_{i=x}^{6*n}f_{n,i} ∑i=x6∗nfn,i
状态转移: f i , j + = ∑ k = 1 6 f i − 1 , j − k f_{i,j}+=\sum_{k=1}^{6}f_{i-1,j-k} fi,j+=∑k=16fi−1,j−k,相当于就是枚举这一位投的点数,然后根据加法原理累加
由于总方案数为
6
n
6^n
6n,所以最后除个
6
n
6^n
6n,再用欧几里得算法化一下简就可以了
【代码】
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
using namespace std;
long long f[N][N];
int main()
{
int n,i,j,k,x;
scanf("%d%d",&n,&x);
f[0][0]=1;
for(i=0;i<n;++i)
for(j=0;j<=6*i;++j)
for(k=1;k<=6;++k)
f[i+1][j+k]+=f[i][j];
long long ans=0,sum=1;
for(i=1;i<=n;++i) sum*=6;
for(i=x;i<=6*n;++i) ans+=f[n][i];
long long t=__gcd(ans,sum);
if(sum/t==1) printf("%lld",ans/t);
else printf("%lld/%lld",ans/t,sum/t);
return 0;
}