逆元求组合数

long long pow_mod(long long x, long long n, long long mod)
{
    long long res = 1;
    while (n) {
        if (n & 1) res = res * x % mod;
        x = x * x % mod;
        n >>= 1;
    }
    return res;
}
long long fac[Max];
long long n, m, p;
int main() {
    while (~scanf("%lld %lld", &n, &m)) {
        p=1e9+7;
        fac[0] = 1;
        for (int i = 1; i <= n; i++) {fac[i] = fac[i - 1] * i % p;}
        //组合数 = n!*(m!%p的逆元)*((n-m)!%p的逆元)%p
        printf("%lld\n", fac[n] * pow_mod(fac[m], p - 2, p) % p * pow_mod(fac[n - m], p - 2, p) % p);
    }
}

原文地址:https://www.cnblogs.com/ww123/p/9402302.html

时间: 2024-10-17 08:39:04

逆元求组合数的相关文章

知识点:逆元求组合数取模

目标:求出C(n,m)%p 这里p是一个素数! 方法:费马小定理求逆元 因为膜的性质并不对除法适用,比如(a/b)%c; 但是,当我们知道了b%c的逆元d时,问题可以转化为:(a*d)%c=((a%c)*(b%c))%c; 考虑费马小定理: a^p-1=1(mod p) 显然有: a*a^p-2=1(mod p) 那么a^p-2就是a膜p意义下的逆元 利用快速幂即可求出! 然后我们需要预处理出m!,(n-m)!,分别对他们求出逆元,再与n!乘并取模即可! 原文地址:https://www.cnb

求组合数小结

今天学了一天数学,觉得自己都要转竞了23333 题目链接https://vjudge.net/contest/282927#problem/E 这里说一说求组合数的方法吧 其实就是求阶乘及其逆元的方法: 规定mod为模数,n为数据规模 1.mod为素数 费马小定理:nlogn 线性求逆元(n较小) 2.mod不为素数 欧拉筛出两个阶乘的素数,再计算出每个素数的次幂,最后快速幂乘起来即可 详见 n = read(),p = read(); ans = 1; for(int i = 2;i <= 2

求组合数C(m,n)的多种计算方法

https://ac.nowcoder.com/discuss/187813?type=101&order=0&pos=1&page=0 https://blog.csdn.net/shadandeajian/article/details/82084087 1.简单法---适合n,m很小 #include<bits/stdc++.h> using namespace std; const int MAXN = 1000; int C[MAXN+1][MAXN+1];

动态规划之求组合数

时间:2014.05.29 地点:基地 -------------------------------------------------------- 一.关于动态规划 如果问题是由交叠的子问题构成,则可用动态规划的方式求解.我们在将一个大问题划分为子问题的的过程中,如果递推关系中包含的子问题和大问题具有相同的形式,但由于子问题的交叠性质,我们用递归解决的代价往往很大,这时可考虑动态规划,即对每个交叠的子问题只求解一次,并把结果存储在记录表中,最后得出原始问题的解.直观上,好像动态规划是采取空

求组合数

组合数的计算虽说简单但也不乏有些陷阱,这主要是因为语言中的数据类型在表示范围上是有限的.更何况还有中间结果溢出的现象,所以千万要小心. 输入 求组合数的数据都是成对(M与N)出现的,每对整数M和N满足0<m, n≤20,以EOF结束. 输出 输出该组合数.每个组合数换行. 样例输入 5 2 18 13 样例输出 10 8568 代码 #include<stdio.h> int main(){int isum=1;int m,n,k;while(scanf("%d%d"

POJ 2992 Divisors 求组合数因子个数

题目来源:POJ 2992 Divisors 题意:... 思路:素数分解的唯一性 一个数可以被分解成若干素数相乘 p1^x1*p2^x2*...*pn^xn 根据乘法原理 因子数为 (x1+1)*(x2+1)*...*(xn+1) 不能直接求出组合数 会溢出 也不能把每个乘的数分解因子 这样会超时 C(N,M)=N!/(M!*(N-M)!) 另dp[i][j] 代表为i的阶乘中j因子的个数(j是素数) 那么i素数的个数为dp[n][i]-dp[m][i]-dp[n-m][i] 最后for循环从

[2011山东ACM省赛] Binomial Coeffcients(求组合数)

Binomial Coeffcients nid=24#time" style="padding-bottom:0px; margin:0px; padding-left:0px; padding-right:0px; color:rgb(83,113,197); text-decoration:none; padding-top:0px"> Time Limit: 1000ms   Memory limit: 65536K  有疑问?点这里^_^ 题目描写叙述 输入

hdu 2519 求组合数

求组合数 如果求C5 3 就是5*4*3/3*2*1 也就是(5/3)*(4/2)*(3/1) Sample Input5 //T3 2 //C3 25 34 43 68 0 Sample Output310101 1 # include <iostream> 2 # include <cstdio> 3 # include <cstring> 4 # include <algorithm> 5 # include <cmath> 6 # def

hdu 4927 java求组合数(大数)

import java.util.Scanner; import java.math.BigInteger; public class Main { private static int [] a = new int[3005]; private static int T; private static int n; public static void solve() { BigInteger rs = BigInteger.ZERO; BigInteger temp = BigInteger