bzoj 1996 [Hnoi2010]chorus 合唱队

Description

为了在即将到来的晚会上有更好的演出效果，作为AAA合唱队负责人的小A需要将合唱队的人根据他们的身高排出一个队形。假定合唱队一共有N个人，第i个人的身高为$H_i$毫米(1000<=$H_i$<=2000),并且已知任何两个人的身高都不同。假定最终排出的队形是N个人站成一排，为了简化问题，小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

从文件input.txt中读入数据，输入文件第一行是一个正整数N，表示总人数。输入文件第二行是用空格隔开的N个正整数，从左到右表示小A心中的理想队形：$H_1，H_2，···，H_N$。输入的数据保证1000<=$H_i$<=2000且没有相同的H值，其中30%点数据满足1<=N<=100，100%的数据满足1<=N<=1000。

Output

输出文件output.txt仅包含一个数，表示从多少种初始队形出发能通过上述排队的方式获得输入文件中指定的小A心中的理想队形。因此为满足条件的初始队形数可能很大，所以规定只要输出满足条件的初始队形数mod 19650827的值。

Sample Input

4
1701 1702 1703 1704

Sample Output

8

HINT

样例解释：8种初始队形分别为(1701,1702,1703,1704)，(1704,1703,1702,1701)，(1702,1701,1703,1704)，(1703,1702,1701,1704)，(1702,1703,1701,1704)，(1702,1703,1704,1701)，(1703,1704,1702,1701)，(1703,1702,1704,1701)

Key To Problem

给定一个理想数列，例如：1701,1702,1703,1704，求一共有多少种原队列
通过小A的排队方式很显然可以看出原队列的最后一位应该从1701和1704中选取
如果选择的是1701，那么原问题就产生一个子问题：
    给定一个理想数列，1702,1703,1704求一共有多少种原队列
如果选择的是1704，那么原问题就产生另一个子问题：
    给定一个理想序列，1701,1702,1703求一共有多少种原序列
由此，任何一个大于二的理想序列都可以分成两个子理想序列
一次选择可以选择队首元素，或者队尾元素
于是我想到了vijos-P1037搭建双塔
本题也同样需要用两个二维数组，来回推导
数组f是当本次选择队首元素时一共有多少种原序列
数组g是当本次选择队尾元素时一共有多少种原序列
状态转移方程为：
            if(a[i]<=a[i+1])
                f[i][j]=(f[i][j]+f[i+1][j])%MOD;
            if(a[i]<a[j])
                f[i][j]=(f[i][j]+g[i][j-1])%MOD;
            if(a[i]<a[j])
                g[i][j]=(g[i][j]+f[i-1][j])%MOD;
            if(a[j]>a[j-1])
                g[i][j]=(g[i][j]+g[i][j-1])%MOD;

CODE

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1010
#define MOD 19650827
using namespace std;
int n;
int a[N];
int f[N][N];
int g[N][N];

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
        f[i][i]=1;
    for(int i=n-1;i>=1;i--)
    {
        for(int j=i+1;j<=n;j++)
        {
            if(a[i]<a[i+1])
                f[i][j]=(f[i][j]+f[i+1][j])%MOD;
            if(a[i]<a[j])
                f[i][j]=(f[i][j]+g[i+1][j])%MOD;
            if(a[i]<a[j])
                g[i][j]=(g[i][j]+f[i][j-1])%MOD;
            if(a[j]>a[j-1])
                g[i][j]=(g[i][j]+g[i][j-1])%MOD;
        }
    }
    printf("%d\n",(f[1][n]+g[1][n])%MOD);
    return 0;
}