把一个数的约数个数定义为该数的复杂程度,给出一个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
题解:首先考虑约数。因为是约数,容易想到与素数,又因为每个数a =(p1^b1)*(p2^b2)*......(pn^bn)。所以一个数的约数个数为(b1+1)*( b2 + 1)*......(bn+1)。然后又因为前16个素数的乘积大于1e18,故用搜索,直接搜索。当然这里需要剪枝,因为考虑到约数个数尽量大而数的值尽量小,所以在相同约数个数情况下,较小的约数使用得多则数越小。即b1>=b2>=b3......>=bn。
代码:
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <vector>
#include <set>
#include <stack>
#include <map>
#include <cmath>
#include <ctime>
using namespace std;
typedef long long ll;
typedef pair<ll, int> P;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3f;
const double PI = acos(-1.0);
const double eps = 1e-10;
const int maxn = 1e5+7;
const int mod = 1e9+7;
int T;
ll n;
int p[16] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
ll ans1,ans2;
void dfs(ll now,ll cnt,int pos,int up){
if(pos>16) return;
if(cnt>ans2) ans1 = now,ans2 =cnt;
else if(cnt==ans2&&ans1>now) ans1 = now;
for(int i = 1;i <= up;i++){
if(now>n/p[pos]) break;
dfs(now*p[pos],cnt*(i+1),pos+1,i);
now = now*p[pos];
}
}
int main()
{
scanf("%d",&T);
while(T--){
ans1 = 1,ans2 = 1;
scanf("%lld",&n);
dfs(1,1,0,65);
printf("%lld %lld\n",ans1,ans2);
}
return 0;
}
涉及知识点:
反素数:对于任何正整数n,其约数个数记为f(n),例如f(6) = 4,如果某个正整数n满足:对任意的正整数i(0<i<n),都有f(i)<f(n),那么称n为反素数。
通俗的讲,就是一个数如果其约数个数大于小于它的所有数的约数个数,则这个数称为反素数。
反素数的性质:
1、一个反素数的所有质因子必然是从2开始的连续若干个质数,因为反素数是保证约数个数为x的这个数n尽量小。
2、同样的道理,如果,n = (2^b1)*(3^b2)*(5^b3).....那么必有b1>=b2>=b3>=......。
常见的有关问题有:
1、求约数为n的最小正整数x。
2、求n以内约数个数最多的较小正整数x。
这类问题多用搜索解决,如果TLE,考虑性质进行合理剪枝。
本题注意点:
if(now>n/p[pos]) break;
这个语句中的判断语句如果改为now*p[pos]>n的话会错,因为有可能乘法溢出。所以关于乘法,应该多注意,小心爆int,甚至longlong。