2021杭电多校第7场1004Link with Balls

题目链接：Problem - 7047 (hdu.edu.cn)

题意

给你一个整数n，有2n个桶，编号1~2n；

编号为奇数（2x-1，1,3,5...）的桶只能取k*x个球，k为非负整数；

编号为偶数（2x，2,4,6...）的桶最多取x个球，即可取0~x个；

求，从这些框中取m个球的方案数？（题目要求，答案 mod 1e9+7）。

思路和做法

这一题，当然啦，肯定是要组合数学来做的，就是使用经典的插板法来求方案数。至于为什么？思路如下：

首先有几个要注意的点：

1. 题目中的k为非负整数，且不同编号（奇数）的桶所取k*x个球中的k不一定是相同的；
2. 编号为偶数的桶，可以取的球数也不是固定的，是0~x。

这里有一个巧妙的思路就是：

将可以取0~k-1个球的桶（也就是编号为偶数的桶）与只能取k的倍数个球的桶合并为一个可以取任意个球的桶，这样就得到了n个可以取任意个球的桶（假设为A桶）和一个可以取0~n个球的桶（假设为B桶）。看到这里，估计有人开始迷了，为什么可以得到上面这样的桶呢？不急，先举个栗子：

例：

n=3；

              桶编号：  1        2        3        4        5        6

可以取球的个数：  k      0~1     2k     0~2      3k    0~3

上面思路中，我们需要得到n个A桶，1个B桶

从上面的结果可以得到B桶一个，6号桶，可以取0~3个球；

1号桶是一个A桶，k是非负整数，所以可以取任意个球；

接下来看3号和5号桶：

3号可以取球的个数为，2，4，6...，而2号桶可以取0~1个球，将这两个桶合并就可以得到A桶；

5号可以取球的个数为，3，6，9...，而4号桶可以取0~2个球，将这两个桶合并就可以得到A桶。

从以上分析下来，我们就得到了n个A桶，1个B桶。

得到这些桶后，我们就可以枚举B桶中取0~n个球，剩下的球的选取方案可以使用插板法得到，于是答案为    ，即  。

至于为什么就这么得到答案了，这里的插板法不是直接的插板法，需要转换一下。下面说说是如何得到答案的。

首先我们剩下要取球的个数是m-i，然后我们假设我们得到的n个A桶取球的个数分别为A1，A2，...，An，从而得到一个等式，A1 + A2 + ... + An = m-i，A1~An >= 0，而经典的插板法，是需要每一项都是>0的，即>=1。

这里我们将每一项都加上1，同时等式右边也要加上等值的n，这样就得到了A1 + A2 + ... + An = m-i+n，这时的等式已经满足条件了，那么接下来就是插板法，在m-i+n个数中插n-1个板子，刚好有m-i+n-1个空，所以答案就是 。

哦，对了，忘记解释这个式子是怎么得来的了，我们将求和公式展开，可以得到，，每当我们看到这种上面的数相同的组合数之和时，我们应该要想到合并，根据公式，我们先加上一个和合并为后，继续合并，最后记得把减回来，即可合并成那个式子。

本人第一篇博客，有不足之处，还请谅解以及指出。

AC代码

/*
1004
标签：组合数学，插板法 
*/ 
#include<bits/stdc++.h> 
using namespace std;
const int mod = 1e9+7;
const int N = 2e6+5;
int fac[N],inv[N];//fac为阶乘，inv为逆元
int qmi(int a,int b,int p)//快速幂
{
	int ans = 1;
	while(b)
	{
		if(b&1)
			ans = 1ll*ans*a%mod;
		a = 1ll*a*a%mod;
		b>>=1;
	}
	return ans;
}
void init()//预处理阶乘和逆元
{
	fac[0]=inv[0]=1;
	for (int i=1;i<N;i++)
	{
		fac[i] = 1ll*fac[i-1]*i%mod;
		inv[i] = 1ll*inv[i-1]*qmi(i,mod-2,mod)%mod;
	}
}
int C(int n,int m)//根据公式，求解组合数。
{
	if(n<m)return 0;
	return 1ll*fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main ()
{
	init();
	int t;
	scanf("%d",&t);
	while(t--)
	{
		int n,m;
		scanf("%d%d",&n,&m);
		int ans = C(m+n,n)-C(m-1,n);
		ans = (ans%mod+mod)%mod;
		printf("%d\n",ans);
	}
	return 0;
}