【数论】快速幂+欧拉降幂
2023-02-19 12:20:37 时间
快速幂的原理很简单,下面来看一下
2^6=2^2*2^2*2^2
其实就是通过降幂
long long qpow(int x,int y){
int res=1;
while(y){
if(y&1)res=res*x;
x=x*x;
y>>=1;
}
return res;
}
然后由于MOD的性质
(a*b)%p=(a%p*b%p)%p【当然,下面这个公式可能会更好理解点(a^b)%c=[(a%c)^b]%c】
所以我们可以再优化一下。
long long qpow(int x,int y,int mod){
int res=1;
while(y){
if(y&1)res=res*x%mod;
x=x*x%mod;
y>>=1;
}
return res;
}
再想一想,如果x有点大,那么第一次乘的时候可能会爆(但是一般情况下上面那个就够用了),那么,我们在循环之前先给x取余一下
long long qpow(int x,int y,int mod){
int res=1;
x%=mod;
while(y){
if(y&1)res=res*a%mod;
a=a*a%mod;
y>>=1;
}
return res;
}
上面就是快速幂。
好了,下面来讲一讲欧拉降幂。
在这里,先说说欧拉函数。
f(x)=在小于x的数中,与x互质的数的总个数(更多的若要了解自行百度)
那欧拉降幂就会用到下面的一个公式。
其中
代表的是f(c)的值
所以我们对b还可以进行一次降幂。
int Euler(int n) //1.直接求欧拉函数的值
{
int rea=n;
for(int i=2; i*i<=n; i++)
if(n%i==0)//第一次找到的必为素因子
{
rea=rea-rea/i;
while(n%i==0){
n=n/i;//把该素因子全部约掉
}
}
if(n>1)
rea=rea-rea/n;
return rea;
}
long long qpow(int x,int y,int mod){
int res=1;
x%=mod;
y=y%Euler(y)+Euler(y);
while(y){
if(y&1)res=res*a%mod;
a=a*a%mod;
y>>=1;
}
return res;
}
相关文章
- Jgit的使用笔记
- 利用Github Action实现Tornadofx/JavaFx打包
- 叹息!GitHub Trending 即将成为历史!
- 微软软了?开源社区讨论炸锅,GitHub CEO 亲自来答
- GitHub Trending 列表频现重复项,前后端都没去重?
- Photoshop Elements 2021版本软件安装教程(mac+windows全版本都有)
- (ps全版本)Photoshop 2020的安装与破解教程(mac+windows全版本都有)
- (ps全版本)Photoshop cc2018的安装与破解教程(mac+windows全版本,包括2023
- 环境搭建:Oracle GoldenGate 大数据迁移到 Redshift/Flat file/Flume/Kafka测试流程
- 每个开发人员都要掌握的:最小 Linux 基础课
- 来撸羊毛了!Windows 环境下 Hexo 博客搭建,并部署到 GitHub Pages
- 超实用!手把手入门 MongoDB:这些坑点请一定远离
- 【GitHub日报】22-10-09 zustand、neovim、webtorrent、express 等4款App今日上新
- 【GitHub日报】22-10-10 brew、minio、vite、seaweedfs、dbeaver 等8款App今日上新
- 【GitHub日报】22-10-11 cobra、grafana、vue、ToolJet、redwood 等13款App今日上新
- Photoshop 2018 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2017 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2020 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2023 资源免费下载(mac+windows全版本都有,包括最新的2023)
- 最新版本Photoshop CC2018软件安装教程(mac+windows全版本都有,包括2023