正题
题目链接:https://www.luogu.com.cn/problem/AT2366
题目大意
有一个长度为2n−12n-12n−1的序列aaa,你可以将其重新排列,定义bib_ibi为a1∼2i−1a_{1\sim 2i-1}a1∼2i−1的中位数。
询问有多少种不同的可能的bbb序列。
1≤n≤501\leq n\leq 501≤n≤50
解题思路
先假设aaa中所有数字都不同,我们去考虑一下bbb的限制。
相当于我们每次在aia_iai中加入两个数字x,yx,yx,y,然后如果x,yx,yx,y都小于最后一个bbb则新的bbb取它在aaa中的前驱,如果都大于则取后继,如果一大一小则不动。
基于这个限制我们得到的条件是每次bib_ibi加入数字之后的排名只能变动一位,这是充要的。
考虑去形式化这个条件,将aaa排序后,我们有:
- an−i+1≤bi≤an+i−1a_{n-i+1}\leq b_i\leq a_{n+i-1}an−i+1≤bi≤an+i−1
- 不存在j<ij<ij<i满足bi<bj<bi+1b_i<b_j<b_{i+1}bi<bj<bi+1或bi+1<bj<bib_{i+1}<b_j<b_ibi+1<bj<bi
这样我们其实并不需要考虑aaa填的情况,如果bbb一次向外跨过了很多个aaa,那么我们可以视这些aaa此时还没有加入,我们只需要限制往内跨的情况出现就好了。
不过这第二个条件依旧不好处理,我们考虑倒着做,那当我们确定一个bib_ibi和bi+1b_{i+1}bi+1之后,相当于(bi,bi+1)(b_i,b_{i+1})(bi,bi+1)之间的数字就都不能选择了。
考虑dpdpdp,设fi,l,rf_{i,l,r}fi,l,r表示目前填了iii个,对于现在的bib_ibi来说,左边还剩下lll个位置,右边还剩下rrr个位置。
然后每次加入两边的aia_iai,如果al=al+1a_{l}=a_{l+1}al=al+1,那么我们将其视为同一个,右边同理。
时间复杂度:O(n4)O(n^4)O(n4)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=110,P=1e9+7;
ll n,m,a[N],f[N][N][N];
signed main()
{scanf("%lld",&n);m=2*n-1;for(ll i=1;i<=m;i++)scanf("%lld",&a[i]);sort(a+1,a+m+1);f[n][0][0]=1;for(ll i=n-1;i>=1;i--){ll l=i,r=2*n-i;ll dl=(a[l]!=a[l+1]),dr=(a[r]!=a[r-1]);for(ll l=0;l<=m;l++)for(ll r=0;r<=m;r++){if(!f[i+1][l][r])continue;(f[i][l+dl][r+dr]+=f[i+1][l][r])%=P;for(ll k=0;k<l+dl;k++)(f[i][k][r+dr+1]+=f[i+1][l][r])%=P;for(ll k=0;k<r+dr;k++)(f[i][l+dl+1][k]+=f[i+1][l][r])%=P;}}ll ans=0;for(ll i=0;i<=m;i++)for(ll j=0;j<=m;j++)(ans+=f[1][i][j])%=P;printf("%lld\n",ans);return 0;
}