bzoj3028 食物

Description 明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！
我们暂且不讨论他有多么NC，他又幻想了他应该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。
他这次又准备带一些受欢迎的食物，如：蜜桃多啦，鸡块啦，承德汉堡等等 当然，他又有一些稀奇古怪的限制： 每种食物的限制如下：
承德汉堡：偶数个
可乐：0个或1个
鸡腿：0个，1个或2个
蜜桃多：奇数个
鸡块：4的倍数个
包子：0个，1个，2个或3个
土豆片炒肉：不超过一个。
面包：3的倍数个

注意，这里我们懒得考虑明明对于带的食物该怎么搭配着吃，也认为每种食物都是以‘个’为单位（反正是幻想嘛），只要总数加起来是N就算一种方案。因此，对于给出的N,你需要计算出方案数，并对10007取模。

生成函数。推导如下：

$$\begin{align} & (1+x^2+x^4+\cdots)(1+x)(1+x+x^2)(x+x^3+x^5+\cdots )(1+x^4+x^8+\cdots)(1+x+x^2+x^3)(1+x)(1+x^3+x^6+\cdots)\\ = & \frac{(1+x)^2(1+x+x^2)x(1+x+x^2+x^3)}{(1-x^2)^2(1-x^4)(1-x^3)(1-x)(x^2+x+1)} \\ = &\frac x{(1-x)^4} \\ = & x(1-x)^{-4} \\ \end{align}$$
$n$次项系数为
$$\begin{align} & (-1)^{n-1}C_{-4}^{n-1}\\ = & (-1)^{n-1} \frac{(-4-n+1)(-4-n+2)\cdots (-4)}{n!} \\ = &\frac{4\times 5\times \cdots \times(n+3)}{n!}\\ = & C_{n+2} ^n\\ = & \frac{(n+2)(n+1)n} 6\\ \end{align}$$

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int p=10007,inv=1668;
int rd()
{
    int x=0;
    char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9')
    {
        x=(x*10+c-'0')%p;
        c=getchar();
    }
    return x;
}
int main()
{
    int n=rd();
    printf("%d\n",n*(n+1)%p*(n+2)%p*inv%p);
}