首页 > 代码库 > bzoj 3028: 食物 -- 母函数

bzoj 3028: 食物 -- 母函数

3028: 食物

Time Limit: 3 Sec  Memory Limit: 128 MB

Description

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

Input

输入样例1
  1
输出样例1
  1
 
输入样例2
  5
输出样例2
  35
 数据范围
   对于40%的数据,1<=N<=100000;
   对于所有数据,1<=n<=10^500;

HINT

关于母函数,要膜拜一下百度百科,讲解的非常详细,在这里我就不赘述了  http://baike.baidu.com/link?url=ks_0J2hQkdV4Bx-3BHgEDyAL7WSJjHcW86sdmMPfZfvaljHA7_-0IQbXhcap6kdOLyYyO8Rn1amc6wgVvR9EZuJl56bOCKpCtVNim8e0ig9L6XpaS-y8wj-FfBnj_qeK

首先我们先求出每个东西的母函数

技术分享

相乘,化简得

                 技术分享

我们知道

      技术分享

所以

      技术分享

所以答案

      技术分享

 

 

 >_< 不知道为什么用 gets 输入就WA了。。

#include<cstdio>#define P 10007int n;char s[510];int main(){    scanf("%s",s+1);     for(int i=1;s[i];i++)(n=(n<<1)+(n<<3)+(s[i]-0))%=P;    printf("%d\n",n*(n+1)%P*(n+2)%P*1668%P);}

 

bzoj 3028: 食物 -- 母函数