bzoj3028 食物 (普通型生成函数)

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

思路：

因为是求方案数，所以是普通型生成函数。

求出每个物品的生成函数:
汉堡：1+x²+x⁴…=1/(1-x²) //等比数列求和，指数n趋于无穷，所以分子为1，下同
可乐：1+x
鸡腿：1+x+x²=(1-x³)/(1-x)
蜜桃多：x+x³+x⁵…=x/(1-x²)
鸡块：1+x⁴+x⁸+x¹⁶…=1/(1-x⁴)
包子：1+x+x²+x³=(1-x⁴)/(1-x)
土豆片炒肉：1+x
面包：x³+x⁶+x⁹…=1/(1-x³)
全部乘起来化简一下得：x/(1-x)⁴

然后要做的就是求xⁿ的系数：

下面的没学过

方法1：广义二项式定理：

广义二项式定理可以求出指数为负整数的二项式系数

当-1<=x<=1，且n为正整数时，有式子

原式子为x(1-x)^-4，因为有一个一次项x，所以只需要带入n=4和k=n-1，求得系数为C(n+2,n-1)
又因为C(n,m)=C(n,n-m)
可推出C(n+2,n-1)=C(n+2,3)
只有三项，直接算就行了

方法2：麦克劳林展开：
有空再补

由上得答案为C(n+2,3)
((n+2)(n+1)n)/6
由取模性质可知每个n都要取模，题目的n到达了10⁵⁰⁰
读入的时候直接模拟取模就行了

code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxm=505;
const int mod=1e4+7;
const int inv=1668;//6的逆元
char s[maxm];
signed main(){
    scanf("%s",s);
    int x=0;
    int len=strlen(s);
    for(int i=0;i<len;i++){
        x=(x*10+s[i]-'0')%mod;
    }
    int ans=1;
    for(int i=x+2;i>=x;i--){
        ans=ans*i%mod;
    }
    ans=ans*inv%mod;//记得乘上逆元
    printf("%d\n",ans);
    return 0;
}