NYOJ 571 整数划分（三）

题目链接：http://acm.nyist.net/JudgeOnline/problem.php?pid=571

转载出处：http://www.cnblogs.com/xiaoxian1369/archive/2011/09/12/2174212.html

整数划分 －－－ 一个老生长谈的问题:

　　1) 练练组合数学能力.
　　2) 练练递归思想
　　3) 练练DP
　　总之是一道经典的不能再经典的题目:
　　这道好题求:
　　1. 将n划分成若干正整数之和的划分数。
　　2. 将n划分成k个正整数之和的划分数。
　　3. 将n划分成最大数不超过k的划分数。
　　4. 将n划分成若干奇正整数之和的划分数。
　　5. 将n划分成若干不同整数之和的划分数。
 
1.将n划分成不大于m的划分法： 
 　　1).若是划分多个整数可以存在相同的：
 　　 dp[n][m]= dp[n][m-1]+ dp[n-m][m]  dp[n][m]表示整数 n 的划分中，每个数不大于 m 的划分数。
     　　则划分数可以分为两种情况:
     　　a.划分中每个数都小于 m，相当于每个数不大于 m- 1, 故划分数为 dp[n][m-1].
    　　 b.划分中有一个数为 m. 那就在 n中减去 m ,剩下的就相当于把 n-m 进行划分， 故划分数为 dp[n-m][m];
　　2).若是划分多个不同的整数：
　　dp[n][m]= dp[n][m-1]+ dp[n-m][m-1]   dp[n][m]表示整数 n 的划分中，每个数不大于 m 的划分数。
   　　 同样划分情况分为两种情况：
    　　a.划分中每个数都小于m,相当于每个数不大于 m-1,划分数为 dp[n][m-1].
    　　b.划分中有一个数为 m.在n中减去m,剩下相当对n-m进行划分，
　　　并且每一个数不大于m-1，故划分数为 dp[n-m][m-1]
　　2.将n划分成k个数的划分法：
　　　　dp[n][k]= dp[n-k][k]+ dp[n-1][k-1];
   　　方法可以分为两类：
   　　　　第一类: n 份中不包含 1 的分法，为保证每份都 >= 2，可以先拿出 k 个 1 分
   　　到每一份，然后再把剩下的 n- k 分成 k 份即可，分法有: dp[n-k][k]
  　　 　　第二类: n 份中至少有一份为 1 的分法，可以先那出一个 1 作为单独的1份，剩
   　　下的 n- 1 再分成 k- 1 份即可，分法有：dp[n-1][k-1]
　　
　　3.将n划分成若干奇数的划分法:
　　　　g[i][j]:将i划分为j个偶数
　　　　f[i][j]:将i划分为j个奇数
　　   g[i][j] = f[i - j][j];
   　　f[i][j] = f[i - 1][j - 1] + g[i - j][j];
　　　　方法可以分为两类：
　　　　　　a.对于偶数g[i][j]: j份 中不包含 1 的分法，为保证每份都 >= 2，可以先拿出 j 个 1 分
　　　　到每一份，然后再把剩下的 i - j 分成 j 个奇数即可，分法有: f[i-j][j]
　　 　　　 b.对于奇数f[i][j]: j 份中至少有一份为 1 的分法，可以先那出一个 1 作为单独的1份，剩
　　　　下的 i-1 再分成 j- 1 份即f[i - 1][j - 1],还有就是没有一份为 1 的分法可以先拿出 j 个 1 分

　　　　到每一份，然后再把剩下的i - j 分成 j 个偶数即g[i-j][j]，分法有：g[i-1][j-1] + f[i-j][j]

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn = 50 + 5;
int dp1[maxn][maxn]; //表示整数i化成不大于整数j的种数有多少
int dp2[maxn][maxn]; //表示整数i化成不大于整数j且每个数都不相等的种数有多少
int dp3[maxn][maxn]; //表示整数i化成j个整数的种类有多少
int f[maxn][maxn];  //表示整数i化成j个奇数的种数有多少
int g[maxn][maxn];  //表示整数i化成j个偶数的种数有多少
void inin()
{
    memset(dp1, 0, sizeof(dp1));
    memset(dp2, 0, sizeof(dp2));
    memset(dp3, 0, sizeof(dp3));
    memset(f, 0, sizeof(f));
    memset(g, 0, sizeof(g));
    for(int i = 1; i <= 50; i++)
        for(int j = 1; j <= 50; j++)
            if(i > j)
            {
                dp1[i][j] = dp1[i][j-1] + dp1[i-j][j];
                dp2[i][j] = dp2[i][j-1] + dp2[i-j][j-1];
                dp3[i][j] = dp3[i-j][j] + dp3[i-1][j-1];
            }
            else if(i == j)
            {
                dp1[i][j] = dp1[i][j-1] + 1;//当i==j时， dp1[i][j] 就比dp1[i][j-1] 多了一个单独j的一种情况， 所以为dp1[i][j-1]+1；
                dp2[i][j] = dp2[i][j-1] + 1;<span style="font-family: Arial, Helvetica, sans-serif;">//当i==j时， dp2[i][j] 就比dp2[i][j-1] 多了一个单独j的一种情况， 所以为dp[2i][j-1]+1；</span>
                dp3[i][j] = 1;  //只有一种，且全部都为1
            }
            else
            {
                dp1[i][j] = dp1[i][i];   //把i分成不大于j的种数最多时i == j; j再增加i也不能再划分了， 故当j>i时， dp1[i][j] = dp[i][i]
                dp2[i][j] = dp2[i][i];  //把i分成不大于j的种数最多时i == j; j再增加i也不能再划分了， 故当j>i时， dp1[i][j] = dp[i][i]
                dp3[i][j] = 0; //总维i，分成j个是不可能的，所以为零
            }
    for(int i = 1; i <= 50; i++)
        for(int j = 1; j <= i; j++)
            if(i > j)
            {
                g[i][j] = f[i-j][j];
                f[i][j] = f[i-1][j-1] + g[i-j][j];
            }
            else
            {
                f[i][j] = 1; //当i与j相等时 组合也就只有一种，即全部都为以1
                g[i][j] = 0; //不可能，所以为0
            }
}
int main()
{
    inin();
    int n, k;
    while(~scanf("%d%d", &n, &k))
    {
        int ans1 = dp1[n][n];
        int ans2 = dp3[n][k];
        int ans3 = dp1[n][k];
        int ans4 = 0;
        for(int i = 1; i <= n; i++)
            ans4 += f[n][i];
        int ans5 = dp2[n][n];
        printf("%d\n%d\n%d\n%d\n%d\n\n", ans1, ans2, ans3, ans4, ans5);
    }
    return 0;
}