hdu 2082 母函数模板

转:http://blog.csdn.net/dellaserss/article/details/7757056#

终于算弄明白母函数是什么东西了 = = , 呼~

感觉它利用了高中学的二项式定理(就是用到组合数的那个)

完完全全的利用啊!以前用二项式求出最后那超长的x的表达式( 1 + C(n,1)*x^1 + C(n,2)*x^2 +~~~+ C(n,n-1)*x^(n-1) + C(n,n)*x^n )之后好歹还会给x赋个值算个答案
到母函数这连赋值都省了……x就代表一个实体,⊙﹏⊙b汗

它只汲取了二项式定理的精华,从每个因子中的取一个系数进行组合的思想。至于怎么实现的 ,给两个传送门 hdu 2082 母函数的应用分析 Tanky Woo 建议先看后面这个,看不懂了再看前面的,

下面给出hdu2082的AC代码
根据自己对代码的理解添加修改了部分注释
凑活着看吧 - -
~~

#include<iostream>  
using namespace std;  
int main()  
{  
    int i,j,t,num,sum;     
    int X[27];                     //保存Xi   

    int c1[51],c2[51];     /*c1[51] 存储指数为1到 50 的项的系数  例如c1[5]表示  ax^5这一项的系数a*/  


    cin>>t;  
    while(t--)  
    {  
        for(i=1;i<=26;i++)  
            cin>>X[i];  
        memset(c1,0,sizeof(c1));  
        memset(c2,0,sizeof(c2));  
        for(i=1;i<=26;i++)                       //首先对第一个因子的各个系数初始化(称上面的式子中用一个括号括起来的为一个因子)  
            if(X[i]!=0)                                         //Xi为0表示第i个因子为0(不存在) ;用第一个因子初始化c1   
            {  
                for(num=0;num<=X[i];num++)  
                {  
                    if(num*i>50)  
                        break;                                        //指数相乘如果大于50 则初始化完毕跳出  
                    c1[num*i]=1;                          //第一个因子里面的各个项的系数当然为1   
                }  
                break;  
            }  
        for(i=i+1;i<=26;i++)                                  //从第二个乘数因子开始 字母开始做接下去的26-i次乘法  
        {  
        //每次循环模拟两个因子相乘, c1储存了前i个因子相乘得到的结果(每一项的系数)   

            for(j=0;j<=50;j++)                                //每一次模拟的计算范围控制在指数小于50  
                if(c1[j]!=0)                              //c1[j]==0 ,说明指数为j的那一项的系数为零,即该项不存在   
                    for(num=0;num<=X[i];num++)              
                    {  
                        if(j+i*num>50) break;          //两个相乘的数的指数相加,第i个因子的第k项的指数为num*i  
                        c2[j+i*num]+=c1[j]*1;            /* 1.c1[i*num+j]+=c1[j] 是错误的写法 ~~ ,因为c1还有利用价值(34行的判断)  
                                                            暂时保存至c2;  
                                                            2.两项相乘 指数相加(j+i*num) ,因数相乘 c1[j]*1;  
                                                         */   
                    }  

            for(j=0;j<=50;j++)                         //把两个因子相乘之后的结果当做一个因子,再重新保存至c1  
            {  
                c1[j]=c2[j];  
                c2[j]=0;  
            }  
        }  

        for(sum=0,i=1;i<=50;i++)                            //前50项之和  
            sum+=c1[i];  

        printf("%d\n",sum);  
    }  
    return 0;  
}  
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值