本笔记主要参考ACWing基础课数论和繁凡初等数论笔记,笔记不会一次性全部写完,会持续更新。本笔记的内容可能并不完整,只会记录做题时的高频知识点,如有遗漏欢迎私信我来补充。
整除
- b/a的余数为0,记作a|b
- a|b,b|c => a|c
- a|b => a|bc
- a|b,a|c=>a|bK+cL,K,L为任意整数
质数(素数)
定义:只有1和它本身两个因数的数。
质数的判定
可以根据定义,通过for循环暴力从2到其本身枚举是否i|n。
- 试除法
inline bool is_prime(int x){
if(x < 2)
return false;
for( int i = 2; i <= x / i;i++)
if(x % i == 0)
return false;
return true;
}
这里注意到并没有循环到n
,而是循环到了n/i
,因为我们惊奇的发现每个数字的因数有一大一小,如果找到小的因子时候就已经可以证明该数为非质数了。例如:21的因子有3和7只需要找到3就可以了。
bool v[N];
void get_prime(int n)
{
for (int i = 2; i <= n; i++)
{
if (v[i]) continue;
cout<< i <<endl;
for (int j = i ; j <= n / i ; j ++)
v[j] = true;
}
}
思想:如果x是合数,那么x的倍数也一定是合数。
- 线性筛(欧拉筛)
const int N = 1e5+5;
int n,primes[N],cnt;
bool vis[N];
inline void get_prime()
{
for(int i = 2,i <= n;i++)
{
if(!vis[i])
primes[cnt++] = i;
for(int j = 1;j <= cnt && i*primes[j] <= n;j++)
{
vis[i * primes[j]] = 1;
if(i % primes[j] == 0)
break;
}
}
}
思想:保证每个合数只都被自己最小的质因子筛掉一遍。
质数筛相关例题:
ACWing 筛质数
AcWing 1293. 夏洛克和他的女朋友
质数判定二
哥德巴赫猜想
质数距离
//质数距离题解
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int N = 1000010;
int primes[N], cnt;
bool st[N];
void init(int n)
{
memset(st, 0, sizeof st);
cnt = 0;
for (int i = 2; i <= n; i ++ )
{
if (!st[i]) primes[cnt ++ ] = i;
for (int j = 0; primes[j] * i <= n; j ++ )//欧拉筛模板
{
st[i * primes[j]] = true;
if (i % primes[j] == 0) break;
}
}
}
int main()
{
int l, r;
while (cin >> l >> r)
{
init(50000);
memset(st, 0, sizeof st);
for (int i = 0; i < cnt; i ++ )
{
LL p = primes[i];
for (LL j = max(p * 2, (l + p - 1) / p * p); j <= r; j += p)
//j初始化为KP,使得KP>L,即找到最少P需要多少倍可以大于等于左边界L,并且至少需要是二倍。
st[j - l] = true;
}
cnt = 0;
for (int i = 0; i <= r - l; i ++ )
if (!st[i] && i + l >= 2) //特判,因为P可能为1
primes[cnt ++ ] = i + l;
if (cnt < 2) puts("There are no adjacent primes.");
else
{
int minp = 0, maxp = 0;
for (int i = 0; i + 1 < cnt; i ++ )// 筛完之后暴力跑
{
int d = primes[i + 1] - primes[i];
if (d < primes[minp + 1] - primes[minp]) minp = i;
if (d > primes[maxp + 1] - primes[maxp]) maxp = i;
}
printf("%d,%d are closest, %d,%d are most distant.\n",
primes[minp], primes[minp + 1],
primes[maxp], primes[maxp + 1]);
}
}
return 0;
}
先记这么多,下次更新约数orz