hdu 5155

动态规划里比较常见的类型,按行进行动态规划,dp[i][j],i表示当前是的行,j是当前已经有放宝石的列,

然后之后每一行至少放一个珠宝,然后利用排列组合的乘法原则得到如下的转移方程

dp[i][j] += dp[i-1][k] *c[m-k][j-k]*c[k][t-j+k]

基本思想就是每次从选过的列和没选过的列根据状态分别选取,进行转移;

代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define MAX 57
#define MOD 1000000007

using namespace std;

typedef long long LL;

LL c[MAX][MAX];

void init ( )
{
    c[0][0]=1;
    for ( int i = 1 ; i <= 50 ; i++ )
        for ( int j = 0 ; j <= i ; j++ )
            if ( j==0 || j==i ) c[i][j] = 1;
            else c[i][j] = (c[i-1][j-1] + c[i-1][j])%MOD ;
}

int n , m;

LL dp[MAX][MAX];

int main ( )
{
    init ();
    while ( ~scanf ( "%d%d" , &n , &m ) )
    {
        memset ( dp , 0 , sizeof ( dp ) );
        dp[0][0] = 1;
        for ( int i = 1 ; i <= n ; i++ )
            for ( int j = 1 ; j <= m ; j++ )
                for ( int k = 0 ; k <= j ; k++ )
                    for ( int t = j-k>0?j-k:1 ; t <= j ; t++ )
                    {
                        dp[i][j] += dp[i-1][k]*c[m-k][j-k]%MOD*c[k][t-j+k]%MOD;
                        dp[i][j] %= MOD;
                    }
        printf ( "%lld\n" , dp[n][m] );
    } 
}