素数、最大公因数、快速幂的函数

一、求最大公因数的函数
int gcd(int x,int y)//辗转相除法求最大公因式数 
{
	int t;
	do{
		if(x<y)
		{
			t=x;
			x=y;
			y=t;
		}
		x=x%y;
	}while(x);
	return y;
}


int gcd2(int x,int y)
{
	return y?gcd2(y,x%y) : x;
}


int gcd3(int a,int b)
{
	if(b==0)
	return a;
	return gcd3(b,a%b);
}


int gcd4(int n,int m)
{
	int i=0;
	while(m)
	{
		i=n%m;
		n=m;
		m=i;
	}	
	return n;
	
} 
如果求最小公倍数 lcm=a*b/gcd(a,b);
早c++中有内置求最大公因数的函数  __gcd(a,b);
直接可以调用

二、素数的判定
//判断a是否为素数 
bool prim(int a)
{
	if(a<2)
	return 0;
	int len = sqrt(a+1); 
	for(i=2;i<=len;i++)
	{
		if(a%i==0)
		return 0;
	}
	return 1;
}

 
 //素数打表法
const int maxn = 1e6+5;
int prim[maxn];
 
void init()
{
	memset(prim,0,sizeof(prim));
	prim[0] = prim[1] = 1;
	int len = sqrt(maxn + 1);
	for(int i = 2;i < len;i++)
	{
		if(prim[i])
		{
			for(int j = i + i;j<maxn;j += i)
			prim[j]=1;	
		}	
	}	
}

 三、快速幂的函数
 //求快速幂的函数 
int poww(int a,int b,int mod)

 - List item

{
	int ans=1;
	while(b)
	{
		if(b & 1)
		ans = ans * a % mod;
		
		a = a * a % mod;
		b>>=1;//b/=2
			
	}	
	return ans;
}