母函数的基本模板讲解

母函数的基本代码模板

自己理解:对于(#式)  (1+x+x^2+x^3+x^4+x^5+....)*(1+x^2+x^4+x^6+x^8+x^10+....)*(1+x^3+x^6+x^9+x^12....).....

第一个for给c1 和 c2 赋值 , 把上面#式的第一个括号(1+x+x^2+x^3+x^4+x^5+....)的系数给放在c1中,

从而再次计算从 # 的 第二个括号开始 , 所以 i 就是代表的# 式第几个括号,

而 用程序模拟手工计算 , 就是 先计算第一个括号 与 第二个括号 计算 , 把结果放到c2中,

在把结果与第三个括号计算 , 把结果放到c2中 , 在和第四个括号计算,........

所以j 就是指的 已经计算出 的 各项的系数 ,比如第一次之后 1+x+x^2+x^3+x^4+x^5+... , j=0指向1 ,

j=2 指向x , .... , 而 k 就是指 将要计算的那个括号中的项 , 因为第i个括号 , 中的指数为0 , i , 2i....所以 k要 + i ;

而结果 c2[j+k] += c1[j]; 就是把 以计算出的 j项的系数 和 现在正在计算的括号的k项相乘 , 所以指数为j+k , 所以结果放到c2[j+k] 中 , 这就是这几个for的作用!

最后刷新下结果 , 下一组数据计算!

c1[n]保存着Xn项前的系数,也就是n个(a1,...,an)组合数

//整数拆分模板
#include <iostream>
using namespace std;
const int lmax=10000;
//c1是用来存放展开式的系数的,而c2则是用来计算时保存的,
//他是用下标来控制每一项的位置,比如 c2[3] 就是 x^3 的系数。
//用c1保存,然后在计算时用c2来保存变化的值。
int c1[lmax+1],c2[lmax+1];
int main()
{
            int n, i, j, k ;
           // 计算的方法还是模拟手动运算,一个括号一个括号的计算,
           // 从前往后
           while ( cin>>n )

          {
                     //对于 1+x+x^2+x^3+ 他们所有的系数都是 1
                     // 而 c2全部被初始化为0是因为以后要用到 c2[i] += x ;
                     for ( i=0; i<=n; i++ )

                     {
                                c1[i]=1;
                                c2[i]=0;
                     }
                      //第一层循环是一共有 n 个小括号,而刚才已经算过一个了
                      //所以是从2 到 n
                     for (i=2; i<=n; i++)

                   {
                                 // 第二层循环是把每一个小括号里面的每一项,都要与前一个
                                 //小括号里面的每一项计算。
                                for ( j=0; j<=n; j++ )
                                 //第三层小括号是要控制每一项里面 X 增加的比例
                                 // 这就是为什么要用 k+= i ;
                                         for ( k=0; k+j<=n; k+=i )

                                        {
                                                 // 合并同类项,他们的系数要加在一起,所以是加法,呵呵。
                                                 // 刚开始看的时候就卡在这里了。
                                                 c2[ j+k] += c1[ j];//因为指数为k的项的系数恒为1,所以与指数为j的项相乘之后得到的指数为k+j的项的系数就自加j指数的项的系数,看了半天
                                         }
                               // 刷新一下数据,继续下一次计算,就是下一个括号里面的每一项。
                              for ( j=0; j<=n; j++ )

                              {
                                          c1[j] = c2[j] ;
                                          c2[j] = 0 ;
                              }
                   }
                    cout<<c1[n]<<endl;
        }
         return 0;
}
//这句 c2[j+k] += c1[j];的理解还要自己好好的体会体会啊!


母函数的基本模板讲解

时间: 2024-10-08 07:30:58

母函数的基本模板讲解的相关文章

组合数学 - 母函数的运用 --- 模板题

Holding Bin-Laden Captive! Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 15064    Accepted Submission(s): 6750 Problem Description We all know that Bin-Laden is a notorious terrorist, and he h

组合数学 - 母函数的运用 + 模板 --- hdu : 2082

找单词 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 4093    Accepted Submission(s): 2933 Problem Description 假设有x1个字母A, x2个字母B,..... x26个字母Z,同时假设字母A的价值为1,字母B的价值为2,..... 字母Z的价值为26.那么,对于给定的字母,可以找到

母函数及其应用+模板

部分摘自这位大佬的博客https://www.cnblogs.com/linyujun/p/5207730.html 生成函数即母函数,是组合数学中尤其是计数方面的一个重要理论和工具. 最早提出母函数的人是法国数学家LaplaceP.S.在其1812年出版的<概率的分析理论>中明确提出“生成函数的计算”. 生成函数的应用简单来说在于研究未知(通项)数列规律,用这种方法在给出递推式的情况下求出数列的通项. 在这里我们不去高深地研究数学上母函数,而是讲讲简单的母函数应用. 1.母函数引入 就是把一

母函数入门【模板】

正整数拆分 hdu1028 解: 对于正整数 $n$ 的拆分,其母函数为 $$f(x) = (1+x+x^2+...)(1+x^2+x^4+...)(1+x^3+x^6+x^9+...)...$$ 答案就是多项式展开后 $x^n$ 项的系数. Code: //其实就是模拟,从前往后一一合并 #include<bits/stdc++.h> using namespace std; const int _max = 10001; int c1[_max], c2[_max]; //c1存放前面项计

快速幂模板及讲解

(这篇其实是我用来练习公式编辑器滴,所以讲的内容略水,大佬们也赏脸看看吧) 定义 快速幂即快速求幂(下文为求a的x次幂模m的结果),但我们一般只在要求对一个数的幂取模时才使用,因为有可能结果很大,有可能long long都存不下,但是因为我们有: \((ab)\%m=(a\%m)(b\%m)\) 通过转换,可得: \((a^x)\%m\) \(= (a\times a\times a\times -\times a(共x个a相乘))\%m\) \(= (a\%m)\times (a\times

二分答案模板

[模板+讲解]二分答案 !阅读须知||阅读本博文前笔者认为读者已经学会(或了解)了: 1.基础语言与算法 2.标准二分法(二分思想) 3.二分查找 定义 二分答案与二分查找类似,即对有着单调性的答案进行二分,大多数情况下用于求解满足某种条件下的最大(小)值. 答案单调性 答案的单调性大多数情况下可以转化为一个函数,其单调性证明多种多样,如下: 移动石头的个数越多,答案越大(NOIP2015跳石头). 前i天的条件一定比前 i + 1 天条件更容易(NOIP2012借教室). 满足更少分配要求比满

字符串部分模板总结

KMP kmp处理题型总结 Manacher POJ - 3974 Palindrome (Manacher模板+讲解) 最大最小表示法 HDU-3374 String Problem (最小最大表示法) AC自动机 AC自动机总结 原文地址:https://www.cnblogs.com/Tianwell/p/11523410.html

【网络流-二分图最大匹配】poj3041Asteroids

/* 这道题将每行x看成是结点x,没列y看成是结点y,而障碍物的坐标xy看成是从x到y的 一条边.建图后问题就变成了,找最少的点,使得这些点与所有的边相邻,即最小 点覆盖,用匈牙利算法解决. ------------------------------- 定理:最小点覆盖数 = 最大匹配数,即求图的最大匹配即可,匈牙利算法 ------------------------------- 模板讲解: bool find(int v) { for(int i=1; i<=n; i++) { if(g

Static Simple Graphics

两种五角星画法 homework of Computer Graphic 花花的三角形, 点和线的像素宽度, 写字 同时出现两个窗口 描边 交互式画图 国民党旗 共青团徽 自定义边数的花 五角星 可以和下面的五角星对比一下 参考模板讲解 // #pragma comment (lib, "opengl32.lib") // #pragma comment (lib, "glu32.lib") // #pragma comment (lib, "glut32