博客
关于我
[bzoj5093][第二类斯特林数][NTT]图的价值
阅读量:86 次
发布时间:2019-02-26

本文共 2338 字,大约阅读时间需要 7 分钟。

Description

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

给定n和k,请计算所有n个点的带标号的简单无向图的价值之和。 因为答案很大,请对998244353取模输出。

Input

第一行包含两个正整数n,k(1<=n<=10^9,1<=k<=200000)。

Output

输出一行一个整数,即答案对998244353取模的结果。

Sample Input

6 5

Sample Output

67584000

题解

一般都要单独考虑一个点的贡献

显然这个点的贡献是独立的,最后乘一个n就可以
那单独一个点的贡献就是
∑ i = 1 n − 1 C n − 1 i ∗ 2 ( n − 1 ) ( n − 2 ) 2 ∗ i k \sum_{i=1}^{n-1}C_{n-1}^i*2^{\frac{(n-1)(n-2)}{2}}*i^k i=1n1Cn1i22(n1)(n2)ik
2的次幂是常数,提到前面去,把n减一,显然我们只需要算
∑ i = 1 n C n i ∗ i k \sum_{i=1}^nC_{n}^{i}*i^k i=1nCniik
根据第二类斯特林数的性质,我们可以知道
∑ i = 1 n C n i ∑ j = 0 m i n ( i , k ) S ( k , j ) ∗ j ! ∗ C i j \sum_{i=1}^nC_{n}^{i}\sum_{j=0}^{min(i,k)}S(k,j)*j!*C_{i}^{j} i=1nCnij=0min(i,k)S(k,j)j!Cij
乘法分配律一下,显然可以把斯特林数提前
∑ i = 0 m i n ( n , k ) S ( k , i ) ∗ i ! ∗ ∑ j = i n C n j ∗ C j i \sum_{i=0}^{min(n,k)}S(k,i)*i!*\sum_{j=i}^{n}C_n^j*C_{j}^{i} i=0min(n,k)S(k,i)i!j=inCnjCji
后面的式子的意义就是在n中选i个数,其它可选可不选,所以就是
∑ i = 0 m i n ( n , k ) S ( k , i ) ∗ i ! ∗ C n i ∗ 2 n − i \sum_{i=0}^{min(n,k)}S(k,i)*i!*C_{n}^{i}*2^{n-i} i=0min(n,k)S(k,i)i!Cni2ni
现在就是要算前面的斯特林数了
可以容斥然后NTT n l o g n nlogn nlogn求,点一下就知道啦

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define LL long long#define mp(x,y) make_pair(x,y)#define pll pair
#define pii pair
using namespace std;inline LL read(){ LL f=1,x=0;char ch=getchar(); while(ch<'0'||ch>'9'){ if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){ x=x*10+ch-'0';ch=getchar();} return x*f;}int stack[20];inline void write(int x){ if(x<0){ putchar('-');x=-x;} if(!x){ putchar('0');return;} int top=0; while(x)stack[++top]=x%10,x/=10; while(top)putchar(stack[top--]+'0');}inline void pr1(int x){ write(x);putchar(' ');}inline void pr2(int x){ write(x);putchar('\n');}const LL mod=998244353;const int MAXN=200005;LL pow_mod(LL a,LL b){ LL ret=1; while(b) { if(b&1)ret=ret*a%mod; a=a*a%mod;b>>=1; } return ret;}LL A[MAXN*4],B[MAXN*4];int R[MAXN*4],L;void NTT(LL *y,int len,int on){ for(int i=0;i
=0;i--)inv[i]=inv[i+1]*(i+1)%mod; n=read();K=read();n--; init(); int ln; for(ln=1;ln<=2*K;ln<<=1)L++; for(int i=0;i
>1]>>1)|(i&1)<<(L-1); int temp=-1; for(int i=0;i<=K;i++) { temp*=-1; A[i]=temp*inv[i]%mod; LL u1=pow_mod(i,K),u2=inv[i]; B[i]=u1*u2%mod; } NTT(A,ln,1);NTT(B,ln,1); for(int i=0;i

转载地址:http://immu.baihongyu.com/

你可能感兴趣的文章
MySQL5.7.37windows解压版的安装使用
查看>>
mysql5.7免费下载地址
查看>>
mysql5.7命令总结
查看>>
mysql5.7安装
查看>>
mysql5.7性能调优my.ini
查看>>
MySQL5.7新增Performance Schema表
查看>>
Mysql5.7深入学习 1.MySQL 5.7 中的新增功能
查看>>
Webpack 之 basic chunk graph
查看>>
Mysql5.7版本单机版my.cnf配置文件
查看>>
mysql5.7的安装和Navicat的安装
查看>>
mysql5.7示例数据库_Linux MySQL5.7多实例数据库配置
查看>>
Mysql8 数据库安装及主从配置 | Spring Cloud 2
查看>>
mysql8 配置文件配置group 问题 sql语句group不能使用报错解决 mysql8.X版本的my.cnf配置文件 my.cnf文件 能够使用的my.cnf配置文件
查看>>
MySQL8.0.29启动报错Different lower_case_table_names settings for server (‘0‘) and data dictionary (‘1‘)
查看>>
MYSQL8.0以上忘记root密码
查看>>
Mysql8.0以上重置初始密码的方法
查看>>
mysql8.0新特性-自增变量的持久化
查看>>
Mysql8.0注意url变更写法
查看>>
Mysql8.0的特性
查看>>
MySQL8修改密码报错ERROR 1819 (HY000): Your password does not satisfy the current policy requirements
查看>>