【SDOI 2017】龙与地下城(组合)
概率论太难了,不会。但这不能阻止我们过题。
相信大家都会一个基于背包的暴力做法,我们可以将其看成是卷积的形式就可以用fft优化了。形式化讲,就是求幂级数$ (\sum\limits_{i = 0}^{x - 1} \frac{1}{x} z^i)^y $在$[z^A, z^B]$之间的系数和。
不在模意义下的做法
直接将上述幂级数暴力倍增卷积求出来复杂度是$O(x*y*log(xy))$,不太能过的。但如果不在模意义下做我们就可以尝试爆精度爆过去。很容易发现最后求得的点数和很大概率就是在均值附近的,过大过小的概率几乎为0。也就是我们中间在做卷积和有很多项几乎为0,我们把他们忽略掉是在精度的承受范围之内的。于是我们在做倍增求出上述幂级数的时候我们可以每次只保留中间的一小段有值的部分,其余的扔掉就行了。也就是说我们设定一个阈值$\epsilon$,每次卷积后把多项式中值$\le \epsilon$的都扔掉。这样复杂度就是$O(len*log(len + y))$,其中$len$是最终幂级数中值$> \epsilon$的个数。实践证明,当$\epsilon$取$10^{-9}$时,$len$大概是两三万,并且此时的精度可以达到小数点后4位。
在模意义下的做法
如果答案可以对某个质数取模,这题就更好做了。
考虑最开始的那个幂级数:
$(\sum\limits_{i = 0}^{x - 1} \frac{1}{x} z^i)^y = (\frac{1}{x} \frac{1 - z^x}{1 - z})^y = (\frac{1}{x})^y (1 - z^x)^y(\frac{1}{1 - z})^y$
我们想要求其在第$[A, B]$之间的系数和,我们可以乘上一个$\frac{1}{1 - x}$来做一次前缀和,这样我们可以转化成两次形如求一个$z^n$的系数的问题。也就是:
$[z^n] (\frac{1}{x})^y (1 - z^x)^y(\frac{1}{1 - z})^{y + 1}$
手动展开后面的两个:
$=[z^n] (\frac{1}{x})^y(\sum\limits_{i = 0}^{y}\binom{y}{i}(-1)^iz^{ix})(\sum\limits_{i = 0}\binom{y + i}{y}z^i)$
$=(\frac{1}{x})^y\sum\limits_{i = 0}^{y}\binom{y}{i}(-1)^i\binom{y + n - ix}{y}$
预处理组和数之后随便做一下就好了,复杂度$O(x*y)$。
$\bigodot$套路与技巧:
- 隔板法用于展开形如$(\frac{1}{1 - z})^y$的幂级数。
做法一:
#include <bits/stdc++.h> using namespace std; namespace PO {
const int N = 4e5 + ;
const double PI = acos(-); struct Com {
double x, y;
friend Com operator + (Com a, Com b) {
return (Com){ a.x + b.x, a.y + b.y };
}
friend Com operator - (Com a, Com b) {
return (Com){ a.x - b.x, a.y - b.y };
}
friend Com operator * (Com a, Com b) {
return (Com){ a.x * b.x - a.y * b.y, a.x * b.y + a.y * b.x };
}
} ta[N], tb[N], w[N];
int rev[N], L = ;
typedef vector<Com> Poly; void Init(int l) {
for (; L < l; L <<= ) {
for (int i = L; i < (L << ); ++i) {
w[i] = (Com){ cos(PI / L * (i - L)), sin(PI / L * (i - L)) };
}
}
for (int i = ; i < l; ++i) {
rev[i] = (rev[i >> ] >> ) | (i & ? l >> : );
}
} void Dft(Com *a, int l) {
for (int i = ; i < l; ++i)
if (i < rev[i]) swap(a[i], a[rev[i]]);
for (int i = ; i < l; i <<= ) {
for (int j = ; j < l; j += i << ) {
Com *l = a + j, *r = l + i, *wx = w + i, y;
for (int k = ; k < i; ++k, ++l, ++r, ++wx) {
y = (*r) * (*wx);
*r = (*l) - y;
*l = (*l) + y;
}
}
}
}
void Idft(Com *a, int l) {
reverse(a + , a + l);
Dft(a, l);
for (int i = ; i < l; ++i) {
a[i].x /= l;
a[i].y /= l;
}
} vector<double> Mul(vector<double> a, vector<double> b) {
int n = a.size(), m = b.size(), l;
for (l = ; l < n + m - ; l <<= );
Init(l);
a.resize(l), b.resize(l);
for (int i = ; i < l; ++i) {
ta[i] = (Com){ a[i], };
tb[i] = (Com){ b[i], };
}
Dft(ta, l), Dft(tb, l);
for (int i = ; i < l; ++i) {
ta[i] = ta[i] * tb[i];
}
Idft(ta, l);
for (int i = ; i < l; ++i) {
a[i] = ta[i].x;
}
a.resize(n + m - );
return a;
} } void Cut(vector<double> &a, int &base) {
const double GAMA = 1e-;
int id = ;
while (a[id] < GAMA) ++id;
for (int i = ; i + id < a.size(); ++i) {
a[i] = a[i + id];
}
while (a.back() < GAMA) {
a.pop_back();
}
base += id;
} int main() {
int tc;
for (cin >> tc; tc--; ) {
int x, y;
cin >> x >> y; vector<double> ans(x * y);
vector<double> F(x, 1.0 / x), G(, );
int base_f = , base_g = ;
for (int ex = y; ex; ex >>= ) {
if (ex & ) {
G = PO::Mul(G, F);
base_g += base_f;
Cut(G, base_g);
}
F = PO::Mul(F, F);
base_f *= ;
Cut(F, base_f);
}
for (int i = ; i < G.size(); ++i) {
ans[i + base_g] = G[i];
}
for (int i = ; i < x * y; ++i) {
ans[i] += ans[i - ];
} for (int cas = ; cas--; ) {
int l, r;
cin >> l >> r;
printf("%.12f\n", ans[r] - (l? ans[l - ] : ));
}
} return ;
}
【SDOI 2017】龙与地下城(组合)的更多相关文章
- SDOI 2017 Day1
日期:2017-04-10 题解: 第一题: 题目大意:求fi(gcd(i,j))的乘积 i,j属于[1,1e6],数据组数1000组. 类别:套路题. 第二题:BZOJ原题. 题解:LCT套线段树 ...
- [SDOI 2017]新生舞会
Description 题库链接 给你个 \(2\times N\) 的带权二分图,两个权值 \(a,b\) ,让你做匹配使得 \[\frac{\sum a}{\sum b}\] 最大. \(1\le ...
- [SDOI 2017]数字表格
Description 题库链接 记 \(f_i\) 为 \(fibonacci\) 数列的第 \(i\) 项. 求 \[\prod_{i=1}^n\prod_{j=1}^mf_{gcd(i,j)}\ ...
- [BZOJ 4817] [SDOI 2017] 树点涂色
Description Bob有一棵 \(n\) 个点的有根树,其中 \(1\) 号点是根节点.Bob在每个点上涂了颜色,并且每个点上的颜色不同. 定义一条路径的权值是:这条路径上的点(包括起点和终点 ...
- [BZOJ 4819] [SDOI 2017] 新生舞会
Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴. 有 \(n\) 个男生和 \(n\) 个女生参加舞会买一个男生和一个女生一起跳舞,互为舞伴. C ...
- [BZOJ 4818] [SDOI 2017] 序列计数
Description Alice想要得到一个长度为 \(n\) 的序列,序列中的数都是不超过 \(m\) 的正整数,而且这 \(n\) 个数的和是 \(p\) 的倍数. Alice还希望,这 \(n ...
- 解题:SDOI 2017 硬币游戏
题面 板板的生成函数做法太神仙了,我跑了 朴素的做法是建立AC自动机变成图上的随机游走问题 来仔细考虑一下转移,把状态分成非结尾状态和结尾状态.在一个非结尾状态后补一个串是一定能到达目标串的,但是如果 ...
- BZOJ.4909.[SDOI2017]龙与地下城(正态分布 中心极限定理 FFT Simpson积分)
BZOJ 洛谷 https://www.luogu.org/blog/ShadowassIIXVIIIIV/solution-p3779# 正态分布 正态分布是随机变量\(X\)的一种概率分布形式.它 ...
- SDOI 2017 天才黑客
/* 根据claris的 博客以及 beginend 的博客来写的 首先考虑如何求出最短路 可以从样例看出 路径是从边走到边的, 所以我们将边看作点 有共同端点的两边之间 互相连边, 边权为lcp. ...
随机推荐
- 2018 Multi-University Training Contest 2 部分简单题解析
Preface 多校第二场,依靠罚时优势打到了校内的Rank 2 暴力分块碾标算系列 T4 Game 题目大意:在一个数集\([1,n]\)中两个人轮流选择其中的一个数,并从数集中删去这个数所有约数. ...
- asp.net core部署时自定义监听端口,提高部署的灵活性
另一种方式 https://www.cnblogs.com/stulzq/p/9039836.html 代码截图: 贴一下代码,方便复制: //默认端口号5000 string port = &quo ...
- Python零基础入门(安装步骤,验证方式, 简单操作)
本篇文章适合新人小白初步了解Python,涵盖Python的介绍.安装以及简单的基础操作. 1.Python简介 Python 是一个高层次的结合了解释性.编译性.互动性和面向对象的脚本语言.它的设 ...
- Centos下MooseFS(MFS)分布式存储共享环境部署记录
分布式文件系统(Distributed File System)是指文件系统管理的物理存储资源不一定直接连接在本地节点上,而是通过计算机网络与节点相连,分布式文件系统的实际基于客户机/服务器模式.目前 ...
- react/React Native 在 import 导入时,有的带花括号{},有的不带原理解析
在使用import引用模块时,如何正确使用{} 例如:有两个文件,home.js.user.js 一:不使用{}: 当需要在home.js中引入user.js的时候 //home.js 文件中impo ...
- 关于QQ的NABCD模型
关于QQ的NABCD模型 N--Need 随着电脑的普及,人们在网络上进行交流的时间越来越多,由于现有的交流工具还不是那么的完善,还不能够完全满足人们在交流时的需求.因此为了满足人们更多的需求,我们设 ...
- 小学生四则运算APP核心代码公布
Mainactivity类: package com.example.XXSCYS; import java.io.ByteArrayOutputStream; import java.io.File ...
- Ajax cross domain
xhrFields:{ withCredentials:true}, https://stackoverflow.com/questions/2054316/sending-credentials-w ...
- MySQL服务器监控注意事项
当开发,测试,生产的数据库环境配置不一致(比如:配置字符集不同)时而导致特殊现象时,可以从Navicat的<工具>-<服务器监控>-<变量>里查找原因.
- eclipse webproject activiti
https://stackoverflow.com/questions/42858723/activiti-eclipse-maven-project-to-dynamic-web-project-a ...