Two Divisors[欧筛,埃筛]

思路参考:

对于gcd有
gcd(a,b)=gcd(a+b,b)
且若gcd(a,b)=1,
则gcd(a,c)=gcd(a,c*b)=gcd(a,c * b^n)
现有gcd(d1,d2)=1,->gcd(d1+d2,d1)=1,gcd(d1+d2,d2)=1
->gcd(d1+d2,d1 * d2)=1,->gcd(d1+d2,d1^n * d2^n)

现在再来看题目要我们求什么，我们只需找到ai的两个互质因数d1,d2 , ai=d1^ n * d2可得gcd(d1,d2)=1,由如上性质gcd(d1+d2，d1^n*d2)=1,即gcd(d1+d2,ai)=1。

在此我们用埃拉托色尼筛法在 O(NloglogN)求出1e7内每个数的第一个质因数d1，第二个与d1互质因数只需用ai当ai%d1==0时不断除d1即可求出，如果最终d2=1表示找不到两个互质因数，反则可以。
————————————————
版权声明：本文为CSDN博主「_Rikka_」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/asbbv/article/details/113750050

问题:用埃筛求出来的并不是第一个质因数,而用欧拉筛求出来的才是第一个质因数,但两种写法都可以ac.

原因:不用必须是第一个质因数,因为根据唯一分解定理,一个和数只能被唯一分解为一对质因数的乘积的形式,所以只要找到一个质因数就可以,并根据这一个质因数找到这一对中的另一个质因数.

上文的链接中有用埃筛的写法,鄙人在原有的基础上改了一个欧筛的写法,时间更优.

#include<iostream>
using namespace std;
typedef long long ll;
const int Max = 5e5 + 5;
int prim[Max*20], vis[Max*20], h[Max*20], x = 0;//h存储合数找到的第一个质因数
int a[Max], b[Max];
void eratos(int n)
{
	for (int i = 2;i <= n;i++)
	{
		if (!vis[i])prim[++x] = i;
		for (int j = 1;j * i <= n;j++)
		{
			if(prim[j] * i > n){
				break;
			}
			vis[prim[j] * i] = 1;
			h[prim[j] * i] = prim[j];
			if(i % prim[j] == 0){
				break;
			}
		//	if(j!=1)vis[i * j] = 1;
		//	h[i * j] = i;
		}
	}
}

int main()
{
	int n;cin >> n;
	eratos(1e7+3);
//	printf("%d",h[12]);
	for (int i = 1;i <= n;i++)
	{
		int p;cin >> p;
		if (vis[p] == 0)a[i] = -1, b[i] = -1;//素数本身就是-1 -1 
		else
		{
			//printf("%d",h[p]); 
			int k = h[p];
			while (p % k == 0)p /= k;
			if (p == 1)a[i] = -1, b[i] = -1;
			else a[i] = k, b[i] = p;
		}
	}
	for (int i = 1;i <= n;i++)cout << a[i] << " ";cout << endl;
	for (int i = 1;i <= n;i++)cout << b[i] << " ";
}