题面

Words

有 \(n\) 天,每天插入一个字符集大小为 \(c\) 长度为 \(l\) 的字符串,求每一天建立 \(\tt Trie\) 树的期望节点数(根节点不算)模 \(998244353\)。

数据范围:\(1\le n\le 10^5\),\(1\le c,l\le 10^9\)。

  • Input
5 4 3
  • Output
4
911976330
792083550
276733174
815453946

正解

转化问题:有一颗高度为 \(l+1\) 的 \(c\) 叉完全 \(\tt Trie\),从根到叶走 \(n\) 遍,求期望经过节点数(根节点不算)。

一层一层考虑,答案为(因为不考虑根节点,所以 \(i\) 从 \(1\) 开始):\(Ans(n)=\sum_{i=1}^l a(i,n)\)。

考虑新的路径会不会走新节点,可以递推:

\[a(i,n)=a(i,n-1)+\frac{c^i-a(i,n-1)}{c^i}\\
a(i,n)=1+a(i,n-1)(1-\frac{1}{c^i})\\
\]

特征方程:\(t=1+(1-\frac{1}{c^i})t\Longrightarrow t=c^i\)。

\[a(i,n)-c^i=(1-\frac{1}{c^i})(a(i,n-1)-c^i)\\
\begin{split}
a(i,n)=&(1-\frac{1}{c^i})(a(i,n-1)-c^i)+c^i\\
=&(1-\frac{1}{c^i})^n(a(i,0)-c^i)+c^i\\
=&c^i-c^i(1-\frac{1}{c^i})^n\\
\end{split}
\]

带回到上面求每一天答案的式子,展开二项式 \(^{\color{#bcbcee}{[1]}}\),交换枚举顺序 \(^{\color{#bceebc}{[2]}}\),抵消 \(^{\color{#eebcbc}{[3]}}\):

\[\begin{split}
Ans(n)=&\sum_{i=1}^l a(i,n)\\
=&\sum_{i=1}^l\left(c^i-c^i(1-\frac{1}{c^i})^n\right)\\
=&\sum_{i=1}^l c^i-\sum_{i=1}^lc^i(1-\frac{1}{c^i})^n\\
=&\sum_{i=1}^l c^i-\sum_{i=1}^lc^i\sum_{j=0}^n\left(\frac{1}{c^i}\right)^j(-1)^j{n\choose j}&^{\color{#bcbcee}{[1]}}\\
=&\sum_{i=1}^l c^i-\sum_{i=0}^n(-1)^i{n\choose i}\sum_{j=1}^l c^j\left(\frac{1}{c^j}\right)^i&^{\color{#bceebc}{[2]}}\\
=&-\sum_{i=1}^n(-1)^i{n\choose i}\sum_{j=1}^l c^j\left(\frac{1}{c^j}\right)^i&^{\color{#eebcbc}{[3]}}\\
=&\sum_{i=1}^n(-1)^{i+1}{n\choose i}\sum_{j=1}^l c^j\left(\frac{1}{c^j}\right)^i\\
=&\sum_{i=1}^n\frac{n!}{i!(n-i)!}(-1)^{i+1}\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(i-1)}\\
=&n!\sum_{i=1}^n\frac{1}{i!(n-i)!}(-1)^{i+1}\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(i-1)}\\
\end{split}
\]

发现有 \(i!\) 和 \((n-i)!\),并且每个 \(Ans(x)\) 都要求,所以想到 \(\tt NTT\)。

\[\begin{split}
Ans(n)=&n!\sum_{i=1}^n\frac{1}{i!(n-i)!}(-1)^{i+1}\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(i-1)}\\
=&n!\sum_{i+j=n}\left(\frac{(-1)^{i+1}\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(i-1)}}{i!}\right)\left(\frac{1}{j!}\right)\\
\end{split}
\]

只需让 \(f(n)=\frac{(-1)^{n+1}\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(n-1)}}{n!}\),\(g(n)=\frac{1}{n!}\) 卷积即可。

其中 \(\sum_{j=1}^l \left(\frac{1}{c}\right)^{j(n-1)}\) 可以用等比数列的公式算。

时间复杂度 \(\Theta(n\log)\)。


代码

#include <bits/stdc++.h>
using namespace std; //Start
typedef long long ll;
typedef double db;
#define mp(a,b) make_pair(a,b)
#define x first
#define y second
#define be(a) a.begin()
#define en(a) a.end()
#define sz(a) int((a).size())
#define pb(a) push_back(a)
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3f; //Data
const int N=1e5;
const int mod=998244353;
int n,l,c; //Pow
int Pow(int a,int x){
if(!a) return 0; int res=1;
for(;x;a=(ll)a*a%mod,x>>=1)if(x&1) res=(ll)res*a%mod;
return res;
} //NTT
int iv[N+1],f[N<<2],g[N<<2];
int up(int len){return 1<<int(ceil(log2(len)));}
void Getpoly(){
for(int i=iv[0]=1;i<=n;i++) iv[i]=(ll)iv[i-1]*i%mod;
for(int i=0;i<=n;i++) iv[i]=Pow(iv[i],mod-2);
for(int i=0;i<=n;i++) g[i]=iv[i];
for(int i=1;i<=n;i++){
int ic=Pow(Pow(c,i-1),mod-2);
int q=(ic==1)?l:(ll)ic*(Pow(ic,l)-1+mod)%mod*Pow(ic-1,mod-2)%mod;
f[i]=(ll)iv[i]*q%mod*((i&1)?1:mod-1)%mod;
}
}
const int G=3,iG=332748118;
int lim,r[N<<2];
void NTT(int a[],int t){
for(int i=0;i<lim;i++)if(i<r[i]) swap(a[i],a[r[i]]);
for(int mid=1;mid<lim;mid<<=1){
int wn=Pow(t==1?G:iG,(mod-1)/(mid<<1));
for(int j=0;j<lim;j+=(mid<<1))
for(int w=1,k=j;k<mid+j;w=(ll)w*wn%mod,k++){
int x=a[k],y=(ll)w*a[mid+k]%mod;
a[k]=(x+y)%mod,a[mid+k]=(x-y+mod)%mod;
}
}
if(t==-1){
int ilim=Pow(lim,mod-2);
for(int i=0;i<lim;i++) a[i]=(ll)a[i]*ilim%mod;
}
}
void Mulpoly(){
lim=up(n+n);
for(int i=0;i<lim;i++) r[i]=(r[i>>1]>>1)|((i&1)*lim>>1);
NTT(f,1),NTT(g,1);
for(int i=0;i<lim;i++) f[i]=(ll)f[i]*g[i]%mod;
NTT(f,-1);
} //Main
int main(){
ios::sync_with_stdio(0);
cin.tie(0),cout.tie(0);
cin>>n>>l>>c;
Getpoly(),Mulpoly();
for(int i=1,fa=1;i<=n;i++) fa=(ll)fa*i%mod,f[i]=(ll)f[i]*fa%mod;
for(int i=1;i<=n;i++) cout<<f[i]<<'\n';
return 0;
}

祝大家学习愉快!

题解-Words的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

  10. JSOI2016R3 瞎BB题解

    题意请看absi大爷的blog http://absi2011.is-programmer.com/posts/200920.html http://absi2011.is-programmer.co ...

随机推荐

  1. quic 2 ietf-transport-draft-ietf-quic-transport-09

    quic的优点 Low-latency connection establishment Multiplexing without head-of-line blocking Authenticate ...

  2. 一篇文章了解_unittest

    1. 基本概念 2018年10月7日 星期日 11:39 unittest是python自带的单元测试框架,有时候又被称为"PyUnit",是python版本的JUint实现. 该 ...

  3. 统计数字问题(Java)

    Description 一本书的页码从自然数1 开始顺序编码直到自然数n.书的页码按照通常的习惯编排,每个页码都不含多余的前导数字0.例如,第6 页用数字6 表示,而不是06 或006 等.数字计数问 ...

  4. ClickHouse数据库数据定义手记之数据类型

    前提 前边一篇文章详细分析了如何在Windows10系统下搭建ClickHouse的开发环境,接着需要详细学习一下此数据库的数据定义,包括数据类型.DDL和DML.ClickHouse作为一款完备的D ...

  5. webug第十三关:XSS

    第十三关:XSS 点击链接 xss,弹框框

  6. Java的BigDecimal,对运算封装

    添加maven依赖 <dependency> <groupId>com.google.guava</groupId> <artifactId>guava ...

  7. Mac下载工具folx如何下载常用的软件

    最近,多档综艺节目都开展得如火如荼,比如<中国新说唱>.<这就是街舞>等深受年轻人喜欢的综艺节目.虽然手机端也可观看,但可以的话,当然是使用电脑屏幕观看节目比较过瘾. 接下来, ...

  8. 2个快速制作完成一幅思维导图的iMindMap思维导图用法

    随着思维导图的流行,与其相关的思维导图制作软件如雨后春笋,纷纷进入我们的视野中,更让人难以选择.那想要入门的萌新该如何开始这个新的旅途呢? 各式各样的思维导图制作软件当中,有一个软件得到了大家一致的好 ...

  9. Vegas视频的音频叠加效果怎么实现,可以用其他软件吗

    有时我们会用Vegas为某段影片配音,我们要怎么把配音和背景声融合在一起呢?想必马上会有人反应过来:让配音和背景声分别置于两条轨道上就好了.这当然是一个相当好的方式. 可是,如果我想要把两段音频合成一 ...

  10. 循序渐进VUE+Element 前端应用开发(29)--- 高级查询条件的界面设计

    在系统模块中的业务列表展示里面,一般我们都会在列表中放置一些查询条件,如果是表字段不多,大多数情况下,放置的条件有十个八个就可以了,如果是字段很多,而这些条件信息也很关键的时候,就可能放置很多条件,但 ...