HDU:1028 Ignatius and the Princess III

HDU:1028 Ignatius and the Princess III


Problem Description

"Well,it seems the first problem is too easy. I will let you know how foolish you arelater." feng5166 says.

"The second problem is, given an positive integer N, we define an equationlike 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

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

 

 

Output

Foreach test case, you have to output a line contains an integer P which indicatethe different equations you have found.

 

 

SampleInput

4

10

20

 

 

SampleOutput

5

42

627

 

题目思路:

就如题目中所说的那样,给出测试数据4,

4 = 4

4 = 1+3

4 = 1+1+2

4 = 2+2

4 = 1+1+1+1

可以构造母函数。想想为什么呢?

我们在数学中学过这样的公式:x^m*x^n= x^(m+n),由此我们知道可以由x的m次方乘上x的n次方,得到x的(m+n)次方,这样可以将乘积转化为和,则这就是m和n的组合。

因此这个题目可以构造母函数:

G(x)=(x^0+x^1+x^2+x^3+x^4)*(x^0+x^2+x^4)*(x^0+x^3)*(x^0+x^4);

=(1+ x^1+x^2+x^3+x^4)*(1+ x^2+x^4)*(1+ x^3)*(1+x^4)

其中x代表数,不管这个数是几,要记住x只是数.

则式子(x^0+x^1+x^2+x^3+x^4)的解释为:x^0代表1取0个,则和为0,则x的幂为0;x^1代表1取了一个,和为1,则x的幂为1,x^2代表1取了两个,和为2,则x的幂为2;x^3代表取了3个1,和为3,则x的幂为3;x^4代表取了4个1,和为4,则x的幂为4;因为要求组成4的组合,因此最多就有4个1,则1的情况已经考虑完了。

同理式子(x^0+x^2+x^4)的解释为,x^0代表数字2取0个,和为0,则x的幂为0;x^2代表数字2取1个,和为2,则x的幂为2,x^4代表数字2取两个,和为4,则x的幂为4;

因为4最多由两个2构成,因此2的个数只能取0,1,2.

同理式子(x^0+x^3)的解释为,x^0代表数字3取0个,和为0,则x的幂为0;x^3代表数字3取1个,和为3,则x的幂为3;因为4中最多包含一个3,则3的个数只取0个或1个。

 同理式子(x^0+x^4)的解释为,x^0代表数字4取0个,和为0,则x的幂为0;x^4代表数字4取1个,和为4,则x的幂为4;因为4中最多包含一个4,则4的个数只取0个或1个。

G(x)=(x^0+x^1+x^2+x^3+x^4)*(x^0+x^2+x^4)*(x^0+x^3)*(x^0+x^4);

=(1+ x^1+x^2+x^3+x^4)*(1+ x^2+x^4)*(1+ x^3)*(1+x^4)

代表0个1与2个2的组合

代表0个1与1个2的组合

代表0个1与0个2的组合

=(1+ x^1+x^2+x^3+x^4)*(1+ x^2+x^4)*( …………………)

 =[(1+1*X^2+1*x^4)+(x^1*1+x^1*x^2+x^1*x^4)+(x^2*1+x^2*x^2+x^2*x^4)+……(x^4*1+x^4*x^2+x^4*x^4)]*(……………………)

=(1+x+2x^2+2x^3+3x^4+2x^5+2x^6+x^7+x^8)* *(……………………)

 

要求的是组成4的方案数,幂为5,6,7,8的相超界直接不予考虑

 =(1+x+2x^2+2x^3+3x^4)*(……………………)

 

这些项代表当前由(0,1,2,3,4)个1与(0,1,2)个2所构成的组合,其系数等于当前的方案数。

                                                             例如x^1系数为1,构成1的只有1

                                                             这一种,x^2的系数为2,构成2的

                                                            方案有两种:1+1,2;构成x^3的方案也有两种,1+1+1,1+2X^4系数为3,代表构成4的方案数有3种。1+1+1+1,1+1+2,2+2这三种。

继续向下求。

=(1+x+2x^2+2x^3+3x^4)*(1+x^3)*()      1,2组合后产生的结果在和3进行组合。

=1+x+2x^2+3x^3+4x^4+2^x5+2x^6+3x^7*() 幂大于4的项不予考虑

=1+x+2x^2+3x^3+4x^4*1+x^4

=1+x+2^2+3x^3+5^x4;

x^4前的系数就是就是构成4的方案数。答案就是5.

明白了母函数的构造以及这样构造的道理,接下来就是怎么将这个思路用代码实现。


int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int i,j,k;
        int a[130]; ///a[i]存放当前组成数字i的方案数
        int b[130]; ///b[i]用来前面已经求出的x^i的系数
        for(i = 0; i <= n; i++)
        {
            a[i]=0;   //初始化,刚开始构成任何数i的方案数都是0
            b[i]=1;  //(x^0+x^1+x^2+x^3+x^4)的系数都是1
        }
        for(i = 2; i <= n; i++)  //从2开始考虑
        {
            for(j = 0; j <= n; j++)
                for(k = 0; k+j <= n; k=k+i)
                    a[j+k] += b[j];  
            for(j = 0; j <= n; j++)
            {
                b[j]=a[j];        ///更新x^j前的系数
                a[j]=0;
            }
        }
        printf("%d\n",b[n]);
    }
    return 0;
}


 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值