2021ICPC济南区域赛题解CDJKL

2021ICPC济南区域赛题解

题目链接:
https://pintia.cn/problem-sets/1459829212832296960

讲义、代码链接:
https://hytidel.lanzoub.com/b031c9dti
密码:cj1j

视频讲解链接:
https://www.bilibili.com/video/BV1Qe4y1D7c4

---

---

K. Search For Mafuyu

题意

有一棵包含编号为$1\sim n$的$n$个节点的树.初始时你在$1$号节点,你想找的人M等概率地在$2\sim n$号节点.每一秒,你可走到与当前节点相邻的节点处.若你采取最优策略,求你找到M所需的时间的期望.

有$t(1\le t\le 1000)$组测试数据.每组测试数据输入一个整数$n(2\le n\le 100)$,表示节点数.接下来$(n-1)$行每行输入两个数$a,b(1\le a,b\le n)$,表示节点$a$与$b$间有边.数据保证输入的信息构成一棵树.

对每组测试数据,输出你采取最优策略时找到M所需的时间的期望,误差不超过$1\mathrm{e}-9$.

思路

最优策略是走树的Euler序,即不走到底不回头.

[证] 若不然,设(遍历时,即深度增加时)你从节点$u$走到节点$v$(非叶子节点)后又走回节点$u$.

(1)走到节点$u$时游戏未结束,说明M不在节点$u$处.走到节点$v$时游戏未结束,说明M不在节点$v$处.

​ 此时若从$v$返回$u$,则需多花$1\mathrm{s}$检查一个已经确认没有M的节点,显然这不是最优的.

(2)若M在节点$v$的子树中,显然继续遍历$v$的子树比先回到$u$再遍历$u$的子树更优.

(3)若M不在节点$v$的子树中(这需遍历完$v$的子树才能确定),又M不在已遍历过的节点中,则需从$v$回溯到$u$,继续遍历其他子树.

​ 显然这种遍历顺序是Euler序.

[注] 由(2)(3)知:交换两子树的遍历顺序,所需的时间的期望不变.

对M在$i(2\le i\le n)$节点的情况,只需在DFS过程中更新节点出入栈的时间即可得到Euler序,将Euler序累加到答案上,最后除以$(n-1)$即可.

代码 -> 2021ICPC济南-K(Euler序+期望)

const int MAXN = 105;
int n;  // 节点数
vi edges[MAXN];
double ans = 0;
int tim;  // 时间

void dfs(int u, int fa) {
	tim++;  // u入栈
	for (auto v : edges[u]) {
		if (v == fa) continue;

		ans += tim;  // 将Euler序累加到答案上
		dfs(v, u);
	}
	tim++;  // u出栈
}

int main() {
	CaseT{
		cin >> n;
		ans = tim = 0;  // 清空
		for (int i = 1; i <= n; i++) edges[i].clear();

		for (int i = 0; i < n - 1; i++) {
			int u, v; cin >> u >> v;
			edges[u].push_back(v), edges[v].push_back(u);
		}

		dfs(1, -1);  // 从起点开始搜,起点无前驱节点
		cout << fixed << setprecision(18) << ans / (n - 1) << endl;
	}
}

---

---

C. Optimal Strategy

题意

有编号$1\sim n$的$n(1\le n\le 1\mathrm{e}6)$个物品,其中第$i(1\le i\le n)$个物品的价值为$a_i(1\le a_i\le n)$.A和B两人轮流操作,A先手.当前轮到的人在剩余的物品中取走一个物品.要求当所有物品取完时,两人各自手中的物品的价值之和都最大.两人都采取最优策略,求有多少种游戏进行序列(每轮玩家取的物品的编号构成的序列),答案对$998244353$取模.

思路

常见误区:最优策略是每次都拿剩下的物品中价值最大的物品.对样例$1$,若应拿价值最大的物品,则先手不应拿$1$,与样例解释不符.事实上,要使得游戏结束时两人各自手中的物品价值之和都最大,只需保证两人都拿到取得自己物品价值之和的最大值应拿的那些物品即可,与什么时候拿到那个物品无关.

显然游戏的结果只有先手必胜和平局两种情况,但这题不是博弈论,而是纯纯的组合计数.

具体地考察最优策略中两人分别应拿到哪些数.设$cnt[i]$表示价值为$i$的数的个数.①$cnt[i]$为偶数时,显然A和B各拿一半;②$cnt[i]$为奇数时,A先手会多拿走一个,转化为B先手且$cnt[i]$为偶数的情况.综上,对$cnt[i]$个数$i$,只需考虑其中$\lfloor \frac{cnt[i]}{2}\rfloor$个的分配情况.

注意到$1\le a_i\le 1\mathrm{e}6$,不妨按顺序枚举每个数值的分配情况,而两人何时拿到自己应拿的物品不影响最终两人各自手中的物品的价值之和,故从小到大或从大到小枚举每个数值的分配情况都可以,下面以从小到大枚举为例.设当前准备分配数$i$,则前面的$1\sim (i-1)$共$sum=\sum _{j=1}^{i-1}cnt[j]$个数都已分配好,只需在前$\left(sum+\lfloor \frac{cnt[i]}{2}\rfloor \right)$个数中选$\lfloor \frac{cnt[i]}{2}\rfloor$个给先手即可,方案数为$C_{sum+\lfloor cnt[i]/2\rfloor }^{\lfloor cnt[i]/2\rfloor }$.因以不同顺序拿到相同数值的数视为不同方案,故$ans=\sum _{i=1}^{n}cnt[i]!\cdot C_{sum+\lfloor cnt[i]/2\rfloor }^{\lfloor cnt[i]/2\rfloor }$.直接计算时间复杂度为$O(n^2)$,注意到$sum$是$cnt[j]$的前缀和,只需在转移的同时维护即可,时间复杂度$O(n)$.

代码 -> 2021ICPC济南-C(组合计数+推公式+前缀和)

const int MAXN = 1e6 + 5;
const int MOD = 998244353;
int n;
int cnt[MAXN];  // cnt[i]表示值为i的数的个数
int fac[MAXN], ifac[MAXN];

void init() {  // 预处理阶乘及其逆元
	fac[0] = 1;
	for (int i = 1; i < MAXN; i++) fac[i] = (ll)fac[i - 1] * i % MOD;

	ifac[MAXN - 1] = qpow(fac[MAXN - 1], MOD - 2, MOD);
	for (int i = MAXN - 1; i; i--) ifac[i - 1] = (ll)ifac[i] * i % MOD;
}

int C(int n, int m) {  // 组合数C(n,m)
	return (ll)fac[n] * ifac[m] % MOD * ifac[n - m] % MOD;
}

int main() {
	init();

	cin >> n;
	for (int i = 0; i < n; i++) {
		int a; cin >> a;
		cnt[a]++;
	}

	int ans = 1;  // 注意初始值为1而不是0
	int pre = 0;  // cnt[i]的前缀和
	for (int i = 1; i <= n; i++) {
		if (!cnt[i]) continue;

		ans = (ll)ans * fac[cnt[i]] % MOD * C(pre + cnt[i] / 2, cnt[i] / 2) % MOD;
		pre += cnt[i];
	}
	cout << ans;
}

---

---

D. Arithmetic Sequence

题意

给定一个长度为$n(1\le n\le 2\mathrm{e}5)$的数列$a_1,\cdots ,a_n(0\le |a_i|\le 1\mathrm{e}13)$.现有两种操作:①选择数列中的一个数$+1$;②选择数列中的一个数$-1$.求将该数列变为等差数列的最小操作次数.

思路 by : keguaiguai

设最优解得到的数列为$b_1,\cdots ,b_n$,其公差为$d$,则$b_i=b_1+(i-1)d$,所需的操作次数为$\sum _{i=1}^n|a_i-b_i|=\sum _{i=1}^n\left|[a_i-(i-1)d]-b_1\right|\stackrel{c_i=a_i-(i-1)d}{}\sum _{i=1}^n|c_i-b_1|$,易想到货舱选址问题.注意到此时未知$c_i$的值,而$c_i$的值可由$a_i$和$d$确定,考虑先求出$d$.

设将序列变为以$d$为公差的等差数列所需的最小操作次数为$f(d)$,它是一个下凸函数.

[证] 注意到$\left|[a_i-(i-1)(d+1)]-b_1\right|+\left|[a_i-(i-1)(d-1)]-b_1\right|$

​ $=|[a_i-(i-1)d-b_1]-(i-1)|+|[a_i-(i-1)d-b_1]+(i-1)|\ge 2|[a_i-(i-1)d]-b_1|$,

两边求和得$f(d+1)+f(d-1)\ge 2f(d)$,故证.

注意到下凸函数是单谷的,可三分出其最小值点.注意每次三分都要求出序列的中位数,用sort可能会TLE,故改用nth_element函数在$O(n)$的时间复杂度内求出序列的中位数.

注意到$n$最大$2\mathrm{e}5,a_i$最大$1\mathrm{e}13$,则累加后可能会爆ll,故要开__int128.注意对__int128要重载abs函数.

代码 -> 2021ICPC济南-D(三分+nth_element)

#define ll __int128

const int MAXN = 2e5 + 5;
int n;
ll a[MAXN];  // 原数列
ll b[MAXN];  // 目标数列

ll abs(ll x) { return x < 0 ? -x : x; }  // 对__int128重载abs()

ll get_ans(ll d) {
	for (int i = 1; i <= n; i++) b[i] = a[i] - (i - 1) * d;

	ll res = 0;
	nth_element(b + 1, b + (n + 1) / 2, b + n + 1);  // 将中位数调整到下标(n+1)/2处
	ll mid = b[(n + 1) / 2];  // 中位数

	for (int i = 1; i <= n; i++) res += abs(a[i] - mid - (i - 1) * d);
	return res;
}

int main() {
	read(n);
	for (int i = 1; i <= n; i++) read(a[i]);

	ll l = -2e13, r = 2e13;
	while (l < r) {
		ll lmid = l + (r - l) / 3, rmid = r - (r - l) / 3;
		if (get_ans(lmid) < get_ans(rmid)) r = rmid - 1;
		else l = lmid + 1;
	}
	write(get_ans(l));
}

---

---

J. Determinant

题意

$n$阶矩阵$A$有性质:$\det (A^{T}A)=(\det A{)}^2$.给定$n$阶矩阵$A$和$|\det A|$,判断$\det A$的正负.

有$t(1\le t\le 100)$组测试数据.每组测试数据第一行输入整数$n(1\le n\le 100)$.第二行输入一个长度不超过$1\mathrm{e}4$的大数表示$|\det A|$.接下来输入$n\times n$的矩阵$A$,其中元素的绝对值不超过$1\mathrm{e}9$.

对每组测试数据,若$\det A>0$,输出’+‘;否则输出’-'.

思路 by : Paranoid5

[引理] 对$\forall $非零整数$x_1$及其相反数$x_2$,对奇模数$p$,有$x_1\not\equiv x_2\ \ (\mathrm{mod}\ p)$,且$x_1\ \mathrm{mod}\ p\neq 0,x_2\ \mathrm{mod}\ p\neq 0$.

[证] $x_1+x_2=0$,则$(x_1+x_2)\mathrm{mod}p=0$,进而$(x_1\mathrm{mod}p+x_2\mathrm{mod}p)\mathrm{mod}p=0$.

设$t_i=x_i\mathrm{mod}p(i=1,2)$.若$x_1\equiv x_2(\mathrm{mod}p)$,则$t_1=t_2$,且 t 1 , t 2 < p t_1,t_2

t1​,t2​<p.

​ 此时$t_1+t_2$为偶数,且$0<t_1+t_2<2p$.

因$(t_1+t_2)\mathrm{mod}p=0$,则$t_1+t_2=kp(k\in \mathbb{Z})$,结合$0<t_1+t_2<2p$知$k\in (0,2)$,即$k=1$.

​ 此时$t_1+t_2=p$,$\mathrm{LHS}$为偶数,$\mathrm{RHS}$为奇数,矛盾.

故$x_1\not\equiv x_2(\mathrm{mod}p)$.

$|\det A|$对大素数$MOD$取模,用Gauss消元求$\det A$在模$MOD$意义下的值,比较两者是否相等即可.常见的模数可取$1\mathrm{e}9+7$.

代码 -> 2021ICPC济南-J(Gauss消元+大数取模)

const int MAXN = 605;
const int MOD = 1e9 + 7;
int n;
int a[MAXN][MAXN];

int gauss() {
	int sgn = 1;  // 行列式的值的符号
	for (int i = 1; i <= n; i++) {
		for (int j = i + 1; j <= n; j++) {
			while (a[i][i]) {
				int d = a[j][i] / a[i][i];
				for (int k = i; k <= n; k++)
					a[j][k] = ((ll)a[j][k] - (ll)d * a[i][k] % MOD + MOD) % MOD;
				swap(a[i], a[j]);
				sgn *= -1;
			}
			swap(a[i], a[j]);
			sgn *= -1;
		}
	}
	
	int res = 1;
	for (int i = 1; i <= n; i++) res = (ll)res * a[i][i] % MOD;
	return (res * sgn + MOD) % MOD;
}

int main() {
	CaseT{
		cin >> n;
		string s; cin >> s;
		int ans = 0;
		for (auto ch : s) ans = ((ll)ans * 10 % MOD + ch - '0') % MOD;

		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= n; j++) cin >> a[i][j];

		if (ans == gauss()) cout << '+' << endl;
		else cout << '-' << endl;
	}
}

---

---

L. Strange Series

题意 ($2\mathrm{s}$)

给定一个$n$次多项式$P(x)=a_0+a_{1}x+\cdots +a_n{x}^n$,其中$a_i\ge 0(i=0,\cdots ,n)$.考察级数$S=\sum _{m=0}^{\infty }\frac{P(m)}{m!}$.可以证明$S$是良好定义的,且$S=p\mathrm{e}$,其中$p$是一个非负整数,$\mathrm{e}=\sum _{m=0}^{\infty }\frac{1}{m!}\approx 2.71828$.求$p$.

有$t(1\le t\le 100)$组测试数据.每组测试数据输入一个整数$n(0\le n\le 1\mathrm{e}5)$.第二行输入$(n+1)$个整数$a_0,a_1,\cdots ,a_n(0\le a_i<998244353)$.数据保证至多有$4$组测试数据满足$n>1000$.

对每组测试数据,输出$p$对$998244353$取模的结果.

思路

[引理] $\sum _{m=0}^{\infty }\frac{m^n}{m!}{x}^m={\mathrm{e}}^x\cdot Q(x)$,其中$Q(x)$是关于$x$的$n$次多项式.

[证] 用数学归纳法证明.不妨设$n-1$时结论成立.

$\sum _{m=0}^{\infty }\frac{1}{m!}\cdot m^n{x}^m=x\sum _{m=0}^{\infty }\frac{1}{m!}\cdot m^{n-1}\cdot m{x}^{m-1}=x\sum _{m=0}^{\infty }\frac{1}{m!}\cdot m^{n-1}\cdot \frac{\mathrm{d}}{\mathrm{d}x}(x^m)=x\cdot \frac{\mathrm{d}}{\mathrm{d}x}\sum _{m=0}^{\infty }\frac{1}{m!}\cdot m^{n-1}{x}^m$

​ $=x\cdot \frac{\mathrm{d}}{\mathrm{d}x}[{\mathrm{e}}^{x}Q(x)]$,其中$Q(x)$是关于$x$的$(n-1)$次多项式.

注意到对${\mathrm{e}}^{x}Q(x)$求导后多项式部分的最高次项仍为$x^{n-1}$,乘上前面的$x$即证.

$P(x)=a_0+a_{1}x+\cdots +a_n{x}^n(a_i\ge 0),S=\sum _{m=0}^{\infty }\frac{P(m)}{m!}=p\mathrm{e}$,则$p=\frac{1}{\mathrm{e}}S=\frac{1}{\mathrm{e}}\underset{s\to \infty }{\lim }\sum _{m=0}^s\frac{P(m)}{m!}$.

因$\frac{P(m)}{m!}=\frac{1}{m!}(a_0+a_{1}m+\cdots +a_n{m}^n)=\sum _{i=0}^n{a}_i\frac{k^i}{k!}$,

​ 则$p=\frac{1}{\mathrm{e}}\underset{s\to \infty }{\lim }\sum _{m=0}^s\sum _{i=0}^n{a}_i\frac{m^i}{m!}=\sum _{i=0}^n{a}_i\left(\frac{1}{\mathrm{e}}\underset{s\to \infty }{\lim }\sum _{m=0}^s\frac{m^i}{m!}\right)=\sum _{i=0}^n{a}_i\left[\frac{1}{\mathrm{e}}\sum _{m=0}^{\infty }\frac{m^i}{m!}\right]$.

由引理知:$\sum _{m=0}^{\infty }\frac{m^i}{m!}{x}^m={\mathrm{e}}^x\cdot Q(x)$,其中$Q(x)$是关于$x$的$i$次多项式,进而$p=\sum _{i=0}^n{a}_i\cdot Q(1)$,同时也证明了$p\in \mathbb{N}$.

下面讨论如何求$Q(x)$的系数.

$i=0$时,$\sum _{m=0}^{\infty }\frac{1}{m!}{x}^m={\mathrm{e}}^x\cdot 1$,此时$Q(x)=1,Q(1)=1$.

$i=1$时,$\sum _{m=1}^{\infty }\frac{m}{m!}{x}^m={\mathrm{e}}^x\cdot x$,此时$Q(x)=x,Q(1)=1$.

$i=2$时,$\sum _{m=1}^{\infty }\frac{m^2}{m!}{x}^m={\mathrm{e}}^x\cdot x(x+1)$,此时$Q(x)=x(x+1),Q(1)=2$.

$i=3$时,$\sum _{m=1}^{\infty }\frac{m^3}{m!}{x}^m={\mathrm{e}}^x\cdot x(x^2+3x+1)$,此时$Q(x)=x(x^2+3x+1),Q(1)=5$.

$i=4$时,$\sum _{m=1}^{\infty }\frac{m^4}{m!}{x}^m={\mathrm{e}}^x\cdot x(x^3+6x^2+7x+1)$,此时$Q(x)=x(x^3+6x^2+7x+1),Q(1)=15$.

$i=5$时,$\sum _{m=1}^{\infty }\frac{m^5}{m!}{x}^m={\mathrm{e}}^x\cdot x(x^4+10x^3+25x^2+15x+1)$,

​ 此时$Q(x)=x(x^4+10x^3+25x^2+15x+1),Q(1)=52$.

易联想到Bell数:$B_0=1,B_1=1,B_2=2,B_3=5,B_4=15,B_5=52,\cdots$,它可用递推求得.

注意到$n$最大为$1\mathrm{e}5$,暴力递推出Bell数会TLE.考虑用Bell数的生成函数求Bell数.

Bell数的定义:将$n$个不同的元素的集合划分为任意多个非空子集的方案数.

设$n$个元素的方案数为$f_n$.固定最后一个元素$x$,枚举最后一个元素所在的集合的元素个数$i\in [1,n]$,从剩下的$(n-1)$个元素中选$(i-1)$个元素与$x$在同一集合中,其他$(n-i)$个元素任意划分,故$f_n=\sum _{i=1}^n{C}_{n-1}^{i-1}\cdot f_{n-i}$,调整下标得$f_{n+1}=\sum _{i=0}^n{C}_n^i\cdot f_{n-i}$.

设$f$的EGF为$F(x)$,则$\frac{f_{n+1}}{(n+1)!}=\sum _{i=0}^n\frac{f_{n-i}}{(n-i)!}\cdot \frac{1}{i!(n+1)}$,进而$(n+1)\cdot \frac{f_{n+1}}{(n+1)!}=\sum _{i=0}^n\frac{f_{n-i}}{(n-i)!}\cdot \frac{1}{i!}$,

​ 即$[x^n]F’(x)=\sum _{i=0}^n[x^{n-i}]F(x)\cdot [x^i]{\mathrm{e}}^x$,亦即$F^{\prime }(x)=F(x){\mathrm{e}}^x$.

令$y=F(x)$,则$\frac{\mathrm{d}y}{\mathrm{d}x}={\mathrm{e}}^{x}y$,移项得$\frac{\mathrm{d}y}{y}={\mathrm{e}}^x\mathrm{d}x$,两边积分得$\ln y={\mathrm{e}}^x+C$,故$y=\exp ({\mathrm{e}}^x+C)$.

$n=1$时只有一种划分方案,代入得$C=-1$,故$F(x)=\exp ({\mathrm{e}}^x-1)$.多项式exp求即可,时间复杂度$O(n\log n)$.

代码 -> 2021ICPC济南-L(Bell数+生成函数)

const int MAXN = 6e5 + 5;
const int MOD = 998244353;

int rev[MAXN];  // rev[x]表示二进制数x镜像后的结果

void NTT(int* y, int lim, int op) {  // op=1时为NTT,op=-1时为INTT
	for (int i = 0; i < lim; i++)  // 做一遍BRP,将系数位置调整为递归后的位置
		if (rev[i] < i) swap(y[i], y[rev[i]]);

	for (int l = 2; l <= lim; l <<= 1) {  // 模拟合并过程,l为当前区间长度
		int g = qpow(3, (MOD - 1) / l, MOD);
		for (int i = 0; i < lim; i += l) {
			int cur = 1;
			int k = l >> 1;
			for (int j = 0; j < k; j++, cur = (ll)cur * g % MOD) {
				int tmp = (ll)y[i + j + k] * cur % MOD;
				y[i + j + k] = (y[i + j] - tmp + MOD) % MOD;  // 奇次项
				y[i + j] = ((ll)y[i + j] + tmp) % MOD;  // 偶次项
			}
		}
	}

	if (op == -1) {  // INTT
		reverse(y + 1, y + lim);
		int inv = qpow(lim, MOD - 2, MOD);
		for (int i = 0; i < lim; i++) y[i] = (ll)y[i] * inv % MOD;
	}
}

void PolyMul(int* a, int* b, int n, int m) {  // n次多项式a乘m次多项式b,结果放在a中
	int lim = 1, len = 0;
	while (lim < (n + m << 1)) lim <<= 1, len++;
	for (int i = 1; i < lim; i++) rev[i] = (rev[i >> 1] >> 1 | ((i & 1) << (len - 1)));

	NTT(a, lim, 1), NTT(b, lim, 1);
	for (int i = 0; i < lim; i++) a[i] = (ll)a[i] * b[i] % MOD;
	NTT(a, lim, -1);
}

int n;
int f[MAXN], derf[MAXN], invf[MAXN];  // f、f的导数、f的逆
int bell[MAXN];

int tmp[MAXN];  // 递归过程中临时存放系数的数组

void PolyInv(int deg, int* a, int* res) {
	if (deg == 1) {
		res[0] = qpow(a[0], MOD - 2, MOD);
		return;
	}

	PolyInv((deg + 1) >> 1, a, res);  // 递归求解

	int lim = 1, len = 0;
	while (lim < (deg << 1)) lim <<= 1, len++;

	for (int i = 1; i < lim; i++) rev[i] = (rev[i >> 1] >> 1 | ((i & 1) << (len - 1)));

	for (int i = 0; i < lim; i++) tmp[i] = a[i] * (i < deg);  // i≥deg的项在模意义下为0

	NTT(tmp, lim, 1), NTT(res, lim, 1);
	for (int i = 0; i < lim; i++)
		res[i] = (ll)(2 - (ll)tmp[i] * res[i] % MOD + MOD) % MOD * res[i] % MOD;
	NTT(res, lim, -1);

	for (int i = deg; i < lim; i++) res[i] = 0;  // i≥deg的项在模意义下为0
}

void PolyDer(int deg, int* a, int* res) {  // 对deg次多项式a求导,结果放在res中
	for (int i = 1; i < deg; i++) res[i - 1] = (ll)i * a[i] % MOD;
	res[deg - 1] = 0;
}

void PolyInt(int deg, int* a, int* res) {  // 对deg次多项式a求不定积分(C=0),结果放在res中
	res[0] = 0;
	for (int i = 1; i < deg; i++) res[i] = (ll)a[i - 1] * qpow(i, MOD - 2, MOD) % MOD;
}

void PolyLn(int n, int* a, int* res) {  // 对n次多项式a求ln,结果放在res中
	memset(derf, 0, so(derf));
	memset(invf, 0, so(invf));

	int deg = 1;
	while (deg <= n) deg <<= 1;

	PolyDer(deg, a, derf);  // a的导数
	PolyInv(deg, a, invf);  // a的逆
	PolyMul(derf, invf, deg, deg);
	PolyInt(deg, derf, res);
}

int lnres[MAXN];  // res的ln
int tmpa[MAXN];  // 递归过程中临时存放系数的数组

void PolyExp(int deg, int* a, int* res) {  // 对deg次多项式a求exp,结果放在res中
	if (deg == 1) {
		res[0] = 1;
		return;
	}

	PolyExp(deg + 1 >> 1, a, res);

	PolyLn(deg, res, lnres);

	int lim = 1, len = 0;
	while (lim < (deg << 1)) lim <<= 1, len++;
	for (int i = 1; i < lim; i++) rev[i] = (rev[i >> 1] >> 1 | ((i & 1) << (len - 1)));

	for (int i = 0; i < deg; i++) tmpa[i] = a[i];
	for (int i = deg; i < lim; i++) tmpa[i] = lnres[i] = 0;
	for (int i = 0; i < lim; i++) tmpa[i] = (((ll)tmpa[i] - lnres[i]) % MOD + MOD) % MOD;
	tmpa[0]++;

	NTT(tmpa, lim, 1), NTT(res, lim, 1);
	for (int i = 0; i < lim; i++) res[i] = (ll)res[i] * tmpa[i] % MOD;
	NTT(res, lim, -1);
	for (int i = deg; i < lim; i++) res[i] = 0;
}

int fac[MAXN], ifac[MAXN];

void init() {  // 预处理Bell数
	fac[0] = 1;
	for (int i = 1; i < MAXN; i++) fac[i] = (ll)fac[i - 1] * i % MOD;

	ifac[MAXN - 1] = qpow(fac[MAXN - 1], MOD - 2, MOD);
	for (int i = MAXN - 1; i; i--) ifac[i - 1] = (ll)ifac[i] * i % MOD;

	n = 1e5 + 1;
	for (int i = 1; i < n; i++) f[i] = ifac[i];

	PolyExp(n, f, bell);

	for (int i = 0; i < n; i++) bell[i] = (ll)bell[i] * fac[i] % MOD;
}

int main() {
	init();  // 预处理Bell数
	
	CaseT{
		int n; cin >> n;
		int ans = 0;
		for (int i = 0; i <= n; i++) {
			int a; cin >> a;
			ans = ((ll)ans + (ll)a * bell[i]) % MOD;
		}
		cout << ans << endl;
	}
}

---

---