Educational Codeforces Round 89 (Rated for Div. 2) D. Two Divisors(线性筛)

题意：

给出$n$个询问，每次询问给出一个数$x$，问是否能找到$x$的两个因子$d_1$,$d_2$ ,满足$gcd(d_1+d_2,x)=1$

题解：

需要先知道$gcd$的一个性质：

如果$gcd(x,y)=1$ ,那么$gcd(x,y)=gcd(x+y,x\cdot y)=1$

先引出两个定理：

定理$1.$ $gcd(x,y)=gcd(x-y,y)$ ,$gcd(x,y)=gcd(x+y,y)$

定理$2.$ 若$gcd(x,y)=1$,则 $gcd(x,ky)=gcd(x,k)$

证明如下：

因为 $gcd(x,y)=gcd(x+y,y)=1$ ,所以 $gcd(x+y,xy)=gcd(x+y,x)=gcd(y,x)=1$

再来看这道题：

我们只需要使得 $d_1\cdot d_2=x$ 且 $gcd(d_1,d_2)=1$ 即可满足条件

那么将$x$ 唯一分解成 $a^{t_1}\cdot b^{t_2}\cdot c^{t_3}\cdot d^{t_4}....$ ，然后取$a^{t_1}$ 为$d_1$ ,剩下的为 $d_2$ 就行。

最后就是简单的套线性筛的板子预处理即可。

代码：

#pragma GCC diagnostic error "-std=c++11"
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<set>
#include<ctime>
#define iss ios::sync_with_stdio(false)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+7;
const int MAXN=1e7+5;
const int inf=0x3f3f3f3f;
int prime[MAXN],cnt;
int isprime[MAXN];
int dp[500005][2];
void getprime(int n)
{
    for(ll i=2;i<=n;i++)
    {
        if(!isprime[i])
        {
            prime[++cnt]=i;
            isprime[i]=i;
        }
        for(int j=1;j<=cnt&&i*prime[j]<=n;j++)
        {
            isprime[i*prime[j]]=prime[j];
            if(i%prime[j]==0) break;
        }
    }
}
int main()
{
    getprime(1e7);
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        int x;
        scanf("%d",&x);
        ll ans=1;
        ll temp=x;
        while(temp!=1&&temp%isprime[x]==0)
        {
            temp/=isprime[x];
            ans=ans*isprime[x];
        }
        if(temp==1)
        {
            dp[i][0]=-1;
            dp[i][1]=-1;
        }
        else 
        {
            dp[i][0]=ans;
            dp[i][1]=x/ans;
        }
    }
    for(int i=1;i<=n;i++) printf("%d ",dp[i][0]);
    printf("\n");
    for(int i=1;i<=n;i++) printf("%d ",dp[i][1]);
}