题目链接: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;
}