[自用]C语言函数

/* gcd 求最大公约数 */
unsigned int gcd(unsigned int a,unsigned int b){
	while(b^=a^=b^=a%=b);
	return a;
} 

/*     扩展欧几里德算法是用来在已知a, b求解一组p，q使得ax+by = Gcd(a, b) =d(根据数论中的相关定理,解一定存在)。
扩展欧几里德常用在求解模线性方程及方程组中。*/
//C #include <math.h>
int exGcd(int a, int b, int &x, int &y){
	int r,t;
	if(b == 0) 
	{x = 1;y = 0;return a;} 
	r = exGcd(b, a % b, x, y); 
	t = x; 
	x = y; 
	y = t - a / b * y; 
	return r; 
} 


/*
欧拉φ函数：φ(n)是所有小于n的正整数里，和n互素的整数的个数。n是一个正整数。 　　
欧拉证明了下面这个式子： 　　
如果n的标准素因子分解式是p1^a1*p2^a2*……*pm^am，其中众pj(j=1,2,……,m)都是素数，而且两两不等。
则有 　　φ(n)=n(1-1/p1)(1-1/p2)……(1-1/pm) 　　
利用容斥原理可以证明它。
*/
int euler(int x)
{
    int i, res=x;
    for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)
        if(x%i==0)
        {
            res = res / i * (i - 1);
            while (x % i == 0) x /= i; // 保证i一定是素数
        }
    if (x > 1) res = res / x * (x - 1);
    return res;
}
/*
求逆序数
n为数组长度 a为数组地址
*/
#include    <stdio.h>  
#include    <string.h>  

int inv(int n,char *a)  
 { 
 	int   l=n> > 1,r=n-l,i,j; 
 	char   *v=new   char(n); 
 	int   ret=(r> 1?(inv(l,a)+inv(r,a+l)):0); 
 	for   (i=j=0;i <=l;v[i+j]=a[i],i++) 
 		for   (ret+=j;j <r&&(i==l||a[i]> a[l+j]);v[i+j]=a[l+j],j++); 
 	memcpy(a,v,sizeof(char)*n); 
 	return   ret; 
 }