HDU 4465 数值计算，避免溢出

数学,数值计算,求期望题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=4465题目描述:有两个盒子,每个中有n个糖果,(n<10^5)每次任选一个盒子,如果有糖就吃掉,没糖就去开另一个盒子.选中盒子1的概率为p,选中盒子2的概率为1-p.问当发现一个盒子里没有糖时,另一个盒子中糖果的个数的数学期望.解法:利用数学期望的定义,结果一共为x = 0,1,2,```,n.如果知道p(x),求sum(x*p(x))即可.为方便计算,设吃掉了i个糖果时发现…

Hdu 4465 Candy (快速排列组合+概率)

题目链接: Hdu 4465 Candy 题目描述: 有两个箱子,每个箱子有n颗糖果,抽中第一个箱子的概率为p,抽中另一个箱子的概率为1-p.每次选择一个箱子,有糖果就拿走一颗,没有就换另外一个箱子.问换箱子的时候,另外一个箱子中剩下糖果的期望值. 解题思路: 注意题目描述,其中任意一个箱子没有糖果,另一个箱子中剩下糖果个数的期望,而不是第一个箱子没有糖果.不是把其中一个箱子取空时,另一个箱子剩下糖果的期望,而是其中一个箱子取空再换另外一个箱子时,这个箱子的期望. 可以根据期望性质画出公式:an…

hdu 4465 Candy(二次项概率)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4465 参考博客:http://www.cnblogs.com/goagain/archive/2012/11/20/2778633.html 看他的分析足够了下面的代码也是他写的,觉得优美就贴下来: #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #includ…

hdu 4465 概率称号

http://acm.hdu.edu.cn/showproblem.php?pid=4465 第一直觉概率DP但很快被否定,发现只有一个简单的二项分布,但感情的表达,没有对生命和死亡的例子.然后找到准确的问题,将不被处理, 事实上与思考C递归成为O(1).每次乘以p 代码看这里http://fszxwfy.blog.163.com/blog/static/44019308201338114456115/ 贴一个我没看懂的代码 //#pragma comment(linker, "/STACK:1…

HDU 4465 Candy （数学期望）

题意:有两个盒子各有n个糖(n<=2*105),每天随机选1个(概率分别为p,1-p),然后吃掉一颗糖.直到有一天打开盒子一看,这个盒子没有糖了.输入n,p,求此时另一个盒子里糖的个数的数学期望. 思路:假设没糖的是A盒子,而B盒子还有0~n个糖.由于B盒子还有0个糖的情况的期望必为0,所以省略,只需要计算1~n的. (1)当A盒没有糖时,B盒就可能有1~n个糖,概率为C(n+i,i)*(pn+1)*(1-p)n-i.为啥还带个大C?这是情况的种数(想象取糖时还有个顺序,有C种可能的顺序),不然…

hdu 4465 Candy 数学

思路:易知结果为 ∑(n-k)*C(n+k,k)*(p^(n+1)*q^k+q^(n+1)*p^k). 注意不能直接算,注意点技巧!!!看代码代码如下: #include<iostream> #include<stdio.h> #include<algorithm> #include<iomanip> #include<cmath> using namespace std; int main(){ ; double p,q,ans1,ans2,…

hdu 4465 Candy

题解: 由题意得需要运用: C(m,n)=exp(logC(m,n)) f[]=; ; i<=; i++) f[i]=f[i-]+log(i*1.0); double logC(int m,int n) { return f[m]-f[n]-f[m-n]; } #include<stdio.h> #include<math.h> ]; double logC(int m,int n) { return f[m]-f[n]-f[m-n]; } int main() { f[]…

HDU 4465 - Candy（概率与数学优化）

2012成都Regional的B题,花了一个小时推出了式子,但是搞了好久发现都控制不了精度,后来突然想到组合数可以用log优化,改了之后就AC了比较水的概率题 #include <stdio.h> #include <math.h> #define maxn 200005 double f[2*maxn]; double c(int a,int b) { return f[a]-f[a-b]-f[b]; } int main() { f[0]=0; for(int i=1;i&l…

hdu 4465 Candy 2012 成都现场赛

/** 对于大数的很好的应用,,缩小放大,,保持精度 **/ #include <iostream> #include <cmath> #include <algorithm> #include <cstdio> using namespace std; int main() { double n,p; ; while(cin>>n>>p){ double p1 = log(p+0.0); -p+0.0); )*p1; )*p2; ,…

hdu 4465 Candy（2012 ACM-ICPC 成都现场赛）

简单概率题,可以直接由剩余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的形式,这样可以保存精度. #include<algorithm> #include<iostream&g…

HDU - 4465 期望 + 取log优化

思路:这个求期望的公式很容易得到,但是在算的时候我们会遇到一个问题,就是组合数太大了根本存不下, 这时候就可以在计算的时候都取log,最后复原... 以前没遇到过.. #include<bits/stdc++.h> #define LL long long #define fi first #define se second #define mk make_pair #define PII pair<int, int> using namespace std; ; const i…

HDU 4465 递推与double的精确性

题目大意不多说了这里用dp[i][0] 代表取完第一个盒子后第二个盒子剩 i 个的概率,对应期望就是dp[i][0] *i dp[i][1] 就代表取完第二个盒子后第一个盒子剩 i 个的概率 dp[i][0] = p^(n+1) * (1-p)^(n-i) * C(2*n-i , n-i) = p^(n+1) * (1-p)^(n-i) * (2*n-i)! / (n-i)! / n! dp[i+1][0] = p^(n+1) * (1-p)^(n-i-1) * C(2*n-i-1 ,…

hdu 4465 求期望（C（m,n）太大用log优化）

/* 坑啊数学函数的运用log处理,exp还原 tle好长时间,一直用g++交,最后把别人正确的代码交上也是tle,用c++交一遍ac 题意:有两个数量为n的糖果,一个人开始吃,吃到最后有一堆剩余为0的时候不吃,把另一堆留给另外一个另一个人,求另一个人得到糖果的期望, 这个人吃第一堆糖果的概率是p,第二堆糖果的概率是1-p,并且当他把一堆糖果吃完,再次吃的时候发现没有了才结束. 解:分别求出当第一堆剩余的期望,第二堆剩余的期望加起来 (n-i)*C(n+i,i)*pow(p,n+1)*pow(…

hdu 4465 Candy （非原创）

LazyChild is a lazy child who likes candy very much. Despite being very young, he has two large candy boxes, each contains n candies initially. Everyday he chooses one box and open it. He chooses the first box with probability p and the second box wi…

2012Chengdu B （快速组合数）

B - Candy Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit Status Practice HDU 4465 Appoint description: System Crawler (2014-10-17) Description LazyChild is a lazy child who likes candy very much. Despite being v…

pandas函数get_dummies的坑

转载:https://blog.csdn.net/mvpboss1004/article/details/79188190 pandas中的get_dummies得到的one-hot编码数据类型是uint8,进行数值计算时会溢出!!! import pandas as pd import numpy as np a = [1, 2, 3, 1] one_hot = pd.get_dummies(a) print(one_hot.dtypes) print(one_hot) print(-one_…

Java中的基本数据类型语法补充

变量要先赋值后使用不给变量赋值代表什么不赋值就使用会怎样 (会报错) 计算并赋值运算符作用是为了让代码更加简洁.比如 a = a + 10,可以简化为 a+=10 += -= *= /= %= &= ^= |= <<= >>= >>>= public class CalcAndAssign { public static void main(String[] args) { int a = 16; a >>= 2; System.out.…

程序员之---C语言细节20（符号和有符号之间转换、两数相加溢出后数值计算）

主要内容:无符号和有符号之间转换.两数相加溢出后数值计算 #include <stdio.h> /* 这个函数存在潜在漏洞 */ float sum_elements(float a[], unsigned length) { int i; float result = 0; for(i = 0; i <= length - 1; i++) { result += a[i]; printf("a[%d] = %f \n",i,a[i]); } return resul…

hdu 2519 新生晚会（求排列组合时容易溢出）

#include<stdio.h> #include<algorithm> using namespace std; void cal(int n,int m) { ; m=min(m,n-m); int j=m; ;m--,i--) { ans*=i; } ;i<=j;i++) { ans/=i; } printf("%I64d\n",ans); } int main() { int _case; int n,m; scanf("%d"…

HDU 2227 Find the nondecreasing subsequences (DP+树状数组+离散化)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2227 Find the nondecreasing subsequences Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) …

HDU 5047 Sawtooth（大数优化+递推公式）

http://acm.hdu.edu.cn/showproblem.php?pid=5047 题目大意: 给n条样子像“m”的折线,求它们能把二维平面分成的面最多是多少. 解题思路: 我们发现直线1条:2平面:2直线:4平面:3直线:7平面......因为第n条直线要与前面n-1条直线都相交,才能使分的平面最多,则添加第n条直线,平面增加n个: 所以公式是面F = 2 + 2 + 3 + ......+ n = (1+n)*n/2 + 1 因为题目的是“M”的折线,一个“M”有4条线将平面分成2…