基准时间限制:1 秒 空间限制:131072 KB 分值: 40
把一个数的约数个数定义为该数的复杂程度,给出一个n,求1-n中复杂程度最高的那个数。
例如:12的约数为:1 2 3 4 6 12,共6个数,所以12的复杂程度是6。如果有多个数复杂度相等,输出最小的。
Input
第1行:一个数T,表示后面用作输入测试的数的数量。(1 <= T <= 100)
第2 - T + 1行:T个数,表示需要计算的n。(1 <= n <= 10^18)
Output
共T行,每行2个数用空格分开,第1个数是答案,第2个数是约数的数量。
Input示例
5
1
10
100
1000
10000
Output示例
1 1
6 4
60 12
840 32
7560 64
思路:假设我们知道在范围中x是我们要的答案,那么我们唯一分解定理展开一下:
x=pa11pa22...pakk
x
=
p
1
a
1
p
2
a
2
.
.
.
p
k
a
k
假设p1到pk从小到大,我们观察一下a1到ak的性质,发现:设 i<j i < j ,若 pi<pj p i < p j 则必然 ai>=aj a i >= a j ,
证明是反证法,假设 ai<aj a i < a j ,那么我们把ai和aj对调,对因子总数不会发生变化,但是x的大小反而会跟小,也就是说答案不会更差,得证。算是贪心的思想。
代码:
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define ull unsigned long long
#define ll long long
#define ul unsigned int
#define maxn 60
using namespace std;
int prime[maxn];
bool isP[maxn];
int cnt=0;
void init()
{
for(int i=2;i<maxn;i++)
{
if(!isP[i]){prime[cnt++]=i;}
for(int j=0;j<cnt&&i*prime[j]<maxn;j++)
{
isP[i*prime[j]]=true;
if(i%prime[j]==0)break;
}
}
}
ll n;
ll ans;
ll _max;
void dfs(int cur,int limit,ll temp,ll d)
{
if(temp>n)return ;
if(_max<d)
{
_max=d;
ans=temp;
}
if(_max==d&&ans>temp)
{
_max=d;
ans=temp;
}
ll now=temp;
for(int i=1;i<=limit;i++)
{
if(prime[cur]>n/now)return;
now*=prime[cur];
//cout<<now<<endl;
dfs(cur+1,i,now,d*(i+1));
}
}
int main()
{
init();
int t;
cin>>t;
while(t--)
{
cin>>n;
_max=1;
ans=1;
dfs(0,70,1,1);
cout<<ans<<" "<<_max<<endl;
}
return 0;
}