数论扫盲
非专业分类,仅供参考。
素数与合数
约数
筛素数
欧拉筛:
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);
}
欧拉定理
费马小定理
威尔逊定理
逆元
莫乌比斯反演
原根
整数的阶
求原根
离散对数