Description

“简单无向图”是指无重边、无自环的无向图(不一定连通)。 一个带标号的图的价值定义为每个点度数的k次方的和。 给定n和k,请计算所有n个点的带标号的简单无向图的价值之和。

因为答案很大,请对998244353取模输出。$n \le 10^9,k \le 200000$

化学学考时含义推式子+手动打表找规律得到了一个$O(nlogn)$的式子开心的很我以为我要AC了回来看数据范围就升天了。

问NC大神这题用到了什么:斯特林数/伯努利数。然后就自闭了学了一天的知识点还去做了点伯努利数的别的题。

然而最后没想到伯努利数的做法,是用第二类斯特林数解决的。

先列出来答案的式子,分别考虑每个点的贡献:

$ans=2^{\frac{(n-1)(n-2)}{2}} \times n \times \sum\limits_{i=0}^{n} C_{n-1}^i \times i^k$

就是首先在n个点里选出你当前考虑贡献答案的是哪个点,然后在剩余的n-1个点的子图里联通情况随意就是$2^{\frac{(n-1)(n-2)}{2}}$

然后从n-1个点里选出i个与当前贡献答案的点连边,选法是组合数,贡献的价值是$i^k$。

可以发现,这里有个求自然数幂和乘个组合数的形式,而自然数幂和的范围很大不能直接求解。

动用学过的(和刚学的)所有数学知识点看有什么能把幂拆成可以接受的形式。

1,伯努利数:$\sum\limits_{i=1}^{n} i^m = \frac{\sum\limits_{i=1}^{m} B_i \times C_{m+1}^{i} \times n^{m+1-i} }{m+1}$

这不妨是个好主意,但是这道题里要求的是自然数幂乘组合数后再求和,并不能直接适用,所以至少我现在还不会做。

等会上网学习一下再来fix这里。

学习完了,网上并没有这样的做法,但是伯努利数和第二类斯特林数可以互相转化。。。

对这题并没有什么帮助。。。

2,第二类斯特林数:$x^n=\sum\limits_{i=0}^{n} \left \{ ^n_k \right \} \times \frac{x!}{(x-i)!}$

这样的话虽然不能直接求和,但是好歹可以把难以解决的高次幂转化为经常处理的阶乘什么的了。

首先说说关于斯特林数怎么求:先弄明白含义,第二类斯特林数$\left \{ ^n_k \right \}$表示把n个元素放进k个无标号非空集合的方案数。

然后我们再弄几个类似的定义:设f(n,k)表示把n个元素放进k个有编号非空集合的方案数,g(n,k)表示n个元素放进k个有标号可空集合的方案数。

对于是否有编号的问题,其实就是多了个内部顺序,所以$\left \{ ^n_k \right \}=\frac{f(n,k)}{k!}$

然后g其实是可以直接求解的:$g(n,k)=k^n$

然后至于f和g之间的关系,我们可以通过枚举它具体有几个集合是空的,得到$g(n,k)=\sum\limits_{i=0}^{k} C_k^i \times f(n,i)$

就是选出到底哪几个不是空的,在这些非空的里面进行分配。

然后发现这个形式就是二项式反演的基本形式,所以反演它:$f(n,k)=\sum\limits_{i=0}^{k} (-1)^{k-i} \times C_k^i \times g(n,i)$

这样我们就能得到f,进而求解出第二类斯特林数。

总的式子就是$\left \{ ^n_k \right \}=\sum\limits_{i=0}^{k} \frac{(-1)^{k-i} \times i^k}{(k-i)! \times i!}$

这样的话单点求第二类斯特林数就是$O(n)$的了。而上式是个卷积形式,所以求一行第二类斯特林数就可以用NTT在$O(nlogn)$之内解决。

回到这道题来,求解斯特林数的复杂度是$O(klogk)$可以接受,那么现在答案的式子是什么?

$ans=2^{\frac{(n-1)(n-2)}{2}} \times n \times \sum\limits_{i=0}^{n} C_{n-1}^{i} \sum\limits_{j=0}^{k} \left \{ ^k_j \right \} \times \frac{i!}{(i-j)!}$

把各种组合数拆一拆放在一起消一消,得到:$ans=2^{\frac{(n-1)(n-2)}{2}} \times n! \times \sum\limits_{i=0}^{n} \sum\limits_{j=0}^{k} \left \{ ^k_j \right \} \times \frac{1}{(i-j)! \times (n-i-1)!}$

换一下ij的枚举顺序,发现最后那个分母上两个阶乘的底之和是一个与i无关的定值,上下都乘一个$(n-j-1)!$合成一个组合数。它等于$\frac{C_{n-j-1}^{i-j}}{(n-j-1)!}$

然后再套上外面那个求和。。。一行组合数求和得到的是2的次幂,所以这个和式就可以扔掉了得到下面形式:

$ans=2^{\frac{(n-1)(n-2)}{2}} \times n! \times \sum\limits_{j=0}^{k} \frac{\left \{ ^k_j \right \} \times 2^{n-j-1}}{(n-j-1)!} $

现在的瓶颈就在于式子中的两个阶乘不能在合理复杂度内求出。

然而它们一个在分子上,一个在分母上,而且它们底数的差值在k范围内。。。

抵消啊!就剩下一个下降幂了!

所以最后可做的式子就是$2^{\frac{(n-1)(n-2)}{2}} \sum\limits_{j=0}^{k} \left \{ ^k_j \right \} \times 2^{n-j-1} \times n^{\underline{j+1}}$

 #include<cstdio>
#define int long long
#define mod 998244353
#define S 555555
int n,k,a[S],s[S],len=,fac[S],inv[S],rev[S],dw[S],ans;
int pow(int b,int t,int a=){for(;t;t>>=,b=b*b%mod)if(t&)a=a*b%mod;return a;}
void NTT(int *a,int opt){
for(int i=;i<len;++i)if(rev[i]>i)a[i]^=a[rev[i]]^=a[i]^=a[rev[i]];
for(int mid=;mid<len;mid<<=)
for(int i=,t=pow(,(mod-)//mid*opt+mod-);i<len;i+=mid<<)
for(int j=i,x,y,w=;j<mid+i;++j,w=w*t%mod)
x=a[j],y=a[j+mid]*w%mod,a[j]=(x+y)%mod,a[j+mid]=(x-y+mod)%mod;
if(opt==-)for(int i=,iv=pow(len,mod-);i<len;++i)a[i]=a[i]*iv%mod;
}
signed main(){
scanf("%lld%lld",&n,&k);fac[]=dw[]=;
while(len<=k<<)len<<=;
for(int i=;i<len;++i)rev[i]=rev[i>>]>>|(i&?len>>:);
for(int i=;i<=k;++i)fac[i]=fac[i-]*i%mod;
inv[k]=pow(fac[k],mod-);
for(int i=k-;~i;--i)inv[i]=inv[i+]*(i+)%mod;
for(int i=;i<=k;++i)a[i]=pow(i,k)*inv[i]%mod,s[i]=(i&?mod-:)*inv[i]%mod;
NTT(a,);NTT(s,);for(int i=;i<len;++i)s[i]=s[i]*a[i]%mod;NTT(s,-);
for(int i=;i<=k+;++i)dw[i]=dw[i-]*(n+-i)%mod;
for(int i=;i<=k;++i)ans=(ans+dw[i+]*pow(,n-i-)%mod*s[i])%mod;
printf("%lld\n",ans*pow(,(n-)*(n-)/)%mod);
}

真累啊。。。

bzoj5093图的价值:多项式,斯特林数(二项式反演)的更多相关文章

  1. 【bzoj5093】 [Lydsy1711月赛]图的价值 组合数+斯特林数+NTT

    Description "简单无向图"是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向 ...

  2. BZOJ.5093.[Lydsy1711月赛]图的价值(NTT 斯特林数)

    题目链接 对于单独一个点,我们枚举它的度数(有多少条边)来计算它的贡献:\[\sum_{i=0}^{n-1}i^kC_{n-1}^i2^{\frac{(n-2)(n-1)}{2}}\] 每个点是一样的 ...

  3. 【题解】BZOJ5093图的价值(二项式+NTT)

    [题解]BZOJ5093图的价值(二项式+NTT) 今天才做这道题,是我太弱了 强烈吐槽c++这种垃圾语言tmd数组越界不re反倒去别的数组里搞事情我只想说QAQ 推了一张A4纸的式子 考虑每个点的度 ...

  4. [CF932E]Team Work & [BZOJ5093]图的价值

    CF题面 题意:求\(\sum_{i=0}^{n}\binom{n}{i}i^k\) \(n\le10^9,k\le5000\) 模\(10^9+7\) BZOJ题面 题意:求\(n*2^{\frac ...

  5. BZOJ5093图的价值(斯特林数)

    题目描述 “简单无向图”是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向图的价值之和. 因为答案很大,请对 ...

  6. BZOJ5093 图的价值(NTT+斯特林数)

    显然每个点会提供相同的贡献.于是现在只考虑1号点的贡献.若其度数为i,则在2~n号点选i个连上,剩下的边随便连,这样可以算出答案为 这个式子可以O(n)计算.发现k比较小,于是考虑如何将这个式子化为与 ...

  7. bzoj5093:图的价值(第二类斯特林数+NTT)

    传送门 首先,题目所求为\[n\times 2^{C_{n-1}^2}\sum_{i=0}^{n-1}C_{n-1}^ii^k\] 即对于每个点\(i\),枚举它的度数,然后计算方案.因为有\(n\) ...

  8. BZOJ5093 图的价值——推式子+第二类斯特林数

    原题链接 题解 题目等价于求这个式子 \[ans=n2^{\frac{(n-1)(n-2)}{2}}\sum\limits_{i=0}^{n-1}\binom{n-1}{i}i^k\] 有这么一个式子 ...

  9. [BZOJ5093]图的价值(NTT+第二类Stirling数)

    5093: [Lydsy1711月赛]图的价值 Time Limit: 30 Sec  Memory Limit: 256 MBSubmit: 250  Solved: 130[Submit][Sta ...

随机推荐

  1. 打开TXT文件并显示

    <!DOCTYPE html> <html> <head> <title></title> <script> function ...

  2. 记录一次VMware与xshell远程链接的总结

    VMware 与xshell链接 说在前面 用了一年多的virtualbox,虚拟机环境的配置算是理清了,可最近参加红帽培训,大家用的都是VMware,想想我也试试吧,刚好再熟悉一下,于是就捣鼓了很久 ...

  3. eclipse 创建 springboot项目

    file --> new --> project --> Spring Boot  --> Spring  start  project Group:公司域名倒置,一般是com ...

  4. 《Java知识应用》Java压缩文件和解压缩

    今天通过Java实现一下:文件的压缩和解压缩. 1. 压缩 准备文件: 准备代码:(压缩) import java.io.BufferedInputStream; import java.io.Buf ...

  5. 守护线程,需要通过调用线程方法:setDaemon(boolean on)来进行设置

    package seday08.thread;/*** @author xingsir * 守护线程又称为后台线程,默认创建出来的线程都是普通线程, 守护线程需要通过调用线程方法:setDaemon( ...

  6. 动态代理模式_应用(Redis工具类)

    本次使用动态代理的初衷是学习Redis,使用Java操作Redis时用到Jedis的JedisPool,而后对Jedis的方法进一步封装完善成为一个工具类.因为直接使用Jedis对象时,为了保证性能, ...

  7. HTTP认知(请求与响应)

    web的工作是:浏览器发送请求报文 + 服务端返回响应报文 通俗的说一下web工作的一个流程: 浏览器向服务端发送HTTP请求报文:这条请求报文组成由请求行.请求头.请求体三大部分组成: 1.请求行 ...

  8. Vim 基本的使用

    三种模式 按 ESC 进入命令模式 命令模式下输入 Shift + : 进入末行模式 命令模式下输入插入命令,如(i,a,o) 进入输入模式 进入 vim 文件名 vim直接编辑一个文件,如果是已经存 ...

  9. docker拉取慢(net/http: request canceled while waiting for connection (Client.Timeout exceeded while awaiting head)

    systemctl restart docker在docker-hub拉取慢,因为服务器在外网 直接配置阿里云镜像就可以 首先: vim /etc/docker/daemon.json加入下面的那句 ...

  10. STL常用结构与方法简明总结

    C++常用的数据结构 序列式容器 vector(向量.有序数列),list(双向链表),deque(双端队列) 适配器容器 stack(栈),queue(队列) 关联式容器 map(映射.键值对二叉树 ...