[洛谷]P3205 [HNOI2010]合唱队 (#区间dp)

题目描述

为了在即将到来的晚会上有更好的演出效果，作为AAA合唱队负责人的小A需要将合唱队的人根据他们的身高排出一个队形。假定合唱队一共N个人，第i个人的身高为Hi米(1000<=Hi<=2000),并已知任何两个人的身高都不同。假定最终排出的队形是A 个人站成一排，为了简化问题，小A想出了如下排队的方式：他让所有的人先按任意顺序站成一个初始队形，然后从左到右按以下原则依次将每个人插入最终棑排出的队形中：

-第一个人直接插入空的当前队形中。

-对从第二个人开始的每个人，如果他比前面那个人高(H较大)，那么将他插入当前队形的最右边。如果他比前面那个人矮(H较小)，那么将他插入当前队形的最左边。

当N个人全部插入当前队形后便获得最终排出的队形。

例如，有6个人站成一个初始队形，身高依次为1850、1900、1700、1650、1800和1750,

那么小A会按以下步骤获得最终排出的队形：

1850

• 1850 , 1900 因为 1900 > 1850

• 1700, 1850, 1900 因为 1700 < 1900

• 1650 . 1700, 1850, 1900 因为 1650 < 1700

• 1650 , 1700, 1850, 1900, 1800 因为 1800 > 1650

• 1750， 1650, 1700，1850, 1900, 1800 因为 1750 < 1800

因此，最终排出的队形是 1750，1650，1700，1850, 1900，1800

小A心中有一个理想队形，他想知道多少种初始队形可以获得理想的队形

输入输出格式

输入格式：

输出格式：

注意要mod19650827

输入输出样例

输入样例#1

4
1701 1702 1703 1704

输出样例#1

8

说明

30%的数据：n<=100

100%的数据：n<=1000

---

思路

区间dp好题！！！！！！！

这题是和P1220有点相似的。给定一个理想数列，例如：1701,1702,1703,1704，求一共有多少种原队列。

通过排队方式很显然可以看出原队列的最后一位应该从1701和1704中选取。如果选择的是1701，那么原问题就产生一个子问题： 给定一个理想数列，1702,1703,1704求一共有多少种原队列；如果选择的是1704，那么原问题就产生另一个子问题： 给定一个理想序列，1701,1702,1703求一共有多少种原序列。对于一个最终序列区间[i,j]，最后一个添加的只可能是第i个人或者第j个人。

所以如果最后一个是i，则对于这个区间[i,j]，a[i]<a[i+1]或a[i]<a[j]

如果最后一个是j，则对于这个区间[i,j]，a[j]>a[j-1]或a[j]>a[i]

令dp[i][j][0]为在最终序列的区间[i,j]中，最后一个添加i的方案数；dp[i][j][1]为在最终序列的区间[i,j]中，最后一个添加j的方案数。[i,j]这个区间可以从[i-1,j]或[i,j+1]转移过来。

1.当r（右端点）的数大于上一个且上一个插在右端时

dp[i][j][1]=dp[i][j][1]+dp[i][j-1][1] (a[j]>a[j-1])

2.当右端点的数大于左端点且上一个插入在左端时

dp[i][j][1]=dp[i][j][1]+dp[i][j-1][0] (a[j]>a[i])

3.当左端点的数小于后一个且上一个插入在左端时

dp[i][j][0]=dp[i][j][0]+dp[i+1][j][0] (a[i]<a[i+1])

4.当左端点的数小于右端点且上一个插入在右端时

dp[i][j][0]=dp[i][j][0]+dp[i+1][j][1] (a[i]<a[j])

答案为最后一个插入的为右端点和最后一个插入的为左端点之和。

#include <stdio.h>
#include <iostream>
using namespace std;
int dp[1001][1001][2],n,s,a[1001]; 
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	register int i,j,l;
	cin>>n;
	for(i=1;i<=n;i++)
	{
		cin>>a[i];
	}
	for(i=1;i<=n;i++)
	{
		dp[i][i][0]=1;//临界条件：i==j 
	}
	for(l=2;l<=n;l++)
	{
		for(i=1;i+l-1<=n;i++)
		{
			int j=i+l-1;
			if(a[j]>a[j-1])
			{
				dp[i][j][1]+=dp[i][j-1][1];
			}
			if(a[j]>a[i])
			{
				dp[i][j][1]+=dp[i][j-1][0];
			}
			if(a[i]<a[i+1])
			{
				dp[i][j][0]+=dp[i+1][j][0];
			}
			if(a[i]<a[j])
			{
				dp[i][j][0]+=dp[i+1][j][1];
			}
			dp[i][j][0]%=19650827;
			dp[i][j][1]%=19650827;
		}
	}
	cout<<(dp[1][n][0]+dp[1][n][1])%19650827;
	return 0;
}