简 单 数 论

1.质数判断（试除法）

#include <bits/stdc++.h>
using namespace std;
int isprime(int x); 
int main()
{
	for(int i=2;i<=100;i++)
	{
		if(isprime(i)==1)
		{
			cout<<i<<" ";
		}
	 } 
	return 0;
}
int isprime(int x)
{
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{
			return 0;//表示x是合数 
		} 
	}
	return 1;//如果return 0 没发生 表示x是质数 ； 
}

1.2 埃氏筛法 

//埃氏筛法 （以空间换时间） 
#include <bits/stdc++.h>
using namespace std;
#define N 10000000
int isprime[N+1];//0表示为质数 
int main()
{
	for(int i=2;i<=N;i++)
	{
		if(isprime[i]==0)
		{		
			for(int j=i+i;j<=N;j+=i)
			{
				isprime[j]=1;	
			}
		}
	}
	for(int i=2;i<=N;i++)
	{
		if(isprime[i]==0)
		{
			cout<<i<<" ";
		}
	}
	return 0;
}

2.最大公约数（辗转相除法，欧几里德算法）

#include <bits/stdc++.h>
using namespace std;
int gcd (int x,int y)
{
	return y==0?x:gcd(y,x%y);
}
int main()
{
	cout<<gcd(10,15);
	return 0;
}

3.最小公倍数（辗转相除法，大数倍增法）

#include <bits/stdc++.h>
using namespace std;
int gcd (int x,int y)
{
	return y==0?x:gcd(y,x%y);
}
int lcm (int x,int y)
{
	return x*y/gcd(x,y);
}

int main()
{
	cout<<" "<<lcm(10,15);
	return 0;
}

4.同余方程

#include <bits/stdc++.h>
using namespace std;
int m1(int a,int b,int c)//计算(a+b)%c 
{
	return (a%c+b%c)%c;
}
int m2(int a,int b,int c)//计算(a*b)%c 
{
	return (a%c*b%c)%c;
}
int main()
{
	cout<<m1(3,5,2)<<endl;
	cout<<m2(3,5,2);
	return 0;
}

5.杨辉三角

#include <bits/stdc++.h>
using namespace std;
long long a[100][100];
void tranyang(int n)
{
	a[1][1]=1;
	for(int i=2;i<=n;i++)
	{
		for(int j=1;j<=i;j++)
		{
			a[i][j]=a[i-1][j-1]+a[i-1][j];//左上加上 
		}
	}	
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=i;j++)
		{
			cout<<a[i][j]<<" "; 
		}
		cout<<endl;
	}
}
int main()
{
	tranyang(5);
	return 0;
}

6.二项式定理

6.1 排列

#include <bits/stdc++.h>
using namespace std;
long long A(int m,int n)
{
	long long f=1;
	for(int i=m;i>=m-n+1;i--)
	{
		f*=i;
	}
	return f;
}
int main()
{
	cout<<A(6,3);
	return 0;
}

 

 6.2 组合

#include <bits/stdc++.h>
using namespace std;
long long A(int m,int n)
{
	long long f=1;
	for(int i=m;i>=m-n+1;i--)
	{
		f*=i;
	}
	return f;
}
long long C(int m,int n)
{
	return A(m,n)/A(n,n);
}
int main()
{
	cout<<C(5,2);
	return 0;
}

 6.3 二项式定理

 

#include <bits/stdc++.h>
using namespace std;
long long A(int m,int n)
{
	long long f=1;
	for(int i=m;i>=m-n+1;i--)
	{
		f*=i;
	}
	return f;
}
long long C(int m,int n)
{
	return A(m,n)/A(n,n);
}
void exs(char a,char b,int n)
{
	for(int i=0;i<=n;i++)
	{
		cout<<'+'<<C(n,i)<<a<<'^'<<n-i<<'*'<<b<<'^'<<i;
	}
}
int main()
{
	exs('a','b',2);
	return 0;
}

7.卡特兰数

#include <bits/stdc++.h>
using namespace std;
long long A(int m,int n)
{
	long long f=1;
	for(int i=m;i>=m-n+1;i--)
	{
		f*=i;
	}
	return f;
}
long long C(int m,int n)
{
	return A(m,n)/A(n,n);
}
long long catalan(int n)
{
//	return C(2*n,n)-C(2*n,n+1);//公式1 
	return C(2*n,n)/(n+1);//公式2
}
int main()
{
	for(int i=0;i<=10;i++)
	{
		cout<<catalan(i)<<" ";
	}
	return 0;
}

8.容斥原理