CF1349F Slime and Sequences （欧拉数，扩展拉格朗日反演）

最新推荐文章于 2022-12-26 21:24:25 发布

Freopen

最新推荐文章于 2022-12-26 21:24:25 发布

阅读量709

点赞数

分类专栏：组合计数 FFT 多项式

本文链接：https://blog.csdn.net/qq_35950004/article/details/106372885

版权

这篇博客介绍了如何将好序列与排列进行映射，并探讨了具体数学中的欧拉数及其应用。通过一系列数学推导，解释了如何利用拉格朗日反演解决序列计数问题，涉及公式展开、二项式反演和扩展拉格朗日反演等概念，最终给出了解决问题的算法思路。

摘要由CSDN通过智能技术生成

这篇博客大部分抄自此人

好序列为对于一个长度为 $n$ 的序列，其中若有 $a_i= k(k>1)$ ，则必存在 $a_{j\in[1,i-1]} = k-1$
求对于所有长度为 $n$ 的好序列，分别求出 $k = 1 . . . n$ 的出现次数和。

第一步神仙映射：
把所有的好序列映射到排列中，发现可以一一对应。
对应方法如下：
排列到好序列：
对于一个排列 $p_i$ ，如果 $p_1,p_2...p_i$ 之间有 $k$ 个升高， $a_{p_i} = k+1$ 。
$p_{i-1}$ 和 $p_i$ 之间有一个升高指的是 $p_i > p_{i-1}$ 。
（这里的定义沿用具体数学欧拉数板块的定义）。
好序列到排列：
考虑好序列中的 $1$ ，一定是放在最前面，所以可以知道对于都是 $a_{p_i} = 1$ 的数，假设有 $k$ 个，他们一定是占据了 $p_{1...k}$ ，并且都满足 $p_{i} > p_{i+1}$ 。
再考虑 $2$ ，可以知道 $2$ 一定是摆在 $1$ 后面并且在 $a_{p_i} = 2$ 中都满足 $p_i > p{i+1}$ 。
那么我们就得到了一个好序列到排列，排列到好序列的映射，不同的排列都能映射到不同的好序列，不同的好序列都能映射到不同的排列，他们数量相同，他们是一一对应。

翻阅具体数学可以得到：
有 $k$ 个升高的长度为 $n$ 的序列数表示为 $\left \langle\begin{aligned} n \\ k \end{aligned} \right\rangle$
那么根据具体数学的描述有
$\left \langle\begin{aligned} n \\ k \end{aligned} \right\rangle = (k+1)\left \langle\begin{aligned} n&-1 \\ & \ k \end{aligned} \right\rangle + (n-k)\left \langle\begin{aligned} n-1 \\ k-1 \end{aligned} \right\rangle$
具体来说就是把 $n$ 插入一个 $1 . . n - 1$ 的排列然后分类讨论即可得到该式子。

所以 $Ans_{i+1} = \sum \left \langle\begin{aligned} i \\ k \end{aligned} \right\rangle \binom{n}{i} (n-i)!$
所以可以先写个 $e a s y$ 版的 $O(n^2)$ 压压惊。
$\mathcal AC\ Code$

#include<bits/stdc++.h>
#define maxn 5005
#define mod 998244353
using namespace std;

int n;
int C[maxn][maxn],E[maxn][maxn],fac[maxn];

int main(){
   
	scanf("%d",&n);
	for(int i=C[0][0]=E[0][0]=1;i<=n;i++)
		for(int j=C[i][0]=E[i][0]=1;j<=i;j++){
   
			C[i][j] = (C[i-1][j-1] + C[i-1][j]) % mod;
			E[i][j] = (1ll * (j+1) * E[i-1][j] + 1ll * (i-j) * E[i-1][j-1]) % mod;
		}
	fac[0] = 1;
	for(int i=1;i<=n;i++) fac[i] = 1ll * fac[i-1] * i % mod;
	for(int i=1;i<=n;i++){
   
		int ans = 0;
		for(int j=i;j<=n;j++) 
			ans = (ans + 1ll * E[j][i-1] * C[n][j] % mod * fac[n-j]) % mod;
		printf("%d%c",ans," \n"[n==i]);
	}
}

接下来是令人目不接暇的炫酷推式子环节了。
~~有点长~~

~~先去掉那我打不来latex的欧拉数~~
设
$\sum_{k=v}^{i-1}\left \langle\begin{aligned} i \\ k \end{aligned} \right\rangle \binom kv$
表示钦定 $v$ 个间隙中是升高的方案数对每种钦定方案的和。

那么可以通过二项式反演得到:
$\left \langle\begin{aligned} i \\ v \end{aligned} \right\rangle = \sum_{k=v}^{i-1} F(i,k)\binom kv(-1)^{k-v}$
然后带入之前的式子：
$Ans_ k = \sum_{i=k+1} \binom ni (n-i)!\sum_{j=k}^{i-1} F(i,j)\binom jk(-1)^{j-k}$
考虑

最低0.47元/天解锁文章