食物分配，数学推导题

 

新版食物（shiwu）

明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！

我们暂且不讨论他有多么NC，他又幻想了他应该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。

他这次又准备带一些受欢迎的食物，如：蜜桃多啦，鸡块啦，承德汉堡等等

当然，他又有一些稀奇古怪的限制：

每种食物的限制如下：

       承德汉堡：偶数个

       可乐：0个或1个

            鸡腿：0个，1个或2个

            蜜桃多：奇数个

            鸡块：4的倍数个

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

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

            面包：3的倍数个

 

 

 

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

 

输入样例1

  1

输出样例1

  1

 

输入样例2

  5

输出样例2

  35

 数据范围

   对于40%的数据，1<=N<=100000;

   对于所有数据，1<=n<=10^500;

 

 

 

我们试着从母函数角度解释（聪明的同学应该还会发现还可以用逐层作差法解决，即首先找小数据的规律）。

首先，

汉堡是2的倍数。

H(X)=1+X^2+X^4+X^6+……=1/(1-X^2);

可乐0个或1个。

G(X)=1+X=(1-X^2)/(1-X);

依此类推

鸡腿：x/(1-x^3)

蜜桃多：......

鸡块：......

包子：......

土豆片炒肉：......

面包：......

利用卷积规则：最后方案数母函数表示为：1/(X^4+1)

而它的第项系数刚好为C3 N+2

这个，同学们不太理解的可以参考毛杰明的论文，只需要前面一点点知识就行了。

 

 

最后记得加高精

 

 

program s; var i,j,na,nb,nc,ans,t:longint; st:ansistring; a,b,c:array[0..10000]of longint; procedure cheng; var i,j:longint; begin for i:=1 to na do for j:=1 to nb do c[i+j-1]:=a[i]*b[j]+c[i+j-1]; nc:=na+nb; for i:=1 to nc do if c[i]>=10 then begin c[i+1]:=c[i+1]+c[i] div 10; c[i]:=c[i] mod 10; end; if c[nc]=0 then dec(nc); end; begin assign(input,'s.in'); reset(input); assign(output,'s.out'); rewrite(output); readln(st); for i:=1 to length(st) do a[length(st)-i+1]:=ord(st[i])-48; b:=a; na:=length(st); nb:=na; if na=0 then begin writeln(0);close(input);close(output);exit;end; inc(b[1]); cheng; na:=nc;a:=c; inc(b[1]); fillchar(c,sizeof(c),0); cheng; for i:=nc downto 1 do begin c[i-1]:=10*(c[i]mod 6)+c[i-1]; c[i]:=c[i] div 6; end; if c[nc]=0 then dec(nc); i:=nc; while i>5 do begin t:=100000*c[i]+10000*c[i-1]+1000*c[i-2]+100*c[i-3]+10*c[i-4]+c[i-5]; c[i]:=0;c[i-1]:=0;c[i-2]:=0;c[i-3]:=0;c[i-4]:=0;c[i-5]:=0; t:=t mod 10007; i:=i-5; c[i]:=t; while c[i]>=10 do begin c[i+1]:=c[i] div 10; c[i]:=c[i] mod 10; inc(i); end; end; ans:=10000*c[5]+1000*c[4]+100*c[3]+10*c[2]+c[1]; write(ans mod 10007); close(input); close(output); end.