HDU 5698 瞬间移动 [数论] [逆元] [组合数取模]

最新推荐文章于 2021-11-13 22:32:49 发布

JacquesdeH

最新推荐文章于 2021-11-13 22:32:49 发布

阅读量414

点赞数

分类专栏：数论逆元费马小定理组合数文章标签：数论费马小定理逆元组合数

本文链接：https://blog.csdn.net/ourfutr2330/article/details/52489925

版权

数论同时被 3 个专栏收录

20 篇文章 0 订阅

订阅专栏

组合数

4 篇文章 0 订阅

订阅专栏

逆元

2 篇文章 0 订阅

订阅专栏

该博客讨论了HDU 5698瞬间移动问题，这是一个涉及数论和组合数学的算法题。作者指出，问题可以转化为求解组合数C(m+n-4, m-2)对1000000007取模的值。提出了三种解决方案，包括预处理阶乘、预处理逆元以及直接求阶乘和递归求逆元，并给出了采用预处理阶乘方法的实现思路，强调阶乘预处理需要到2N的规模以应对m+n-4的规模。" 112910060,8473821,国科大高级人工智能2018期末试卷解析,"['人工智能', '机器学习', '深度学习', '自然语言处理', '神经网络']

摘要由CSDN通过智能技术生成

瞬间移动
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 893 Accepted Submission(s): 476

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

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

Output
一个整数表示答案

Sample Input
4 5

Sample Output
10

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

首先找出是一个杨辉三角。。。
那么很容易得到 ans = C (m+n-4,m-2)
然后的话组合数取模，但是mod太大，这时候有三种方法：

预处理阶乘，然后直接定义法求解
预处理N以内的逆元，然后用组合数线性递推公式边走边求
直接求出阶乘，然后用递归求逆元，此时的数可能超过N

下面给出第一种方法，注意阶乘要预处理到2N，因为是m+n-4数量级的。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
typedef long long LL;
const int INF=0x3f3f3f3f;
const int mod = 1000000007;
const int maxn = 200005;
const int N = 200000;
int fac[maxn];
int quick_exp(int a,int p)
{
    if(!p) return 1;
    if(p==1) return a%mod;
    int tmp=quick_exp(a,p>>1);
    tmp=(LL)tmp*tmp%mod;
    if(p&1) return (LL)tmp*a%mod;
    else return tmp;
}
int C(int m,int n)
{
    int tmp=(LL)fac[m]*quick_exp(fac[n],mod-2)%mod;
    tmp=(LL)tmp*quick_exp(fac[m-n],mod-2)%mod;
    return tmp;
}
int n,m;
void init()
{
    fac[0]=1;
    for(int i=1;i<=N;i++) fac[i]=(LL)fac[i-1]*i%mod;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("rapid.in","r",stdin);
    freopen("rapid.out","w",stdout);
#endif
    init();
    while(~scanf("%d%d",&n,&m))
    {
        if(n==1||m==1)
        {
            if(n==1&&m==1) printf("1\n");
            else printf("0\n");
            continue;
        }
        int ans = C(m+n-4,m-2);
        printf("%d\n",ans);
    }
    return 0;
}