【POJ1322】Chocolate (生成函数)
Description
有 c c c种颜色的巧克力, 每种颜色有无限个. 现在每次取出一个巧克力, 其颜色等概率为 1 … c 1\dots c 1…c中的一种.
问最终有 m m m种颜色的巧克力个数为奇数的概率.
n ≤ 1 0 6 , c ≤ 100 n\le 10^6, c\le 100 n≤106,c≤100
Solution
睿智dp题。
设 f i , j f_{i,j} fi,j表示当前取了 i i i个、有 j j j种颜色是奇数个的概率,直接dp即可,时间复杂度 O ( n c ) O(nc) O(nc)。
Code
/*************************************** Au: Hany01* Prob: poj1322* Date: Nov, 18th, 2018* Email: hany01@foxmail.com
**************************************/#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>using namespace std;typedef long long LL;
typedef long double LD;
typedef pair<int, int> PII;
typedef vector<int> VI;
#define File(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)
#define Rep(i, j) for (register int i = 0, i##_end_ = j; i < i##_end_; ++ i)
#define For(i, j ,k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define SZ(a) ((int)(a.size()))
#define ALL(a) a.begin(), a.end()
#define PB(a) push_back(a)
#define MP(a, b) make_pair(a, b)
#define X first
#define Y second
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define y1 WoXiHuanNiA
#ifdef hany01
#define debug(...) fprintf(stderr, __VA_ARGS__)
#else
#define debug(...)
#endiftemplate<typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template<typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }
template<typename T> inline T read() {register char c_; register T _, __;for (_ = 0, __ = 1, c_ = getchar(); !isdigit(c_); c_ = getchar()) if (c_ == '-') __ = -1;for ( ; isdigit(c_); c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);return _ * __;
}
//EOTconst int maxn = 105;double fac[maxn];
inline void Init(int n) {fac[0] = 1;For(i, 1, n) fac[i] = fac[i - 1] * i;
}
inline double C(int n, int m) { return fac[n] / fac[m] / fac[n - m]; }inline double fpm(double a, int b) {register double ans = 1;for ( ; b; b >>= 1, a *= a)if (b & 1) ans *= a;return ans;
}double p[maxn], p_[maxn];int main() {#ifdef hany01File("poj1322");
#endifInit(100);for (register int n, m, c; scanf("%d", &c) != EOF && c; ) {n = read<int>(), m = read<int>();if (m > c || m > n || ((n - m) & 1)) { puts("0.000"); continue; }if (!n && !m) { puts("1.000"); continue; }register double ans = 0;Set(p, 0), Set(p_, 0);For(i, 0, m) For(j, 0, c - m) {int t = c - 2 * i - 2 * j;if (t > 0) p[t] += C(m, i) * C(c - m, j) * ((i & 1) ? -1 : 1) / fpm(2, c);else p_[-t] += C(m, i) * C(c - m, j) * ((i & 1) ? -1 : 1) / fpm(2, c);}For(k, 1, c) ans += fpm(1. * k / c, n) * (p[k] + ((n & 1) ? -1 : 1) * p_[k]) * C(c, m);printf("%.3f\n", ans);}return 0;
}
True-Solution
这题的时间复杂度可以做到 O ( c log ) O(c\log ) O(clog)。
首先要保证 m ≤ c , m ≤ n m\le c, m\le n m≤c,m≤n且 n , m n,m n,m同奇偶,否则概率为 0 0 0。
由于是排列,考虑指数型生成函数。
对于出现奇数次的巧克力,其指数型生成函数为:
f 1 ( x ) = e x − e − x 2 f_1(x) = \frac{e^x - e^{-x}}{2} f1(x)=2ex−e−x
出现偶数次的巧克力的生成函数为:
f 2 ( x ) = e x + e − x 2 f_2(x) = \frac{e^x + e^{-x}}{2} f2(x)=2ex+e−x
方案数即为 ( c m ) × [ x n ] ( f 1 ( x ) m × f 2 ( x ) c − m ) × n ! \binom{c}{m}\times [x^n](f_1(x)^m\times f_2(x)^{c-m})\times n! (mc)×[xn](f1(x)m×f2(x)c−m)×n!。
(在 c c c中颜色中选择 m m m个取奇数次, n ! n! n!是指数型生成函数需要乘的)
这时,时间复杂度的瓶颈在于 n n n,考虑怎么把 n n n去掉。
我们可以将 f 1 ( x ) m × f 2 ( x ) c − m = ( e x − e − x 2 ) m × ( e x + e − x 2 ) c − m f_1(x)^m\times f_2(x)^{c-m}=\left(\frac{e^x-e^{-x}}{2}\right)^m\times \left(\frac{e^x+e^{-x}}{2}\right)^{c-m} f1(x)m×f2(x)c−m=(2ex−e−x)m×(2ex+e−x)c−m看作一个关于 e x e^x ex的多项式。 e k x e^{kx} ekx的系数对第 n n n项的贡献是 k n n ! \frac{k^n}{n!} n!kn,就可以把 n ! n! n!消掉了!
用FFT计算系数做到 c log c\log clog(上面的代码是 O ( c 2 ) O(c^2) O(c2)计算的)。
【POJ1322】Chocolate (生成函数)相关推荐
- POJ-1322 Chocolate 动态规划
这题当M=N=0的时候要输出1.000 刚写的时候默认从第二次开始取了. 详见代码: #include <cstdlib> #include <cstring> #inclu ...
- [指数型生成函数专练]chocolate,红色病毒问题,排列组合,字串数
文章目录 T1:chocolate 题目 题解 code T2:"红色病毒"问题 题目 题解 code T3:排列组合 题目 题解 code T4:字串数 题解 code T1:c ...
- ACM里的生成函数初探
生成函数 前情提要:由于本人不太喜欢写博客(懒),所以写的时候一般是由于某些原因写的专题内容.这次是2020CCPC长春热身赛的D题没做出来. 大佬论文放前面:推荐阅读:毛杰明 母函数的性质及应用 基 ...
- ASPJPEG缩略图生成函数
好久没有发文章,贴一段代码出来晒晒!一段aspjpeg组件生成缩略图的代码,有4种生成方式,建议用最后一种,生成的缩略图最清晰而且不会拉伸.变形!做图片生成最好不过! // 缩略图生成函数 Code ...
- 2018 ICPC 青岛 L. Sub-cycle Graph(生成函数)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 L. Sub-cycle Graph Weblink https://zoj.pintia.cn/pr ...
- (每日一题)P4841 [集训队作业2013]城市规划 (无向连通图计数)(普通生成函数 + 多项式求逆)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 每日一题(莫反 / 多项式 / 母函数 / 群论) 2021.4.14 生成函数 + 多项式求逆 Pr ...
- 2019 ACM - ICPC 上海网络赛 E. Counting Sequences II (指数型生成函数)
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
- 解题报告(四)生成函数(ACM/ OI)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量的题解和代码,题目难度不一 ...
- 解题报告(二)E、(BZOJ3513) [MUTC2013] idiots(生成函数 + FFT + 组合计数)
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
最新文章
- Java基础——组合关系
- Tomcat根目录下静态html怎样访问,避免走坑
- 【服务器】一次对Close_Wait 状态故障的排查经历
- 重磅!英特尔终于挤出10nm芯片 六大技术战略,震动芯片届
- pytorch: 在训练中保存模型,加载模型
- 服务器不删档的设置_不删档预捏脸开启,快来体验次世代黑科技!
- IT行业的职员加班到底有没有价值?
- Tensorflow——Dropout(解决过拟合问题)
- 微信小程序云开发教程-小程序代码发布和上线运行
- 人人商城小程序消息服务器配置,微擎系统人人商城小程序前端配置教程
- 删除下拉框只找23火星软件_下拉框软件找28火星下拉
- excel从只有省市县的地址中分别提取省、市、县
- 手把手教你实现echarts3的折线图下钻drilldown功能系列篇二
- 边缘计算在物联网领域的发展前景
- win7打印机服务器修改ip,win7系统电脑更换IP地址后打印机不能打印文件了的解决方法...
- 零基础:21天搞定Python分布爬虫
- 网线插座板上网络模块接线
- [人工智能-深度学习-29]:卷积神经网络CNN - 全连接网络与卷积网络结构的互为等效与性能比较
- Python pywin32(一)
- chrome调试与设置 remote设备调试
热门文章
- 擦窗机器人毕业设计_自动擦窗机器人系统设计
- mysql rrdtool_RRDtool(开源高性能数据库)
- 转子接地保护原理_转子接地保护调试
- 【无标题】龙迅LT6711GX
- el-date-picker使用中遇到无法点击叉号清除,且不清除状态下无法重新选择时间的问题
- 科猫网:20万元创业,15年后市值170亿,小熊电器的增长逻辑是什么|领军企业
- 不要过分信赖别人,即便是你的影子,也会在黑暗中离开你
- react+antd的Upload 通过证书 上传大文件视频至oss(上传视频)
- 正态分布与威布尔分布图-matplotlib
- 使用Google object_detection API与opencv实现简单的动态目标实时检测