算法的数学知识

1.欧几里得定理

又称辗转相除法，即gcd(a,b)=gcd(b,a%b); gcd(x,y)表示x与y的最大公约数,在求最小公倍数的时候可以直接用ab/gcd(a,b)；
扩展欧几里得：存在一组数x、y，使得ax1+by1=gcd(a,b);因此延伸出对x、y的求解。
首先 gcd(a,b)=gcd(b,a%b);将右式按照上式进行代换，得到bx2+[a-a/bb]y2=gcd(b,a%b);(注：此处除法为整除)
由于两式右端相等，因此可以得到ax1+by1=bx2+[a-a/bb]y2;按照a,b进行整理，
得到ax1+by1=ay2+b*(x2-a/by2);
因此返回的时候，写法如下，int t=x;x=y;y=t-a/by;递归函数的定义形式为exgcd(int a,int b,int &x,int &y);
3.a b交换
1).三元交换法
int t=a;
a=b;
b=t
2).作差法
int t=a+b;
a=t-a;
b=t-a;
3).快速写法
b^=a=b^=a;
3.a^=b=a^{=b;这是基于二进制的方法，理解方法如下：a}=b得到的结果是a与b在二进制下不同的位为1，将这个结果和b抑或之后，得到的就是a,因为a^{b得到的数中的某些1是一定是b提供的，与b抑或之后，这些位就不存在了。还有一个疑问就是当1与1抑或成0怎么办，也很简单，如果是这种情况，假设是ci=ai}bi=0，aibi1.那么当bi或者ai与这个0抑或时又会变回1，所以没有影响。
于是欧几里得算法可以用位运算进行加速 while(b^=a=b^=a%=b); 解析：gcd(int a,int b);终止条件为 gcd(x,0); 因此终止时b为0，赋值运算符是从右到左结合的，所以while里的表达式的最终返回值为b的值，因此是可以在b为0时停止的。
b^=a=b^=a%=b;
第一步：a%=b;
第二步到最后：就是在交换a,b的值。
4.斯特林公式

用来取n阶乘近似值的数学公式。n!=sqrt(2Pin)*(n/e)^n；
5.求数的位数

(int)floor（log10(x)+0.5）;
注意：c\c++里面的log是log10(X),而不是log2(X)。

如何在c\c++里使用lnx？

1).exp(x)函数 表示自然对数e^x。

2).换底公式 loga(b)/loga©=logc(b); 所以 lnx=log(X)/log(exp(1));
6.数的奇偶判断
x&1 举个栗子，假设目前类型是8位short类型,那么1的写法就是00000001,也就是x&1之后得到的只是最后一位。奇偶判断原理是，奇数=小于他的最大偶数+1 偶数=他自己+0
7.矩阵乘法
for (int i = 0; i <= 1; i++)
{
for (int j = 0; j <= 1; j++)
{
for (int k = 0; k <= 1; k++)
{
A[i][j] += B[i][k] * C[k][j];
}
}