【AGC030F】Permutation and Minimum DP

题目大意

  有一个长度为序列 (a),其中某些位置的值是 (-1)

  你要把 (a) 补成一个排列。

  定义 (b_i=min(a_{2i-1},a_{2i})),求有多少种可能的 (b)

  (nleq 300)

题解

  如果 (a_{2i-1})(a_{2i}) 都有值,就把这两个位置扔掉。

  记 (c_i) 表示 (i) 这个值是否在初始的 (a) 中。

  从后往前DP。记 (f_{i,j,k}) 表示已经处理完了 (i) 后面的数,有多少个 (j>i,c_j=1) 的数匹配的是 (leq i) 的数,有多少个 (j>i,c_j=0) 的数匹配的是 (leq i) 的数。

  如果 (c_i=1) 且往后匹配的是 (c_j=0),那么方案数为 (1)。(因为 (min=i)

  如果 (c_i=0) 且往后匹配的是 (c_j=0),那么先暂定方案数为 (1)。(因为暂时不能确定 (i) 填在哪个位置。)记这种匹配对数为 (cnt)

  如果 (c_i=0) 且往后匹配的是 (c_j=1),那么方案数为 (j)。(因为可以确定 (i) 填在哪个位置。)

  最后方案数要乘上 (cnt!),因为这些位置的 (b) 可以随便交换。

  时间复杂度:(O(n^3))

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<functional>
#include<cmath>
#include<vector>
#include<assert.h>
//using namespace std;
using std::min;
using std::max;
using std::swap;
using std::sort;
using std::reverse;
using std::random_shuffle;
using std::lower_bound;
using std::upper_bound;
using std::unique;
using std::vector;
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
typedef std::pair<int,int> pii;
typedef std::pair<ll,ll> pll;
void open(const char *s){
#ifndef ONLINE_JUDGE
	char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
void open2(const char *s){
#ifdef DEBUG
	char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
int rd(){int s=0,c,b=0;while(((c=getchar())<'0'||c>'9')&&c!='-');if(c=='-'){c=getchar();b=1;}do{s=s*10+c-'0';}while((c=getchar())>='0'&&c<='9');return b?-s:s;}
void put(int x){if(!x){putchar('0');return;}static int c[20];int t=0;while(x){c[++t]=x%10;x/=10;}while(t)putchar(c[t--]+'0');}
int upmin(int &a,int b){if(b<a){a=b;return 1;}return 0;}
int upmax(int &a,int b){if(b>a){a=b;return 1;}return 0;}
const ll p=1000000007;
const int N=310;
void add(ll &a,ll b)
{
	a=(a+b)%p;
}
int n;
int a[2*N];
ll f[2*N][N][N];
int b[2*N];
int c[2*N];
int t;
int main()
{
	open2("f");
	scanf("%d",&n);
	for(int i=1;i<=2*n;i++)
		scanf("%d",&a[i]);
	int cnt=0;
	for(int i=1;i<=2*n;i+=2)
		if(a[i]!=-1&&a[i+1]!=-1)
			b[a[i]]=b[a[i+1]]=2;
		else if(a[i]!=-1)
			b[a[i]]=1;
		else if(a[i+1]!=-1)
			b[a[i+1]]=1;
		else
			cnt++;
	for(int i=1;i<=2*n;i++)
		if(b[i]==1)
			c[++t]=1;
		else if(b[i]==0)
			c[++t]=2;
	f[t][0][0]=1;
	for(int i=t;i>=1;i--)
		for(int j=0;j<=t&&j<=n;j++)
			for(int k=0;k<=t&&k<=n;k++)
				if(c[i]==1)
				{
					add(f[i-1][j+1][k],f[i][j][k]);
					if(k)
					 	add(f[i-1][j][k-1],f[i][j][k]);
				}
				else
				{
					add(f[i-1][j][k+1],f[i][j][k]);
					if(k)
						add(f[i-1][j][k-1],f[i][j][k]);
					if(j)
						add(f[i-1][j-1][k],f[i][j][k]*j);
				}
	ll ans=f[0][0][0];
	for(int i=1;i<=cnt;i++)
		ans=ans*i%p;
	printf("%lld
",ans);
	return 0;
}
原文地址:https://www.cnblogs.com/ywwyww/p/10198176.html