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;
Output
Sample Input
Sample Output
HINT
Source
母函数好题 答案是所有母函数积x^n的系数
根据公式∑(i=0 to inf)x^i=1/(1-x) 和∑(i=0 to n)x^i=(1-x^(n+1))/(1-x)
故所有函数的母函数
1:f(x)=∑(i=0 to inf)x^(2i)=1/(1-x^2)
2:f(x)=1+x=(1-x^2)/(1-x)
3:f(x)=1+x+x^2=(1-x^3)/(1-x)
4:f(x)=∑(i=0 to inf)x^(2i+1)=x/(1-x^2)
5:f(x)=∑(i=0 to inf)x^(4i)=1/(1-x^4)
6:f(x)=1+x+x^2+x^3=(1-x^4)/(1-x)
7:f(x)=1+x=(1-x^2)/(1-x)
8:f(x)=∑(i=0 to inf)x^(3i)=x/(1-x^3)
全部乘起来 约分得f(x)=x/((1-x)^4)=x*(1/(1-x))^4=x*(∑(i=0 to inf)x^i)^4
我们考虑(∑(i=0 to inf)x^i)^2=∑(i=0 to inf)(i+1)x^i(按次幂算贡献很容易得到) ①
把上式平方 考虑多项式乘法
A*B=C
记F(A,i)为多项式A中x^i的系数
有F(C,x)=∑(i=0 to x)(F(A,i)*F(B,x-i))
计算①的平方的值 记为C 要求的就是F(C,n-1)(因为外层还有一个x)
F(C,n-1)=∑(i=0 to n-1)(F(A,i)*F(A,n-1-i))
由于F(A,i)=i+1 故答案为∑(i=0 to n-1)(i+1)*(n-i)=∑(i=0 to n-1)(i*n+n-i-i^2)
分别计算 得答案为C(n+2,3)
#include<bits/stdc++.h>
using namespace std;
const int mod=10007;
int n;
char s[502];
int main()
{
scanf("%s",s+1);
for(int i=1;i<=strlen(s+1);i++) n=(n*10+s[i]-'0')%mod;
printf("%d\n",n*(n+1)%mod*(n+2)%mod*1668%mod);
}