您的位置:首页 > 其它

洛谷.3807.[模板]卢卡斯定理(Lucas)

2018-02-15 21:34 519 查看
题目链接

Lucas定理

日常水题...sublime和C++字体死活不同步怎么办...

//想错int范围了...不要被longlong坑
//这个范围现算阶乘比预处理快得多
#include <cstdio>
typedef long long LL;
const int N=1e5+5;

LL n,m,p;//,fac[N+3];

LL FP(LL x,LL k,LL p)
{
LL t=1;
for(; k; k>>=1,x=x*x%p)
if(k&1) t=t*x%p;
return t;
}
inline LL inv(LL x,LL p)
{
return FP(x,p-2,p);
}
//LL C(LL n,LL m)
//{
//  if(n<m) return 0ll;
//  return fac
*inv(fac[m],p)%p*inv(fac[n-m],p)%p;
//}
LL C(LL n,LL m)
{
if(n<m) return 0ll;
LL up=1ll,down=1ll;
for(LL i=n-m+1; i<=n; ++i) (up*=i)%=p;
for(LL i=2; i<=m; ++i) (down*=i)%=p;
return up*inv(down,p)%p;
}
LL Lucas(LL n,LL m,LL p)
{
LL ans=1;
for(; m && ans; n/=p, m/=p)
(ans*=C(n%p,m%p))%=p;
return ans;
}

int main()
{
//  fac[0]=fac[1]=1;
LL t; scanf("%lld",&t);
while(t--)
{
scanf("%lld%lld%lld",&n,&m,&p);
//      for(LL i=2; i<=p; ++i) fac[i]=i*fac[i-1]%p;
printf("%lld\n",Lucas(n+m,m,p));
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: