[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{n(n-1))}{2}-(n-1)}*\sum_{i=0}^{n-1}\binom{n-1}{i}i^k\)
\(n\le10^9,k\le2*10^5\)
模\(998244353\)
第二类斯特林数
赶紧去学第二类斯特林数啊
第二类斯特林数:\(S(n,m)\),表示把\(n\)个不同的的球放到\(m\)个相同的盒子且不允许空盒的方案数。
首先有一个\(O(n^2)\)的递推。考虑第\(n\)个球的放置,如果单独放置方案数就等于\(S(n-1,m-1)\),否则就是\(m*S(n-1,m)\),表示枚举放在那个盒子里。写出来就是
\]
边界条件:\(S(0,0)=1,S(i,0)=0(i>0)\)
然后还有各种考虑组和意义得出来的柿子
第二类斯特林数存在通项公式。因为\(m^n\)表示把\(n\)个不同的球放到\(m\)个不同的盒子里且允许空盒的方案数,我们想办法用这个东西构造出第二类斯特林数的组和意义。
首先盒子是否相同很好办,直接乘/除以\(m!\)即可。关键问题在于空盒的处理。
考虑容斥,计算“至少有\(k\)个空盒然后剩下的随便是不是空盒的方案数”,乘上一个容斥系数
\]
还有一个,反过来用第二类斯特林数来表示\(m^n\)
其实就是一个反演,但如果也同样考虑组和意义的话岂不妙哉。
枚举有几个盒子不是空的,可以直接得到如下柿子:
\]
sol
所以说这道题怎么做呢?
开始推柿子辣~
(k,j)*j!*\binom{n}{j}*2^{n-j}\\=\sum_{j=0}^{k}S
(k,j)*j!*\binom{n}{j}*2^{n-j}\]
注意\(\sum_{i=j}^{n}\binom{n}{i}\binom{i}{j}\)的组合意义:从\(n\)个里面选出任意多个(大于等于\(j\))再从中选出\(j\)个,相当于从\(n\)个中选\(j\)个然后剩下的\(n-j\)个随便选。
柿子里面的阶乘\(j!\)、组合数\(\binom{n}{j}\)、2的次幂\(2^{n-j}\)都可以\(O(k)\)处理出来。复杂度瓶颈在于第二类斯特林数的处理。
\(k\le5000\)可以直接用递推式\(O(k^2)\)处理出来,而观察其通项式会发现其实就是一个卷积形式:
\]
所以直接上\(NTT\),复杂度\(O(k\log{k})\)
但是\(CF\)上的那道题不能用\(NTT\)诶,毕竟\(10^9+7\)怎么\(NTT\)(我没说不可以啊,只是我还不会呀。。。)
update 3.28 放心\(MTT\)这个坑已经填了。
code
[CF932E]Team Work
这里我强行用通项公式然后写了个\(O(k^2)\)的多项式乘法
#include<cstdio>
#include<algorithm>
using namespace std;
const int mod = 1e9+7;
int gi()
{
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
int fastpow(int a,int b)
{
int res=1;
while (b) {if (b&1) res=1ll*res*a%mod;a=1ll*a*a%mod;b>>=1;}
return res;
}
const int N = 5005;
int jc[N],C[N],two[N],inv,a[N],b[N],s[N],ans;
int main()
{
int n=gi(),k=gi();
jc[0]=1;
for (int i=1;i<=k;++i) jc[i]=1ll*jc[i-1]*i%mod;
C[0]=1;
for (int i=1;i<=k;++i) C[i]=1ll*C[i-1]*(n-i+1)%mod*fastpow(i,mod-2)%mod;
two[0]=fastpow(2,n);inv=fastpow(2,mod-2);
for (int i=1;i<=k;++i) two[i]=1ll*two[i-1]*inv%mod;
for (int i=0;i<=k;++i) a[i]=i&1?mod-fastpow(jc[i],mod-2):fastpow(jc[i],mod-2);
for (int i=0;i<=k;++i) b[i]=1ll*fastpow(i,k)*fastpow(jc[i],mod-2)%mod;
for (int i=0;i<=k;++i)
for (int j=0;j<=i;++j)
(s[i]+=1ll*a[j]*b[i-j]%mod)%=mod;
for (int i=1;i<=k;++i) (ans+=1ll*s[i]*jc[i]%mod*C[i]%mod*two[i]%mod)%=mod;
printf("%d\n",ans);
return 0;
}
[BZOJ5093]图的价值
#include<cstdio>
#include<algorithm>
using namespace std;
const int _ = 800005;
const int mod = 998244353;
int gi()
{
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
int fastpow(int a,int b)
{
int res=1;
while (b) {if (b&1) res=1ll*res*a%mod;a=1ll*a*a%mod;b>>=1;}
return res;
}
int n,N,k,l;
int jc[_],C[_],two[_],inv,a[_],b[_],rev[_],ans;
void NTT(int *P,int opt)
{
for (int i=0;i<N;++i) if (i>rev[i]) swap(P[i],P[rev[i]]);
for (int i=1;i<N;i<<=1)
{
int W=fastpow(3,(mod-1)/(i<<1));
if (opt==-1) W=fastpow(W,mod-2);
for (int j=0,p=i<<1;j<N;j+=p)
{
int w=1;
for (int k=0;k<i;++k,w=1ll*w*W%mod)
{
int x=P[j+k],y=1ll*P[j+k+i]*w%mod;
P[j+k]=(x+y)%mod;P[j+k+i]=(x-y+mod)%mod;
}
}
}
if (opt==-1)
{
int Inv=fastpow(N,mod-2);
for (int i=0;i<N;++i) P[i]=1ll*P[i]*Inv%mod;
}
}
int main()
{
n=gi()-1;k=gi();
jc[0]=1;
for (int i=1;i<=k;++i) jc[i]=1ll*jc[i-1]*i%mod;
C[0]=1;
for (int i=1;i<=k;++i) C[i]=1ll*C[i-1]*(n-i+1)%mod*fastpow(i,mod-2)%mod;
two[0]=fastpow(2,n);inv=fastpow(2,mod-2);
for (int i=1;i<=k;++i) two[i]=1ll*two[i-1]*inv%mod;
for (int i=0;i<=k;++i) a[i]=i&1?mod-fastpow(jc[i],mod-2):fastpow(jc[i],mod-2);
for (int i=0;i<=k;++i) b[i]=1ll*fastpow(i,k)*fastpow(jc[i],mod-2)%mod;
for (N=1;N<=k*2;N<<=1) ++l;--l;
for (int i=0;i<N;++i) rev[i]=(rev[i>>1]>>1)|((i&1)<<l);
NTT(a,1);NTT(b,1);
for (int i=0;i<N;++i) a[i]=1ll*a[i]*b[i]%mod;
NTT(a,-1);
for (int i=0;i<=k;++i) (ans+=1ll*a[i]*jc[i]%mod*C[i]%mod*two[i]%mod)%=mod;
int Pow=(1ll*n*(n+1)/2-n)%(mod-1);
printf("%d\n",1ll*ans*fastpow(2,Pow)%mod*(n+1)%mod);
return 0;
}
[CF932E]Team Work & [BZOJ5093]图的价值的更多相关文章
- 【学术篇】CF932E Team Work && bzoj5093 图的价值
两个题的传送门 对于CF这道题, 分别考虑每种可能的集合大小, 每个大小为\(k\)的集合数量有\(\binom nk\)个, 所以最后的答案就是 \[\sum_{i=0}^n\binom{n}{i} ...
- 【题解】BZOJ5093图的价值(二项式+NTT)
[题解]BZOJ5093图的价值(二项式+NTT) 今天才做这道题,是我太弱了 强烈吐槽c++这种垃圾语言tmd数组越界不re反倒去别的数组里搞事情我只想说QAQ 推了一张A4纸的式子 考虑每个点的度 ...
- [BZOJ5093]图的价值(NTT+第二类Stirling数)
5093: [Lydsy1711月赛]图的价值 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 250 Solved: 130[Submit][Sta ...
- BZOJ5093图的价值(斯特林数)
题目描述 “简单无向图”是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向图的价值之和. 因为答案很大,请对 ...
- bzoj5093图的价值:多项式,斯特林数(二项式反演)
Description “简单无向图”是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向图的价值之和. 因为 ...
- Bzoj5093: 图的价值
题面 Bzoj Sol 一张无向无重边自环的图的边数最多为\(\frac{n(n-1)}{2}\) 考虑每个点的贡献 \[n*2^{\frac{n(n-1)}{2} - (n-1)}\sum_{i=0 ...
- BZOJ5093 图的价值(NTT+斯特林数)
显然每个点会提供相同的贡献.于是现在只考虑1号点的贡献.若其度数为i,则在2~n号点选i个连上,剩下的边随便连,这样可以算出答案为 这个式子可以O(n)计算.发现k比较小,于是考虑如何将这个式子化为与 ...
- bzoj5093:图的价值(第二类斯特林数+NTT)
传送门 首先,题目所求为\[n\times 2^{C_{n-1}^2}\sum_{i=0}^{n-1}C_{n-1}^ii^k\] 即对于每个点\(i\),枚举它的度数,然后计算方案.因为有\(n\) ...
- BZOJ5093 图的价值——推式子+第二类斯特林数
原题链接 题解 题目等价于求这个式子 \[ans=n2^{\frac{(n-1)(n-2)}{2}}\sum\limits_{i=0}^{n-1}\binom{n-1}{i}i^k\] 有这么一个式子 ...
随机推荐
- js压缩上传图片
初学有不当之处,请多多指点, <body> <div class="cc"> <input type="file" id=&quo ...
- shell编程值之shell流程控制(7)
条件判断式 1 按照文件类型判断(常用类型) 测试类型 作用 -d 文件 判断该文件是否存在,并且是否为目录文件(是目录文件为真) -e 文件 判断该文件是否存在(存在为真) -f 文件 判断该文件是 ...
- JavaScript对象创建的几种方式
1 工厂模式 1.1 创建 function createFruit(name,colors) { var o = new Object(); o.name = name; o.colors = co ...
- CodeForces 820C
题意略. 这道题目的出题者竟然被hack了!? 我的思路是:在游戏开始时,为了尽量少地用字母,我应该尽量选取计算机输出的前a个字母中已经使用过的字母.但是为了使电脑也尽量少用字母,我添加的这b个字母应 ...
- css实现多行多列的布局
1.两列多行: HTML: <div class="box1"> box1:实现两列多行布局 <ul> <li>111</li> & ...
- 排序算法(Java实现)
这几天一直在看严蔚敏老师的那本<数据结构>那本书.之前第一次学懵懵逼逼,当再次看的时候,发觉写的是非常详细,非常的好. 那就把相关的排序算法用我熟悉的Java语言记录下来了.以下排序算法是 ...
- java程序员理解js中的闭包
1.闭包概念: 就是函数内部通过某种方式访问一个函数内部的局部变量 再次理解: 闭包产生原因: 1.内部函数引用了外部函数的变量 作用:延长局部变量的生命周期 让函数外部可以调用到函数内部的数据 利用 ...
- PHP 支持加解密的函数
function encrypt($string,$operation,$key=''){ $key=md5($key); $key_length=strlen($key); $string=$ope ...
- PHP判断是手机端还是PC端
function check_wap() { if (isset($_SERVER['HTTP_VIA'])) return true; if (isset($_SERVER['HTTP_X_NOKI ...
- DAY2-JAVA
2018-1-28学习笔记 1.在开发中定义类.方法时也可以先添加文档注释,再用javadoc工具生成自己的API文档.Java文档注释和java工具使用. 2.掌握查阅java JDK和API文档. ...