题面

CTS2019 珍珠

有 \(n\) 个在 \([1,d]\) 内的整数,求使可以拿出 \(2m\) 个整数凑成 \(m\) 个相等的整数对的方案数。

数据范围:\(0\le m\le 10^9\),\(1\le n\le 10^9\),\(1\le d\le 10^5\)。


蒟蒻语

非常巧妙的题,主要要用到二项式反演、指数级生成函数和 NTT

做个广告,这是我读过最好的生成函数讲解:link


蒟蒻解

设 \(c_i\) 表示 \(i\) 这个数的出现次数。

设 \(odd=\sum [c_i\in {\rm odd} ]\),即 \(c_i\) 奇数个数。

很明显最多能凑成 \(\frac{n-odd}{2}\) 对,按题意:

\[\begin{aligned}
\frac{n-odd}{2}&\ge m\\
odd&\le n-2m
\end{aligned}
\]

这里有两个特判,如果 \(n-2m<0\) 答案是 \(0\),如果 \(n-2m\ge d\) 答案是 \(d^n\)。

设 \(g(i)\) 表示 \(odd=i\) 的方案数。

设 \(f(i)\) 表示钦定 \(i\) 个 \(c_i\) 是奇数,剩下随意的方案数(不是 \(odd\ge i\) 的方案数,这里对一些排列会重复统计,但是反演完就没事了)。

\[f(i)=\sum_{x=i}^d{x\choose i}g(x)\Longleftrightarrow g(i)=\sum_{x=i}^d(-1)^{x-i}{x\choose i} f(x)
\]

所以可以先求 \(f(i)\),用到了指数级生成函数,中间把每个 \(e\) 的幂次项展开,最后归成卷积形式:

\[\begin{aligned}
f(i)=&{d \choose i} \left(\frac{e^x-e^{-x}}{2}\right)^i e^{(d-i)x} n![n]\\
=&{d \choose i} (e^x-e^{-x})^i e^{(d-i)x} \frac{n!}{2^i}[n]\\
=&{d \choose i} e^{(d-i)x} \sum_{j=0}^i(-1)^{i-j} {i\choose j}e^{(j-(i-j))x} \frac{n!}{2^i}[n]\\
=&{d \choose i} e^{(d-i)x} \sum_{j=0}^i(-1)^{i-j} {i\choose j}e^{(2j-i)x} \frac{n!}{2^i}[n]\\
=&{d \choose i} \frac{n!}{2^i}\sum_{j=0}^i(-1)^{j} {i\choose j}e^{(d-2j)x} [n]\\
=&{d \choose i} \frac{n!}{2^i}\sum_{j=0}^i(-1)^{j} {i\choose j} \frac{(d-2j)^n}{n!}\\
=&\frac{d!}{i!(d-i)!2^i}\sum_{j=0}^i(-1)^{j} \frac{i!}{j!(i-j)!} (d-2j)^n\\
=&\frac{d!}{(d-i)!2^i}\sum_{j=0}^i\frac{(-1)^{j}(d-2j)^n}{j!}\cdot \frac{1}{(i-j)!} \\
\end{aligned}
\]

最后一个难点是如何求:

\[g(i)=\sum_{x=i}^d(-1)^{x-i}{x\choose i} f(x)
\]

感觉可以凑成卷积形式,但总差一点。尝试把 \(f\) 和 \(g\) 都反过来,即令 \(f'(x)=f(d-x)\),\(g'(x)=g(d-x)\):

\[\begin{aligned}
g(i)=&\sum_{x=i}^d(-1)^{x-i}{x\choose i} f(x)\\
=&\sum_{x=i}^d(-1)^{x-i}{x\choose i} f'(d-x)\\
=&\sum_{x=0}^{d-i}(-1)^{d-x-i}{d-x\choose i} f'(x)\\
\end{aligned}\\
g'(d-i)=\sum_{x=0}^{d-i}(-1)^{d-x-i}{d-x\choose i} f'(x)\\
\begin{aligned}
g'(i)=&\sum_{x=0}^{i}(-1)^{i-x}{d-x\choose d-i} f'(x)\\
=&\sum_{x=0}^{i}(-1)^{i-x}\frac{(d-x)!}{(d-i)!(i-x)!} f'(x)\\
=&\frac{1}{(d-i)!}\sum_{x=0}^{i} (-1)^{i-x}\frac{1}{(i-x)!}\cdot f'(x)(d-x)!\\
\end{aligned}\\
\]

然后就可以求 \(g(i)\) 了。

答案便是 \(\sum_{i=0}^{n-2m} g(i)\)。


代码

#include <bits/stdc++.h>
using namespace std; //Start
typedef long long ll;
typedef double db;
#define mp(a,b) make_pair((a),(b))
#define x first
#define y second
#define bg begin()
#define ed end()
#define sz(a) int((a).size())
#define pb(a) push_back(a)
#define R(i,a,b) for(int i=(a),i##E=(b);i<i##E;i++)
#define L(i,a,b) for(int i=(b)-1,i##E=(a)-1;i>i##E;i--)
const int iinf=0x3f3f3f3f;
const ll linf=0x3f3f3f3f3f3f3f3f; //Data
const int N=1e5;
int n,m,d,ans; //Math
const int mod=998244353;
int Pow(int a,int x){
int res=1; for(;x;x>>=1,a=1ll*a*a%mod)
if(x&1) res=1ll*res*a%mod; return res;
}
const int mN=N+1;
int fac[mN],ifac[mN];
void math_init(){
fac[0]=1;R(i,1,d+1) fac[i]=1ll*fac[i-1]*i%mod;
ifac[d]=Pow(fac[d],mod-2);
L(i,0,d) ifac[i]=1ll*ifac[i+1]*(i+1)%mod;
} //Poly
const int pN=mN<<2;
int f[pN],g[pN];
const int G=3,iG=Pow(3,mod-2);
int rev[pN],pn;
void poly_init(){
pn=1<<int(ceil(log2(d*2+2)));
R(i,0,pn) rev[i]=(rev[i>>1]>>1)|((i&1)*(pn>>1));
}
void NTT(int* a,int t){
R(i,0,pn)if(i<rev[i]) swap(a[i],a[rev[i]]);
for(int mid=1;mid<pn;mid<<=1){
int wn=Pow(~t?G:iG,(mod-1)/(mid<<1));
for(int i=0;i<pn;i+=(mid<<1)){
int w=1;
R(j,i,mid+i){
int x=a[j],y=1ll*a[mid+j]*w%mod;
a[j]=(x+y)%mod,a[mid+j]=(x-y+mod)%mod,w=1ll*w*wn%mod;
}
}
}
if(!~t){
int in=Pow(pn,mod-2);
R(i,0,pn) a[i]=1ll*a[i]*in%mod;
}
} //Main
int main(){
ios::sync_with_stdio(0);
cin.tie(0),cout.tie(0);
cin>>d>>n>>m;
if(n-2*m<0) return cout<<0<<'\n',0;
if(d<=n-2*m) return cout<<Pow(d,n)<<'\n',0;
math_init(),poly_init();
R(i,0,d+1){
f[i]=1ll*Pow((d-2*i+mod)%mod,n)*ifac[i]%mod;
if(i&1) f[i]=(mod-f[i])%mod;
g[i]=ifac[i];
}
R(i,d+1,pn) f[i]=g[i]=0;
NTT(f,1),NTT(g,1);
R(i,0,pn) f[i]=1ll*f[i]*g[i]%mod;
NTT(f,-1);
R(i,d+1,pn) f[i]=g[i]=0;
R(i,0,d+1) f[i]=1ll*f[i]*fac[d]%mod*ifac[d-i]%mod*Pow(2,mod-1-i)%mod;
reverse(f,f+d+1);
R(i,0,d+1) f[i]=1ll*f[i]*fac[d-i]%mod;
R(i,0,d+1){
g[i]=ifac[i];
if(i&1) g[i]=(mod-g[i])%mod;
}
R(i,d+1,pn) f[i]=g[i]=0;
NTT(f,1),NTT(g,1);
R(i,0,pn) g[i]=1ll*f[i]*g[i]%mod;
NTT(g,-1);
R(i,d+1,pn) f[i]=g[i]=0;
R(i,0,d+1) g[i]=1ll*g[i]*ifac[d-i]%mod;
reverse(g,g+d+1);
R(i,0,n-2*m+1) (ans+=g[i])%=mod;
cout<<ans<<'\n';
return 0;
}

祝大家学习愉快!

题解-CTS2019 珍珠的更多相关文章

  1. 【题解】CTS2019珍珠(二项式反演+卷积)

    [题解]CTS2019珍珠 题目就是要满足这样一个条件\(c_i\)代表出现次数 \[ \sum {[\dfrac {c_i } 2]} \ge 2m \] 显然\(\sum c_i=n\)所以,而且 ...

  2. [CTS2019]珍珠——二项式反演

    [CTS2019]珍珠 考虑实际上,统计多少种染色方案,使得出现次数为奇数的颜色数<=n-2*m 其实看起来很像生成函数了 n很大?感觉生成函数会比较整齐,考虑生成函数能否把n放到数值的位置,而 ...

  3. 题解 P5401 [CTS2019]珍珠

    蒟蒻语 这题太玄学了,蒟蒻写篇题解来让之后复习 = = 蒟蒻解 假设第 \(i\) 个颜色有 \(cnt_i\) 个珍珠. \(\sum\limits_{i=1}^{n} \left\lfloor\f ...

  4. Luogu5401 CTS2019珍珠(生成函数+容斥原理+NTT)

    显然相当于求有不超过n-2m种颜色出现奇数次的方案数.由于相当于是对各种颜色选定出现次数后有序排列,可以考虑EGF. 容易构造出EGF(ex-e-x)/2=Σx2k+1/(2k+1)!,即表示该颜色只 ...

  5. LOJ3120 CTS2019 珍珠 生成函数、二项式反演、NTT

    传送门 题目大意:给出一个长度为\(n\)的序列\(a_i\),序列中每一个数可以取\(1\)到\(D\)中的所有数.问共有多少个序列满足:设\(p_i\)表示第\(i\)个数在序列中出现的次数,\( ...

  6. [LOJ#3120][Luogu5401][CTS2019]珍珠(容斥+生成函数)

    https://www.luogu.org/blog/user50971/solution-p5401 #include<cstdio> #include<algorithm> ...

  7. 题解-CTS2019氪金手游

    Problem \(\mathtt {loj-3124}\) 题意概要:给定 \(n\) 个点,\(w_i\) 分别有 \(p_{i,1},p_{i,2},p_{i,3}\) 的概率取 \(1,2,3 ...

  8. 题解-CTS2019随机立方体

    problem \(\mathtt {loj-3119}\) 题意概要:一个 \(n\times m\times l\) 的立方体,立方体中每个格子上都有一个数,如果某个格子上的数比三维坐标中至少有一 ...

  9. [CTS2019]珍珠(NTT+生成函数+组合计数+容斥)

    这题72分做法挺显然的(也是我VP的分): 对于n,D<=5000的数据,可以记录f[i][j]表示到第i次随机有j个数字未匹配的方案,直接O(nD)的DP转移即可. 对于D<=300的数 ...

随机推荐

  1. 从头学起Verilog(三):Verilog逻辑设计

    引言 经过了组合逻辑和时序逻辑的复习,终于到了Verilog部分.这里主要介绍Verilog一些基础内容,包括结构化模型.TestBench编写和仿真.真值表模型. 这部分内容不多,也都十分基础,大家 ...

  2. 使用日志系统graylog获取Ceph集群状态

    前言 在看集群的配置文件的时候看到ceph里面有一个graylog的输出选择,目前看到的是可以收集mon日志和clog,osd单个的日志没有看到,Elasticsearch有整套的日志收集系统,可以很 ...

  3. [web安全原理分析]-XEE漏洞入门

    前言 1 前言 XXE漏洞 XXE漏洞全称(XML External Entity Injection)即xml外部实体注入漏洞,XXE漏洞发生在应用程序解析XML输入时,没有禁止外部实体的加载,导致 ...

  4. 从执行上下文角度重新理解.NET(Core)的多线程编程[1]:基于调用链的”参数”传递

    线程是操作系统能够进行运算调度的最小单位,操作系统线程进一步被封装成托管的Thread对象,手工创建并管理Thread对象已经成为了所能做到的对线程最细粒度的控制了.后来我们有了ThreadPool, ...

  5. 类虚拟机软件CrossOver是什么?它的优势在哪里?

    虚拟机软件对于很多人来说已经不是一个陌生的词汇了.我们可以通过软件来模拟具有完整硬件系统功能的计算机系统.比如我们可以在Mac OS系统上模拟Windows 7 的系统,以此来安装我们想要使用的应用程 ...

  6. FL Studio入门:如何使用Layer插件叠加音色

    Layer控制器也是FL Studio中一个特别有用的插件,主要用来叠加音色,以及通过Layer通道来控制多个打击乐通道. 下面我们一起来看看叠加音色是怎么做出来的. 1.新建一个空白工程,插入3个3 ...

  7. SpringSecurity之整合JWT

    SpringSecurity之整合JWT 目录 SpringSecurity之整合JWT 1. 写在前面的话 2. JWT依赖以及工具类的编写 3. JWT过滤器 4. 登录成功结果处理器 5. Sp ...

  8. CET-6备考丨词组、佳句积累

    一.片段积累 (2020/09/17) 近三年 CET6 - 翻译 Phrases or Expressions Notes drive to commute to and from work 开车上 ...

  9. Java基础教程——线程通信

    线程通信:等待.唤醒 Object方法 这些方法在拥有资源时才能调用 notify 唤醒某个线程.唤醒后不是立马执行,而是等CPU分配 wait 等待,释放锁,不占用CPU资源 notifyAll 唤 ...

  10. 顺序结构(C语言基本结构)

    顺序结构 1.基本概念 语句执行的顺序与顺序程序书写的顺序一致 特点 a.程序执行的顺序和语句书写的顺序一致 b.有一个数据入口,一个数据出口 顺序结构与四则运算 顺序结构是C语言的基本结构 程序由上 ...