题解报告`排列组合 (组合 + 牡牛和牝牛 + 方程的解) 7/30

[1]组合
题目描述
给出组合数 C(n,m) 表示从 n 个元素中选出 m 个元素的方案数。例如 C(5,2)=10,C(4,2)=6。可是当 n,m 比较大的时候，C(n,m) 很大。于是 xiaobo 希望你输出C(n,m)modp 的值。
输入
输入数据第一行是一个正整数 T，表示数据组数；

接下来是 T 组数据，每组数据有 3 个正整数 n,m,p。

对于所有数据，T≤100，1≤m≤n≤109，m≤104，m<p<109，p 是素数。

输出
对于每组数据，输出一个正整数，表示 C(n,m)modp 的结果。
样例输入
2
5 2 3
5 2 61
样例输出
1
10
思路:
Lucas定理模板题,
因为n,m的值太大所以要简化运算, C(m, n) = C(m mod p, n mod p) * C(m / p, n / p) (mod p)

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#define ll long long
using namespace std;

ll fpow(ll a,ll b,ll n){
   
	if(b==1) return a;
	if(b%2 == 1){
   
		ll t = fpow(a,b/2,n);
		t = t * t % n;
		t = t * a % n;
		return t; 
	}
	else{
   
		ll t = fpow(a,b/2,n);
		t = t * t % n;
		return t;
	}
} 

ll comb(ll n,ll m,ll p){
   
	if(m>n)return 0;
	if(m==0)return 1;
	ll up = 1 ,down = 1;
	for(ll i = 1 ; i <= m ; i++){
   
		up = (up * (n-i+1)) % p;
		down = (down * i) % p ; 
	} 
	return up * fpow(down,p-2,p) % p