fzu2020( c(n,m)%p,其中n, m, p (1 <= m <= n <= 10^9, m <= 10^4, m < p < 10^9, p是素数) )

基本的模板题,统计分子分母中p出现的次数,然后求逆元取模.

//
//  main.cpp
//  fzu2020
//
//  Created by 陈加寿 on 15/12/27.
//  Copyright (c) 2015年 chenhuan001. All rights reserved.
//

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <string>
#include <stdlib.h>
using namespace std;

typedef long long ll;

//ax + by = gcd(a,b)
//传入固定值a,b.放回 d=gcd(a,b), x , y
void extendgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(b==0){d=a;x=1;y=0;return;}
    extendgcd(b,a%b,d,y,x);
    y-=x*(a/b);
}

//Ax=1(mod M),gcd(A,M)==1
//输入:10^18>=A,M>=1
//输出:返回x的范围是[1,M-1]
ll GetNi(ll A,ll M)
{
    ll rex=0,rey=0;
    ll td=0;
    extendgcd(A,M,td,rex,rey);
    return (rex%M+M)%M;
}


int main(int argc, const char * argv[]) {
    int T;
    cin>>T;
    while(T--)
    {
        int n,m,p;
        scanf("%d%d%d",&n,&m,&p);
        int cnt=0;//统计有多少个因子p
        long long up=1;
        for(int i=0;i<m;i++)
        {
            int tn = n-i;
            if( tn%p==0 )
            {
                while( tn%p==0 )
                {
                    tn/=p;
                    cnt++;
                }
            }
            up = (up*tn)%p;
        }
        
        long long dn=1;
        for(int i=1;i<=m;i++)
        {
            int tm=i;
            if( tm%p==0 )
            {
                while(tm%p==0)
                {
                    tm /= p;
                    cnt--;
                }
            }
            dn = (dn*tm)%p;
        }
        if(cnt != 0)
        {
            printf("0\n");
            continue;
        }
        
        cout << GetNi(dn, p)*up%p<<endl;
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值