6516. 【GDOI2020模拟03.28】数二数（two）

题目

有一个$[1,n]$的整数，可以询问$[L,R]$，表示整个整数是否在这个区间里。
计算有多少个询问集合，使得这些询问过后，无论整数是$[1,n]$中哪一个都能被唯一确定。
$n\le 300$

思考历程

简化一下题目大意：可以通过询问，将数字分成若干组。
一开始$[1,n]$为一组。
某个询问$[L,R]$之后，将这个东西与当前分成的若干组分别取交，从而进一步分组。
到最后如果每一组如果仅仅包含一个数字，那么就是一种合法的方案。

我只会强行状压DP……然而并没有这档部分分。

---

正解

正解是个玄妙的东西。官方题解在此。

分组的更加严谨的解释：如果两个数字被包含的询问集合完全相同，那么它们就被分到同一组。
给组别分配一个编号，令$V_i$表示数字$i$在哪个组。
如果$V_i$互不相同，那么这就是一个合法的方案；如果有至少一个相同，那么就不合法。

考虑用总方案减去不合法的方案。
有个奇妙的性质：如果存在四个数字$a,b,c,d$满足$a<b<c<d$且$V_a=V_c,V_b=V_d$，那么一定有$V_a=V_b=V_c=V_d$
用人话说就是，每个组别的最左端和最右端形成的区间，不可能真相交。
题解用一种奇妙的方法将$V$缩短。
假如$V=[1,2,3,2,4,2,5,5,10,6,7,8,7,6,9]$
可以通过以下方式构造出$V^{\prime }$：
在$V$中从左往右扫，见到一个组别号$x$，首先将$x$丢到$V^{\prime }$的后面，然后找出$l_x$和$r_x$分别表示类别为$x$的最左端和最右端，将$[l_x,r_x]$从$V$中删除。
这个例子构造出来的$V^{\prime }=[1,2,5,10,6,9]$

至于这个东西有什么用，我先放出题解原文：
You might be wondering why is this helpful? Well, it’s helpful because any set of questions that corresponds to $V$ should respect the following:

1. Some questions are contained by a removed subarray. There is basically no restriction here
2. All the other questions should uniquely identify $V^{\prime }$ - we have the first hint of dynamic programming

本人英语不好，解释得不好也请见谅……
要搞出$V$和$V^{\prime }$之间的联系，先考虑能得到$V$的询问集合。

1. 部分的询问区间是已经被$V$变成$V^{\prime }$的过程中，删去的子区间所包含的。这些询问基本上对$V^{\prime }$没有限制。
2. 除去那些对$V^{\prime }$没有用的询问，剩余询问集合应该唯一地确定$V^{\prime }$。

那么可以从$V^{\prime }$推到$V$，限制只在新增的那些子区间中。
设$f_i$表示$[1,i]$合法的方案数。考虑用总数减去不合法的。
$g_{i,j}$表示$V$总长为$i$，经过缩减处理的$V^{\prime }$总长为$j$的方案数。这个东西就是$V^{\prime }$到$V$新增的限制。
$$f_i=2^{C_{i+1}^2}-\sum _{j=1}^{i-1}{f}_j{g}_{i,j}$$
至于$g_{i,j}$的转移，考虑在$V^{\prime }$后面新增一个组别编号。分为两种情况：这个组别仅仅出现一次；这个组别出现了至少两次，枚举左端和右端之间空隙的长度$k$，在这个长度为$k$的区间里的询问自由组合。
$$g_{i,j}=g_{i-1,j-1}+\sum _{k=0}{g}_{i-k-2,j-1}{2}^{C_{k+1}^1}$$
时间复杂度$O(n^3)$。

在网上查题解的时候还发现这个方法可以用拉格朗日反演来优化，时间复杂度$O(n\lg n)$的。
表示没有看懂，以后心血来潮就来这里看看。

---

代码

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#define N 310
int n,mo;
ll p2[N*N];
ll f[N],g[N][N];
int main(){
//	freopen("in.txt","r",stdin);
	freopen("two.in","r",stdin);
	freopen("two.out","w",stdout);
	scanf("%d%d",&n,&mo);
	p2[0]=1;
	for (int i=1;i<=(n+1)*n>>1;++i)
		p2[i]=p2[i-1]*2%mo;
	g[0][0]=1;
	for (int i=1;i<=n;++i)
		for (int j=1;j<=i;++j){
			ll s=g[i-1][j-1];
			for (int k=0;i-k-2>=0;++k)
				s=(s+g[i-k-2][j-1]*p2[(k+1)*k>>1])%mo;
			g[i][j]=s;
		}
	for (int i=1;i<=n;++i){
		ll s=0;
		for (int j=1;j<i;++j)
			s=(s+f[j]*g[i][j])%mo;
		f[i]=(p2[(i+1)*i>>1]-s%mo+mo)%mo;
	}
	printf("%lld\n",f[n]);
	return 0;
}

---

总结

好思维的一道题目……
如果在比赛的时候给我想，我死也想不出来……
论脑子的重要性……