[母函数 泰勒展开] BZOJ 3028 食物

这个神题

汉堡=$x^0+x^2+x^4+\cdots={1 \over {1-x^2}}$
蜜桃=$x^1+x^3+x^5+\cdots={x \over {1-x^2}}$
面包=$x^0+x^3+x^6+\cdots={1 \over {1-x^3}}$
鸡块=$x^0+x^4+x^8+\cdots={1 \over {1-x^4}}$
土豆=$x^0+x^1={1-x^2 \over {1-x}}$
可乐=$x^0+x^1={1-x^2 \over {1-x}}$
鸡腿=$x^0+x^1+x^2={1-x^3 \over {1-x}}$
包子=$x^0+x^1+x^2+x^3={1-x^4 \over {1-x}}$

乘在一起就是$f(x)={x\over (1-x)^4}$
$[x^n]\ f(x)$就是答案
然后我们在$x=0$处泰勒展开
答案就是$f^n(0) \over n!$

这个n阶导数可以找规律得出

$$f^n(x)=n{(n+2)!\over3!}(1-x)^{-(n+3)}+{(n+3)!\over3!}x(1-x)^{-(n+4)}$$
化简 $f^n(0) \over n!$ 得就是 ${n(n+1)(n+2)\over 6}=C_{n+2}^3$

UPD
具体推导可以看这里
其中有一种简洁的方法

$${x \over (1-x)^4}=x(1+x^1+x^2+x^3+\cdots)^4$$
$[x^n](1+x^1+x^2+x^3+\cdots)^4$就是把n分成4部分的方案数
那么答案就是n-1分成4部分的方案数 由插板法得
$$Ans=C_{n-1+3}^2=C_{n+2}^3$$

#include<cstdio>
using namespace std;

const int P=10007;

int main(){
  char s[501]; int n=0;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%s",s);
  for (int i=0;s[i];i++) n=((n<<1)+(n<<3)+s[i]-'0')%P;
  printf("%d\n",n*(n+1)%P*(n+2)%P*1668%P);
  return 0;
}