2016"百度之星" - 初赛（Astar Round2B）1003 瞬间移动(HDU5698)(逆元+快速幂）

2016”百度之星” - 初赛（Astar Round2B）1003 瞬间移动(HDU5698)(逆元+快速幂）

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

Source
2016”百度之星” - 初赛（Astar Round2B）

题目在这里：题目链接

---

假设到第(n,m)个格子的方法数为s(n,m)

则s(n,m)等于他与起始位置所夹的矩形内的各个格子方法数之和，即 s(n,m)=$\sum_{i=2}^{n- 1}\sum_{j=2}^{m-1} s(i,j)$

不难发现s(n,m)等于s(n-1,m)+s(n,m-1)，即答案是一个斜着的杨辉三角

问题便转换为组合数求模，费马小定理，快速幂求逆元，

预处理阶乘和逆元

---

代码：

#include<stdio.h>
#define max 200005
#define mod 1000000007
int n,m;
long long j[max];           
long long p[max];           

long long ksm(long long a,long long b)
{
    long long ans,base;
    ans=1;
    base=a%mod;
    while(b)
    {
        if(b&1)
            ans=ans*base%mod;
        base=base*base%mod;
        b>>=1;
    }
    return ans;
}

void init()
{
    j[1]=1;
    p[1]=1;
    for(int i=2;i<=max;i++)
    {
        j[i]=(i*j[i-1])%mod;
        p[i]=ksm(j[i],mod-2);
    }
}

int main()
{
    init();
    while(~scanf("%d %d",&n,&m))
    {
        if(n==2||m==2)
        {
            printf("1\n");
            continue;
        }
        printf("%lld\n",j[n+m-4]%mod*p[m-2]%mod*p[n-2]%mod);
    }
    return 0;
}