区间dp--洛谷P3205 [HNOI2010] 合唱队

[HNOI2010] 合唱队

洛谷链接

题目描述

为了在即将到来的晚会上有更好的演出效果，作为 AAA 合唱队负责人的小 A 需要将合唱队的人根据他们的身高排出一个队形。假定合唱队一共 $n$ 个人，第 $i$ 个人的身高为 $h_i$ 米（$1000\le h_i\le 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 心中有一个理想队形，他想知道多少种初始队形可以获得理想的队形。

请求出答案对 $19650827$ 取模的值。

输入格式

第一行一个整数 $n$。
第二行 $n$ 个整数，表示小 A 心中的理想队形。

输出格式

输出一行一个整数，表示答案 $\mathrm{mod}19650827$ 的值。

样例 #1

样例输入 #1

4
1701 1702 1703 1704

样例输出 #1

8

提示

对于 $30\%$ 的数据，$n\le 100$。
对于 $100\%$ 的数据，$n\le 1000$，$1000\le h_i\le 2000$。

题目解答

明显题目应当从局部入手，求出答案队列局部不同方案数累加，且答案区间可以进行划分，明显符合大区间包含小区间的情况，故考虑区间$dp$
考虑当前队列形成情况
要么是从左边进来一个人形成此队列
要么是从右边进来一个人形成队列
故可以划分为以下情况
$f[i,j,0]$表示区间$i$到$j$且是由$i$从左边进来形成的方案数
$i$从左边进来包含两种情况：

1. 上一个进来的数为$a[i+1]$
   可以得出状态转移方程为
   $f[i,j,0]+=f[i+1,j,0]$
   要求$a[i]<a[i+1]$
2. 上一个进来的数为$a[j]$
   可得状态转移方程为
   $f[i,j,0]+=f[i+1,j,1]$
   要求$a[i]<a[j]$

最终可以得到$i$从左边进入得全部答案

---

同理可以得出$j$从右边进入的答案
状态转移方程为
$f[i,j,1]+=f[i,j-1,1]$
要求$a[j]>a[j-1]$
$f[i,j,1]+=f[i,j-1,0]$
要求$a[j]>a[i]$

---

初始化
我们初始化仅初始化一个方向即可，我们可以默认方向为左边，因为我们只有一个人时方案数应该是1，如果我们初始化两边，那么他可以从右边进来，也可以从左边进来，那么会导致方案数被错误记为2

code

---

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=1010;
const int Mod=19650827;
typedef long long ll;
ll f[N][N][2];
int n;
int a[N];
int main()
{   cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    for(int i=1;i<=n;i++)f[i][i][0]=1;//初始化，仅初始化一个方向即可
    for(int len=1;len<=n;len++)//枚举区间长度
        for(int i=1;i+len-1<=n;i++)//枚举左端点
    {
        int j=i+len-1;
        //第一种情况，考虑最后i从左边进入
        if(a[i]<a[i+1])f[i][j][0]+=f[i+1][j][0];
        if(a[i]<a[j])f[i][j][0]+=f[i+1][j][1];
        //第二种情况，考虑j从右边进入
         if(a[j]>a[j-1])f[i][j][1]+=f[i][j-1][1];
         if(a[j]>a[i])f[i][j][1]+=f[i][j-1][0];
         //别忘记取模
         f[i][j][0]%=Mod;
         f[i][j][1]%=Mod;
    }
    cout<<(f[1][n][0]+f[1][n][1])%Mod<<endl;
}