数论基础

素数的判定

bool isprime(int n){
	for(int i = 2; i <= sqrt(n); i++)
		if(n % i == 0) return false;
	return true;
}

素数筛

时间复杂度 $O(n)$

bool isprime[maxn];
int prime[maxn]; //保存的全是素数
int cnt;    //素数的个数
void getprime(int n){
	memset(isprime, true, sizeof isprime);
	isprime[0] = isprime[1] = false;
	cnt=0;
	for(int i = 2; i <= n; i++){
		if(isprime[i])
			prime[cnt++] = i;
		for(int j = 0; j < cnt && i * prime[j] <= n; j++){
			isprime[i * prime[j]] = false;
			if(i % prime[j] == 0) break;
		}
	}
}

素因数分解

$N=p_1^{c_1}{p}_2^{c_2}\dots p_m^{c_m}$
时间复杂度 $O(\sqrt{n})$

int m;
int p[maxn],c[maxn];
void divide(int n){
	m=0;
	for(int i = 2; i <= sqrt(n); i++){
		if(n % i == 0){
			p[++m] = i,c[m] = 0;
			while(n % i == 0) n /= i, c[m]++;
		}
	}
	if(n > 1)  //n是素数
		p[++m] = n, c[m] = 1;
	for(int i = 1; i <= m; i++)
		cout<<p[i]<<'^'<<c[i]<<endl;
}

求 N 的正约数集合

时间复杂度 $O(\sqrt{n})$

int factor[maxn];
void fac(int n){
	int m=0;
	for(int i = 1; i*i <= n; i++){
		if(n % i == 0){
			factor[++m] = i;
			if (i != n/i) factor[++m] = n/i;
		}
	}
	for(int i = 1; i <= m; i++)
		cout<<factor[i]<<endl;
}

扩展欧几里得

$Bezout$ 定理

对于任意整数 $a,b$，存在一对整数 $x,y$，满足 $ax+by=gcd(a,b)$。

//返回a,b的最大公约数，并算出这样的x,y;
int exgcd(int a, int b, int &x, int &y){
	if (b == 0) { x = 1, y = 0; return a; }
	int d = exgcd(b, a % b, x, y);
	int z = x; x = y; y = z - y * (a / b);
	return d;
}

乘法逆元

  若整数 $b,m$ 互质，并且 $b|a$，则存在一个整数 $x$，使得 $a/b=a\ast x(modm)$。称 $x$ 为 $b$ 的模 $m$ 乘法逆元，记为 $b^{-1}(modm)$。
  $b\ast x\equiv 1(modm)$ 的 $x$ 的解就是 $b$ 的逆元，可以通过扩展欧几里得算出。
求解过程：
  由上式子$b\ast x\equiv 1(modm)$ 得 $b\ast x-1$ 是 $m$ 的整数倍，我们假设就为 $-y$ 倍，即 $b\ast x-1=-m\ast y$，所以就有 $b\ast x+m\ast y=1$。
乘法逆元递推
$inv[i]=(p-p/i)\ast inv[p\%i]\%p(p!=2,且为素数)$

线性同余方程

  给定整数 $a,b,m$，求一个整数 $x$ 满足 $a\ast x\equiv b(modm)$，或者给出无解,有解仅当 $gcd(a,m)|b$。因为未知数 $x$ 的指数为 $1$，所以我们称之为一次同余方程，也称线性同余方程。
  $a\ast x\equiv b(modm)$ 等价于 $a\ast x-b$ 是 $m$ 的倍数，不妨设为 $-y$ 倍。于是，该方程可以改写为 $a\ast x+m\ast y=b$，则该 $x$ 就为方程的解。
  根据 $Bezout$ 定理及其证明过程，线性同余方程有解仅当 $gcd(a,m)|b$。
  所以存在整数 $k$ 使得 $k=b/gcd(a,m)$，即 $gcd(a,m)=b/k$。所以有 $a\ast (x/k)+m\ast (y/k)=gcd(a,m)$，令 $x^{{}^{\prime }}=x/k,y^{{}^{\prime }}=y/k$，所以有 $a\ast x^{{}^{\prime }}+m\ast y^{{}^{\prime }}=gcd(a,m)$。可以通过扩展欧几里得算出 $x^{{}^{\prime }}$ 从而算出 $x$。