算术分解定理取根号的注意点

很多时候，对于一个数据范围很大的数，如n$\in$(1, 2e9)分解质因数，若循环n次就会超时，所以通常做法是循环$\sqrt{n}$次，再有因数的性质，算出大于$\sqrt{n}$的情况，然而此方法有些部分需要引起注意，以下代码为例

void find_prime_factor(int n) {
	for (int i = 2; i <= sqrt(n); i++) {
		if (n % i == 0) {
			prime[++cnt] = i;
			while (n % i == 0)n /= i;
		}
	}
	//！！有缺陷
}

则最后n有两种可能，1或大于$\sqrt{n}$的质数(注意n可能改变过)
因此，若题目要求仅“打印所有质因数“，要在！处要加上

if (n > 1)
		prime[++cnt] = n;

以完善最答案

计算一个数的欧拉函数，也有相应的知识点

int phi(int n) {
	int ans = n;
	for (int i = 2; i <= sqrt(n); i++) {
		if (n % i == 0) {
			ans = ans / i * (i - 1); //1式
			while (n % i == 0)n /= i;//2式
		}
		if (n > 1)ans / n * (n - 1);
	}
	return ans;

其中2式好处有2：
1缩小n的值，使循环的次数减少；
2n的最终值可能是1或一个大于$\sqrt{n}$的单个的素因数p，2式排除用过的素因数，从而计算出，补足了缺陷

但有的题目，$\sqrt{n}$的缺陷并不会很容易就能看出，也要注意在最后补上这个。如Hankson 的趣味题
由《算法竞赛进阶指南》p144的解法二，我们将a,b,c,d分别分解质因数，再由各质因数的次数来组成答案

#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int MAX_N = 1e5 + 10;

int prime[MAX_N], cnt_prime = 0, v[MAX_N];
void find_prime() {
	memset(v, 0, sizeof(v));
	for (int i = 2; i <= 50000; i++) {
		if (v[i] == 0) v[i] = i, prime[++cnt_prime] = i;
		for (int j = 1; j <= cnt_prime; j++) {
			if (i * prime[j] > 50000 || prime[j] > v[i])break;
			v[prime[j] * i] = prime[j];
		}
	}
}

int main() {
	find_prime();
	int n;
	scanf("%d", &n);
	while (n--) {
		int a, b, c, d;
		int ans = 1;
		scanf("%d%d%d%d", &a, &b, &c, &d);
		
		int flag = 1;
		for (int i = 1; i <= cnt_prime; i++) {
			if (a == 1 && b == 1 && c == 1 && d == 1) break;
			int ma = 0, mb = 0, mc = 0, md = 0, cntx;
			int p = prime[i];
			while (a % p == 0)	a /= p, ma++;
			while (b % p == 0)	b /= p, mb++;
			while (c % p == 0)	c /= p, mc++;
			while (d % p == 0)	d /= p, md++;
			if (ma > mb && mc < md && mb == md) cntx = 1;
			else if (ma == mb && mc == md && mb <= md) cntx = md - mb + 1;
			else if (ma == mb && mc < md && mb <= md)cntx = 1;
			else if (ma > mb && mc == md && mb <= md)cntx = 1;
			else flag = 0;
			if (flag == 0) {
				ans = 0; break;
			}
			ans = ans * cntx;
		}
		//试除法取根号的缺陷，若最后剩下一个大于根号n的质数，则需特判
		if (!(a == 1 && b == 1 && c == 1 && d == 1)) {
			if (a == b && b != 1) ans <<= 1;
			if (c == d && d != 1) ans <<= 1;
		}
		printf("%d\n", ans);
	}
	return 0;
}

对于这一题，我们先用试除法($\sqrt{n}$)找出几乎所有的质因数，中间隔了一大段操作，在最后输入输出的时候才有关于n的判断，实际竞赛中，可能很容易忘记。故要引起注意

PS：刚学数论部分，例子接触的少，讲解可能又不玩完善处，欢迎指正