bzoj 3028: 食物（生成函数）

最新推荐文章于 2019-03-08 21:38:10 发布

clover_hxy

最新推荐文章于 2019-03-08 21:38:10 发布

阅读量538

点赞数

分类专栏：数论生成函数文章标签：生成函数

本文链接：https://blog.csdn.net/clover_hxy/article/details/70213981

版权

数论同时被 2 个专栏收录

158 篇文章 0 订阅

订阅专栏

生成函数

6 篇文章 0 订阅

订阅专栏

题目描述

传送门

题解

因为没有顺序，所以这是普通的生成函数。

1 ( 1 - a x ) m = 1 + a C m - 1 m x + a 2 C m - 1 m + 1 x 2 + a 3 C m - 1 m + 2 x 3 + . . . . .

${1\over (1-ax)^m}=1+aC_{m}^{m-1}x+a^2C_{m+1}^{m-1}x^2+a^3C_{m+2}^{m-1}x^3+.....$
然后将这题的限制条件全部转换成生成函数中的未知项的指数，可以得到式子

(1 + x 2 + x 4 + x 6 . . . .) (1 + x) (1 + x + x 2) (x 1 + x 3 + x 5 . . . .) (1 + x 4 + x 8 . . . . .) (1 + x + x 2 + x 3) (1 + x) (1 + x 3 + x 6 . . . . .)

$(1+x^2+x^4+x^6....)(1+x)(1+x+x^2)(x^1+x^3+x^5....)(1+x^4+x^8.....)(1+x+x^2+x^3)(1+x)(1+x^3+x^6.....)$
化简约分最后得到

x(1−x)4 $x\over (1-x)^4$ ，相当于

x+C34x2+C34+1x3+C34+2x4+..... $x+C_{4}^{3}x^2+C_{4+1}^{3}x^3+C_{4+2}^{3}x^4+.....$
因为答案是x^n的系数，所以

ans=C3n+2 $ans=C_{n+2}^3$

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define p 10007
#define N 1003
using namespace std;
char s[N]; 
int n,jc[100000],inv[100000];
int calc(int n,int m)
{
    return jc[n]*inv[n-m]%p*inv[m]%p;
}
int main()
{
    freopen("a.in","r",stdin);
    scanf("%s",s+1); int len=strlen(s+1);
    jc[0]=1; inv[0]=inv[1]=1;
    for (int i=1;i<=p;i++) jc[i]=jc[i-1]*i%p;
    for (int i=2;i<=p;i++) inv[i]=(p-p/i)*inv[p%i]%p;
    for (int i=2;i<=p;i++) inv[i]=inv[i]*inv[i-1]%p;
    for (int i=1;i<=len;i++)
     n=n*10+s[i]-'0',n%=p;
    //cout<<n<<endl;
    printf("%d\n",calc(n+2,3));
}