BZOJ 1996 [HNOI2010] chorus 合唱队 区间DP

题目描述 Description
为了在即将到来的晚会上有更好的演出效果，作为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心中有一个理想队形，他想知道多少种初始队形可以获得理想的队形

输入描述 Input Description
第一行一个正整数N,表示总人数
第二行使用空格隔开的N个正整数,从左到右表示小A心中的理想队形:H1,H2,H3…Hn

输出描述 Output Description
仅包含一个数,表示从多少种初始队形出发能通过上述排队的方式获得小A心中的理想队形.因为满足条件的初始队形数可能很大,所以规定只要输出满足条件的初始队形数mod19650827的值

样例输入 Sample Input
4
1701 1702 1703 1704

样例输出 Sample Output
8

数据范围及提示 Data Size & Hint
对于30%的数据,$1<=N<=100$
对于100%的数据,$1<=N<=1000,1000<=Hi<=2000$

Solution

很容易就能看出这是一个区间DP吧(雾
通过手推发现每个序列$[l,r]$只有可能从$[l,r-1]$和$[l+1,r]$转移
同时如何转移还要看之前序列的朝向
所以用
$f[l][r][0]$表示$[l,r]$区间内最后朝向左的序列个数
$f[l][r][1]$表示$[l,r]$区间内最后朝向右的序列个数

显然的通过题目条件转移
$if(a[l]<a[r]) f[l][r][0]+=f[l+1][r][1];$
$if(a[l]<a[l+1]) f[l][r][0]+=f[l+1][r][0];$
$if(a[r]>a[r-1]) f[l][r][1]+=f[l][r-1][1];$
$if(a[r]>a[l]) f[l][r][1]+=f[l][r-1][0];$

代码奇短,好懂,自己YY吧QwQ

代码如下:

#include <bits/stdc++.h>
using namespace std;
#define r l+lenth-1
const int N = 1005;
const int mod = 19650827;
int n;
int a[N],b[N];
int f[N][N][2];
int read() {
    int ans=0,flag=1;
    char ch=getchar();
    while((ch<'0' || ch>'9') && ch!='-') ch=getchar();
    if(ch=='-') { flag=-1;ch=getchar(); }
    while(ch>='0' && ch<='9') { ans=ans*10+ch-'0';ch=getchar(); }
    return ans*flag;
}
int main() {
    n=read();
    for(int i=1;i<=n;++i) a[i]=read(),b[i]=a[i];
    sort(b+1,b+n+1);unique(b+1,b+n+1);
    for(int i=1;i<=n;++i) a[i]=lower_bound(b+1,b+n+1,a[i])-b;
    for(int lenth=1;lenth<=n;++lenth) {
        for(int l=1;l+lenth-1<=n;++l) {
            if(lenth==1) {f[l][r][1]=1;continue;}
            if(a[l]<a[r]) f[l][r][0]+=f[l+1][r][1];
            if(a[l]<a[l+1]) f[l][r][0]+=f[l+1][r][0];
            if(a[r]>a[r-1]) f[l][r][1]+=f[l][r-1][1];
            if(a[r]>a[l]) f[l][r][1]+=f[l][r-1][0];
            f[l][r][0]%=mod;f[l][r][1]%=mod;
        }
    }
    printf("%d\n",(f[1][n][0]+f[1][n][1])%mod);
    return 0;
}