【BZOJ4555】[TJOI&HEOI2016]求和 斯特林数+NTT
Description
在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心。
现在他想计算这样一个函数的值:
S(i, j)表示第二类斯特林数,递推公式为:
S(i, j) = j ∗ S(i − 1, j) + S(i − 1, j − 1), 1 <= j <= i − 1。
边界条件为:S(i, i) = 1(0 <= i), S(i, 0) = 0(1 <= i)
你能帮帮他吗?
Input
输入只有一个正整数
Output
输出f(n)。由于结果会很大,输出f(n)对998244353(7 × 17 × 223 + 1)取模的结果即可。1 ≤ n ≤ 100000
Sample Input
3
Sample Output
87
Sol
本题递推公式没用,因为我们需要在\(nlogn\)时间内求出这个结果。
首先我们根据第二类斯特林数的定义“把i个数字分到j的相同集合的方案数”,得:
\(s(i,j)=\frac{1}{j!}*\sum_{k=0}^{j}(-1)^kC^k_j(j-k)^i\)
意义就是我们假设至少k个集合为空,然后用组合数算出选择空集的方案数,之后剩下\(j-k\)个集合,i个数字,那么方案数就是\((j-k)^i\)。但是这样不保证剩下的严格非空,所以要容斥一波。然后因为上述是有序的,而第二类斯特林数是无序的,所以要除以阶乘。
之后我们把组合数拆开,得到:
\(s(i,j)=\sum_{k=0}^{j}\frac{(-1)^k}{k!}*\frac{(j-k)^i}{(j-k)!}\)
带入所求的式子中,得到:
\(f(n)=\sum_{i=0}^{n}\sum_{j=0}^{i}\sum_{k=0}^{j}\frac{(-1)^k}{k!}*\frac{(j-k)^i}{(j-k)!}\)
因为\(j>i\)的时候右面不会产生贡献,所以我们可以把j的范围写到n:
\(f(n)=\sum_{i=0}^{n}\sum_{j=0}^{n}j!2^j\sum_{k=0}^{j}\frac{(-1)^k}{k!}*\frac{(j-k)^i}{(j-k)!}\)
尽管是卷积,但是还是得\(O(n^2logn)\),但是i只用到了一处,所以我们更换循环顺序,得:
\(f(i)=\sum_{j=0}^{n}j!2^j\sum_{k=0}^{j}\frac{(-1)^k}{k!}*\frac{\sum_{i=0}^{n}(j-k)^i}{(j-k)!}\)
那个带\(\sum\)的是个等比数列,可以\(O(1)\)计算,然后就是\(NTT\)啦。
时间复杂度\(O(nlogn)\)。
Code
#include <cstdio>
int i,j,k,I[100005],IF[100005],a[262145],b[262145],F[262145],B[100005],P=998244353,w,wn,t,n,len,ans;
int ksm(int a,int b){int res=1;for(;b;b>>=1,a=1ll*a*a%P) if(b&1) res=1ll*res*a%P;return res;}
void ntt(int *a,int n,int op)
{
for(i=k=0;i<n;i++){if(i>k) a[i]^=a[k]^=a[i]^=a[k];for(j=(n>>1);(k^=j)<j;j>>=1);}
for(k=2,wn=ksm(3,op==1?(P-1)/k:P-1-(P-1)/k);k<=n;k<<=1,wn=ksm(3,op==1?(P-1)/k:P-1-(P-1)/k))
for(i=0,w=1;i<n;i+=k,w=1) for(j=0;j<(k>>1);j++,w=1ll*w*wn%P)
t=1ll*a[i+j+(k>>1)]*w%P,a[i+j+(k>>1)]=(a[i+j]-t+P)%P,a[i+j]=(a[i+j]+t)%P;
if(op==-1) for(t=ksm(n,P-2),i=0;i<n;i++) a[i]=1ll*a[i]*t%P;
}
int main()
{
scanf("%d",&n);for(len=1;len<=n*2;len<<=1);
F[0]=B[0]=I[1]=F[1]=1,B[1]=2;for(int i=2;i<=n;i++) F[i]=1ll*F[i-1]*i%P,B[i]=1ll*B[i-1]*2%P,I[i]=1ll*I[P%i]*(P-P/i)%P;
IF[n]=ksm(F[n],P-2);for(int i=n-1;~i;i--) IF[i]=1ll*IF[i+1]*(i+1)%P;
a[0]=b[0]=1;for(int i=1;i<=n;i++) b[i]=(((i&1)?-1:1)*IF[i]+P)%P,a[i]=(i==1)?n+1:1ll*IF[i]*(ksm(i,n+1)-1)%P*I[i-1]%P;
ntt(a,len,1);ntt(b,len,1);
for(int i=0;i<len;i++) a[i]=1ll*a[i]*b[i]%P;ntt(a,len,-1);
for(int i=0;i<=n;i++) ans=(ans+1ll*B[i]*F[i]%P*a[i]%P)%P;
printf("%d\n",ans);
}
【BZOJ4555】[TJOI&HEOI2016]求和 斯特林数+NTT的更多相关文章
- [HEOI2016/TJOI2016][bzoj4555] 求和 [斯特林数+NTT]
题面 传送门 思路 首先,我们发现这个式子中大部分的项都和$j$有关(尤其是后面的$2^j\ast j!$),所以我们更换一下枚举方式,把这道题的枚举方式变成先$j$再$i$ $f(n)=\sum_{ ...
- 【BZOJ】4555: [Tjoi2016&Heoi2016]求和 排列组合+多项式求逆 或 斯特林数+NTT
[题意]给定n,求Σi=0~nΣj=1~i s(i,j)*2^j*j!,n<=10^5. [算法]生成函数+排列组合+多项式求逆 [题解]参考: [BZOJ4555][Tjoi2016& ...
- [BZOJ4555 TJOI2016 HEOI2016 求和]
第一篇博客,请大家多多关照.(鞠躬 BZOJ4555 TJOI2016 HEOI2016 求和 题意: 给定一个正整数\(n\)(\(1\leqq n \leqq100000\)),求: \[ ...
- BZOJ4555 [Tjoi2016&Heoi2016]求和 【第二类斯特林数 + NTT】
题目 在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心. 现在他想计算这样一个函数的值: S(i, j)表示第二类斯特林数,递推公式为: S(i, j) = j ∗ S(i − 1, j) + ...
- 【BZOJ4555】【TJOI2016】【HEOI2016】求和 (第二类斯特林数+NTT卷积)
Description 在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心. 现在他想计算这样一个函数的值: $$f(n)=\sum_{i=0}^n\sum_{j=0}^i S(i,j)\tim ...
- bzoj 4555 [Tjoi2016&Heoi2016] 求和 —— 第二类斯特林数+NTT
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4555 关于第二类斯特林数:https://www.cnblogs.com/Wuweizhen ...
- P4091 [HEOI2016/TJOI2016]求和(第二类斯特林数+NTT)
传送门 首先,因为在\(j>i\)的时候有\(S(i,j)=0\),所以原式可以写成\[Ans=\sum_{i=0}^n\sum_{j=0}^nS(i,j)\times 2^j\times j! ...
- BZOJ 4555:[TJOI2016&HEOI2016]求和(第二类斯特林数+NTT)
题目链接 \(Description\) 求 \[\sum_{i=0}^n\sum_{j=0}^iS(i,j)2^jj!\]对998244353取模后的结果. \(n<=10^5\) \(Sol ...
- P4091-[HEOI2016/TJOI2016]求和【斯特林数,NTT】
正题 题目链接:https://www.luogu.com.cn/problem/P4091 题目大意 给出\(n\),求 \[\sum_{i=0}^n\sum_{j=0}^i\begin{Bmatr ...
随机推荐
- 免Oracle客户端程序监听程序配置
Oracle默认安装时,监听程序和tnsnames程序中的监听方式都是默认的localhost,但免客户端的程序是连接不上的.这时需要: 1.将listener中的(HOST = localhost) ...
- Django界面不能添加中文解决办法
Django项目部署好后,界面添加中文会报错,解决办法: 创建数据库时要指定编码格式: CREATE DATABASE blog CHARACTER SET utf8; 如果已经创建完毕则修改: al ...
- CSS——优先级
转自:http://www.planabc.net/2008/05/06/css_specificity/ CSS2.1 中规定了关于 CSS 规则 Specificity(特异性)的计算方式,用一个 ...
- JQuery UI Datepicker中文显示的方法
出自:http://www.aimks.com/method-to-display-the-jquery-ui-datepicker-chinese.html Query UI Datepicker这 ...
- Mycat之日志分析 select * from travelrecord order by id limit 100000,100 的执行过程
解释:mycat在执行分页排序的时候,分成2步走.如果M很大,会改写成 limit 0,m+n, 然后在每个MYSQL分片上排序后还需要在mycat汇总输出,所以会很慢.下面是详细执行计划以及日志输出 ...
- MySQL备份还原之一mydumper
1)源码编译安装 1.下载 mydumper源码 2.解压 [mysql@localhost ~]$ tar -xvf mydumper-0.9.1.tar mydumper-0.9.1/CMakeL ...
- 19-从零玩转JavaWeb-继承关系与方法覆盖
配套视频详解 继承子类与父类的关系 继承子类与父类的关系代码演示 继承方法的覆盖 方法覆盖原则 一.继承编写过程 二.子类可以继承父类当中哪些内容 可以自行把修饰符换下, 然后查看子类能不能访问 ...
- Java多线程-线程的同步(同步代码块)
对于同步,除了同步方法外,还可以使用同步代码块,有时候同步代码块会带来比同步方法更好的效果. 追其同步的根本的目的,是控制竞争资源的正确的访问,因此只要在访问竞争资源的时候保证同一时刻只能一个线程访问 ...
- MAVEN学习总结1
一.安装Maven插件 下载下来的maven插件如下图所示:,插件存放的路径是:E:/MavenProject/Maven2EclipsePlugin
- Nginx静态网站的部署
静态网站的部署 首先先看一下nginx/conf/nginx.conf 配置文件内的信息: #user nobody; worker_processes 1; #error_log logs/erro ...