把一个数的约数个数定义为该数的复杂程度,给出一个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
这个题目是ural上的原题!
题目的意思就是给你一个数n, 让你求1~n范围内因数最多的数,当有多个符合的数时取最小的输出,并且输出因数的个数!
要做这道题首先要先认识一个公式,就是一个数k=p1^(a1)*p1^(a1)*......pn^(an) (唯一分解定理),那么这个数因数的个数等于(a1+1)*(a2+1)*......*(an+1)。
因为这个数十分的大,所以不能枚举,所以我们可以用深搜枚举质因数,因为连乘质数2~53时已经大于10^18所以我们只需要搜到53就行。当然,虽然用dfs方便了很多,但是只由2~53质因数相乘得到的10^18以内的数还是很多,所以可以适当的进行剪枝,不然还是会TLE。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const ll inf = 1e18;
ll n, large, record;
ll prime[50], num = 0;
bool vis[107];
void dfs(ll k, ll s, ll step, ll mm){
// mm是剪枝部分,此次质因子枚举次数比上一个质因数要少,就可以减少深搜的次数。
if(step > num) return;
ll temp = k;
for(ll i = 1; i <= mm; i ++){
if(n/prime[step] < temp) break;
temp *= prime[step];
if(record < s*(i+1)){
record = s*(i+1);
large = temp;
}
else if(record == s*(i+1) && large > temp) large = temp;
dfs(temp, s*(i+1), step+1, i);
}
}
int main(){
int t;
memset( vis, true, sizeof( vis));
for(int i = 2; i <= 60; i ++){
if(!vis[i]) continue;
prime[++ num] = i;
for(int j = i<<1; j <= 60; j += i)
vis[j] = false;
}
//printf("%I64d\n", num);
scanf("%d", &t);
while(t --){
scanf("%I64d",& n);
large = 1;
record = 1;
dfs( 1, 1, 1, 64);
printf("%I64d %I64d\n", large, record);
}
return 0;
}