反素数
#define Y(i) i的因数的个数
我们把n称为反素数如果 Y(n)>Y(∀i)(i<n)Y(n)>Y(∀i)(i<n)
为了让n满足这个式子,我们需要让n的因子数多的前提尽量让n小,让n的因子分布尽量密集
所以,有以下两个结论 :
n的因子从小开始且尽可能的连续(2,3,5,7,11…)
如果可以让一个因子多重复出现几次,那么一定会选择小的因子
即如果n=2q1∗2q2∗2q3∗...n=2q1∗2q2∗2q3∗...,那么q1>=q2>=q3...q1>=q2>=q3...
你会发现你遇到的题目和上面的那句话很相似
Y(i)=n,求min(i)
i<=n,求max(Y(i))
算法思路
n=pq11∗pq22...pqkkn=p1q1∗p2q2...pkqk
n的因子数为(q1+1)∗(q2+1)∗...∗(qk+1)(q1+1)∗(q2+1)∗...∗(qk+1)
假设现在求的是n个因数的最小数
建一棵树,除root所在的层以外,每一层都由一个素数和其i次幂构成
我们暴力跑一遍dfs,跑到一个节点就乘上一个结点上的数
跑的时候记录的有 : 当前因子的个数 , 累乘后数的大小 , 跑了几层(几个素数)
因子的个数可以用上面注释部分的公式来推,假设之前已经有jj个了,现在乘了ikik,因子数便是j∗(k+1)
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
LL pri[16] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
LL ans,n,MAXY;
//求n以内因子最多的数,多个答案输出小的
void dfs(LL Y,LL dep,LL num,LL limit){
//因子数 深度 当前的数
if(dep==16)return;
//if(Y>=n)return;
if(Y>MAXY)MAXY=Y,ans=num;
if(Y==MAXY)ans=min(ans,num);
for(int i=1;i<=limit;i++){
//因为后面的次数不会超过前面,所以这里添加上限就是一个极大的剪枝
if(n/pri[dep]<num)break;//保证不超过n
dfs(Y*(i+1),dep+1,num*=pri[dep],i);
}
}
int main(){
int t;cin>>t;
while(t--){
scanf("%lld",&n);
ans=1e18;
MAXY=0;
dfs(1,0,1,60);
printf("%lld %lld\n",ans,MAXY);
}
}