loj3120 「CTS2019 | CTSC2019」珍珠
.... 感觉自己太颓废了。。。。还是来更题解吧。。。【话说写博客会不会涨 rp 啊 qaq ?
题意:
有 n 个物品,每个都有一个 [1,D] 中随机的颜色,相同颜色的两个物品可以配对。现在要求至少能配 m 对,问方案数? $n,m\leq 10^9,D\leq 10^5$
题解:
配对数量 $\geq m \Longleftrightarrow$ 出现奇数次的权值个数 $\leq n-2m$ 。
一个权值出现偶数次的生成函数: $\frac{e^x +e^{-x}}{2}$
一个权值出现奇数次的生成函数: $\frac{e^x - e^{-x}}{2}$
答案的式子:
$$ \begin{aligned} &n!\sum_{k=0}^{n-2m} \left ( \frac{e^x +e^{-x}}{2} + y \frac{e^x - e^{-x}}{2} \right )^D [x^n][y^k]\\ =&n!\left (\frac12 \right )^D \sum_{k=0}^{n-2m} \left (e^x + e^{-x} + y\left ( e^x - e^{-x} \right )\right )^D[x^n][y^k]\\ =&n!\left (\frac12 \right )^D \sum_{k=0}^{n-2m} \left (e^x (1+y) + e^{-x}(1-y) \right )^D[x^n][y^k]\\ =&n!\left (\frac12 \right )^D \sum_{k=0}^{n-2m} \sum_{i=0}^{D}\binom{D}{i} \left (e^{x}(1+y) \right )^i \left( e^{-x}(1-y)\right )^{D-i} [x^n][y^k]\\ =&n!\left (\frac12 \right )^D \sum_{k=0}^{n-2m} \sum_{i=0}^{D}\binom{D}{i} e^{\left (2i-D\right )x} (1+y)^i (1-y)^{D-i} [x^n][y^k]\\ =&\left (\frac12 \right )^D \sum_{i=0}^{D} \binom{D}{i} (2i-D)^n \sum_{k=0}^{n-2m} (1+y)^i (1-y)^{D-i} [y^k] \end{aligned}$$
这就是最后的式子。
考虑怎么对于每个 i 快速求出后面的 $\sum_{k=0}^{n-2m} (1+y)^i (1-y)^{D-i} [y^k]$ ?
可以发现这是一个生成函数的系数前缀和。
即 $\left ( (1+y)^i (1-y)^{D-i} \right )\left (1+y+y^2+y^3+...\right ) [y^{n-2m}]$ ;
即 $\left ( (1+y)^i (1-y)^{D-i} \right )\left (\frac{1}{1-y}\right ) [y^{n-2m}]$ 。
可以发现当 $i\neq D$ 的时候,上式 $= \left ( (1+y)^i (1-y)^{D-i-1} \right )[y^{n-2m}]$ ;
$i=D$ 的时候上式 $=\sum_{k=0}^{n-2m} (1+y)^D [y^k]$ ,由于 $(1+y)^D$ 中 $y$ 不超过 D 次,所以有效的 $k\leq D$ ,暴力计算即可。
那 $i\neq D$ 的时候怎么办呢?我们记 $k=n-2m$ ,$d=D-1$ ,枚举两边对 k 的贡献,即
$$\sum_{j=0}^k \binom{i}{j} \binom{d-i}{k-j}(-1)^{k-j}$$
这种组合数的形式的转化当然是套路的阶乘展开啦——
$$\begin{aligned} \sum_{j=0}^k \frac{i!}{j!(i-j)!} \cdot \frac{(d-i)!}{(k-j)!(d-i-k+j)!} \cdot (-1)^{k-j}\\ =i!(d-i)!\sum_{j=0}^k \frac{(-1)^{k-j}}{j!(k-j)!} \cdot \frac{1}{(i-j)!(d-i-k+j)!}\\ \end{aligned}$$
就是个卷积形式了, fft 即可。
复杂度 $\mathcal{O}(D \log D).$
code:
#include<bits/stdc++.h>
#define rep(i,x,y) for (int i=(x);i<=(y);i++)
#define ll long long using namespace std; const int N=(<<)+,mod=,inv2=(mod+)>>;
int D,n,m,k,d,fac[N],ifac[N],ans,invn,w[N],r[N],a[N],b[N]; int ksm(int x,int y){
int s=; x=(x+mod)%mod;
for (;y;y>>=,x=(ll)x*x%mod) if (y&) s=(ll)s*x%mod;
return s;
} void init(int n){
fac[]=;
rep (i,,n) fac[i]=(ll)fac[i-]*i%mod;
ifac[n]=ksm(fac[n],mod-);
for (int i=n;i;i--) ifac[i-]=(ll)ifac[i]*i%mod;
} int C(int n,int m){
if (n<m) return ;
return (ll)fac[n]*ifac[m]%mod*ifac[n-m]%mod;
} inline void upd(int &x,int y){x+=y; x-=x>=mod?mod:;} void fft_init(int &n){
int len=; while (<<len<n) len++; n=<<len;
int G=ksm(,(mod-)>>len); invn=ksm(n,mod-);
w[]=w[n]=;
rep (i,,n-) w[i]=(ll)w[i-]*G%mod,r[i]=(r[i>>]>>)|((i&)<<(len-));
} void fft(int *a,int all,int fl){
rep (i,,all-) if (i<r[i]) swap(a[i],a[r[i]]);
int n=,m=,x=all>>;
for (;n<=all;m=n,n<<=,x>>=)
for (int i=;i<all;i+=n)
for (int k=;k<m;k++){
int t=(ll)w[fl?all-x*k:x*k]*a[i+m+k]%mod;
a[i+m+k]=(a[i+k]+mod-t)%mod;
a[i+k]=(a[i+k]+t)%mod;
}
if (fl) rep (i,,all-) a[i]=(ll)a[i]*invn%mod;
} int main(){
scanf("%d%d%d",&D,&n,&m),init(D);
k=min(n-*m,D),d=D-;
rep (i,,k) a[i]=(ll)(((k-i)&)?mod-:)%mod*ifac[i]%mod*ifac[k-i]%mod;
rep (i,,d) if (d-k-i>=) b[i]=(ll)ifac[i]%mod*ifac[d-k-i]%mod;
int _n=k+d+; fft_init(_n);
fft(a,_n,),fft(b,_n,);
rep (i,,_n-) a[i]=(ll)a[i]*b[i]%mod;
fft(a,_n,);
rep (i,,d) a[i]=(ll)a[i]*fac[i]%mod*fac[d-i]%mod;
rep (i,,D){
int v=(ll)C(D,i)*ksm(*i-D,n)%mod,tmp=;
if (i==D) rep (i,,k) upd(tmp,C(D,i));
else tmp=a[i];
upd(ans,(ll)tmp*v%mod);
}
printf("%lld\n",(ll)ans*ksm(inv2,D)%mod);
return ;
}
loj3120 「CTS2019 | CTSC2019」珍珠的更多相关文章
- LOJ 3120: 洛谷 P5401: 「CTS2019 | CTSC2019」珍珠
题目传送门:LOJ #3120. 题意简述: 称一个长度为 \(n\),元素取值为 \([1,D]\) 的整数序列是合法的,当且仅当其中能够选出至少 \(m\) 对相同元素(不能重复选出元素). 问合 ...
- 【LOJ】#3120. 「CTS2019 | CTSC2019」珍珠
LOJ3120 52pts \(N - D >= 2M\)或者\(M = 0\)那么就是\(D^{N}\) 只和数字的奇偶性有关,如果有k个奇数,那么必须满足\(N - k >= 2M\) ...
- @loj - 3120@ 「CTS2019 | CTSC2019」珍珠
目录 @description@ @solution@ @accepted code@ @details@ @description@ 有 \(n\) 个在范围 \([1, D]\) 内的整数均匀随机 ...
- Loj #3124. 「CTS2019 | CTSC2019」氪金手游
Loj #3124. 「CTS2019 | CTSC2019」氪金手游 题目描述 小刘同学是一个喜欢氪金手游的男孩子. 他最近迷上了一个新游戏,游戏的内容就是不断地抽卡.现在已知: - 卡池里总共有 ...
- 「CTS2019 | CTSC2019」氪金手游 解题报告
「CTS2019 | CTSC2019」氪金手游 降 智 好 题 ... 考场上签到失败了,没想容斥就只打了20分暴力... 考虑一个事情,你抽中一个度为0的点,相当于把这个点删掉了(当然你也只能抽中 ...
- 「CTS2019 | CTSC2019」随机立方体 解题报告
「CTS2019 | CTSC2019」随机立方体 据说这是签到题,但是我计数学的实在有点差,这里认真说一说. 我们先考虑一些事实 如果我们在位置\((x_0,y_0,z_0)\)钦定了一个极大数\( ...
- LOJ 3119: 洛谷 P5400: 「CTS2019 | CTSC2019」随机立方体
题目传送门:LOJ #3119. 题意简述: 题目说的很清楚了. 题解: 记恰好有 \(i\) 个极大的数的方案数为 \(\mathrm{cnt}[i]\),则答案为 \(\displaystyle\ ...
- LOJ #3119「CTS2019 | CTSC2019」随机立方体 (容斥)
博客链接 里面有个下降幂应该是上升幂 还有个bk的式子省略了k^3 CODE 蛮短的 #include <bits/stdc++.h> using namespace std; const ...
- 【LOJ】#3123. 「CTS2019 | CTSC2019」重复
LOJ3123 60pts 正难则反,熟练转成总方案数减掉每个片段都大于等于s的字典序的方案 按照一般的套路建出kmp上每个点加一个字符的转移边的图(注意这个图开始字母必须是nxt链中下一个相邻的字符 ...
随机推荐
- PowerBuilder学习笔记之根据新PBL文件替换pbd文件的方法
在程序目录上右键,选择Build Runtime Library 在Build Runtime Library 页面点击OK,注意,会将全部的文件重新编译
- jvm--工具
jps (java process status) == ps / top 作用:显示所有运行中的java进程. jstat 作用:查看类装载,内存,垃圾收集,jit编译的信息. jinfo 作用:实 ...
- Codeforces Round #222 (Div. 1) (ABCDE)
377A Maze 大意: 给定棋盘, 保证初始所有白格连通, 求将$k$个白格变为黑格, 使得白格仍然连通. $dfs$回溯时删除即可. #include <iostream> #inc ...
- Consul作为SpringCloud配置中心
一.背景介绍 在分布式系统中动态配置中,可以避免重复重启服务,动态更改服务参数等.一句话非常重要. 另外一篇文章也是这样说的,哈哈. Consul 作为Spring 推荐的分布式调度系统其也具备配置中 ...
- Java8新特性 - 新时间和日期 API
本地时间和时间戳 主要方法: now:静态方法,根据当前时间创建对象 of:静态方法,根据指定日期/时间创建对象 plusDays,plusWeeks,plusMonths,plusYears:向当前 ...
- Java LowerBound
Java LowerBound /** * <html> * <body> * <P> Copyright 1994-2018 JasonInternational ...
- Java单例最好的设计模式:Spring 之 GlobalAdvisorAdapterRegistry
先上代码: public abstract class GlobalAdvisorAdapterRegistry { /** * Keep track of a single instance so ...
- Suricata Rules
Suricata Rules https://redmine.openinfosecfoundation.org/projects/suricata/wiki/Suricata_Rules https ...
- svn提交时把node_modules忽略掉
空白处右键>选中TortoiseSVN>设置(settings)>常规设置(General)>Subversion>编辑(edit)>在弹出的config文件中找g ...
- React Native 开发豆瓣评分(一)环境搭建&配置模拟器
详细可参考 官方文档,这里进记录一些重要过程. 安装环境 下载 Android Studio 选择 Custom 进行安装: Android SDK Android SDK Platform Perf ...