hdu 3864 D_num

思路：给一个数n，是否只有4个约数(包括1)，也就是找3个大于1的约数。

而任何一个数都可由质数表示，所以对于给定的数，只需要进行质因数分解。这里有

2种情况：如果有3个一样的质因数，则满足条件；否则只需要2个不同的质因子。

代码如下：

 

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<cstring>
#include<vector>
#define ll __int64
#define pi acos(-1.0)
#define MAX 5000001
using namespace std;
ll n,e[5];
int prime[MAX],cnt;
bool f[MAX];
void init()
{
    int i,j;
    cnt=0;
    for(i=2;i<MAX;i++){
        if(f[i]==0) prime[cnt++]=i;
        for(j=0;j<cnt&&i*prime[j]<MAX;j++){
            f[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
void solve()
{
    ll i;
    ll num=n,k=0;
    for(i=0;i<cnt&&prime[i]*prime[i]<=n;i++){
        while(n%prime[i]==0){
            e[k++]=prime[i];
            n/=prime[i];
        }
        if(k>3){
            cout<<"is not a D_num"<<endl;
            return;
        }
    }
    if(n>1) e[k++]=n;
    if(k>3||k<2){
        cout<<"is not a D_num"<<endl;
        return;
    }
    else if(k==3){
        if(e[0]==e[1]&&e[1]==e[2]){
            cout<<e[0]<<' '<<e[0]*e[1]<<' '<<num<<endl;
            return;
        }
        else {
            cout<<"is not a D_num"<<endl;
            return;
        }
    }
    else{
        if(e[0]==e[1]){
            cout<<"is not a D_num"<<endl;
            return;
        }
        else{
            cout<<e[0]<<' '<<e[1]<<' '<<num<<endl;
            return;
        }
    }
}
int main(){
    init();
    while(cin>>n){
        solve();
    }
    return 0;
}