sgu 269 分类： sgu 2015-06...

最新推荐文章于 2015-06-27 13:49:20 发布

转载最新推荐文章于 2015-06-27 13:49:20 发布 · 150 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

原文链接：http://www.cnblogs.com/dashgua/p/4722983.html

本文介绍了解决SGU269问题的一种动态规划(DP)方法，通过排序消除后效性，使用dp[i][j]表示前i行，用j个棋子的方案数，详细阐述了DP状态转移方程，并提供了完整的C++代码实现。

类似 sgu221 的 DP

首先对 b 排序，消除后效性，

dpi,j 表示前 i 行，用 j 个棋子的方案数，

dp0,0=1，dpi,j=dpi−1,j+dpi−1,j−1∗(bi−j+1)

dpn,k 就是问题的答案。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>

const int maxl = 300, maxn = 255, maxk = maxn, Base = 10000;

struct BigNum
{
    int num[maxl], l;

    void digit(int t)
    {
        num[t+1] += num[t]/Base, num[t] %= Base;
    }
    BigNum operator *(const int b)
    {
        BigNum ret = (*this);
        if(!l || !b){ret.l = 0; return ret;}
        for(int i = 1; i <= l; i++) ret.num[i] *= b;
        for(int i = 1; i <= l; i++) ret.digit(i);
        while(ret.num[ret.l+1]) ++ret.l, ret.digit(ret.l);   
        return ret;     
    }
    void operator +=(const BigNum add)
    {
        l = std::max(l, add.l);
        for(int i = 1; i <= l; i++) num[i] += add.num[i];
        for(int i = 1; i <= l; i++) digit(i);
        if(num[l+1]) l++;   
    }
    void prt()
    {
        printf("%d",num[l]);
        for(int i = l-1; i >= 1; i--)
            printf("%04d", num[i]);
    }
}emp, dp[maxk];

int n, k, b[maxn];

int main()
{
#ifndef ONLINE_JUDGE    
    freopen("sgu269.in","r",stdin);
    freopen("sgu269.out","w",stdout);
#endif

    std::cin >> n >> k;

    for(int i = 1; i <= n; i++)
        std::cin >> b[i];

    std::sort(b + 1, b + n + 1);

    dp[0].num[1] = dp[0].l = 1;

    for(int i = 1; i <= n; i++)
    {
        for(int j = k; j >= 1; j--)
            dp[j] += dp[j-1]*(b[i]-j+1);
    }

    dp[k].prt();

#ifndef ONLINE_JUDGE
    fclose(stdin);
    fclose(stdout); 
#endif
    return 0;   
}

转载于:https://www.cnblogs.com/dashgua/p/4722983.html