BZOJ.5093.[Lydsy1711月赛]图的价值(NTT 斯特林数)
对于单独一个点,我们枚举它的度数(有多少条边)来计算它的贡献:$$\sum_{i=0}{n-1}ikC_{n-1}i2{\frac{(n-2)(n-1)}{2}}$$
每个点是一样的,所以$$Ans=n\cdot 2{\frac{(n-2)(n-1)}{2}}\sum_{i=0}{n-1}C_{n-1}iik$$
考虑如何计算\(\sum_{i=0}^{n-1}C_{n-1}^ii^k\)。
然后...dalao看到\(i^k\)就想起了第二类斯特林数:
\(S(n,m)\)即在\(m\)个无区别盒子中放\(n\)个不同小球的方案数(要求盒子非空)。
\(S(n,m)\)的一个公式为$$S(n,m)=\frac{1}{m!}\sum_{k=0}m(-1)kC_mk(m-k)n$$
即利用容斥,枚举空盒子至少有多少个。因为盒子无序所以再除以\(m!\)。
而利用反演,或者是组合意义可以得到:$$mn=\sum_{k=0}mC_m^kS(n,k)k!$$
斯特林数中的盒子是无序的所以再乘个\(k!\)。
(\(\sum\)的上界是\(m\)是\(n\)都可以,看需要)
为了方便先令\(n=n-1\)。
我们把\(m^n=\sum_{k=0}^mC_m^kS(n,k)k!\)代进\(Ans\)的\(\sum\)里:$$\sum_{i=0}nC_niik=\sum_{i=0}nC_ni\sum_{j=0}iC_i^j\cdot S(k,j)\cdot j!$$
然后,还是没法做就把\(j\)放到前面枚举试试:$$\sum_{j=0}^nS(k,j)\cdot j!\cdot\sum_{i=j}nC_niC_i^j$$
考虑一下\(\sum_{i=j}^nC_n^iC_i^j\)的组合意义,即从\(n\)个物品中选任意多个(至少\(j\)个),然后从它们中再选出\(j\)个。也就是从\(n\)个中选出\(j\)个后,其余\(n-j\)个任意选的方案数,即\(C_n^j2^{n-j}\)。
所以式子还可以化成:$$\sum_{j=0}^nS(k,j)\cdot j!\cdot C_n^j\cdot 2^{n-j}$$
后面的三项\(j!\cdot C_n^j\cdot 2^{n-j}\)(\(A_n^j\cdot 2^{n-j
}\))都可以直接算,所以我们只要算\(S(k,j)\)就可以了。同BZOJ4555,把上面的$$S(n,m)=\frac{1}{m!}\sum_{k=0}m(-1)kC_mk(m-k)n$$展开成$$S(n,m)=\sum_{k=0}m\frac{(-1)k}{k!}\cdot\frac{(m-k)^n}{(m-k)!}$$
是卷积形式,就可以用\(NTT\)计算了。(\(n<k\)时\(S(n,k)=0\),所以\(j\)枚举到\(\min(n,k)\)就好了)
//11060kb 11548ms
#include <cstdio>
#include <algorithm>
#define G 3
#define invG 332748118
#define inv2 499122177
#define mod 998244353
#define Mod(x) x>=mod&&(x-=mod)
#define Add(x,v) (x+=v)>=mod&&(x-=mod)
#define Mul(x,y) (1ll*(x)*(y)%mod)
typedef long long LL;
const int N=(1<<19)+5;
int fac[N],ifac[N],f[N],g[N],rev[N];
inline int FP(int x,int k)
{
int t=1;
for(; k; k>>=1,x=Mul(x,x))
if(k&1) t=Mul(t,x);
return t;
}
void NTT(int *a,int lim,int opt)
{
for(int i=1; i<lim; ++i) if(i<rev[i]) std::swap(a[i],a[rev[i]]);
for(int i=2; i<=lim; i<<=1)
{
int mid=i>>1,Wn=FP(~opt?G:invG,(mod-1)/i);
for(int j=0; j<lim; j+=i)
for(int k=j,w=1,t; k<j+mid; ++k,w=Mul(w,Wn))
a[k+mid]=a[k]-(t=Mul(w,a[k+mid]))+mod, Mod(a[k+mid]),
a[k]+=t, Mod(a[k]);
}
if(opt==-1) for(int i=0,inv=FP(lim,mod-2); i<lim; ++i) a[i]=Mul(a[i],inv);
}
int main()
{
int n,K; scanf("%d%d",&n,&K); --n;//!
int m=std::min(n,K);
fac[0]=fac[1]=1, ifac[0]=ifac[1]=1;
for(int i=2; i<=m; ++i) fac[i]=Mul(fac[i-1],i);
ifac[m]=FP(fac[m],mod-2);
for(int i=m; i; --i) ifac[i-1]=Mul(ifac[i],i);
for(int i=0; i<=m; ++i) f[i]=i&1?mod-ifac[i]:ifac[i], g[i]=Mul(FP(i,K),ifac[i]);//x^K/x! (n=K)
int lim=1,l=-1;
while(lim<=m+m) lim<<=1,++l;
for(int i=1; i<lim; ++i) rev[i]=(rev[i>>1]>>1)|((i&1)<<l);
NTT(f,lim,1), NTT(g,lim,1);
for(int i=0; i<lim; ++i) f[i]=Mul(f[i],g[i]);
NTT(f,lim,-1);
LL ans=0; int pw2=FP(2,n),A=1;//2^{n-i} A(n,i)
for(int i=0; i<=m; ++i) ans+=1ll*f[i]*A%mod*pw2%mod, pw2=Mul(pw2,inv2), A=Mul(A,n-i);
printf("%lld\n",ans%mod*(n+1)%mod*FP(2,(1ll*n*(n-1)>>1ll)%(mod-1))%mod);
return 0;
}
BZOJ.5093.[Lydsy1711月赛]图的价值(NTT 斯特林数)的更多相关文章
- bzoj 5093 [Lydsy1711月赛]图的价值 NTT+第二类斯特林数
[Lydsy1711月赛]图的价值 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 245 Solved: 128[Submit][Status][D ...
- BZOJ 5093: [Lydsy1711月赛]图的价值 第二类斯特林数+NTT
定义有向图的价值为图中每一个点的度数的 \(k\) 次方之和. 求:对于 \(n\) 个点的无向图所有可能情况的图的价值之和. 遇到这种题,八成是每个点单独算贡献,然后累加起来. 我们可以枚举一个点的 ...
- BZOJ 5093: [Lydsy1711月赛]图的价值
第二类斯特林数模版题 需要一些组合数的小$ trick$ upd:这里更新了本题巧妙的$ O(k)$做法,虽然常数很大就是了 传送门:here 题意:求所有$ n$个节点的无重边自环图的价值和,定义一 ...
- bzoj 5093 [Lydsy1711月赛]图的价值——第二类斯特林数
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=5093 不要见到组合数就拆! 枚举每个点的度数,则答案为 \( n*\sum\limits_{ ...
- 【bzoj5093】 [Lydsy1711月赛]图的价值 组合数+斯特林数+NTT
Description "简单无向图"是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向 ...
- BZOJ 5093[Lydsy1711月赛]图的价值 线性做法
博主曾更过一篇复杂度为$O( k· \log k)$的多项式做法在这里 惊闻本题有$ O(k)$的神仙做法,说起神仙我就想起了于是就去学习了一波 幂与第二类斯特林数 推导看这里 $$ x^k=\sum ...
- 【bzoj5093】[Lydsy1711月赛]图的价值(NTT+第二类斯特林数)
题意: 给定\(n\)个点,一个图的价值定义为所有点的度数的\(k\)次方之和. 现在计算所有\(n\)个点的简单无向图的价值之和. 思路: 将式子列出来: \[ \sum_{i=1}^n\sum_{ ...
- BZOJ5093 [Lydsy1711月赛]图的价值 【第二类斯特林数 + NTT】
题目链接 BZOJ5093 题解 点之间是没有区别的,所以我们可以计算出一个点的所有贡献,然后乘上\(n\) 一个点可能向剩余的\(n - 1\)个点连边,那么就有 \[ans = 2^{{n - 1 ...
- bzoj5093:[Lydsy1711月赛]图的价值
题目 首先考虑到这是一张有标号的图,每一个点的地位是相等的,因此我们只需要求出一个点的价值和乘上\(n\)就好了 考虑一个点有多少种情况下度数为\(i\) 显然我们可以让除了这个点的剩下的\(n-1\ ...
随机推荐
- Gson解决字段为空是报错的问题
json解析有很多工具,这里说的是最常用也是解析速度最快的Gson,Gson是google家出的,有一个缺点就是无法设置null替换, 我们只能手动的批量替换服务器返回的null了,正常的接口定义的时 ...
- C++ Primer 笔记——固有的不可移植的特性
1.为了支持底层编程,C++定义了一些固有的不可移植的特性,所谓不可移植特性是指因机器而异的特性. 2.一个位域中含有一定数量的二进制位,位域在内存中的布局是机器相关的.位域的类型必须是整型或枚举类型 ...
- windows 7 下用git
参考:http://my.oschina.net/longxuu/blog/141699
- Paddington2
- zeromq的安装,部署(号称最快的消息队列,消息中间件)
1:Storm作为一个实时处理的框架,产生的消息需要快速的进行处理,比如存在消息队列ZeroMQ里面. 由于消息队列ZeroMQ是C++写的,而我们的程序是运行在JVM虚拟机里面的.所以需要jzmq这 ...
- CallContext
1.线程本地存储区的专用集合对象,并提供对每个逻辑执行线程都唯一的数据槽.2.数据槽不在其他逻辑线程上的调用上下文之间共享. class Program { static Jason_TestEnti ...
- 错误 java.lang.ClassCastException: com.xx cannot be cast to ResourceBundle
出现错误: java.lang.ClassCastException: com.xxx cannot be cast to ResourceBundle 百度搜索错误,没有结果.谷歌搜索:http:/ ...
- react-native-video
<Video ref='videoPlayer' source={{uri:'http://www.thinktown.com/video/th.mp4'}} //网络视频 //source={ ...
- cf666 C. Codeword 组合数学
题解: 首先暴力很显然 f[i][j]表示到第i个位置,串匹配到j 这样每次是n^2的 我们假设每个位置匹配的第一个位置 然后从这个到上一个位置一定不能等于这个串的值 ans=simga{i,C(i- ...
- [转]Maven与nexus关系
开始在使用Maven时,总是会听到nexus这个词,一会儿maven,一会儿nexus,当时很是困惑,nexus是什么呢,为什么它总是和maven一起被提到呢? 我们一步一步来了解吧. 一.了解Mav ...