反素数
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;
}