[BZOJ]3028: 食物生成函数

200815147

于 2019-01-08 13:34:02 发布

阅读量189

点赞数

分类专栏：生成函数

本文链接：https://blog.csdn.net/baidu_36797646/article/details/86071065

版权

生成函数专栏收录该内容

2 篇文章 0 订阅

订阅专栏

Description

明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！我们暂且不讨论他有多么NC，他又幻想了他应该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。他这次又准备带一些受欢迎的食物，如：蜜桃多啦，鸡块啦，承德汉堡等等当然，他又有一些稀奇古怪的限制：每种食物的限制如下：
承德汉堡：偶数个
可乐：0个或1个
鸡腿：0个，1个或2个
蜜桃多：奇数个
鸡块：4的倍数个
包子：0个，1个，2个或3个
土豆片炒肉：不超过一个。
面包：3的倍数个
注意，这里我们懒得考虑明明对于带的食物该怎么搭配着吃，也认为每种食物都是以‘个’为单位（反正是幻想嘛），只要总数加起来是N就算一种方案。因此，对于给出的N,你需要计算出方案数，并对10007取模。

Solution

生成函数最简单的题……把所有生成函数乘起来，最后得到 $x\over(1-x)^4$ 。
而 ${x\over(1-x)^4}=\sum_{i>0}C_{i+3}^3x^{i+1}$ 。
所以最后的答案为 $C_{n+2}^3$ 。

Code

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=510;
const int mod=10007;
char s[Maxn];
int Pow(int x,int y)
{
	if(!y)return 1;
	int t=Pow(x,y>>1),re=t*t%mod;
	if(y&1)re=re*x%mod;
	return re;
}
int main()
{
	scanf("%s",s+1);
	int n=0,len=strlen(s+1);
	for(int i=1;i<=len;i++)n=(n*10%mod+s[i]-'0')%mod;
	printf("%d",n*(n+1)%mod*(n+2)%mod*Pow(6,mod-2)%mod);
}