题意:n个建筑高度为1-n,从前能看到f个,从后能看到b个,求可能有多少张排序情况
解析:最高建筑n是一定可以看到的,固定住最高建筑n号楼的位置,将n号楼左边分为f-1组,右边b-1组,且用每组的最高元素代表这一组。而后发现n号楼的左边,组与组之间的元素一定是单调递增的!
且每组中最高元素一定排在该组的最左边,每组的其他元素任意排列(相当于这个组中所有元素的环排列,所谓环排列就是一个元素位置固定其他元素自己动,方案数等于(n-1)!)。
右边同理。
把n-1个数分为分f-1+b-1组,每组做环排列(第一类Stirling数),寻找f-1组放在n的左边
将上面的相乘即是结果;
#include<stdio.h>
#include<iostream>
using namespace std;
#define ll long long
const int maxn=2005;
const int mod=1e9+7;
long long s[maxn][maxn];
long long c[maxn][maxn];
void init()
{
for(int i=0;i<maxn;i++)
{
c[i][0]=1;
c[i][i]=1;
s[i][0]=0;
s[i][i]=1;
for(int j=1;j<i;j++)
{
c[i][j]=(c[i-1][j]%mod+c[i-1][j-1]%mod)%mod;
s[i][j]=(i-1)%mod*s[i-1][j]%mod+s[i-1][j-1]%mod;
}
}
}
int main()
{
int t;
cin>>t;
init();
while(t--)
{
ll n,f,b;
scanf("%lld%lld%lld",&n,&f,&b);
ll ans;
if(b+f-2<=2000)
ans=c[f+b-2][f-1]%mod*s[n-1][f+b-2]%mod;
else
ans=0;
printf("%lld\n",ans);
}
}
PS:杨辉三角居然是求组合数的结果!
第一类Stirling数 s(p,k)
s(p,k)的一个的组合学解释是:将p个物体排成k个非空循环排列的方法数。
s(p,k)的递推公式: s(p,k)=(p-1)*s(p-1,k)+s(p-1,k-1) ,1<=k<=p-1
边界条件:s(p,0)=0 ,p>=1 s(p,p)=1 ,p>=0
递推关系的说明:
考虑第p个物品,p可以单独构成一个非空循环排列,这样前p-1种物品构成k-1个非空循环排列,方法数为s(p-1,k-1);
也可以前p-1种物品构成k个非空循环排列,而第p个物品插入第i个物品的左边,这有(p-1)*s(p-1,k)种方法。
第二类Stirling数 S(p,k)
S(p,k)的一个组合学解释是:将p个物体划分成k个非空的不可辨别的(可以理解为盒子没有编号)集合的方法数。
k!S(p,k)是把p个人分进k间有差别(如:被标有房号)的房间(无空房)的方法数。
S(p,k)的递推公式是:S(p,k)=k*S(p-1,k)+S(p-1,k-1) ,1<= k<=p-1
边界条件:S(p,p)=1 ,p>=0 S(p,0)=0 ,p>=1
递推关系的说明:
考虑第p个物品,p可以单独构成一个非空集合,此时前p-1个物品构成k-1个非空的不可辨别的集合,方法数为S(p-1,k-1);
也可以前p-1种物品构成k个非空的不可辨别的集合,第p个物品放入任意一个中,这样有k*S(p-1,k)种方法。
第一类斯特林数和第二类斯特林数有相同的初始条件,但递推关系不同。
转自:http://blog.csdn.net/acdreamers/article/details/8521134
第一类斯特林数和第二类斯特林数的分组都不是排列,就是说组是没有编号的,如果加上编号,乘上k!即可,或者在递推公式里面乘。
比如: s(p,k)=(p-1)*s(p-1,k)+s(p-1,k-1)*k (表示如果有一个新的组,要考虑它的位置对答案影响)
第一类斯特林数的n个物品循环排列方案事实上等于n-1个物品的排列方案,因为n!/n=(n-1)!
所以如果你要分组,组内一个元素位置要固定,那么就等于循环排列。(如本题)
斯特林数长这样:
第一类斯特林数有以下性质: