D. Two Divisors（素数筛）

题目传送门

题意： 给你n个数，对于每一个a[i]，如果能找到两个a[i]的因子d1>1,d2>1，使得gcd(d1+d2,a[i]）=1，则输出这两个因子，否则输出-1,-1。

思路： 对于一个a[i]，我们把他的最小因子d1除尽，设剩下d2,那么d1和d2肯定互质，那么就有gcd(d1,d2)=1,即gcd(d1+d2,d1)=gcd(d1+d2,d2)=1，即gcd(d1+d2,d1* d2)=1，d1 * d2是a[i]的因子，即a[i]/(d1 * d2)=d1^k,所以gcd（d1+d2,a[i]）=1成立，注意题目条件对d1、d2!=1特判即可。

代码：

#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define pii pair<int,int>
#define ull unsigned long long
#define all(x) x.begin(),x.end()
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read(){int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
return x*f;}
using namespace std;
const int N=1e7+5;
const int inf=0x3f3f3f3f;
const int mod=1e9+7;
const double eps=1e-7;
const double PI=acos(-1);
int a[N/20+5],prim[N/10+5],d[N],ans1[N/20+5],ans2[N/20+5],cnt=0;
bool ok[N];
void euler()
{
    d[1]=1;
    for(int i=2;i<N;i++)
    {
        if(!ok[i])
        {
            prim[++cnt]=i;
            d[i]=i;
        }
        for(int j=1;j<=cnt;j++)
        {
            if(i*prim[j]>=N)
                break;
            ok[i*prim[j]]=1;
            d[i*prim[j]]=prim[j];
            if(i%prim[j]==0)
                break;
        }
    }
}
signed main()
{
    int n;
    cin>>n;
    euler();
    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        int t=d[a[i]];
        if(t==a[i])
        {
            ans1[i]=ans2[i]=-1;
            continue;
        }
        while(a[i]%t==0)
            a[i]/=t;
        if(a[i]==1)
        {
            ans1[i]=ans2[i]=-1;
            continue;
        }
        ans1[i]=t;
        ans2[i]=a[i];
    }
    for(int i=1;i<=n;i++)
        cout<<ans1[i]<<' ';
    cout<<endl;
    for(int i=1;i<=n;i++)
        cout<<ans2[i]<<' ';
}