hdu 5794

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

vvv_557

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

阅读量230

点赞数

分类专栏：数学_lucas

本文链接：https://blog.csdn.net/vvv_557/article/details/75046509

版权

数学_lucas 专栏收录该内容

0 篇文章 0 订阅

订阅专栏

/*
  打表发现路径数与坐标 表现成杨辉三角形关系 
  无障碍下，路径数等于相对坐标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;
}