动态规划DP——整数划分

题目描述
一个正整数n可以表示成若干个正整数之和，形如：n=n1+n2+…+nk，其中n1≥n2≥…≥nk,k≥1

我们将这样的一种表示称为正整数n的一种划分。

现在给定一个正整数n，请你求出n共有多少种不同的划分方法。

输入格式
共一行，包含一个整数n。

输出格式
共一行，包含一个整数，表示总划分数量。

由于答案可能很大，输出结果请对10^9+7取模。

数据范围
1≤n≤1000

输入样例:
5

输出样例:
7

思路：把1,2,3, … n分别看做n个物体的体积，这n个物体均无使用次数限制，问恰好能装满总体积为n的背包的总方案数（完全背包问题变形）

初值问题：
求最大值时，当都不选时，价值显然是 0
而求方案数时，当都不选时，方案数是 1（即前 i 个物品都不选的情况也是一种方案），所以需要初始化为 1
即：for (int i = 0; i <= n; i ++) f[i][0] = 1;
等价变形后： f[0] = 1

状态计算：

f[i][j]f[i][j] 表示前i个整数（1,2…,i）恰好拼成j的方案数
求方案数：

把集合选0个i，1个i，2个i，…全部加起来

f[i][j] = f[i - 1][j] + f[i - 1][j - i] + f[i - 1][j - 2 * i] + ...;
f[i][j - i] =           f[i - 1][j - i] + f[i - 1][j - 2 * i] + ...;

因此

f[i][j]=f[i−1][j]+f[i][j−i];

(这一步类似完全背包的推导）

朴素做法

// f[i][j] = f[i - 1][j] + f[i][j - i]

#include <iostream>

using namespace std;

const int N = 1e3 + 7, mod = 1e9 + 7;

int f[N][N];

int main() {
    int n;
    cin >> n;

    for (int i = 0; i <= n; i ++) {
        f[i][0] = 1; // 容量为0时，前 i 个物品全不选也是一种方案
    }

    for (int i = 1; i <= n; i ++) {
        for (int j = 0; j <= n; j ++) {
            f[i][j] = f[i - 1][j] % mod; // 特殊 f[0][0] = 1
            if (j >= i) f[i][j] = (f[i - 1][j] + f[i][j - i]) % mod;
        }
    }

    cout << f[n][n] << endl;
}

等价变形

// f[i][j] = f[i - 1][j] + f[i][j - i]

#include <iostream>

using namespace std;

const int N = 1e3 + 7, mod = 1e9 + 7;

int f[N];

int main() {
    int n;
    cin >> n;


    f[0] = 1; // 容量为0时，前 i 个物品全不选也是一种方案

    for (int i = 1; i <= n; i ++) {
        for (int j = i; j <= n; j ++) {
            f[j] = (f[j] + f[j - i]) % mod;
        }
    }

    cout << f[n] << endl;
}

第二种做法（非背包问题做法）

最小值1  表示  组成总和是 i 的里面最小的数为1

#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1e4+7,mod = 1e9+7;
int f[N][N];
int main()
{
    int n; scanf("%d",&n);
    f[0][0] = 1;
    for(int i = 1 ; i <= n ;i++)
        for(int j = 1; j <= i; j++)
            f[i][j] = (f[i-1][j-1] + f[i-j][j]) % mod;

    int res = 0;
    for(int i = 1 ; i <= n ; i++) res =(res + f[n][i]) % mod;
    cout << res << endl;
    return 0;
}