一、数论
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)
*/