hdu 5794

最新推荐文章于 2021-10-13 20:56:05 发布

原创最新推荐文章于 2021-10-13 20:56:05 发布 · 298 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数学_lucas 专栏收录该内容

0 篇文章

订阅专栏

本文介绍了一种利用杨辉三角形关系计算无障碍路径数的方法，并通过Lucas定理处理大整数运算，避免溢出。对于有障碍的情况，提出通过减少障碍对路径的影响来求解。

/*
  打表发现路径数与坐标 表现成杨辉三角形关系 
  无障碍下，路径数等于相对坐标C(n,m);
  有障碍下，可以减去障碍坐标对终点，其他障碍的影响。
  计算大整数取mod用lucas定理;
*/
//非标程
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
typedef long long ll;
const double eps=1e-8,pi=acos(-1.0);
const int mod=110119;
ll A[110121];
ll quickm(int a,int n)
{
    ll ans=1;
    ll x=a;
    while(n)
    {
        if(n&1) ans =(ans*x)%mod;
        x =(x*x)%mod;
        n /=2;
    }
    return ans;
}
ll lucas(ll n,ll m)
{
    ll ans=1;
    while(n&&m)
    {
        int n1=n%mod,m1=m%mod;
        if(n1<m1) return 0;
        ans= (ans*A[n1])%mod*quickm(A[m1]*A[n1-m1]%mod,mod-2)%mod;
        cout<<ans<<' '<< A[n1] <<' '<<quickm(A[m1]*A[n1-m1]%mod,mod-2)<<endl;
        n /=mod;m /=mod;
    }
    return ans;
}
void init()
{
    A[0]=1;
    for(int i=1;i<=110120;i++) A[i] =(A[i-1]*i)%mod;
}
int main()
{
  init();
  ll n,m;
  cin>>n>>m;
  cout<<lucas(n,m)<<endl;
  return 0;
}