2011 Dalian Onsite——E Number String

看shi哥的解题报告，应该很清晰能理解方程。

但是其实就是处理'D'的时候有些像我一样的弱菜要理解很久。

所以就写一个报告给没到大牛水平的ACMers共勉一下。

题意：求满足一定升降序列的排列数的方案数

这题最需要理解的一个性质是当加入一个新的数n时，如果想把n放在前n-1个位置的任何一个，

并且不影响升降性质，只需要在之前方案中把当前在n位置的数x(n放前面必然有一个x要放在n的位置）

所有的大于等于x的数都加1即可。比如之前方案是（1，5，3，2，4），现在把3放在第6位。

则对应方案是(1,6,4,2,5,3)。（1,2不变，其他加一）

理解这个性质后，dp转移方程就很简单了。

设dp[i][j]表示长度为i的序列，以j结尾的满足前i个条件的排列数

则当第i个要求为'?'时，很简单，即前面方案的总和。即dp[i][j]=sum(dp[i-1][k])(1<=k<i)

  当第i个要求为'I'时，也容易理解。之前小于j的结尾数的方案都能累加过来。

  即dp[i][j]=sum(dp[i-1][k])(1<=k<j)

但当第i个要求为'D'时，需要特殊理解下，对于某个j，为了保持排列的性质不变，

我们把大于等于j的数字都+1，这样就保证不改变状态了。其实‘I’的时候也是这样处理的。

只是处理后还是(1<=k<i)。但是‘D’的时候，k的取值应该是[j,i)。因为等于j的时候，由于之前

大于等于j的数字都加1了。所以之前的j变成了j+1。一样满足'D'条件。

然后这个方程用部分和可以优化成O(n^2)，得解。

#include <iostream>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <numeric>
#include <algorithm>
#include <functional>
using namespace std;
#define M 1005
#define MOD 1000000007
long long dp[M][M];
char st[M];
int main()
{
    while (scanf("%s",st+1)!=EOF)
    {
        int i,j,k;
        int n=strlen(st+1)+1;
        fill(dp[0],dp[n+1],0LL);
        dp[1][1]=1;
        for (i=2;i<=n;i++)
        {
            if (st[i-1]=='?')
            {
                long long sum=0;
                for (j=1;j<i;j++)
                    sum=(sum + dp[i-1][j]) % MOD;
                for (j=1;j<=i;j++)
                    dp[i][j]=sum;
            }
            else if (st[i-1]=='I')
            {
                dp[i][1]=0;
                for (j=2;j<=i;j++)
                    dp[i][j]=(dp[i][j-1]+dp[i-1][j-1]) % MOD;
            }
            else
            {
                dp[i][i]=0;
                for (j=i-1;j>=1;j--)
                    dp[i][j]=(dp[i][j+1]+dp[i-1][j]) % MOD;
            }
        }
        long long ans=0;
        for (i=1;i<=n;i++)
            ans=(ans + dp[n][i]) % MOD;
        printf("%lld\n",ans);
    }
    return 0;
}