反素数讲解

反素数

Description
对于任何正整数x，其约数的个数记作g(x)。例如g(1)=1、g(6)=4。如果某个正整数x满足：g(x)>g(i) 0<i<x，则称x为反质数。例如，整数1，2，4，6等都是反质数。现在给定一个数N，你能求出不超过N的最大的反质数么？

Input
一个数N（1≤N≤2,000,000,000）。

Output
不超过N的最大的反素数。
Samples
Input
1000
Output
840

本题就是反素数经常考的一种类型–给定n，求n以内因子数最多的数。
一提到因子数，首先想到的就是把素因子分解。把n分解成
n=（p1的k1次方）*（p2的k2次方）**（p3的k3次方）……

那么因子总数 ans=（k1+1）（k2+1）（k3+1）……

但是我们要枚举到多少次幂呢？（好问题）
考虑极端情况，n是某一个素数的k次幂，long long 为2的64次幂，所以我们最多就考虑到64层！！！（emmm）
那什么时候停止呢？
1.当前因子大于我们想要的因子或正好为我们想要的因子
2.已经大于n

下面是已经AC的代码``

ll primes[16] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
int maxd,maxval;
long long n;
//dp 为当前枚举第几个素数。 num为当前因子的数量
//temp为当前因子数位num的数值。up为上一个素数的幂，这次应该小于等于这个幂次
void dfs(int dp,ll temp,int num,int up)
{
    if((num>maxd)||(num==maxd&&maxval>temp)){
    //当前因子数>maxd or 因子数相同，取值小的
        maxd=num;
        maxval=temp;
    }
    if(dp==16) return;
    for(int i=1;i<=up;i++){
        if(primes[dp]*temp>n) break;
        temp*=primes[dp];
        dfs(dp+1,temp,num*(i+1),i);
    }
}
int main()
{
    cin>>n;
    dfs(0,1,1,60);
    cout<<maxval<<endl;
}