首页 > 其他分享 >【CTS2019】珍珠(计数,容斥,NTT)

【CTS2019】珍珠(计数,容斥,NTT)

时间:2022-10-28 20:45:48浏览次数:59  
标签:frac dfrac sum 容斥 NTT leq CTS2019 序列 binom

题意:称一个长度为 \(n\),值域为 \([1,D]\) 整数序列为合法序列,当且仅当序列中能选出 \(m\) 对数相等。问合法序列数。

\(1\leq D\leq 10^5,1\leq n,m\leq 10^9\)。

题解:

设 \(c_i\) 表示数 \(i\) 在序列中的出现次数,那么限制转化成:

\[\sum_{i=1}^D \lfloor\frac{c_i}{2}\rfloor\geq m \]

把 \(\lfloor\dfrac{c_i}{2}\rfloor\) 转化为 \(\dfrac{c_i-c_i\bmod 2}{2}\),这样做的好处是 \(\sum c_i\) 是固定的,于是限制就能转化成:

\[\sum_{i=1}^D c_i\bmod 2\leq n-2m \]

左式的取值范围是 \([0,D]\),所以若右式 \(n-2m\) 不在这个范围的话,可以直接特判掉,现在只考虑 \(n-2m\in [0,D]\) 的情况。

现在我们只关心 \(c_i\) 的奇偶性。

设 \(f_k\) 表示钦定 \(k\) 个 \(c_i\) 为奇数的方案数,求出 \(f\) 之后恰好为 \(k\) 个的方案数可以卷出来。

对于一种 \(c_i\),强制钦定它为奇数的 EGF 为 \(\{0,1,0,1,\cdots\}=\dfrac{e^x-e^{-x}}{2}\),否则任选的 EGF 为 \(e^x\)。

于是:

\[\begin{aligned} f_k&=\binom{D}{k}[x^n]\left(\dfrac{e^x-e^{-x}}{2}\right)^k\left(e^x\right)^{D-k}\\ &=\binom{D}{k}\frac{1}{2^k}[x^n]\sum_{i=0}^k\binom{k}{i}(-1)^{k-i}e^{(D-2k+2i)x}\\ &=\binom{D}{k}\frac{1}{2^k}\sum_{i=0}^k\binom{k}{i}(-1)^{k-i}\frac{(D-2k+2i)^n}{n!}\\ \end{aligned} \]

卷一下即可。

代码咕了。

标签:frac,dfrac,sum,容斥,NTT,leq,CTS2019,序列,binom
From: https://www.cnblogs.com/ez-lcw/p/16837419.html

相关文章

  • 【CFgym102870J】Junction of Orz Pandas(思维,容斥)
    暴力做法就不会做……考虑容斥,用所有数\(\leqa_i\)的方案减去所有数\(<a_i\)的方案得到最大值为\(a_i\)的方案,\(b_i\)同理容斥,时间复杂度\(O(2^{n+m}nm)\)。直......
  • 【AGC005D】_K Perm Counting(容斥,二分图,计数dp)
    首先正面做不太好做,考虑容斥。设\(f(m)\)表示排列中至少有\(m\)处\(|P_i-i|=k\)的方案数。那么答案就是\(\sum\limits_{i=0}^n(-1)^if(i)\)。原题可以看成一个二......
  • 题解 E. NTT 集美大学第九届程序设计竞赛
    传送门现场没推出了,找了个规律,发现是\((n+1)^{n-1}\)就直接冲过了【分析】考虑\(0\leqk<n\),所以\(\min(k,n-1)=k\)因此有:\(\begin{aligned}&\sum_{i=k}^{\mi......
  • 调色盘 (3维k点最小最远点对-容斥原理)
    调色盘(pastele)题目描述Albus得到了一份礼物:来自Polaris的水彩油墨包。Polaris的油墨包里面有N个颜色,现在Albus打算选其中的K种来作一幅风景画。既然是风景画,颜色就不能太......
  • 容斥定理
    用来求解集合计数问题,求解多个集合并的数目,转化为求交,结果等于加上奇数集合交的数目,将去偶数集合交的数目经典题目求错位排列,反求不是错位排列的条件,在交集合的时候可以合......
  • luogu P5339 [TJOI2019]唱、跳、rap和篮球 (容斥,指数型母函数,NTT)
    https://www.luogu.com.cn/problem/P5339要求不含1234的方案,反过来求含至少一个1234的方案。钦定存在i个位置有1234,位置的方案是Cn-3i,i.其他n-4i个位置的方案是多重集......
  • uoj34 多项式乘法 ntt
    ​​http://www.elijahqi.win/2018/03/17/uoj34ntt/​​​这是一道模板题。给你两个多项式,请输出乘起来后的多项式。输入格式第一行两个整数nn和mm,分别表示两个多项......
  • Luogu2167 Bill的挑战 - 容斥 - dp -
    题目链接:https://www.luogu.com.cn/problem/P2167题解:摘录一段描述容斥题目的话:本题中,关于容斥系数,可以先感性理解一下,严格证明可以用即除了自身,自身的超集都计算......
  • AcWing算法提高课 容斥原理
    容斥原理的复杂度是2^n,一般n不会很大形如:  由于容斥原理一共有2^n中选法,可以用二进制枚举,1表示选择某个条件。然后将偶数个1的状态加起来,奇数个1的状态减去例题:ht......
  • Java获取当前系统事件System.currentTimeMillis()方法 ,获取当前时间戳10位 166529114
    Java获取当前系统事件System.currentTimeMillis()方法,获取当前时间戳10位1665291145转为时间字符串yyy-MM-ddSystem.currentTimeMillis()产生一个当前的毫秒,这个毫秒......