ACM学习历程—HDU1028 Ignatius and the Princess（组合数学）

最新推荐文章于 2019-09-20 10:40:35 发布

weixin_30655219

最新推荐文章于 2019-09-20 10:40:35 发布

阅读量49

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/andyqsmart/p/4101747.html

本文探讨了一个整数分解问题，即给定一个正整数N，求出所有可能的N的正整数和形式的组合数量。通过引入动态规划方法，特别是使用插板法和状态转移方程，作者提出了一种高效求解此类问题的算法，并提供了具体实现代码。文章详细解释了算法背后的逻辑，包括如何避免重复计数和优化计算过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Ignatius and the Princess

Description

"Well, it seems the first problem is too easy. I will let you know how foolish you are later." feng5166 says.
"The second problem is, given an positive integer N, we define an equation like this: N=a[1]+a[2]+a[3]+...+a[m]; a[i]>0,1<=m<=N; My question is how many different equations you can find for a given N. For example, assume N is 4, we can find: 4 = 4; 4 = 3 + 1; 4 = 2 + 2; 4 = 2 + 1 + 1; 4 = 1 + 1 + 1 + 1; so the result is 5 when N is 4. Note that "4 = 3 + 1" and "4 = 1 + 3" is the same in this problem. Now, you do it!"

Input

The input contains several test cases. Each test case contains a positive integer N(1<=N<=120) which is mentioned above. The input is terminated by the end of file.

Output

For each test case, you have to output a line contains an integer P which indicate the different equations you have found.

Sample Input

4

10

20 
Sample Output
 5 
42 
627 
 

 刚看到题目是k个数相加为n的题目，第一反应是求1到n元的一次不定方程求解，用插板的方法就能求出通项，但是考虑到这k元未知数是没有先后关系的，也就是说任意(a,a)和(a,a)是相同的。于是不能这样考虑。
然后考虑到每组解必然有个最大值，于是设了这样一个函数（或者数组）f(n,k)，表示和为n且最大数为k的上述元数不定的方程的解的个数。这样一来题目要求求得就是f(n,1)到f(n,n)的和了。
然后我们考虑，对于f(n,k)如果去掉最大值k，那么其子问题就是求和为n-k且最大值小于等于k的上述方程的解的个数。 于是得到递推方程：
但是考虑到n-k或许会比k还要小，于是更新递推方程为：
然后打表求f(n,k)即可，此外用f(n,0)来存所有f(n,i)的和。用数组实现。


代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <vector>
#define inf 0x3fffffff

using namespace std;

int a[125][125];

int main()
{
    //freopen ("test.txt", "r", stdin);
    memset (a, 0, sizeof(a));
    for (int i = 1; i <= 120; ++i)
    {
        for (int j = 1; j <= i; ++j)
        {
            if (i == 1 || i == j)
            {
                a[i][j] = 1;
            }
            else
            {
                for (int k = 1; k <= j && k <= i-j; ++k)
                {
                    a[i][j] += a[i-j][k];
                }
            }
            a[i][0] += a[i][j];
        }
    }
    int n;
    while(scanf("%d", &n) != EOF)
    {
        printf ("%d\n", a[n][0]);
    }
    return 0;
}

 