HDU 5396

一道区间 dp 题目

有 n 个数字，n-1 个操作符，求所有按照不同顺序进行运算的结构的和。。。

举个例子吧。。。

1 2 3 4 

+ * +

结果就是 ( (1+2) * 3 + 4  ) + ( 1+ 2*3 + 4) + ((1+2) * (3+4))

总之就是这样

我们就可以用区间 dp 的思路 dp[ i ] [ j ] 代表 区间 ( i , j ) 之和

那么 dp[ i ] [ j ] =  ( dp [ i ] [ k ] * (j - k - 1) ! +/- dp[ k+1 ] [ j ] * ( k - i) !  ) 是 +/- 的情况

如果是 * 号的话 就不需要乘以阶乘了，因为加法的话 我们理论上 对于左边的i到k来说，假设右边有n个符号，那么右边肯定有n！种方案，然后组合一下 需要乘以阶乘代表全排列，最后还需要乘上 一个 组合数 C[j - i - 1][k - i] 因为 第 k 个运算符是最后取得么假设左边有n个操作符，右边有m个操作符，这么操作符取的顺序不同，最后的方案数也是不同的。

以下是 AC代码

 

 

 

 

 

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define ll long long
const int mod = 1000000007;
char s[106];
ll num[106],dp[106][106];
ll a[106],c[106][106];
void init()
{
    int i,j;
    a[0]=1;
    a[1]=1;
    a[2]=2;
    for(ll i=3;i<=100;i++)
    {
        a[i]=(a[i-1]*i)%mod;
    }
    for(ll i=1;i<=100;i++)c[i][0]=1;

    for(ll i=1;i<=100;i++)
    {
        for(j=1;j<=i;j++)
        {
            if(i==j)c[i][j]=1;
            else if(i>j)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
        }
    }
}

int main()
{
    int n,m;
    ll t;
    init();
    while(scanf("%d",&n)!=EOF)
    {
        for(ll i=1;i<=n;i++)
        {
            scanf("%lld",&num[i]);
            dp[i][i]=num[i];
        }
        scanf("%s",s+1);
        for(ll i=1;i<=n-1;i++)
        {
            if(s[i]=='+')
            {
                dp[i][i+1]=(num[i]+num[i+1])%mod;
            }
            else if(s[i]=='-')
            {
                dp[i][i+1]=(num[i]-num[i+1])%mod;
            }
            else if(s[i]=='*')
            {
                dp[i][i+1]=(num[i]*num[i+1])%mod;
            }
        }
        for(ll len=3;len<=n;len++)
        {
            for(ll i=1;i+len-1<=n;i++)
            {
                ll j=i+len-1;
                dp[i][j]=0;
                for(ll k=i;k<j;k++)
                {
                    t=0;
                    if(s[k]=='*')
                    {
                        t=(dp[i][k]*dp[k+1][j])%mod;
                    }
                    else if(s[k]=='+')
                    {
                        t=(dp[i][k]*a[j-k-1]+dp[k+1][j]*a[k-i])%mod;

                    }
                    else if(s[k]=='-')
                    {
                        t=(dp[i][k]*a[j-k-1]-dp[k+1][j]*a[k-i])%mod;
                    }
                    dp[i][j]=(dp[i][j]+t*c[j-i-1][k-i])%mod;
                }

            }
        }
        printf("%lld\n",(dp[1][n]+mod)%mod);
    }
    return 0;
}