Chess
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 235 Accepted Submission(s): 170
Problem Description
車是中国象棋中的一种棋子,它能攻击同一行或同一列中没有其他棋子阻隔的棋子。一天,小度在棋盘上摆起了许多車……他想知道,在一共N×M个点的矩形棋盘中摆最多个数的車使其互不攻击的方案数。他经过思考,得出了答案。但他仍不满足,想增加一个条件:对于任何一个車A,如果有其他一个車B在它的上方(車B行号小于車A),那么車A必须在車B的右边(車A列号大于車B)。
现在要问问你,满足要求的方案数是多少。
现在要问问你,满足要求的方案数是多少。
Input
第一行一个正整数T,表示数据组数。
对于每组数据:一行,两个正整数N和M(N<=1000,M<=1000)。
对于每组数据:一行,两个正整数N和M(N<=1000,M<=1000)。
Output
对于每组数据输出一行,代表方案数模1000000007(1e9+7)。
Sample Input
1 1 1
Sample Output
1
Source
类似N皇后问题,n*m棋局,最多放,那么就是放n和m中小的那个数,而且是在大的中取小的数量,这样就转成了组合数问题了,数据量较大,用快速幂取模和组合数大数取模计算,这里引用卢卡斯函数,模板带上:
#include <cstdio>
#include <iostream>
#include <algorithm>
#define LL long long
const int mod = 1000000007;
const int maxn = 1010;
using namespace std;
LL Fast_power(LL n,LL k,LL mod)
{
if(k == 0) return 1;
LL temp = Fast_power(n,k / 2,mod);
temp = (temp * temp) % mod;
if(k % 2 == 1) temp = (temp * (n % mod)) % mod;
return temp;
}
LL Comb(LL n,LL m,LL p)
{
if(n < m) return 0;
if(n == m) return 1;
m = min(n - m,m);
LL ans = 1,cn = 1,cm = 1;
for(LL i = 0; i < m; ++ i)
{
cn = (cn * (n - i)) % p;
cm = (cm * (m - i)) % p;
}
ans = (cn * Fast_power(cm,p - 2,p)) % p;
return ans;
}
LL Lucas(LL n,LL m,LL p)
{
LL ans = 1;
while(n && m && ans)
{
ans = (ans * Comb(n % p,m % p,p)) % p;
n = n / p;
m = m / p;
}
return ans;
}
int main()
{
int T;
scanf("%d",&T);
while (T--)
{
int n,m;
scanf("%d%d",&n,&m);
int mm=max(n,m);
int nn=min(n,m);
printf("%I64d\n",Lucas(mm,nn,mod));
}
return 0;
}