模板 2018-01-18 普通型母函数 生成函数

      这个东西在我的博客里面是先有题才有这个模板的, 这里从头讲一遍原理还是有些麻烦, 如果有兴趣的话可以自行搜索, 关键字 普通型母函数 即可.

哦对, 忘了说了, 这个东西叫做生成函数也是可以的.

先是母函数的关键单元, 多项式. 由于母函数名并不会真正用到自变量和因变量的关系, 而是其指数和系数的关系, 所以建立结构体的时候仅需要开一个数字记录系数就行, 下标表示对应的次数, 比如 a[ 7 ]=  8 表示 x^7 的系数是8(你也可以理解成拿7块钱有8种拿法), 然后变量n用于储存最高次项的次数

结构体定义如下:

struct polynomial
{
    long long int a[100];//数组大小视情况决定大小
    int n;
};
typedef polynomial pol;//简写成3个字母的 pol

然后是初始化, 举个例子, 我有三个5元硬币, 应该如何初始化?

void initial_data(pol &data,int value,int n)//第no种硬币,面值为value,有n个
{//学习了C++, 妈妈再也不用担心我的指针和引用啦!
    data.n=n*value;//如果只拿价值value的硬币,最多拿n*value块钱
    for(i=0;i<=n;i++)//中间的跳出循环判断有时候可以用来剪枝,比如最多也只会拿50块可以写一个i*value<=50
    {
        data.a[i*value]=1;
    }
}

当我们把所有的多项式初始化好之后, 接下来就是多项式的相乘和展开了.

pol pol_multi(pol a,pol b)
{
    int i,j;
    pol res;
    memset(res.a,0,sizeof(res.a));
    res.n=0;
    for(i=0;i<=a.n;i++)
    {
        for (j=0;j<=b.n;j++)//这里的跳出循环判断调节可以用来剪枝,如最多只会拿50块钱可以写i+j<=50
        {
            if(a.a[i]*b.a[j]!=0)
                res.n=max(res.n,i+j);//最高次的更新
            res.a[i+j]+=a.a[i]*b.a[j];
        }
    }
    return res;
}