HDU 6114 Chess【Lucas定理】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6114

Chess

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 542 Accepted Submission(s): 333

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个数里面挑m个数有多少种方案

针对这一题的数据量来看有两种方法：
1.组合数递推
组合数递推公式：c(n,m)=c(n-1,m-1)+c(n-1,m)

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL mod=1e9+7;
LL c[1005][800];//数组大小我以为开到1005*505就够了，可是答案竟然错了。开到800就对了
void math()
{
    c[0][0]=1;
    for(int i=1;i<1004;i++)
    {
        c[i][0]=c[i][i]=1;
        for(int j=1;j<i;j++)
        c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }
}
int main()
{
    math();
    int n,m,t;
    scanf("%d",&t);
    while(t--)
    {
        int tt;
        scanf("%d%d",&n,&m);
        if(n<m)
        {
            tt=n;
            n=m;
            m=tt;
        }
        tt=n-m;
        printf("%lld\n",c[n][tt]);
    }
    return 0;
}

2.卢卡斯定理
Lucas定理是用来求 c(n,m) mod p，p为素数的值
题目中所给的1e9+7就是一个素数

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long LL;
const LL mod=1e9+7;
typedef long long LL;
LL n,m;
LL quickPow(LL a,LL k)
{//快速幂求a^k%mod
    LL ans=1;
    while(k)
    {
        if(k&1)//相与运算,如果k是奇数那么结果为真,执行下面操作
        ans=(ans*a)%mod;
        a=(a*a)%mod;
        k>>=1;
    }
    return ans;
}
LL C(LL a,LL b,LL p)
{//求从n个数中挑出m个数的组合数
    if(a<b)
      return 0;
    if(b>a-b)
      b=a-b;
    LL up=1,down=1;
    for(LL i=0;i<b;i++)
    {
        up=up*(a-i)%p;
        down=down*(i+1)%p;
    }
    return up*quickPow(down,p-2)%p;//乘法逆元(费马小定理)
}
LL lucas(LL a,LL b,LL p)
{
    if(b==0)
        return 1;
    return C(a%p,b%p,p)*lucas(a/p,b/p,p)%p;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        LL x;
        scanf("%lld%lld",&n,&m);
        if(n<m)
        {
            x=n;
            n=m;
            m=x;
        }
        printf("%lld\n",lucas(n,m,mod));
    }
    return 0;
}