数论扫盲

数论扫盲

非专业分类，仅供参考。

素数与合数

约数

筛素数
欧拉筛：

const int maxn = 1e7+10;
bool book[maxn];
int prime[maxn],cnt;
void Ola_prime(int n)
{
	cnt=0;
	memset(book,0,sizeof(book));
	book[1] = 1;
	for(int i = 2; i <= n; i++)
	{
		if(book[i]==0)
			prime[cnt++]=i;
		for(int j = 0; j < cnt && prime[j]*i <= n; j++)
		{
			book[i*prime[j]] = 1;
			if(i%prime[j]==0)
				break;
		}
	}
}

区间筛：

typedef long long ll;
const int maxn = 1e6+10;
bool vis[maxn];
bool book[maxn];
int ans;

void region_prime(ll l,ll r)
{
	for(ll i = 2;i*i <= r;i++)
	{
		if(!vis[i])
		{
			for(ll j = i; j*j <= r;j += i)
				vis[j]=1;
			for(ll j=max(2ll,(l+i-1)/i)*i;j <= r;j += i)
				book[j-l]=1;
		}
	}
	ans = 0;
	for(int i = 0;i <= r-l;i++)
		if(book[i]==0)
			ans++;
}

反素数
对于任何正整数x其约数的个数记作g(x).例如g(1)=1,g(6)=4. {1,2,3,6}.如果某个正整数x满足：对于任意的i∈(0,x),都有g(i)<g(x),则称x为反素数。
解决方案：基于因子个数公式，按照质因子数大小递增顺序搜索每一个质因子，枚举每个质因子的个数。
性质一：一个反素数的质因子必然是从2开始连续的质数。
性质二：C₁≥C₂≥C₃≥…≥Cг
n!的质因子分解

大数判素数(不保证正确)
费马小定理：若n是奇数，a是任意整数（1≤a≤n-1），则a^(n-1) ≡1(mod n)
二次探测定理：如果p为质数，且0<x<p则x²≡1(mod p) 的解为x=1||x=p-1;
Miller-Rebin素数测试

pollard-rho大数分解

余数

GCD
辗转相除：

ll gcd(ll a, ll b)
{
	   return !b ? a : gcd(b, a%b);
}

更相减损数

ll gcd(ll a, ll b)
{
    if(a%2==0 &&b%2==0)
    		return 2gcd(a/2,b/2);
    else(a%2==0)
    	return gcd((a-b)/2,b);
}

扩展欧几里得



代码：

int  extend_gcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    int r =  extend_gcd(b,a%b,y,x);
    y-=x*(a/b);
    return r;
}

线性方程

中国剩余定理



二次剩余

https://blog.csdn.net/acdreamers/article/details/10182281

欧拉函数

欧拉函数
欧拉函数是求小于等于n的数中与n互质数的个数，用φ(n)表示。

代码：

void get_eular()//1-N的欧拉函数
{
    for(int i = 1;i < N;i++)
        eular[i]= i;
    for(int i = 2;i < N;i++)
        if(eular[i] == i)
            for(int j = i;j < N;j+=i)
                euler[j] = eular/i*(i-1);
}

欧拉定理

费马小定理

威尔逊定理

逆元

莫乌比斯反演

原根

整数的阶

求原根


离散对数