jzoj5235. 【NOIP2017模拟8.7A组】好的排列

Description

对于一个1->n的排列 ，定义A中的一个位置i是好的，当且仅当Ai-1>Ai 或者Ai+1>Ai。对于一个排列A，假如有不少于k个位置是好的，那么称A是一个好的排列。
现在有q个询问，每个询问给定n,k，问有多少排列是好的。答案对10^9+7取模。

Input

输入文件名为permutation.in。
首先输入q。
接下来输入q个询问n,k 。

Output

输出文件名为permutation.out。
输出q行，每行一个整数代表答案。

Sample Input

8
4 3
6 4
10 7
20 14
50 40
100 72
1000 900
3000 2000

Sample Output

8
448
1433856
868137807
908422882
609421284
150877522
216180189

Data Constraint

对于20%的数据，n<=10,q=1
对于40%的数据，n<=20,q=1
对于60%的数据，n<=100
对于100%的数据，n,k<=3000,q<=10000

题外话

大爷省选加油啊！
重要的事情怎么能不重复能？

题解

这是道套路题。

一开始还以为是什么数位dp，然后仔细读读题就大概知道怎么做了。
以前做过一道题是依次插入一个数，然后就分成很多很多的小块。
大概就是设$f_{i,j,k}$表示当前放第i个数，然后当前分成了j个块，有k个好数。

结果这样分析下来就变成$O(n^3)$ 了。
然后再看看发现似乎贡献其实只有两种：

• 一、就是当前第i个数插在两个原本是好数的数中间，那么好数个数不变。
• 二、当前第i个数插在了一个非好数的旁边，那么好数个数+1（原来非好数变成好数）

然后就可以不用上面中间那维了，直接dp即可。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
const long long mo=1000000007;
const int maxn=3010;

int n,st,en;
long long f[2][maxn][maxn],ans[maxn][maxn],mi[maxn];

int main()
{
	freopen("permutation.in","r",stdin);
	freopen("permutation.out","w",stdout);
//	freopen("data.in","r",stdin);
//	freopen("data.out","w",stdout);
	n=3000;
	f[0][0][1]=1;
	st=0;en=1;
	mi[0]=1;
	for (int i=1;i<=n;i++)
	{
		mi[i]=2*i;
	}
	ans[1][0]=1;
	for (int i=2;i<=n;i++)
	{
		for (int k=1;k<=i;k++)
		{
			ans[i][k]=(ans[i][k]+ans[i-1][k-1]*(i-k)*2%mo)%mo;
			ans[i][k]=(ans[i][k]+ans[i-1][k]*(i-(i-1-k)*2)%mo)%mo;
		}
	}
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		int x;int y;
		scanf("%d%d",&x,&y);
		long long answer=0;
		for (int i=y;i<=x-1;i++) answer=(answer+ans[x][i])%mo;
		printf("%lld\n",answer);
	}
}