初等数论笔记

本笔记主要参考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