2018-2019 ICPC, NEERC, Northern Eurasia Finals 部分题解

$I$
感觉这题非常厉害
为了方便解释，做以下约定
$1.$将全序列$p_1,p_2...p_n$记为$A$
$2.$将$A$的某个严格子区间记作$B_i$

题目把 interval 定义为值域大小和区间大小相等的$B_i$
我们定义maximal interval为不被任意一个interval完整包含的interval

有一个重要结论：如果两个maximal interval $B_x,B_y$相交，那么区间$B_x\cup B_y$一定符合interval的定义，又因为$B_x,B_y$是maximal interval，所以若他们相交，那么并一定就是$A$

考虑算出$n$的全排列里所有不合法的排列
依据这个结论，我们可以将不合法的排列分为两类
类型$1：$存在两个maximal interval $B_u,B_v$，使$B_u\cup B_v=A$
类型$2：$A为大于等于3个maximal interval挨着连在一起组成

我们先考虑类型$1$，它可以分为$B_u$和$B_v$相交或相邻两种情况
若$B_u$和$B_v$相交，记
$B_w=B_u\cap B_v,B_{u^{\prime }}=B_u-B_w,B_{v^{\prime }}=B_v+B_w$，
转化完后得到的$B_{u^{\prime }}和B_{v^{\prime }}$是相邻的，因此我们只要考虑$B_u,B_v$相邻的情况

显然可以发现$B_u,B_v$的值域分别为$[1,k],[k+1,n]$（证明可以用反证）
不妨假设$B_u$在$B_v$的左侧（右侧的情况方案数是相同的）

为了避免重复计数，我们强制令$B_u$的prefix中除了自己外，没有任何一个prefix满足值域为$[1,k]$且为interval （若存在这样的prefix，可以通过上面那个相交变相邻的转化转为不存在）

那么就可以统计数量了，记$f_n$为$n$的所有排列$A$中满足 任何一个prefix都不满足值域为$[1,k]$且为interval 的排列数量
有$$f_n=n!-\sum _{i=1}^{n-1}(f_i\ast (n-i)!)$$

接着考虑类型$2$
对于某个maximal interval $B_u$，他内部是任意排的，我们唯一要满足的限制是不存在某一段连续的maximal interval可以拼成一个interval
这个限制相当于这些maximal interval的排列是interval-free的
由于我们统计的是非法方案，因此至少存在一个maximal interval $B_v$的大小是>=2的，于是这些maximal interval的个数<n，就化成了子问题
为了计算这种的方案数，需要一个辅助的dp

记$g_{n,k}$为将$n$划分成$k$段，每段内部任意排列的方案数
有$g_{n,k}={\sum }_{i=k-1}^{n-1}{g}_{i,k-1}\ast (n-i)!$

然后就可以算答案了，记$an{s}_n$为$n$的所有排列$A$中interval-free的排列数
$$an{s}_n=n!-2\sum _{i=1}^{n-1}{f}_i\ast (n-i)!-\sum _{k=3}^{n-1}{g}_{n,k}\ast an{s}_k$$

maya我的语言表达太差了，写的题解长度可能是原题解的2倍…

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 405;

int n,mod;
inline void add(int &a,const int &b){ a+=b;if(a>=mod)a-=mod; }
inline void dec(int &a,const int &b){ a-=b;if(a<0) a+=mod; }

int s[maxn];
int f[maxn],ans[maxn];
int g[maxn][maxn];

int main()
{
	//freopen("tmp.in","r",stdin);
	//freopen("tmp.out","w",stdout);
	
	int Tcase; scanf("%d%d",&Tcase,&mod);
	
	s[0]=1; for(int i=1;i<maxn;i++) s[i]=(ll)s[i-1]*i%mod;
	f[1]=1;
	for(int i=2;i<maxn;i++)
	{
		f[i]=s[i];
		for(int j=1;j<i;j++) dec(f[i],(ll)f[j]*s[i-j]%mod);
	}
	g[0][0]=1;
	for(int i=1;i<maxn;i++) for(int j=1;j<=i;j++)
	{
		for(int k=j-1;k<i;k++) add(g[i][j],(ll)g[k][j-1]*s[i-k]%mod);
	}
	
	ans[1]=1;
	for(int i=2;i<maxn;i++)
	{
		ans[i]=s[i];
		for(int j=1;j<i;j++) if(i-j>1||j>1)
			dec(ans[i],2ll*f[j]%mod*s[i-j]%mod);
		for(int k=3;k<i;k++) dec(ans[i],(ll)g[i][k]*ans[k]%mod);
	}
	
	while(Tcase--)
	{
		scanf("%d",&n);
		printf("%d\n",ans[n]);
	}
	
	return 0;
}

先写这么多剩下的心情好的时候填坑