bzoj3028（母函数）

Description
明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！我们暂且不讨论他有多么NC，他又幻想了他应
该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。他这次又准备带一些受欢迎的食物，
如：蜜桃多啦，鸡块啦，承德汉堡等等当然，他又有一些稀奇古怪的限制：每种食物的限制如下：
承德汉堡：偶数个
可乐：0个或1个
鸡腿：0个，1个或2个
蜜桃多：奇数个
鸡块：4的倍数个
包子：0个，1个，2个或3个
土豆片炒肉：不超过一个。
面包：3的倍数个
注意，这里我们懒得考虑明明对于带的食物该怎么搭配着吃，也认为每种食物都是以‘个’为单位（反正是幻想嘛
），只要总数加起来是N就算一种方案。因此，对于给出的N,你需要计算出方案数，并对10007取模。
Input
输入一个数字N，1<=n<=10^500
Output
如题
Sample Input
输入样例1
1
输入样例2
5
Sample Output
输出样例1
1
输出样例2
35

---

$solution$:
首先构造每一种食物的生成函数，然后用求和之类的办法得到母函数
把所有母函数相乘得到$\frac{x}{(1-x{)}^4}$
显然展开后的第$n$项就是答案
怎么做？
熟练的$oier$可以直接广义二项式得到这个式子第$n$项就是$C_{n+2}^3$

实际上我们知道${\sum }_{i=0}^{\infty }{x}^i=\frac{1}{1-x}$

故原式为$x\ast (1+x+x^2+...+x^{\infty }{)}^4=x\ast ({\sum }_{i=0}^{\infty }{x}^i{)}^4$
显然这个式子的第$n$项系数可以由插板法得到为$C_{n+2}^3$
直接做即可

#include<bits/stdc++.h>
using namespace std;
#define rep(i,j,k) for(int i = j;i <= k;++i)
#define repp(i,j,k) for(int i = j;i >= k;--i)
#define rept(i,x) for(int i = linkk[x];i;i = e[i].n)
#define P pair<int,int>
#define Pil pair<int,ll>
#define Pli pair<ll,int>
#define Pll pair<ll,ll>
#define pb push_back 
#define pc putchar
#define mp make_pair
#define file(k) memset(k,0,sizeof(k))
#define ll long long
int rd()
{
	int num = 0;char c = getchar();bool flag = true;
	while(c < '0'||c > '9') {if(c == '-') flag = false;c = getchar();}
	while(c >= '0' && c <= '9') num = num*10+c-48,c = getchar();
	if(flag) return num;else return -num;
}
const int p = 10007;
char s[510];
int ksm(int a,int x){int now = 1;for(;x;x>>=1,a=1ll*a*a%p) if(x&1) now=1ll*now*a%p;return now;}
int main()
{ 
    scanf("%s",s+1);int n = strlen(s+1),m=0,inv;
    rep(i,1,n) m = (m*10+s[i]-'0')%p;
    inv = ksm(6,p-2);
    m = 1ll*m*(m+1)*(m+2)%p*inv%p;
    printf("%d\n",m);
	return 0;
}