hdu 1028整数划分

1、母函数

杭州电子科技大学课件

算法就是简单的模拟，

(1+x+x^2+x^3+...)(1+x^2…+x^4+x^6+...)(1+x^3…+x^6+x^9+)...(...)

=[1*(1+x^2…+x^4+x^6+...)+x*(1+x^2…+x^4+x^6+...)+x^2*(1+x^2…+x^4+x^6+...]*(1+x^3…+x^6+x^9+)...(...)

#include <iostream>
#include <cstdio>
#include <memory.h>
using namespace std;
const int maxn=125;
int main()
{
    //freopen("in.txt","r",stdin);
    int n;
    int a[maxn],b[maxn];
    while(scanf("%d",&n)!=EOF){
        memset(b,0,sizeof(b));
        for(int i=0;i<=n;++i){
            a[i]=1;b[i]=0;
        }
        for(int i=2;i<=n;++i){
            for(int j=0;j<=n;++j){
                for(int k=0;k+j<=n;k+=i){
                    b[k+j]+=a[j];//(...)*(x^i+x^2i+x^3i+...),k=i,2i,3i,..
                }
            }
            for(int j=0;j<=n;++j){
                a[j]=b[j];
                b[j]=0;//b[j]注意清空，因为每次从系数为1开始乘//(x^i+x^2i+x^3i+...)

            }
        }
        printf("%d\n",a[n]);
    }
    return 0;
}

2、dp
A.递归，会超时，但比较好理解
设n的一个划分为{m1,m2,...,mi}，其中m=max{m1,m2,...,mi}
dp(n,m)表示n的划分中最大的数不超过m的方案数
则dp(n,m)有以下几种情况
（1）n==1，只有一种划分即{1};
（2）m==1，只有一种划分即n个1，{1,1,1,...,1}; 
（3）n<m，此时dp(n,m)就等于dp(n,n)
（4）n==m，此时dp(n,m)=dp(n,n)=dp(n,n-1)+1
（5）n>m，根据划分中是否包含最大值m，可以分为两种情况： 
         划分中包含m的情况，即{m, {x1,x2,...xi}}, 其中{x1,x2,... xi} 的和为n-m，可能再次出现m，因此是（n-m）的m划分，因此这种划分个数为f(n-m, m); 
        划分中不包含m的情况，则划分中所有值都比m小，即n的(m-1)划分，个数为f(n,m-1);因此 f(n, m) = f(n-m, m)+f(n,m-1);
#include <iostream>
#include <cstdio>
using namespace std;

const int maxn=125;
int dp[maxn][maxn];
int f(int n,int m)//递归会超时
{
    if(n==1||m==1)return 1;
    if(n==m)return f(n,n-1)+1;
    if(n<m)return f(n,n);
    if(n>m)return f(n,m-1)+f(n-m,m);
}

int main()
{
    //freopen("in.txt","r",stdin);
    int n;
    while(scanf("%d",&n)!=EOF){
        //printf("%d\n",f(n,n));
        for(int i=1;i<=n;++i){
            dp[1][i]=1;
            dp[i][1]=1;
        }
        for(int i=1;i<=n;++i){
            for(int j=1;j<=n;++j){
                if(i<j)dp[i][j]=dp[i][i];
                else if(i==j)dp[i][j]=dp[i][j-1]+1;
                else if(i>j)dp[i][j]=dp[i][j-1]+dp[i-j][j];
            }
        }
        printf("%d\n",dp[n][n]);
    }
    return 0;
}