母函数的分析与总结
转自:
关于什么是母函数 , 以及在现实生活中的应用 , 请大家详看 或者 HDU 母函数 PPT:
http://www.cppblog.com/MiYu/archive/2010/08/05/122290.html
对于给出的母函数模板 , 让人理解起来比较费劲的!以下给出几种解释 , 和自己理解!
//made by syx
//time 2010年9月11日 10:17:27
//母函数例题
/*//整数拆分模板
#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];
}
// 刷新一下数据,继续下一次计算,就是下一个括号里面的每一项。
for ( j=0; j<=n; j++ )
{
c1[j] = c2[j] ;
c2[j] = 0 ;
}
}
cout<<c1[n]<<endl;
}
return 0;
}
这句 c2[j+k] += c1[j];的理解还要自己好好的体会体会啊!
*/
自己理解:对于(#式) (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的作用!
最后刷新下结果 , 下一组数据计算
1 //整数拆分母函数模板 2 #include <iostream> 3 using namespace std; 4 const int lmax=10000; 5 int c1[lmax+1],c2[lmax+1]; 6 int main() 7 { int n,i,j,k; 8 while (cin>>n) 9 { 10 //首先对c1初始化,由第一个表达式(1+x+x2+..xn)初始化, 11 //把质量从0到n的所有砝码都初始化为1. 12 for (i=0;i<=n;i++) 13 { 14 c1[i]=1; 15 c2[i]=0; 16 } 17 18 //i从2到n遍历,这里i就是指第i个表达式, 19 //上面给出的第二种母函数关系式里,每一个括号括起来的就是一个表达式。 20 for (i=2;i<=n;i++) 21 { 22 //j 从0到n遍历,这里j就是只一个表达式 23 //里第j个变量,比如在第二个表达式里:(1+x2+x4….)里,第j个就是x2*j. 24 for ( j=0;j<=n;j++) 25 for (k=0;k+j<=n;k+=i)//k表示的是第j个指数,所以k每次增i(因为第i个表达式的增量是i)。 26 c2[ j+k]+=c1[j]; 27 //把c2的值赋给c1,而把c2初始化为0,因为c2每次是从一个表达式中开始的 28 for (j=0;j<=n;j++) 29 { 30 c1[j]=c2[j]; 31 c2[j]=0; 32 } 33 } 34 cout<<c1[n]<<endl; 35 } 36 37 return 0; 38 } 39 40 /*加个例题! 41 //HDU 1085 42 #include <iostream> 43 using namespace std; 44 45 int c1[10000], c2[10000]; 46 int num[4]; 47 int main() 48 { 49 int nNum; 50 while(scanf("%d %d %d", &num[1], &num[2], &num[3]) && (num[1]||num[2]||num[3])) 51 { 52 int _max = num[1]*1+num[2]*2+num[3]*5; 53 // 初始化 54 for(int i=0; i<=_max; ++i) 55 { 56 c1[i] = 0; 57 c2[i] = 0; 58 } 59 for(int i=0; i<=num[1]; ++i) 60 c1[i] = 1; 61 for(int i=0; i<=num[1]; ++i) 62 for(int j=0; j<=num[2]*2; j+=2) 63 c2[j+i] += c1[i]; 64 for(int i=0; i<=num[2]*2+num[1]*1; ++i) // 看到范围的变化了吗? 65 { 66 c1[i] = c2[i]; 67 c2[i] = 0; 68 } 69 70 for(int i=0; i<=num[1]*1+num[2]*2; ++i) 71 for(int j=0; j<=num[3]*5; j+=5) 72 c2[j+i] += c1[i]; 73 for(int i=0; i<=num[2]*2+num[1]*1+num[3]*5; ++i) 74 { 75 c1[i] = c2[i]; 76 c2[i] = 0; 77 } 78 int i; 79 80 for(i=0; i<=_max; ++i) 81 if(c1[i] == 0) 82 { 83 printf("%d\n", i); 84 break; 85 } 86 if(i == _max+1) 87 printf("%d\n", i); 88 } 89 return 0; 90 } 91 */ 92 93 /* 94 //HDOJ_1398 Square Coins 95 //G(x)=(1+x+x2+x3+x4+…)(1+x4+x8+x12+…)(1+x9+x18+x27+…)… 96 97 #include <iostream> 98 using namespace std; 99 const int lmax=300; 100 int c1[lmax+1],c2[lmax+1]; 101 int main(void) 102 { 103 int n,i,j,k; 104 while (cin>>n && n!=0) 105 { 106 for (i=0;i<=n;i++) 107 { 108 c1[i]=1; 109 c2[i]=0; 110 } 111 for (i=2;i<=17;i++) 112 { 113 for (j=0;j<=n;j++) 114 for (k=0;k+j<=n;k+=i*i) 115 { 116 c2[j+k]+=c1[j]; 117 } 118 for (j=0;j<=n;j++) 119 { 120 c1[j]=c2[j]; 121 c2[j]=0; 122 } 123 } 124 cout<<c1[n]<<endl; 125 } 126 return 0; 127 } 128 129 */
有模板好……