数论小集（1）
话说当初我不选择数学竞赛与其他竞赛的原因，就是我以为信息学竞赛有关数学的方面不怎么较难。结果到数论时，我就GG了，所以这次借着复习，巩固数论。
（注：本文以实用和总结为主，大部分定理与公式是没有证明的）
数论小集（基础篇）
1.质数
既然质数的英文名叫做“Prime number”，而Prime又有“基础的”的意思，那么它是很重要的。
判断一个数是不是质数可以用O(sqrt（n）)的复杂度（枚举1~sqrt(n)的数
），但是若要求1~n的质数，则可以用欧拉筛以O（n）的算法求出。代码：
void Prime(int n)//
{
cnt=0;//cnt从零开始
memset(vis,0,sizeof(vis));
for(int i=2;ivoid Prime(int num)
{
cnt=0;
memset(minn,0,sizeof(minn));
for(int i=2;i<=num;i++)
{
if(!minn[i])
{
prime[cnt++]=i;
minn[i]=i;
}
for(int j=0;j2.gcd与lcm
gcd是最大公因数，lcm是最小公倍数。
gcd用辗转相除法，复杂度为O（logn）的，而lcm则建立在gcd之上：lcm（a，b）=a/gcd（a，b）*b
代码：
int gcd(int a,int b)//最大公约数
{
return b==0?a:gcd(b,a%b);
}
int lcm(int a,int b)//最小公倍数
{
return a/gcd(a,b)*b;//注意：这里应该先除后乘，减少溢出的可能
}3.模运算基础
有：
（a+b）%mod=（（a%mod）+（b%mod
））%mod
（a*b）%mod=（（a%mod）*（b%mod））%mod
（a-b）%mod=（（a%mod）-（b%mod）+mod）%mod
以及快速幂（O（logn））取模：
int pow(int base,int b)//建议使用循环，而不用递归（容易崩栈！血的教训！）
{
int res=1;
while(b)
{
if(b&1) res=res*base%mod;
base=base*base%mod;
b>>=1;
}
return res;
}4.（重点）扩展欧几里得算法
扩展欧几里得能够找出一对（x，y），使得ax+by=gdc(a,b)，并有以下简易结论：
结论1： 若方程ax+by=c的一组整数解为(x0,y0),则它的任意整数解都可以写成(x0+kb',y0-ka'),其中a'=a/gcd(a,b),b'=b/gcd(a,b),k取任意整数。
结论2： g=gcd(a,b),方程ax+by=g的一组解是(x0,y0),则当c是g的倍数时ax+by=c的一组解是(x0*c/g,y0*c/g);否则无解。
结论2： g=gcd(a,b),方程ax+by=g的一组解是(x0,y0),则当c是g的倍数时ax+by=c的一组解是(x0*c/g,y0*c/g);否则无解。
它的作用很多，以后应该会讲到，其中最基本的就是他可以解两元一次方程组（利用以上两个结论）。
代码：
void gcd(int a,int b,int g,int &x,int &y)
{
if(!b) {g=a;x=1;y=0;}
else {gcd(b,a%b,g,y,x);y-=x*(a/b);}
}
5.计数问题
基础的加法原理，乘法原理，以及容斥原理应该都很简单（很多dp都会用）
主要是C（n，m），其实P（n，m）都用得很少。
关于C（n，m）简单的有以下几种算法
void C(int n,int m)
{
/*1.*/for(int i=1;i<=n;i++)
{
f[i][i]=f[i][0]=1;
for(int j=1;j<=i;j++)
{
f[i][j]=f[i-1][j]+f[i-1][j-1];//排列数的递推公式
}
}//递推式
/*2.*/f[n][m]=n!/m!(n-m)!//公式，可以预处理阶乘，但由于要用除法，n和m太大后要溢出，
//取模又没有除法运算（但有逆元）,总之就是比较蛋疼的算法
/*3.*/c[0]=1;利用C(n,k)=(n+1-k)/k*C(n,k-1)
for(int i=1;i<=n;i++)c[i]=c[i-1]*(n-i+1)/i;实际意义不明显，但。。。很好玩
} 以上就是一部分笔者觉得重要的基础，虽然应用感觉很少，但是接下来大多数都建立在这些之上，至少把代码记住还是很重要的。