【51nod】最复杂的数 （简单搜索）

基准时间限制：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=p_1^{a_1}p_2^{a_2}...p_k^{a_k}$$
假设p1到pk从小到大，我们观察一下a1到ak的性质，发现：设 $i<j$,若 $p_i<p_j$则必然 $a_i>=aj$,
证明是反证法，假设 $a_i<aj$,那么我们把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;
}