题意:有两个盒子各有n个糖,每次随机选一个(概率分别为p,1-p),然后吃掉,直到有一次,你打开盒子发现,没糖了!

输入n,p,求另一个盒子里糖的个数的数学期望。

析:先不说这个题多坑,首先要用long double来实现高精度,我先用的double一直WA,后来看了题解是用long double,

改了,可一直改不对,怎么输出结果都是-2.00000,搞了一晚上,真是无语,因为我输入输出数据类型是long double,

结果一直不对 ,可能是我的编译器是C89的吧,和C语言,输入输出格式不同,但是我改成C99也不对。。。

最后一看题解明白了,输入输出不需要那么高的精度,用double就足够了,可算是废了一晚上的时间。。。悲惨

下面分析这个题:

首先我们不知道最后打开的是哪个盒子,所以我们需要都考虑,当最后打开是第一个时,此时另一个有i个,那么在些之前一定是打开过

n+(n-i)次盒子,其中n次是第一个盒子,n-i次是第二个盒子,这是明显是一个二项分布,所以概率为C(2n-i, n)pn+1(1-p)n-i,这个地方

一定要注意是n+1,不是n,因为我们也要把最后一次的算上。

这个式子很明显是正确的,但是有一个大问题,那就是C(2n-i, n)这数据实在是太大了(因为n太大了),肯定无法用整形保存,

甚至无法用浮点保存,所以我们对它进行高精度处理,大数选的是对数法,令v1(i) = ln(C(2n-i, n)) + (n+1)ln(p) + (n-i)ln(1-p).

这个地方一定是加减,不是乘除,这是对数,刚开始又忘了,结果一直不对。。。

那么最后打开的是第一个盒子数学期望是ev1(i)

同理最后一个打开的是第二个盒子时,对数为v2(i) = ln(C(2n-i, n)) + (n+1)ln(1-p) + (n-i)ln(p),概率为ev2(i)

那么总的数学期望是sum{i*(ev1(i)+ev2(i))}((ev1(i)+ev2(i))这个是概率)。

这个题需要高精度,所以用long double来存储,真是坑。。。。很难想到。

代码如下:

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath> using namespace std;
typedef long double LD;
const int maxn = 400000 + 10;
LD log_ld[maxn]; void init(){
log_ld[1] = 0.0;
for(int i = 2; i < maxn; ++i)
log_ld[i] = log_ld[i-1] + (LD)log((LD)i);
} double solve(int n, LD p){
LD ans = 0.0;
for(int i = 1; i <= n; ++i)
ans += (LD)i * exp(log_ld[2*n-i] - log_ld[n] - log_ld[n-i] + (LD)(n+1)*log(p) + (LD)(n-i)*log(1.0-p)); for(int i = 1; i <= n; ++i)
ans += i * exp(log_ld[2*n-i] - log_ld[n] - log_ld[n-i] + (LD)(n+1)*log(1.0-p) + (LD)(n-i)*log(p));
return ans;
} int main(){
// freopen("in.txt", "r", stdin);
int n, kase = 0;
double p;
init(); while(~scanf("%d %lf", &n, &p)){
double ans = solve(n, p);
printf("Case %d: %.6f\n", ++kase, ans);
}
return 0;
}

UVa 1639 Candy (数学期望+组合数学+高精度存储)的更多相关文章

  1. 紫书 例题 10-17 UVa 1639(数学期望+分数处理+处理溢出)

    设当前有k个,那么也就是说拿到其他图案的可能是(n-k)/n 那么要拿到一个就要拿n/(n-k)次 所以答案就是n(1/n + 1/(n-1) ......1/2 + 1 / 1) 看起来很简单,但是 ...

  2. 紫书 例题 10-17 UVa 1639(数学期望+对数保存精度)

    设置最后打开的是盒子1, 另外一个盒子剩下i个 那么在这之前打开了n + n - i次盒子 那么这个时候的概率是C(2 * n - i, n) p ^ (n+1) (1-p)^ (n - i) 那么反 ...

  3. UVa 1639 - Candy(数学期望 + 精度处理)

    链接: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem& ...

  4. UVA - 1639 -Candy

    题目链接:https://vjudge.net/problem/UVA-1639 题目大意: 有两个糖果盒,每个盒子里面有n个糖果,每天随机选一个(概率分别为p,1-p),然后吃一颗糖.直到有一天,打 ...

  5. uva 1639 Candy (对数处理精度)

    https://vjudge.net/problem/UVA-1639 有两个盒子各有n(n≤2*10 5 )个糖,每天随机选一个(概率分别为p,1-p),然后吃一颗糖. 直到有一天,打开盒子一看,没 ...

  6. 紫书 例题 10-16 UVa 12230(数学期望)

    感觉数学期望的和化学里面求元素的相对原子质量的算法是一样的 就是同位素的含量乘上质量然后求和得出 这道题因为等待时机是0到2*l/v均匀分配的,所以平均时间就是l/v 再加上过河的l/v, 最后加上步 ...

  7. UVA 1639 Candy (组合数+精度)

    题意:两个箱子,每个箱子有n颗糖,每次有p的概率拿1号箱子的一颗糖出来(有1-p的概率拿2号箱子的一颗糖出来),问当打开某个箱子为空的时候,另一个箱子的期望糖的数量是多少 题解:枚举另一个箱子的糖的数 ...

  8. UVA - 1639 Candy (概率,精度)

    X表示剩下的糖数量,如果最后打开的是p对应的盒子.划分:Xi表示剩下i个糖,最后一次选的概率为p, 前面的服从二项分布.根据全概率公式和期望的线性性,求和就好了. 精度处理要小心,n很大,组合数会很大 ...

  9. 数学期望和概率DP题目泛做(为了对应AD的课件)

    题1: Uva 1636 Headshot 题目大意: 给出一个000111序列,注意实际上是环状的.问是0出现的概率大,还是当前是0,下一个还是0的概率大. 问题比较简单,注意比较大小: A/C & ...

随机推荐

  1. UI5-文档-4.22-Expression Binding

    有时预定义的SAPUI5类型不够灵活,您希望在视图中执行简单的计算或格式化——这正是表达式真正有用的地方.我们使用它们根据数据模型中的当前数字格式化价格. Preview The price is n ...

  2. PostgreSql别名区分大小写的问题

    PostgreSql是区分大小写的,如果别名的大小不一致就会提示错误: SELECT * FROM ( SELECT cpi."product_item_id" "PRO ...

  3. EL表达式获取请求

    1EL 表达式语法: ${ EL 表达式}    它的特点: 1自动转换类型,  2使用简单. 2 EL表达式具有类型无关性,可以使用“.”或者“[]”操作符在相应的作用域(page  ,  requ ...

  4. one by one 项目 part 2

    在网上百度了一下mySQL常用语句 ,整理如下: 1.说明:创建数据库 CREATE DATABASE database-name 2.说明:删除数据库 drop database dbname 3. ...

  5. js中的Event对象

    event代表事件的状态,例如触发event对象的元素,鼠标的位置及状态,按下的键等等 event对象只在事件发生的过程中才有效. <!DOCTYPE html><html lang ...

  6. window 下复制 cmd内容

    https://www.microsoft.com/resources/documentation/windows/xp/all/proddocs/en-us/windows_dos_copy.msp ...

  7. Anaconda中python加入环境变量

    1.我的电脑---高级系统设置 2.选中环境变量,保存. 3.在系统环境变量PATH中,加入Anaconda3及Script路径加入其中 4.测试python

  8. 通过jquery,从json中读取数据追加到html中

    1.下载安装jquery   可通过下面的方法引入在线版本的js: <script src="https://apps.bdimg.com/libs/jquery/2.1.4/jque ...

  9. vs code 插件推荐

    通用插件 HTML Snippets 超级实用且初级的 H5代码片段以及提示 HTML CSS Support 让 html 标签上写class 智能提示当前项目所支持的样式新版已经支持scss文件检 ...

  10. 第八章 高级搜索树 (xa1)红黑树:动机