NKOJ-4052 糖果盒

P4052糖果盒
时间限制 : - MS 空间限制 : 65536 KB
评测说明 : 1s
问题描述

何老板有很多颗相同的糖果，同时还有n个不同的糖果盒。何老板想把其中一些糖果放入这些盒子里，要求所有盒子里总的糖果数不超过m，问总共有多少种不同的方案？
每个盒子的容量无限大，可以选择装糖，也可以不装。
答案可能很大，mod p后在输出，p是一个素数。

输入格式

第一行，三个空格间隔的整数 n,m,p （1 <= n, m <= 100000000, 1 < p < 100000）

输出格式

一行，一个整数，表示计算结果

样例输入 1

2 1 5

样例输出 1

3

样例输入 2

3 8 7

样例输出 2

4
——————-说在前面———————
这是一道描述何老板如何闲的**装糖的题，表面上看是要你算组合数学，实际上，它是在告诉我们，一个人闲起来到底有多可怕
———————正文————————
按照老板的大板子+假装有球法
把m个糖分成n堆，而且还可以分到的数为0，其实就是相当于

在（m个实球和n个虚球（就是只有聪明人才看到的球））的m+n-1个间隙里插n-1个板子

例如：

    我有3个实球（x)，分3堆，我就加两个虚球(y)
    那么方案如下;
    3+2个球 插 3-1个板子
        x|x|xyy
        x|xx|yy
        x|xxy|y
        xx|x|yy
        xx|xy|y
        xxx|y|y
        =C(3+2-1,3-1)

结论：
m个球分为n堆，一堆的个数可以为0时，方案数为：

        C(m+n-1,n)     (肯定是对的，不要问我为什么)(因为这是老板说的呀哈哈哈)

但你以为这就完了么

稚嫩

我们最终要算的是0-m个球分到n个盒子里头的方案数
    C(n+0-1,0)+C(n+1-1,1)+(n+2-1,2)+...+C(n+m-1,m)
=   C(n-1,0)+C(n,1)+C(n+1,2)+...+C(n+m-1,n-1)
=   C(n,0)+C(n,1)+C(n+1,2)+...+C(n+m-1,n-1)
(根据C(n,m)=C(n-1,m)+C(n-1,m-1),就有:)
=   C(n+1,1)     +C(n+1,2)+...+C(n+m-1,n-1)
=   C(n+2,2)     +...+C(n+m-1,n-1)
...
=   C(n+m-1,n)+C(n+m-1,n-1)
=   C(n+m,n);

附上一段半抄半写的代码

#include <iostream>
#include <cstdio>
using namespace std;

long long p;

long long KSM(long long a,long long b,long long c)
{
    long long ans=1;
    a=a%c;
    while(b>0)
    {
        if(b&1)ans=(ans*a)%c;
        a=(a*a)%c;
        b=b/2;
    }
    return ans;
}
long long C(long long a,long long b)
{
    if(a<b)return 0;
    if(a==b)return 1;
    if(b>a-b)b=a-b;
    long long A=1,B=1;
    for(long long i=0;i<b;++i)
    {
        A=(A*(a-i))%p;
        B=(B*(b-i))%p;
    }
    long long add=KSM(B,p-2,p);
    return A*add%p;
}
long long Pro(long long n,long long m)
{
    if(m==0)return 1;
    return C(n%p,m%p)*Pro(n/p,m/p)%p;
}



int main()
{
    long long n,m;
    scanf("%lld%lld%lld",&n,&m,&p);
    long long ans=Pro(n+m,n);
    cout<<ans;
}