【组合数学】 放小球之隔板法

题目描述

N个盒子排成一行(1<=N<=20)。你有A个红球和B个蓝球。0 <= A <= 15, 0 <= B <= 15。球除了颜色没有任何区别。你可以将球放进盒子。一个盒子可以同时放进两种球,也可以只放一种,也可以空着。球不必全部放入盒子中。编程计算有多少种放置球的方法。

输入

一行,N,A,B,用空格分开。

输出

一行,输出放置方案总数。

样例输入

 (如果复制到控制台无换行,可以先粘贴到文本编辑器,再复制)

2 1 1

样例输出

9

分析:

主要是想通过这道题整理一下组合数学的一些基本概念。
题目说球不必全部放入盒子中,我们先假定它要全部放进去,这时应该怎么算呢?
先考虑红球,A个红球放进N个盒子里,因为盒子是排成一行的,即盒子是有顺序的,所以我们可以设:
第一个盒子放x1个.....第n个盒子放xn个,即有:x1+x2+...+xn=A。
显然,我们需要把A拆成N份,但是每份可以为空,所以我们把每个x都加上一,即转化为:
(x1+1)+(x2+1)+....+(xn+1)=A+N。
这时候就简单了,把A+N分成N份,每份不能为空,这时候就用隔板法。
如图,在A+N-1个空中选N-1个放隔板:C(a+n-1,n-1).
但这是球必须放完的情况,那么想,如果球不放完,我们加一个盒子,把剩下的球放到第N+1个盒子中,是不是就是C(a+n,n)了呢?

至此,求出放A的方案和放B的方案,两者相乘,即为答案。

那么我又想bibi两句了,若何求C(n,k)呢?若是直接用n!/k!(k-n)!,肯定要炸,所以要边乘边除,
即:for(int i=1;i<=k;i++,n--) 
s=s*n/i;
为什么要这么写,而不是s*n/k呢,因为这样可能出现除不尽的情况,但是为何那样写就可以除尽呢?

因为连续的n个正整数的乘积一定能被n!整除

证明:设a为非负数,
则n个正整数的乘积可表示为(a+1)(a+2)....(a+n) = (a+n)! / a! = (a+n)!/(a!n!) * n! 
因为(a+n)!/a!n! = C(a+n,a)
所以(a+1)....(a+n)/n! = C(a+n,a),而C(a+n,a)一定为正整数,证毕。

题解至此,可以上代码了吧。

#include<iostream>
using namespace std;
#define LL unsigned long long
int n,a,b;
int main()
{
    cin>>n>>a>>b;
    LL s=1,s0=1;
    for(int i=1,j=n+a;i<=a;i++,j--)
        s=s*j/i;
    for(int i=1,j=n+b;i<=b;i++,j--)
        s0=s0*j/i;
    cout<<s*s0;
}



这道题可以使用递归来解决。假设现在需要将 $A$ 个红球和 $B$ 个蓝球进 $N$ 个盒子中,可以分为以下三种情况: 1. 当前盒子球,继续考虑下一个盒子; 2. 当前盒子红球蓝球,继续考虑下一个盒子,并将 $A$ 或 $B$ 的值减去 1; 3. 当前盒子红球蓝球,继续考虑下一个盒子,并将 $A$ 和 $B$ 的值均减去 1。 使用递归的方式,将问题规模逐渐缩小,直到所有球都被盒子中,此时方案数加 1。代码如下: ```c #include <stdio.h> int count = 0; void dfs(int n, int a, int b) { if (n == 0) { // 所有盒子都已考虑完毕,方案数加 1 count++; return; } // 当前盒子球 dfs(n - 1, a, b); // 当前盒子红球 if (a > 0) { dfs(n - 1, a - 1, b); } // 当前盒子蓝球 if (b > 0) { dfs(n - 1, a, b - 1); } // 当前盒子红球蓝球 if (a > 0 && b > 0) { dfs(n - 1, a - 1, b - 1); } } int main() { int n, a, b; scanf("%d%d%d", &n, &a, &b); dfs(n, a, b); printf("%d", count); return 0; } ``` 具体实现中,使用全局变量 `count` 记录方案数,使用函数 `dfs` 进行递归。在函数中,首先判断是否所有盒子都已考虑完毕,如果是,则方案数加 1 并返回;否则,分别考虑当前盒子球、只红球、只蓝球红球蓝球的四种情况,并继续递归考虑下一个盒子。最后,在主函数中读入输入数据,调用 `dfs` 函数计算方案数,并输出结果。 需要注意的是,当 $A$ 和 $B$ 的值较大时,递归深度可能会非常大,导致栈溢出。可以使用尾递归优化或改为迭代实现来解决这个问题。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值