【HDU5698】瞬间移动（杨辉三角+逆元）

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=5698

题目：

---

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)

Problem Description

有一个无限大的矩形，初始时你在左上角（即第一行第一列），每次你都可以选择一个右下方格子，并瞬移过去（如从下图中的红色格子能直接瞬移到蓝色格子），求到第n行第m列的格子有几种方案，答案对1000000007取模。

Input

多组测试数据。

两个整数n,m(2≤n,m≤100000)

Output

一个整数表示答案

Sample Input

4  5

Sample Output

10

 

解题思路：

---

手动列一些项得到如下内容，可以看出这是一个斜置的杨辉三角，且远同一行的杨辉三角数在这里都位于一条斜线上，且n+m是个定值。

根据杨辉三角规律：第n行第m个数=，所以最终可以推出当n≥2,m≥2时，到大第n行第m列位置的方法数=

先把阶乘预处理出来，因为预处理的时候要取模，所以最终的分子未必能整除分母（也可能比分母小），如果再取模的话就可以转换成求逆元

 

ac代码：

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 200005;
const ll mod = 1000000007;
ll c[maxn];
ll qkpow(ll a, ll p, ll mod)//a^p
{
    ll t = 1, tt = a % mod;
    while(p)
    {
        if(p&1) t = (t*tt) % mod;
        tt = (tt*tt) % mod;
        p >>= 1;
    }
    return t;
}
ll getInv(ll a, ll mod)
{
    return qkpow(a, mod-2, mod);
}
int main()
{
    //freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    int n, m;
    c[0] = 1;
    for(int i = 1; i < maxn; i++)
        c[i] = (c[i-1]*i) % mod;
    while(~scanf("%d", &n))
    {
        scanf("%d", &m);
        ll N = n + m - 4, M = n - 2;
        ll ans = ((c[N] * getInv(c[M], mod)) % mod * getInv(c[N-M], mod)) % mod;
        cout << ans << endl;
    }
    return 0;
}