Chess

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）。

 

Output

对于每组数据输出一行，代表方案数模1000000007（1e9+7）。

 

Sample Input

  
  

   
   1
1 1
  
  

 

Sample Output

  
  

   
   1
  
  

 

Source

2017"百度之星"程序设计大赛 - 初赛（B）

类似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;
}