bzoj3028 食物（生成函数+Lucas）

Description

明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！我们暂且不讨论他有多么NC，他又幻想了他应该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。他这次又准备带一些受欢迎的食物，如：蜜桃多啦，鸡块啦，承德汉堡等等当然，他又有一些稀奇古怪的限制：每种食物的限制如下：

承德汉堡：偶数个

可乐：0个或1个

鸡腿：0个，1个或2个

蜜桃多：奇数个

鸡块：4的倍数个

包子：0个，1个，2个或3个

土豆片炒肉：不超过一个。

面包：3的倍数个

注意，这里我们懒得考虑明明对于带的食物该怎么搭配着吃，也认为每种食物都是以‘个’为单位（反正是幻想嘛），只要总数加起来是N就算一种方案。因此，对于给出的N,你需要计算出方案数，并对10007取模。

Input

输入一个数字N，1<=n<=10^500

Output

如题 

Sample Input

输入样例1

1

输入样例2

5

Sample Output

输出样例1

1

输出样例2

35

[ Submit][ Status][ Discuss] 

分析：
这道题只要是看过生成函数内容的同学
一定一眼就可以看出这是生成函数的题：

我们可以得到以下的式子：

我们看一下每个式子是怎么化简的：

都乘起来：

我们之所以化出这个式子，是因为有这样一个关系：
设e1,e2,e3,…,e8为各种食物的数量，则答案就是一下方程的解：

e1+e2+e3+….+e8=n

我们化出这个式子后，就有：

x^e1*x^e2*x^e3*….*x^e8=x^n

因此x^n的系数就是n的答案

所以这道题的答案就是：

C(n+2,3)

top

花了好长时间化出了式子，我一看n的数据范围
我的ma呀，这是要高精度的节奏啊
但是模数格外的小
所以可以用Lucas定理来完成

C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p
(C(n%p, m%p)*Lucas(n/p, m/p))%p

//这里写代码片
#include<cstdio>
#include<iostream>
#include<cstring>
#define ll long long

using namespace std;

const int p=10007;    //质数 
int inv[p+10],jc[p+10],n;
char s[510];

int C(int n,int m)
{
    return (jc[n]%p*inv[m]%p*inv[n-m]%p)%p;
}

int main()
{
    scanf("%s",s);

    inv[0]=1; inv[1]=1;
    for (int i=2;i<=p;i++) inv[i]=((p-p/i)%p*inv[p%i]%p)%p; 
    for (int i=2;i<=p;i++) inv[i]=(inv[i]*inv[i-1])%p;
    jc[0]=1; jc[1]=1;
    for (int i=2;i<=p;i++) jc[i]=(i*jc[i-1])%p;

    n=0;
    for (int i=0;i<strlen(s);i++)
        n=((n*10)+s[i]-'0')%p;
    printf("%d\n",C((n+2)%p,3));
    return 0;
}