bzoj3028:食物

传送门

生成函数的模板题

---

前置知识：

\[\sum_{i=0}^{+\infty}x^i=\frac{1}{1-x} \]

其实就是等比数列求和公式，这就是公比为\(x\)的等比数列，然后取\(x\in(-1,1)\)

也就是你只要会等比数列求和就行了

也就有

\[(1+x+x^2+x^3+x^4+...)^k=(1-x)^k \]

然后还有一个结论，就是\((1+x+x^2+x^3+x^4+...)^k\)的第\(i\)项系数就是\(\binom{i+k-1}{k-1}\)

证明的话插板法就好了

---

承德汉堡：\(\frac{1}{1-x^2}\)

可乐：\(\frac{1-x^2}{1-x}\)

鸡腿：\(\frac{1-x^3}{1-x}\)

蜜桃多：\(\frac{x}{1-x^2}\)

鸡块：\(\frac{1}{1-x^4}\)

包子：\(\frac{1-x^4}{1-x}\)

土豆片炒肉：\(\frac{1-x^2}{1-x}\)

面包：\(\frac{1}{1-x^3}\)

然后乘起来就是\(\frac{x}{(1-x)^4}\)

设

\[f(x)=\frac{x}{(1-x)^4}\\ f(x)=x(1+x+x^2+x^3+...)^4 \]

那么我们要求的也就是\(f(x)\)的第\(n-1\)项系数，也就是\(\binom{n+2}{3}\)

实现的话其实也没必要写高精度，读入的时候按位读，边读边模就好了

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
void read(int &x) {
	char ch; bool ok;
	for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
	for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define rg register
const int maxn=1e3+10,mod=10007,inv=1668;
int n,len,m=1;char a[maxn];
int main()
{
	scanf("%s",a+1),len=strlen(a+1);
	for(rg int i=len;i;i--){
		n=(n+m*(a[i]-'0'))%mod;
		m=m*10%mod;
	}
	printf("%d\n",n*(n+1)%mod*(n+2)%mod*inv%mod);
}