hdu 4465 Candy(2012 ACM-ICPC 成都现场赛)
2013-09-01 20:37
351 查看
简单概率题,可以直接由剩余n个递推到剩余0个。现在考虑剩余x个概率为(1-p)的candy时,概率为C(2 * n - x, x) * pow(p, n + 1) *pow(1 - p, n - x);
在写出x - 1的情况,就可以发现组合数可以直接递推,所以可以直接求。但是考虑到p可能很小,n可能很大,这样的话直接用pow函数会丢失精度,我们可以把double类型写成log10的形式,这样可以保存精度。
在写出x - 1的情况,就可以发现组合数可以直接递推,所以可以直接求。但是考虑到p可能很小,n可能很大,这样的话直接用pow函数会丢失精度,我们可以把double类型写成log10的形式,这样可以保存精度。
#include<algorithm> #include<iostream> #include<cstring> #include<cstdlib> #include<fstream> #include<sstream> #include<bitset> #include<vector> #include<string> #include<cstdio> #include<cmath> #include<stack> #include<queue> #include<stack> #include<map> #include<set> #define FF(i, a, b) for(int i=a; i<b; i++) #define FD(i, a, b) for(int i=a; i>=b; i--) #define REP(i, n) for(int i=0; i<n; i++) #define CLR(a, b) memset(a, b, sizeof(a)) #define debug puts("**debug**") #define LL long long #define PB push_back #define SL(a) strlen(a) using namespace std; const int N = 11; const int MOD = 1e9 + 7; int main() { int n, cas = 1, i, j; double p, ans, now; while(scanf("%d%lf", &n, &p) != EOF) { now = log10(p) * (n + 1); ans = pow(10.0, now) * n; for(i = n - 1; i > 0; i --) { now = now + log10(2 * n - i + 0.0) - log10(n - i + 0.0) + log10(1 - p); ans += pow(10.0, now) * i; } p = 1 - p; now = log10(p) * (n + 1); ans += pow(10.0, now) * n; for(i = n - 1; i > 0; i --) { now = now + log10(2 * n - i + 0.0) - log10(n - i + 0.0) + log10(1 - p); ans += pow(10.0, now) * i; } printf("Case %d: %.6lf\n", cas ++, ans); } }
相关文章推荐
- hdu 4465 Candy(2012 ACM-ICPC 成都现场赛)
- 2012 ICPC/ACM 成都现场赛 Candy
- hdu 4472 Count (2012 ACM-ICPC 成都现场赛)
- hdu 4472 Count (2012 ACM-ICPC 成都现场赛)
- 【数学期望】【2012 ACM/ICPC 成都赛区现场赛】【B.Candy】
- HDU-4472-Count (2012 ACM/ICPC成都现场赛)
- HDU-4464-Browsing History (2012 ACM/ICPC成都现场赛!)
- HDU 4790 2013 ACM/ICPC 成都赛区现场赛 J. Just Random
- 2012ACM/ICPC成都赛区现场赛A题---Browsing History(hdu4464)
- 2012 ACM/ICPC 成都赛区网络赛 HDU 4296
- 2012ACM/ICPC成都赛区现场赛I题---Count(hdu4472)
- HDU 4790 Just Random (2013ACM/ICPC成都赛区现场赛J题)
- HDU 4465 Candy 概率 (2012成都区域赛B题)
- 【DP】【2012 ACM/ICPC 成都赛区现场赛】【I.Count】
- hdu 4473 Exam (2012ACM成都赛区现场赛J题)
- HDU 5536 2015ACM-ICPC长春赛区现场赛J题
- 2014ACM/ICPC亚洲区鞍山赛区现场赛I (hdu 5078)
- HDU 5073 Galaxy (2014 ACM/ICPC 鞍山赛区现场赛D题)
- Hdu 4035 Maze (dp求期望) - 2011 ACM/ICPC 成都赛区网络预选赛 1005
- HDU 5115 Dire Wolf(2014ACM/ICPC北京赛区现场赛D)