a^b
题目链接
求 a的 b次方对 p 取模的值。
输入格式
三个整数 a,b,p ,在同一行用空格隔开。
输出格式
输出一个整数,表示a^b mod p的值。
样例输入
3 2 7
样例输出
2
思路
- 首先我们会想到直接使用 pow 得到次方后结果,然后取余直接得到结果
- 但是数据范围过大,pow 的方法不可用
- 然后我们可以使用 快速幂
坑点
数据范围过大,pow 方法只能过部分样例,所以我们需要将 pow 方法 pass掉
快速幂模板
int st(int a,int b,int p)
{
int res=1,t=a;
while(b)
{
if(b&1)//==>b%2=1
{
res=res*t%p;
}
t=t*t%p
b>>=1;//==>b/=2
}
return res;
}
算法:快速幂
代码
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
long long int st(long long int a,long long int b,long long int c)
{
long long int ans=1%c;//当b==0,p==1时
while(b)
{
if(b&1)// ==> b%2=1
{
ans=(ans*a)%c;
}
a=(a*a)%c;
b/=2;
}
return ans;
}
int main()
{
long long int n,m,p;
cin>>n>>m>>p;
cout<<st(n,m,p)<<endl;
return 0;
}
总结
一道简单的快速幂模板应用题
蓝桥杯算法提高VIP-欧拉函数
题目链接
给定一个大于1,不超过2000000的正整数n,输出欧拉函数,phi(n)的值。如果你并不了解欧拉函数,那么请参阅提示。
提示
欧拉函数phi(n)是数论中非常重要的一个函数,其表示1到n-1之间,与n互质的数的个数。显然的,我们可以通过定义直接计算phi(n)。
当然,phi(n)还有这么一种计算方法。
首先我们对n进行质因数分解,不妨设n=p1^a1 * p2^a2 * … * pk^ak (这里a^b表示a的b次幂,p1到pk为k个互不相同的质数,a1到ak均为正整数),那么
phi(n)=n(1-(1/p1))(1-(1/p2))…(1-(1/pk))
稍稍化简一下就是
phi(n)=n(p1-1)(p2-1)…(pk-1)/(p1p2…*pk)
计算的时候小心中间计算结果超过int类型上界,可通过调整公式各项的计算顺序避免(比如先做除法)!
输入格式
在给定的输入文件中进行读入:
一行一个正整数n。 不超过2000000的正整数n
输出格式
将输出信息输出到指定的文件中:
一行一个整数表示phi(n)。
样例输入
17
样例输出
16
题意
给定一个数 n ,在1~ n 的范围内,寻找与 n 只存在最大公约数为 1 的个数;
思路
涉及最大公约数使用 gcd
gcd 模板
int gcd(int a,int b)
{
if(a%b==0)
{
return b;
}
else
{
return gcd(b,a%b);
}
}
算法:gcd(最大公约数)
代码
#include<stdio.h>
int gcd(int a,int b)
{