[自用]C语言函数
来源:互联网 发布:dsp广告平台源码 编辑:程序博客网 时间:2024/06/07 07:01
/* 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; }