我们忽略掉
A
2
i
−
1
A_{2i-1}
A2i−1和
A
2
i
A_{2i}
A2i均不为
−
1
-1
−1的
i
i
i,称最终序列中相邻两个数为一对匹配。
考虑从大到小分配权值,每对匹配只有在较小的权值出现的时候才确定相对顺序。设
F
[
i
]
[
j
]
[
k
]
F[i][j][k]
F[i][j][k]表示考虑了不小于
i
i
i的权值,剩余没匹配的权值中恰有
j
j
j个出现在
A
A
A序列中,恰有
k
k
k个没出现在
A
A
A序列中。
转移的时候讨论一下
i
−
1
i-1
i−1是否出现在
A
A
A序列中,如果出现的话,我们可以让它作为最小值(匹配一个
k
k
k)或者不匹配,如果没出现的话也类似,但需要留意一下,如果它的选择是跟后面的某个
j
j
j匹配,对应的都是不相同的方案,如果是跟
k
k
k匹配的话,具体匹配哪个是不会影响答案的,只需要确定最终位置即可。
时间复杂度
O
(
n
3
)
\mathcal O(n^3)
O(n3)。
#include <bits/stdc++.h>
#define MOD 1000000007
using namespace std;
typedef long long ll;
inline void update(int &x,ll y) {
x=(x+y)%MOD;
}
int f[2][305][305];
int num[605],kind[605];
int main() {
int n;
scanf("%d",&n);
for(int i=1;i<=2*n;i++) scanf("%d",&num[i]);
int cnt=0;
for(int i=1;i<=n;i++)
if (num[2*i-1]==-1||num[2*i]==-1) {
if (num[2*i-1]!=-1) kind[num[2*i-1]]=1;
if (num[2*i]!=-1) kind[num[2*i]]=1;
if (num[2*i-1]==-1&&num[2*i]==-1) cnt++;
}
else kind[num[2*i-1]]=kind[num[2*i]]=-1;
int cur=0;
f[cur][0][0]=1;
int s1=0,s2=0;
for(int i=2*n;i>0;i--)
if (kind[i]!=-1) {
cur^=1;
memset(f[cur],0,sizeof(f[cur]));
for(int j=0;j<=min(s1,n);j++)
for(int k=0;k<=min(s2,n);k++)
if (f[cur^1][j][k]) {
if (kind[i]) {
if (j<n) update(f[cur][j+1][k],f[cur^1][j][k]);
if (k) update(f[cur][j][k-1],f[cur^1][j][k]);
}
else {
if (k<n) update(f[cur][j][k+1],f[cur^1][j][k]);
if (j) update(f[cur][j-1][k],(ll)j*f[cur^1][j][k]);
if (k) update(f[cur][j][k-1],(ll)(cnt-((s2-s1+j-k)>>1))*f[cur^1][j][k]);
}
}
if (kind[i]) s1++; else s2++;
}
printf("%d\n",f[cur][0][0]);
return 0;
}