【数学】

一、数论

1.质数

质数的判定——试除法

//复杂度：O（√n）
bool is_prime(int n){
	if(n < 2) return 0;
	for(int i = 2; i <= n/i; ++i)
		if(n%i==0) return 0;
	return 1;
}

2. 分解质因数——试除法

/*
分解质因数，暴力遍历所有<=n的数 
复杂度：O（n）
*/
void divide(int n){
	for(int i = 2; i <= n; ++i){
		if(n%i == 0){
			int s = 0;
			while(n%i == 0){ // 36/2=18;18/2=9;9/2除不开 
				n/=i;
				s++;
			}
			cout << "因数："<<i<<",个数："<<s<<endl;
		}
	}
}

/*
优化：因为n中只包含一个大于√n的质因子，所以for循环跑到√n即可。
证明：若存在两个大于√n的因子，则相乘大于n

复杂度：O（√n）
当x=2^n时，复杂度为O(logn)
*/
void divide(int n){
	for(int i = 2; i <= n/i; ++i){
		if(n%i == 0){
			int s = 0;
			while(n%i == 0){ // 36/2=18;18/2=9;9/2除不开 
				n/=i;
				s++;
			}
			cout << "因数："<<i<<",个数："<<s<<endl;
		}
	}
	if(n > 1) //如果最后一个数大于1，那么他也是n的因子	
		cout << "因数："<<n<<",个数："<<1<<endl;
}

3. 素数筛

素数定理：
1~n中，最多有n/lnn个素数

/*
线性筛

复杂度：O（n）
*/
const int N = 10010;
int cnt,prine[N];
bool vis[N];
void Prime(int n){
	for(int i = 2; i <= n; ++i){
		if(!vis[i]) prime[cnt ++] = i;
		for(int j = 0; prime[j] <= n/i; ++j){
			vis[prime[j]*i] = 1;
			if(i % prime[j] == 0) break; 
		}	
	}
}

约数

/*
(1)求一个数的约数
复杂度：O（√n）

*/
vector<int> get_divisors(int n){
	vector<int>res;
	for(int = 1; i <= n/i; ++i){//n=16 
		if(n%i == 0){
			res.push_back(i);// 存入：1，2，4 
			if(i != n/i) res.push_back(n/i); //存入：16，8;去充根号n 
		}
	}
	sort(res.begin(),res.end());
	return res;
}

/*
如果 N = p1^c1 * p2^c2 * ... *pk^ck
约数个数： (c1 + 1) * (c2 + 1) * ... * (ck + 1)
约数之和： (p1^0 + p1^1 + ... + p1^c1) * ... * (pk^0 + pk^1 + ... + pk^ck)
*/

二、组合计数

三、高斯消元

四、简单博弈论