BZOJ 3028: 食物

\(\color{#0066ff}{ 题目描述 }\)

明明这次又要出去旅游了,和上次不同的是,他这次要去宇宙探险!我们暂且不讨论他有多么NC,他又幻想了他应

该带一些什么东西。理所当然的,你当然要帮他计算携带N件物品的方案数。他这次又准备带一些受欢迎的食物,

如:蜜桃多啦,鸡块啦,承德汉堡等等当然,他又有一些稀奇古怪的限制:每种食物的限制如下:

承德汉堡:偶数个

可乐:0个或1个

鸡腿:0个,1个或2个

蜜桃多:奇数个

鸡块:4的倍数个

包子:0个,1个,2个或3个

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

面包:3的倍数个

注意,这里我们懒得考虑明明对于带的食物该怎么搭配着吃,也认为每种食物都是以‘个’为单位(反正是幻想嘛

),只要总数加起来是N就算一种方案。因此,对于给出的N,你需要计算出方案数,并对10007取模。

\(\color{#0066ff}{输入格式}\)

输入一个数字N,1<=n<=\(10^{500}\)

\(\color{#0066ff}{输出格式}\)

如题

\(\color{#0066ff}{输入样例}\)

1
    
5

\(\color{#0066ff}{输出样例}\)

1
    
    
    
35

\(\color{#0066ff}{数据范围与提示}\)

none

\(\color{#0066ff}{ 题解 }\)

生成函数

\(ans = (x^0+x^2+x^4+x^6...)*(x^0+x^1)*(x^0+x^1+x^2)*(x^1+x^3+x^5...)*(x^0+x^4+x^8...)*(x^0+x^1+x^2+x^3)*(x^0+x^1)*(x^0+x^3+x^6...)\)的第n项系数

因为有\(x^0+x^1+x^2+x^3+...=\frac{1}{1-x}\)

所以原式为\(ans = (\frac{1}{1-x^2})*(1+x)*(1+x+x^2)*(\frac{x}{1-x^2})*(\frac{1}{1-x^4})*(1+x+x^2+x^3)*(1+x)*(\frac{1}{1-x^3})\)

还有一些项可以等比数列求和

\(ans = \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}\)

大力约分

\(ans = \frac{x}{(1-x^4)}\)的n次方项系数

这个怎么搞呢?

因为\(\begin{aligned}\frac{1}{(1-x)^{n + 1}}=\sum_{i\ge 0}C_{n+i}^i x^i\end{aligned}\)

因此可以化为\(ans=x*\begin{aligned}\sum_{i\ge 0}C_{3+i}^i x^i\end{aligned}\)

把x弄进去\(ans=\begin{aligned}\sum_{i\ge 0}C_{3+i}^i x^{i +1}\end{aligned}\)

改变一下i的枚举\(ans=\begin{aligned}\sum_{i\ge 1}C_{2+i}^{i - 1} x^{i}\end{aligned}\)

则$ans=\begin{aligned}\sum_{i\ge 1}C_{2+i}^{i - 1} x{i}\end{aligned}的第n项系数=C_{2+n} =C_{2+n}^{3} $

然后就没有然后了,暴力组合数展开。。。qwq

#include<bits/stdc++.h>
#define LL long long
const int mod = 10007;
LL in() {
    char ch; LL x = 0, f = 1;
    while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
    for(x = ch ^ 48; isdigit(ch = getchar()); x = ((x << 1) + (x << 3) + (ch ^ 48)) % mod);
    return x * f;
}
int main() {
    LL n = in();
    printf("%lld\n", (n * (n + 1) * (n + 2) / 6) % mod);
    return 0;
}
posted @ 2019-01-16 07:20  olinr  阅读(171)  评论(0编辑  收藏  举报