算法基础课 第四章 数学知识

算法基础课 第四章 数学知识

数论

质数

1. 定义：对于所有大于1的自然数，如果只包含1和本身这两个约数，就被称为质数或者素数
2. 质数的判定：试除法（枚举的时候只要枚举到d<=n/d）

#include<iostream>
#include<algorithm>
using namespace std;
bool is_prime(int n)
{
   
	if (n < 2)return flase;
	for (int i = 2;i <= n;i++)
		if (n % i == 0)
			return false;
	return true;
}

3. 分解质因数：试除法：从小到大尝试n的所有因数

#include<iostream>
#include<algorithm>
using namespace std;
int n;//n中最多只包含一个大于根号n的质因子
void divide(int x)
{
   
	for(int i=2;i<= n;i++)
		if (n % i == 0)
		{
   
			int s = 0;
			while (n % i == 0)//i一定是质数
			{
   
				n /= i;
				s++;
			}
			printf("%d %d\n", i, s);
		}
	if (n > 1)printf("%d %d\n", n, 1);
	puts("");

}
int main()
{
   
	int n;
	scanf_s("%d", &n);
	while (n--)
	{
   
		int x;
		scanf_s("%d", &x);
		divide(x);
		
	}
	return 0;
}

筛选质数

从前往后看，将每一个数的所有倍数删除，删除完后，所有剩下的数都是质数
朴素筛选法

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1000010;
int primes[N], cnt;
bool st[N];
void get_primes(int n)
{
   
	for (int i = 2;i <= n;i++)
	{
   
		if (!st[i])
		{
   
			primes[cnt++] = n;
		}
		for (int j = i + i;j <= n;j += i)st[j] = true;
	}
}
int main()
{
   
	int n;
	cin >> n;
	get_primes(n);
	cout << cnt << endl;
	return 0;
}

改进（埃式筛法）

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1000010;
int primes[N], cnt;
bool st[N];
void get_primes(int n)
{
   
	for (int i = 2;i <= n;i++)
	{
   
		if (!st[i])
		{
   
			primes[cnt++] = n;
			for (int j = i + i;j <= n;j += i)st[j] = true;
		}
	}
}
int main()
{
   
	int n;
	cin >> n;
	get_primes(n);
	cout << cnt << endl;
	return 0;
}

线性筛法
n只会被它的最小质因子筛选掉
![

](https://img-blog.csdnimg.cn/94c8b8ded4b24e2ab139ed03edb04777.png)

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1000010;
int primes[N], cnt;
bool st[N];
void get_primes(int n)
{
   
	for (int i = 2;i <= n;i++)
	{
   
		if (!st[i])primes[cnt++] = i;
		for (int j = 0;primes[j] <= n / i;j++)
		{
   
			st[primes[j] * i] = true;
			if (i % primes[j] == 0)break;//primes一定是i的最小质因子
		}
	}
}
int main()
{
   
	int n;
	cin >> n;
	get_primes(n);
	cout << cnt << endl;
	return 0;
	return 0;
}

约数

1. 试除法求所有约数：从小到大判断，如果这个数能整出当前数，说明是其约数。优化思路：如果d能整除n，那么n/d也一定能整除n。即一个数的约数一定是成对出现的，在枚举时可以只枚举到较小的约数。

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
vector<int>get_divisors(int n)
{
   
	vector<int>res;
	for(int i=0;i<=n/i;i++)
	{
   
		if(n%i==0)
		{
   
			res.push_back(i);
			if(i!=n/i)res.push_back(n/i);
		}
	}
	sort(res.begin(),res.end());
	return res;
		
}
int main()
{
   
	int n;
	cin>>n;
	while(n--)
	{
   
		int x;
		cin>>x;
		auto res=get_divisors(x);
		for(auto t:res)cout<<t<<' ';
		cout<<endl;
	}
	return 0;
}

2. 约数个数：

3. 约数之和：
   

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<unordered_map>//哈希表
using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;
int main()
{
   
	int n;
	cin >> n;
	unordered_map<int,int>primes;
	while (n--)
	{
   
		int x;
		cin >> x;
		for (int i = 2;i <= x / i;i++)
		{
   
			while (x % i == 0)
			{
   
				x /= i;
				primes[i]++;//i因子的指数加加
			}
			if (x > 1)primes