组合数专题

最新推荐文章于 2022-11-12 11:10:12 发布
最新推荐文章于 2022-11-12 11:10:12 发布
阅读量1.1k 收藏
点赞数
分类专栏： shuoj 快速幂 组合数 文章标签： 快速幂 shuoj 组合数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sinat_30062549/article/details/47301767
版权
shuoj 同时被 3 个专栏收录
11 篇文章 0 订阅
订阅专栏
快速幂
4 篇文章 0 订阅
订阅专栏
组合数
2 篇文章 0 订阅
订阅专栏

 

shuoj1937-组合数level 0

 

Description
求C(n,m)对1000000007取模后的值，已知公式如下

 

Input
第一行有一个整数T，表示有T组数据(T≤1000)

接下来有T行，每行有两个整数n，m，用空格隔开。

0≤n≤20,0≤m≤n

Output
对于每组数据，输出答案并换行


Sample Input
2
4 2
5 2

Sample Output
6
10

题解：：通过看题本题的n,m<=20，可以应用组合数的知识直接求解。但有点要注意mod不满足除法，不能在求阶乘时mod。

代码：：

#include <iostream>
#include <sstream>
#include <ios>
#include <iomanip>
#include <functional>
#include <algorithm>
#include <vector>
#include <string>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <climits>
#include <cctype>
using namespace std;
#define XINF INT_MAX
#define INF 0x3FFFFFFF
#define MP(X,Y) make_pair(X,Y)
#define PB(X) push_back(X)
#define REP(X,N) for(int X=0;X<N;X++)
#define REP2(X,L,R) for(int X=L;X<=R;X++)
#define DEP(X,R,L) for(int X=R;X>=L;X--)
#define CLR(A,X) memset(A,X,sizeof(A))
#define IT iterator
typedef long long ll;
typedef pair<int,int> PII;
typedef vector<PII> VII;
typedef vector<int> VI;
const int MOD = 1e9+7;
ll factorial(ll a)
{
	ll fac = 1,cnt = a;
	for(int i = 1;i<cnt;i++){
		fac = fac*a;
		a -=1;
	}
	return fac;	
}
int main(){
	int T;
	cin>>T;
	while(T--)
	{
		ll a,b;
		cin>>a>>b;
		ll fa,fb,fc;
		fa = factorial(a);
		fb = factorial(b);
		fc = factorial(a - b);
		cout<<(fa/fb/fc)%MOD<<endl;
	}
	return 0;
}

shuoj1938-组合数level 1

Description
求C(n,m)对1000000007取模后的值，已知公式如下
Input
第一行有一个整数T，表示有T组数据(T≤10^5)

接下来有T行，每行有两个整数n，m，用空格隔开。0≤n≤1000,0≤m≤n

Output
对于每组数据，输出答案并换行
Sample Input
2
4 2
5 2
Sample Output
6
10
HINT
输入和输出的数据规模较大，请使用scanf和printf
Source
xyiyy
题解：： 这道题很明显的不同之处（1）数据量变大了（2）n,m的值变大了。这时候取模的作用真正的用上了。我们的问题转化成组合数取模而不是单单的求组合数，应用上面的方法求解不再满足要求（会超时）。
       但我们在初中都学过C(i,j) = C(i-1,j) + C(i-1,j-1).很自然的我们想到用dp的思想或者记忆化搜索的方法将0<=n,m<=1000的所有值存到数组中，每次询问直接调用。

dp代码：：

#include <iostream> 
#include <sstream> 
#include <ios> 
#include <iomanip> 
#include <functional> 
#include <algorithm> 
#include <vector> 
#include <string> 
#include <list> 
#include <queue> 
#include <deque> 
#include <stack> 
#include <set> 
#include <map> 
#include <cstdio> 
#include <cstdlib> 
#include <cmath> 
#include <cstring> 
#include <climits> 
#include <cctype> 
using namespace std; 
#define XINF INT_MAX 
#define INF 0x3FFFFFFF 
#define MP(X,Y) make_pair(X,Y) 
#define PB(X) push_back(X) 
#define REP(X,N) for(int X=0;X<N;X++) 
#define REP2(X,L,R) for(int X=L;X<=R;X++) 
#define DEP(X,R,L) for(int X=R;X>=L;X--) 
#define CLR(A,X) memset(A,X,sizeof(A)) 
#define IT iterator 
typedef long long ll; 
typedef pair<int,int> PII; 
typedef vector<PII> VII; 
typedef vector<int> VI; 
const ll MOD = 1000000007; 
ll dp[1005][1005]; 

void com_num(int t,int mod)//生成1000以内的所有组合数
{
	REP2(i,0,t) { 
		REP2(j,0,i) { 
            if(j == 0||j==i)dp[i][j] = 1; 
            else dp[i][j] = (dp[i-1][j]%mod+dp[i-1][j-1]%mod)%mod; 
        } 
    }
}
int main() 
{ 
     com_num(1000,MOD);
    int T; 
    scanf("%d",&T); 
    while(T--) 
    { 
        int a,b; 
        scanf("%d%d",&a,&b); 
        printf("%lld\n",dp[a][b]); 
    } 
    return 0; 
  
}

记忆化搜索代码：： 

#include <sstream> 
#include <ios> 
#include <iomanip> 
#include <functional> 
#include <algorithm> 
#include <vector> 
#include <string> 
#include <list> 
#include <queue> 
#include <deque> 
#include <stack> 
#include <set> 
#include <map> 
#include <cstdio> 
#include <cstdlib> 
#include <cmath> 
#include <cstring> 
#include <climits> 
#include <cctype> 
using namespace std; 
#define XINF INT_MAX 
#define INF 0x3FFFFFFF 
#define MP(X,Y) make_pair(X,Y) 
#define PB(X) push_back(X) 
#define REP(X,N) for(int X=0;X<N;X++) 
#define REP2(X,L,R) for(int X=L;X<=R;X++) 
#define DEP(X,R,L) for(int X=R;X>=L;X--) 
#define CLR(A,X) memset(A,X,sizeof(A)) 
#define IT iterator 
typedef long long ll; 
typedef pair<int,int> PII; 
typedef vector<PII> VII; 
typedef vector<int> VI; 
const ll MOD = 1000000007; 
ll dp[1005][1005]; 

//计算C(m,n)对mod取模，每求一次就要调用一次;A,B是记忆化搜索的关键
ll com_num(int m,int n,int mod){
	if(dp[m][n]) return dp[m][n];//已经搜索过的避免重复
	if(m == n||n == 0)return 1;//终止条件
	dp[m][n] = (com_num(m-1,n-1,mod)+com_num(m-1,n,mod))%mod;
	return dp[m][n];
}
int main() 
{ 
    int T; 
    scanf("%d",&T); 
    while(T--) 
    { 
        int a,b; 
        scanf("%d%d",&a,&b); 
        printf("%lld\n",com_num(a,b,MOD));//直接输出com_num()的返回值即可 
    } 
    return 0; 
  
}

shuoj1939-组合数level2

Description

求C(n,m)对k取模后的值，已知公式如下

Input

第一行有一个整数T，表示有T组数据(T≤10)

接下来有T行，每行有三个整数n，m，k，用空格隔开。

0≤n≤10^5，0≤m≤n，1≤k≤10^9+7

Output

对于每组数据，输出答案并换行

Sample Input

2
4 2 100000000 7
5 2 1000000007

Sample Output



6
10

Source

xyiyy

题解：：可以发现m，n从1000涨到100000。我们可以想到的一种方法通过把n,m,(n-m)分解素因子将因子存放在数组中，让n中的减掉m,(n-m)中的然后乘起来对mod取模。这是一种解题的好思路，但是会遇到两个问题：：（1）分解素因子（2）把多个因子乘起来（如果说有1e5个2,1e5个3……就很有可能超时）所以我们还需要一个求快速幂的方法。这两个问题解决了就ok了。

（1）我们可以通过枚举[2 , 根号n]的方式求出n的素因子分解

（2）我们发现看k^13 = k^8*k^4*k^1用二进制表示k^(1101) = k^(1000)*k^(100)*k^(1)，求快速幂的诀窍就在这里。思路代码里见或者点下面的链接：快速幂

代码：：

#include <iostream>
#include <sstream>
#include <ios>
#include <iomanip>
#include <functional>
#include <algorithm>
#include <vector>
#include <string>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <climits>
#include <cctype>
using namespace std;
#define XINF INT_MAX
#define INF 0x3FFFFFFF
#define MP(X,Y) make_pair(X,Y)
#define PB(X) push_back(X)
#define REP(X,N) for(int X=0;X<N;X++)
#define REP2(X,L,R) for(int X=L;X<=R;X++)
#define DEP(X,R,L) for(int X=R;X>=L;X--)
#define CLR(A,X) memset(A,X,sizeof(A))
#define IT iterator
typedef long long ll;
typedef pair<int,int> PII;
typedef vector<PII> VII;
typedef vector<int> VI;
//const int MAXN = 10010;
//#define INF 0x3FFFFFFF
//const int MOD = 1e9+7;
int a[100005];//此处数组不能小于100000
int b[100005];
int c[100005];
void fenjie1(int x){//分解竭诚
    int k;
    for(k = x;k>1;k--){
        int j = k;int i;
        for(i = 2;i*i<=j;i++){
            while(j%i==0){
                a[i]++;
                j /=i;
            }
        }
        if(j!=1)
        a[j]++;
    }
}
void fenjie2(int x){
    int k;
    for(k = x;k>1;k--){
        int j = k;int i;
        for(i = 2;i*i<=j;i++){
            while(j%i==0){
                b[i]++;
                j /=i;
            }
        }
        if(j!=1)
        b[j]++;
    }
}
void fenjie3(int x){
    int k;
    for(k = x;k>1;k--){
        int j = k;int i;
        for(i = 2;i*i<=j;i++){
            while(j%i==0){
                c[i]++;
                j /=i;
            }
        }
        if(j!=1)
        c[j]++;
    }
}
ll fast_mod(int n,int m,int k){
    ll ans = 1;
    while(m){
        if(m&1)ans = ans * n % k;
        m>>=1;
        n = n*n%k;
    }
    return ans;
}
int main(){
    int T ;
    int m,n,k;
    ll answer;
    cin>>T;
    while(T--){
        cin>>m>>n>>k;
        answer = 1;
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        memset(c,0,sizeof(c));
        fenjie1(m);
        fenjie2(n);
        fenjie3(m-n);
        REP(i,100005){
            a[i] = a[i] - b[i] - c[i];
            if(a[i]!= 0)
            answer = answer * fast_mod(i,a[i],k)%k;
        }
        cout<<answer<<endl;
    }
    return 0;
}


1940: 组合数level3

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 26   Solved: 13
[ Submit][ Status][ Web Board]

Description

求C(n,m)对k取模后的值，已知公式如下

Input

第一行有一个整数T，表示有T组数据(T≤100)

接下来有T行，每行有三个整数n，m，k，用空格隔开。

0≤n≤10^5，0≤m≤n，1≤k≤10^9+7，并且保证k为素数

Output

对于每组数据，输出答案并换行

Sample Input

2
4 2 1000000007
5 2 1000000007

Sample Output

6
10

Source

xyiyy



题解：：由题知，组数T从10 变成100，这样会造成TLE，同时k也有了变化，从任意数变成素数。这是就要转换思路求解，我们知道通过费马小定理知道：


a*x =1 mod k 

>>

 x = a ^ -1 mod k 

>>

 (当k为素数时，a^k mod k = a mod k:费马小定理)

>>

a^(k-1)  modk = 1 mod k

>>

x = a^-1 mod k = a^(k-2) mod k  (a的逆元)

所以设inv(a,k)为求a在mod为k时的逆元

C(n,m) = fac[n]/(fan[m]*fan[n-m] mod k) = fac[n] * inv (fan[m]*fan[n-m] mod k, k) mod k;

当且仅当gcd(a,k) == 1时逆元存在

#include <iostream>
#include <sstream>
#include <ios>
#include <iomanip>
#include <functional>
#include <algorithm>
#include <vector>
#include <string>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <climits>
#include <cctype>
using namespace std;
#define XINF INT_MAX
#define INF 0x3FFFFFFF
#define MP(X,Y) make_pair(X,Y)
#define PB(X) push_back(X)
#define REP(X,N) for(int X=0;X<N;X++)
#define REP2(X,L,R) for(int X=L;X<=R;X++)
#define DEP(X,R,L) for(int X=R;X>=L;X--)
#define CLR(A,X) memset(A,X,sizeof(A))
#define IT iterator
typedef long long ll;
typedef pair<int,int> PII;
typedef vector<PII> VII;
typedef vector<int> VI;

ll fac[100005];
void qfac(int n,int k){
    fac[1] = 1;
    REP2(i,2,n)
    fac[i] = fac[i-1]*i%k;
}
ll fast_mod(ll n,int m,int k){
    ll ans = 1;
    while(m){
        if(m&1) ans = ans*n%k;
        m>>=1;
        n = n*n%k;
    }
    return ans;
}
ll inv(ll n,int k){
    ll ans = fast_mod(n,k-2,k);
    return ans;
}

int main(){
    int m,n,k;
    int T;
    cin>>T;
    while(T--){
        cin>>m>>n>>k;
        memset(fac,0,sizeof(fac));
        qfac(m,k);
        //REP(i,10)cout<<fac[i]<<endl;
        cout<<fac[m]*inv(fac[n],k)%k*inv(fac[m-n],k)%k<<endl;
    }
    return 0;
}



6
10

Source

xiaohuihuicb
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
为什么要对1000000007取模（取余）
skr？
04-30 883
大数阶乘，大数的排列组合等，一般都要求将输出结果对1000000007取模（取余） 为什么总是1000000007呢= = 大概≖‿≖✧是因为：（我猜的，不服你打我呀~） 1000000007是一个质数（素数），对质数取余能最大程度避免冲突～ int32位的最大值为2147483647，所以对于int32位来说1000000007足够大 int64位的最大值为2^63-1，对于1000000007来说它的平方不会在int64中溢出 所以在大数相乘的时候，因为(a∗b)%c=((a%c)∗(b%c))%c，
【贪心】【大数比较大小】组合数问题
weixin_44111457的博客
09-13 1303
【描述】 给你两个数 n 和k ，希望你找到 k 个不同的组合数使得这 k个组合数的和最大。所谓不同的组 合数，即对于组合数 和 ，若a1!=a2 或者 b1!=b2，则我们认为这两个组合数是不同 的。现在小葱希望你找到这样 k 个不同的组合数，使得它们互不相同且对于其中任何一个组 合数 有0<=a<=b<=n 。问这 k个组合数的和最大是多少？ 【输入】 第一行两个整数 n,k...
求组合数 计蒜客T1984
最新发布
m0_54149707的博客
11-12 150
从 n 个不同元素中，任取 m(m \le n)m(m≤n) 个元素并成一组，叫做从 nn 个不同元素中取出 mm 个元素的一个组合；从 nn 个不同元素中取出 m(m \le n)m(m≤n) 个元素的所有组合的个数，叫做从 nn 个不同元素中取出 mm 个元素的组合数，记为 C_n^mCnm​。组合数的计算方法为 C_n^m = \frac{n!}{m!(n-m)!}Cnm​=m!(n−m)!n!​，其中!!表示阶乘。
牛客小白月赛8-C（组合数公式）
Alibaba_lhl的博客
10-21 281
Description clccle一行
组合数（卢卡斯定理）
★探梦少年☆的博客
02-24 4650
组合数 从n个不同元素中，任取m(m≤n)个元素并成一组，叫做从n个不同元素中取出m个元素的一个组合；从n个不同元素中取出m(m≤n)个元素的所有组合的个数，叫做从n个不同元素中取出m个元素的组合数。 公式：C(n,m)=n!/((n-m)!*m!)（m≤n） 性质1：C(n,m)= C(n,n-m) 性质2：C(n,m)=C(n-1,m-1)+C(n-1,m) 定义 组合是数学...
高中数学排列组合专题.doc
10-07
）和组合数(C(n,r)=n!/r!(n-r)!)，并学会在实际问题中合理使用这些公式。此外，还要注意排除法、插空法、捆绑法、倍缩法等解决排列组合问题的策略。通过这些题目的练习，可以加深对排列组合概念的理解，并提高解决...
专题资料（2021-2022年）产品策略产品组合策略DOC 49页.doc
10-08
产品组合策略涉及到产品宽度（产品线的数量）、深度（每条产品线内的品种数）、关联度（产品线之间的关联性）和产品生命周期管理。企业需要根据市场需求和自身资源，决定产品组合的结构，例如： 1. 附带服务的产品...
排列与组合专题教案设计.doc
10-21
在本专题教案中，教学目标聚焦于理解和掌握排列的定义以及排列数的计算公式，旨在提升学生的数学计算技能、计算工具使用技能和数学思维能力。 【排列】指的是从给定的不同元素中，按照特定顺序选取一部分元素形成...
IT通用技术专题-存储
08-03
磁盘阵列，如RAID（冗余磁盘阵列），通过组合多个硬盘来提高性能和数据冗余，以防止单个硬盘故障导致的数据丢失。控制器负责管理数据流，而磁盘柜则用于容纳多个硬盘。磁带库则用于长期存储和归档，它们经济且适合大...
组合数学课件及试题——哈工程
04-16
这些专题进一步丰富了组合数学的学习，并为学生提供了更深层次的理论基础和问题解决技巧。 总的来说，这个资源对于想要深入理解组合数学的学生来说是一份宝贵的资料，它不仅提供理论知识，还包含了实践应用，有助于...
组合数详解
a_forever_dream的博客
10-17 3794
介绍 组合数，也就是形如 CnmC_n^mCnm​ 这样的东西。 意义有很多，最经典的是：表示在 nnn 个不同的物品里面选 mmm 个物品的方案数。 下面就以这个意义来讲递推式和通项式。 递推式 Cnm=Cn−1m−1+Cn−1mC_n^m=C_{n-1}^{m-1}+C_{n-1}^mCnm​=Cn−1m−1​+Cn−1m​ 从 nnn 个里面选 mmm 个，那么对于第 nnn 个物品，有两种...
组合数
Blogs of zcy1221
07-22 612
本文将介绍组合数的定义、性质和实际计算方法.
组合数与排列数
我思故我在
03-09 3万+
1 组合数与排列数是什么 组合：有一个袋子，里面有10个标有1-10数字的球，问如果每次拿出3个球，一共有多少不同的组合？（1，2，3）和（3，1，2）算一种，不考虑次序。 排序：有一个袋子，里面有10个标有1-10数字的球，问如果每次拿出3个球，一共有多少不同的排序？（1，2，3）和（3，1，2）算二种，考虑次序。 2 如何表示 ...
计算组合数的三种方式
GD_ONE的博客
03-19 1万+
本文将详细介绍计算组合数的三种方式，这三种方式分别适用于不同的场景。
求组合数-算法学习
dalao_whs的博客
12-05 235
求组合数
常用组合数计算公式及推算
热门推荐
ccsu_deer
04-27 7万+
参考：博客1 博客2 更多更详细请看博客2 组合数的通项公式： 公式1： 证明： n个不同的数选择m个，第m个的选择方案为： 1、选第m个： 2、不选第m个: 公式2： 证明： 性质3： 证明： 性质4： 证明： 性质5： ...
组合数的各种性质和定理
litble的成(tui)长(fei)史
07-23 5万+
从m个物品里选出n个的方案数，记作CnmCmnC_m^n，即为组合数 组合数有很多很多的性质和定理。。。 注意由于本人沉迷玩梗无法自拔，如果看见您看不懂的梗请随意跳过。 组合数通项公式 Cnm=m!n!∗(m−n)!Cmn=m!n!∗(m−n)!C_m^n=\frac{m!}{n!*(m-n)!} 证明：现在我们从m个不同的数里选出n个数组成一个排列，第一个位子上的数有m种取法，第二...
组合数（阶乘的质因子的个数，组合数的计算）
qq_51825761的博客
07-22 209
的Pi的个数，个数分别为下，x,y,z.则C(n,m)中质因子Pi的个数为(x-y-z);将C(n,m)进行质因子分解，对遍历不超过n的质数Pi分别计算。通过定义式计算将C(n,m)进行质因子分解；支持n
组合c(m,n)的计算方法
烟杀流年
07-10 1万+
问题：求解组合数C(n,m)，即从n个相同物品中取出m个的方案数，由于结果可能非常大，对结果模10007即可。       共四种方案。ps：注意使用限制。 方案1: 暴力求解，C(n,m)=n*(n-1)*...*(n-m+1)/m!， n ； int Combination(int n, int m) { const int M = 10007;
方法技巧专题21 排列组合与二项式定理（解析版）.docx
10-15
方法技巧专题21 排列组合与二项式定理（解析版）.docx

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值