Sit sit sit HDU - 5151 （区间dp +组合计数）

这道这道区间DP，我也开始觉得其实区间DP是一种应用型的思想，做这类题目一个重要的点是在于题目情景的把握,这道题的一个情景就是数学的排列组合问题. 首先应用根据小区间推出大区间的思路,我们可以先固定一个位置k,k位置是最后做的位置,那么我们要算出在这种情况下符合的方法数,假如k是头或者尾,那不用说, 此时的方法为dp[i+1][j]或者是dp[i][j-1]; 假如k是有邻居的，那我们就要考虑邻居的颜色，如果颜色不同，那么方法便是0，否则的话，可以根据排列组合的方法，方法数自然等于dp[i][k-1]*dp[k+1][j]*X; 这个X便是考你对于情景的理解，i-k-1以及k+1-j这两个区间之内已经是各自相对有序做好的位置，那在形成一组的时候我们也要保证相对有序 这里我们用组合的思想进行分析： 首先总共有j-i个位置（因为不考虑k位置），那么我们有两个各自相对有序的数组要合二为一，问最终有多少种不同的序列。 这边我一开始的想法便是用插空法，但一直过不了，小的数据过的去，大的数据过不去，我猜应该是在大数据的时候数据溢出了，因为插空的思想准没错。 所以这边退而求其次，还记得组合当中有所谓的什么元素法还有位置法，你可以一开始固定元素去插，也可以固定位置去排。 这边我们选择固定位置，一开始为k-i个元素选择好位置，那另外一组数的位置自然也固定了。 这边也牵扯到一个十分重要的思想：便是我们在考虑两个变量的时候，可能会由于各自运算的时间可能刚好符合给的限定时间，假如合起来的话会超过。 这时候我们不妨寻找两个变量之间的关系，从而只做一个数，从关系当中，推出另外一个数，这是一种十分有效，并且能够突破思维的方法。（蓝桥杯的教训） 给出状态转移方程：dp[i][j]=(dp[i][j]+dp[i][k-1]*dp[k+1][j]*c[j-i][k-i]) --------------------- 本文来自 jason_star 的CSDN 博客 ，全文地址请点击：https://blog.csdn.net/u013611908/article/details/43193651?utm_source=

前面的基本他已经分析的非常清楚了，我补充一下对那个组合数的理解：

区间[i,j]中，我们已经知道了区间[i,k-1],[k+1,j],他们都是已经计算好的，有序的。现在我们考虑进行区间合并，那么现在的问题就是有多少种不同的合并方法，总人数有j-i+1,但是要去掉第k个，所以还有j-i个，[i,k-1]区间有k-i个人，我们表面上需要进行一次排列才能够计算出方法数，但是问题是[i,k-1]已经是有序的了，我们不要在考虑排列问题。

啰嗦了半天，它就是[i,k-1]内部类似一个自动排列机一样，只要你把人分给它，那它就拥有了不同的排列方法，所以我们只要考虑将j-i个人，分k-i个出来，就可以解决问题了（对于后面的j-k个，只要前面的确定了，那么它也就确定了）。

总结：这类dp+组合数学的问题，不能只考虑dp或组合数学的内容，就像上面的分析，合并区间的方法数，看似应该是排列数，但是加上我们的dp定义以后，我们只要求出组合数就行了。

代码：

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+7;
const int inf=0x3f3f3f3f;
typedef long long ll;
ll C[120][102],A[120];
const int mod=1e9+7;
int n;
ll dp[120][120];
int a[120];
void init()
{
    A[0] = 1;
    for(int i=1; i<=100; i++)
        A[i] = (A[i-1] * i)%mod;
    C[0][0] = 1;
    for(int i=1; i<=100; i++)
    {
        C[i][0] = 1;
        C[i][i]=1;
        for(int j=1; j<=i; j++)
            C[i][j] = (C[i-1][j-1] + C[i-1][j])%mod;
    }
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    init();
    while(scanf("%d",&n)!=EOF)
    {
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            dp[i][i]=1;
            if(i<n)
            {
                dp[i][i+1]=2;
            }
        }   
        for(int l=2;l<n;l++)
        {
            for(int i=1;i<=n-l;i++)
            {
                int j=l+i;
                for(int k=i;k<=j;k++)
                {
                    if(i==k)//两边直接放
                    {
                        dp[i][j]+=dp[i+1][j];
                    }
                    else if(j==k)
                    {
                        dp[i][j]=(dp[i][j]+dp[i][j-1])%mod;
                    }
                    else
                    {
                        if(a[k-1]==a[k+1])
                        {
                            dp[i][j]=(dp[i][j]+dp[i][k-1]*dp[k+1][j]%mod*C[j-i][k-i]%mod)%mod;
                        }
                    }
                }
            }
        }
        printf("%lld\n",dp[1][n]);
    }
    return 0;
}