bzoj3028 食物

题目

母函数，数学题，大概就是推公式吧。

承德汉堡:$1+x^2+x^4+\dots=\frac{1}{1-x^2}$
可乐:$1+x$
鸡腿:$1+x+x^2$
蜜桃多:$x+x^3+x^5+\dots=\frac{x}{1-x^2}$
鸡块:$1+x^4+x^8+\dots=\frac{1}{1-x^4}$
包子:$1+x+x^2+x^3$
土豆片炒肉:$1+x$
面包:$1+x^3+x^6+\dots=\frac{1}{1-x^3}$

然后全部乘起来
可得$\frac{x}{(1-x)^4}$是总的函数。
接下来就看它是什么数列的母函数就好了。
我们先对它求n阶导。
$f(x)=\frac{x}{(1-x)^4}$
莱布尼兹公式展开得
$f^{(n)}(x)=\frac{\frac{(n+3)!}{3!}x}{(1-x)^{n+4}}+n\frac{\frac{(n+2)!}{3!}}{(1-x)^{n+3}}$
则用泰勒展开后$x^n$项系数为$\frac{f^{(n)}(0)}{n!}$
代入得：
$g_{n}=n\frac{(n+2)!}{3!n!}=C_{n+2}^{3}=\frac{(n+2)(n+1)n}{6}$
最后，$g_{n}$为答案。

#include<bits/stdc++.h>
using namespace std;
const int mod=10007,inv=1668;
char s[1005];
int x,ans;
int main()
{
    scanf("%s",s);
    for(int i=0;i<strlen(s);i++)x=(x*10+s[i]-'0')%mod;
    ans=x;
    ans=ans*(x+1)%mod;
    ans=ans*(x+2)%mod;
    ans=ans*inv%mod;
    cout<<ans;
    return 0;
}