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

关注

原创最新推荐文章于 2023-11-28 19:49:43 发布 1.1k 阅读 AI 写同款· GEO 优化 ›

限时加码！20+主流AI编程工具免费用购周边加赠Coding Plan Lite，Claude Code、Cursor等即刻畅享，学习进阶更高效！阅读详情

这个神题

汉堡= $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 ) ! 3 ! (1 - x) - (n + 3) + ( n + 3 ) ! 3 ! x (1 - x) - (n + 4)

$f^n(x)=n{(n+2)!\over3!}(1-x)^{-(n+3)}+{(n+3)!\over3!}x(1-x)^{-(n+4)}$
化简

fn(0)n! $f^n(0) \over n!$ 得就是

n(n+1)(n+2)6=C3n+2 ${n(n+1)(n+2)\over 6}=C_{n+2}^3$

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

x ( 1 - x ) 4 = x (1 + x 1 + x 2 + x 3 + \dots) 4

${x \over (1-x)^4}=x(1+x^1+x^2+x^3+\cdots)^4$

[xn](1+x1+x2+x3+⋯)4 $[x^n](1+x^1+x^2+x^3+\cdots)^4$ 就是把n分成4部分的方案数
那么答案就是n-1分成4部分的方案数由插板法得

A n s = C 2 n - 1 + 3 = C 3 n + 2

$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;
}