组合数取模(费马小定理)

费马小定理:   

 

p为质数) 

我们由此可以推出:

m是质数

所以a在mod m 意义下的逆元等价于a^m-2,这样我们就会求一个数的逆元了。

注意:费马小定理只能在a,m互质时使用!

这里我们会想可不可以直接带公式求解呢?

因为涉及到除法的取模,所以要用到逆元求解,自然想到使用费马小定理来求一个数的逆元

思路 

这样我们就可以快速的求出组合数取模后的结果了

题目链接:组合数取模

我的代码:

### Codeforces 上与费马小定理相关的题目及其实现 #### 费马小定理简介 费马小定理指出,如果 $p$ 是质数且 $a$ 不是 $p$ 的倍数,则有: $$ a^{p-1} \equiv 1 \ (\text{mod}\ p) $$ 由此可以推导出,在 $p$ 下求逆元的方法:对于任意整数 $a$ (满足 $\gcd(a, p) = 1$),其在 $p$ 意义下的逆元为 $a^{p-2}$。 --- #### 经典题目解析与实现 以下是几道经典的 Codeforces 题目以及它们基于费马小定理的解决方法: --- ##### **题目一:Codeforces 1380G** 该题要求计算一系列操作后的结果,并涉及大量运算。其中的关键在于利用费马小定理快速求解逆元[^3]。 ###### 解决方案 通过预处理阶乘及其逆元数组,可以在常数时间内完成组合数的查询。下面是完整的代码实现: ```cpp #include <bits/stdc++.h> using namespace std; typedef long long LL; const int MAXN = 3e5 + 5; const LL MOD = 998244353; // 快速幂板 LL qpow(LL a, LL b){ LL res = 1; while(b){ if(b & 1) res = res * a % MOD; a = a * a % MOD; b >>= 1; } return res; } LL fact[MAXN], inv_fact[MAXN]; void preprocess(){ fact[0] = 1; for(int i=1;i<MAXN;i++) fact[i] = fact[i-1]*i%MOD; inv_fact[MAXN-1] = qpow(fact[MAXN-1], MOD-2); for(int i=MAXN-2;i>=0;i--){ inv_fact[i] = inv_fact[i+1]*(i+1)%MOD; } } // 计算组合数 C(n,k) mod MOD LL comb(LL n, LL k){ if(k > n || k < 0) return 0; return fact[n] * inv_fact[k] % MOD * inv_fact[n-k] % MOD; } int main(){ ios::sync_with_stdio(false); cin.tie(0); preprocess(); // 输入处理... } ``` 此代码片段实现了对阶乘和逆元的预处理功能,从而能够高效地支持后续的大规组合数计算需求[^3]。 --- ##### **题目二:Advertising Agency (Codeforces Round #697 Div. 3)** 这是一道关于广告公司分配预算的问题,核心思路也是围绕着组合数展开讨论[^4]。 ###### 解决方案 由于直接计算可能会超出计算机所能表达的最大数值范围,因此采用动态规划配合费马小定理间接获答案成为必然选择之一。下面给出了一种可行的做法: ```cpp #include<bits/stdc++.h> using namespace std; const int N=2e5+5,MOD=1e9+7; long long fac[N],inv[N]; vector<int>v; long long power(long long x,long long y){ long long ret=1; while(y){ if(y&1)ret=(ret*x)%MOD; x=x*x%MOD;y>>=1; } return ret; } void precompute(){ fac[0]=fac[1]=1; for(int i=2;i<N;i++)fac[i]=(fac[i-1]*i)%MOD; inv[N-1]=power(fac[N-1],MOD-2); for(int i=N-2;i>=0;i--)inv[i]=(inv[i+1]*(i+1))%MOD; } long long C(int n,int r){ if(r>n||r<0)return 0; return ((fac[n]*inv[r]%MOD)*inv[n-r])%MOD; } int main(){ precompute(); int T,n,m,x,y; scanf("%d",&T); while(T--){ v.clear(); scanf("%d%d%d%d",&n,&m,&x,&y); for(int i=1;i<=n;i++){ int tmp;scanf("%d",&tmp); if(tmp<=y)v.push_back(tmp); } sort(v.begin(),v.end()); long long ans=0; for(int i=v.size();i>0;i--){ if((long long)i*m>x)break; ans+=C(v.size()-i+i-1,i-1); ans%=MOD; } printf("%lld\n",ans); } } ``` 这里不仅包含了基本的幂次方函数定义,还额外加入了针对特定条件筛选有效元素列表的操作步骤。 --- #### 总结 上述两例充分体现了费马小定理在实际竞赛编程中的广泛应用价值——无论是用来加速单步除法操作还是辅助构建更大规的数据结构型都极为重要! ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值