快速幂
int qupow(int a,int b,int m)
{
int ans=1;
while(b!=0)
{
if(b%2==1)
ans=ans*a%m;
a=a*a%m;
b=b/2;
}
return ans;
}
位运算代码为
int qupow(int a,int b,int m)
{
int ans=1;
while(b)
{
if(b&1)
ans=ans*a%m;
a=a*a%m;
b>>=1;
}
return ans;
}
快速幂将原来0(b)的时间复杂度降为了O(logb)
素数
首先我们在一般方式求素数时,可以采用以下方式
bool ifprime(int x)
{
if(x==1) return false;
for(int i=2;i*i<=x;i++)
{
if(x%i==0) return false;
}
return true;
}
为何用ii<=x; 因假如一个数x=nm,那么n和m必定有一个小于x/2;所以在x/2前不存在能整除x的整数,那么x就一定是素数。
在做题中如果需要判断大量的数,如还用上述方法将使效率大大降低,这种情况下可以选择使用素数筛,素数筛也就是在查找数之前将所有的数都判断好是否为素数,存于数组,等输入要判断的数时就可直接通过数组记录来判断出是否为素数。
普通素数筛
时间复杂度为 O(nloglogn)
bool v[maxn+5];
memset(v,true,sizeof(v));
for(int i=2;i<=n;i++)
{
if(v[i])
{
for(int j=2;i*j<=n;j++)
v[i*j]=false;
}
}
普通素数筛有很多数组都被重复遍历,所以可以采用线性筛,也称欧拉筛。
欧拉筛的时间复杂度为O (n)
欧拉筛
bool v[maxn+5];
int prime[maxn+5];
void isprime()
{
memset(v,true,sizeof(v));
v[1]=false;
int c=0;
for(int i=2;i<=n;i++)
{
if(v[i])
prime[c++]=i;
for(int j=0;j<c&&i*prime[j]<=n;j++)
{
v[i*prime[j]]=false;
if(i%prime[j]==0) break; //保证每个合数只会被它的最小质因数筛去,因此每个数只会被标记一次
}
}
}