题目大意:n*n方格放m个车,求m个车不会 互相攻击(在同一行或同一列)的方案数。
当m>n时,直接输出0然后。。没有然后了。
原因很显然:只有n行,每行最多放一个,但有超过n个,显然不可能(或者说列也一样的)
否则,分两步:
n行里面选m行放车,方案数C(n,m)(从n个里选出m个)
n列里面选m列放车,方案数A(n,m)(因为要分到m列)=C(n,m)* m!
所以结果就为:
C(n,m) * A(n,m) = C(n,m)^2*m!
于是我们可以:
先预处理出m!
预处理出C(n,m)(=C(n-1,m)+C(n-1,m-1)),同时计算出C(n,m)*A(n,m) = C(n,m)^2*m!
接着就是O(1)的询问啦!
代码比较丑:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Up(a,b,c) for(ll a(b),end##a(c);a<=end##a;++a)
#define Down(a,c,b) for(ll a(b),end##a(c);a>=end##a;--a)
#define sUp(a,c,b) for(ll a(b),end##a(c);a<=end##a;++a)
#define sDown(a,b,c) for(ll a(b),end##a(c);a>=end##a;--a)
typedef long long ll;
using namespace std;
ll c[33][33],f[33][33],jc[33];
int main(){
jc[0]=1;
Up(i,1,30)jc[i]=jc[i-1]*i;
c[0][0]=f[0][0]=1;
Up(i,1,30){
c[i][0]=c[i][i]=f[i][0]=1;
f[i][i]=jc[i];
Up(j,1,i-1){
c[i][j]=c[i-1][j]+c[i-1][j-1];
f[i][j]=c[i][j]*c[i][j]*jc[j];
}
}
ll T,n,m;
scanf("%lld",&T);
Up(t,1,T){
scanf("%lld%lld",&n,&m);
printf("Case %lld: %lld\n",t,m<=n?f[n][m]:0);
}
return 0;
}