uva 11605 - Lights inside a 3d Grid(概率)
2016-02-25 21:10
316 查看
题目链接:uva 11605 - Lights inside a 3d Grid
题目大意:给定一个三维坐标系大小,每一个位置有一个灯。初始状态为关。每次随机选中两个点,以这两点为对角线的长方体内全部灯转变状态。操作K次。问说平均情况下。最后会有多少栈灯亮着。
解题思路:枚举坐标系上的点。计算单个点亮着的概率,然后累加即使总体的期望。
对于一个点x,y,z,分别考虑每维坐标系,比如x。选中的概率为px=2∗(n−x+1)∗x−1n∗n,三维坐标均选中的概率p即为该点被选中的概率。
可是对于一点来说,由于操作K次,所以最后灯为亮的话。操作到灯的次数一定要为奇数才行,所以有∑C(iK)pi(1−p)K−i(i为奇数)
===》(1−p+p)K−(1−p−p)K2
===》1−(1−2p)K2
题目大意:给定一个三维坐标系大小,每一个位置有一个灯。初始状态为关。每次随机选中两个点,以这两点为对角线的长方体内全部灯转变状态。操作K次。问说平均情况下。最后会有多少栈灯亮着。
解题思路:枚举坐标系上的点。计算单个点亮着的概率,然后累加即使总体的期望。
对于一个点x,y,z,分别考虑每维坐标系,比如x。选中的概率为px=2∗(n−x+1)∗x−1n∗n,三维坐标均选中的概率p即为该点被选中的概率。
可是对于一点来说,由于操作K次,所以最后灯为亮的话。操作到灯的次数一定要为奇数才行,所以有∑C(iK)pi(1−p)K−i(i为奇数)
===》(1−p+p)K−(1−p−p)K2
===》1−(1−2p)K2
#include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; int N, M, P, K; inline double getp (double n, double x) { double s = n * n; double t = 2 * (n - x + 1) * x - 1; return t / s; } inline double handle (double p) { return (1 - pow(1 - 2 * p, K)) / 2; } double solve () { double ret = 0; for (int x = 1; x <= N; x++) { double px = getp(N, x); for (int y = 1; y <= M; y++) { double py = getp(M, y); for (int z = 1; z <= P; z++) { double pz = getp(P, z); ret += handle(px * py * pz); } } } return ret; } int main () { int cas; scanf("%d", &cas); for (int kcas = 1; kcas <= cas; kcas++) { scanf("%d%d%d%d", &N, &M, &P, &K); printf("Case %d: %.10lf\n", kcas, solve()); } return 0; }
相关文章推荐
- 大型网站系统与Java中间件读书笔记
- HTML,CSS和Javascript综合练习-1
- JSP中的数据库使用
- 性能测试工具 nGrinder 项目剖析及二次开发
- CramFS文件系统
- C语言静态函数静态变量
- 02—策略模式
- 数据结构——二叉查找树
- 往文件中写入内容(代码示例)
- new与malloc的区别
- python logging 替代print 输出内容到控制台和重定向到文件
- 字符串的定义、替换、比较、截取
- 谷歌推荐的图片加载框架--Glide的简单介绍
- Java 常见面试题整理
- HDU1532 - Drainage Ditches(网络流)
- 我们编程吧 之 Linux 学习手册0.41
- 网络编程_UDP编程原理_发送类型_cs与bs区别JAVA187-188
- 【POJ 2750】 Potted Flower(线段树套dp)
- TWinControl.DefaultHandler里的CallWindowProc(FDefWndProc)还挺有深意的,TButton对WM_PAINT消息的处理就是靠它来处理的(以前不明白为什么总是要调用inherited,其实就是没有明白TWinControl.DefaultHandler的真正用处,而且还很有用)
- Log4j详解