牛客练习赛4反素数

链接： https://www.nowcoder.com/acm/contest/82/A
来源：牛客网

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld

题目描述

t次询问，每次给你一个数n，求在[1,n]内约数个数最多的数的约数个数

输入描述:

第一行一个正整数t
之后t行，每行一个正整数n

输出描述:

输出t行，每行一个整数，表示答案

示例1

输入

5
13
9
1
13
16

输出

6
4
1
6
6

备注:

对于100%的数据，t <= 500 , 1 <= n <= 1000000000000000000

题解：反素数的应用。就是找一个数P  = p1^a1*p2^a2*...pk^ak  <= N , 且max {(a1+1) * (a2 + 1) * ... (ak+1)}。我们直接暴力DFS，但有一个强剪枝就是指数a1,a2,...ak非递增。

代码：

#include <stdio.h>
#include <math.h>
#include <algorithm>
#include <queue>
#include <vector>
#include <deque>
#include <stack>
#include <string>
#include <string.h>
#include <bitset>
#include <map>
#include <set>
#include <assert.h>
#include <iostream>
#include <unordered_map>

using namespace std;

typedef long long ll;
typedef long double ld;

#define x0 x0___
#define y0 y0___
#define pb push_back
#define SZ(X) ((int)X.size())
#define mp make_pair
#define fi first
#define se second
#define pii pair<int,int>
#define pll pair<ll,ll>
#define pli pair<ll,int>
#define pil pair<int,ll>
#define ALL(X) X.begin(),X.end()
#define RALL(X) X.rbegin(),X.rend()
#define rep(i,j,k) for(int i = j;i <= k;i ++)
#define per(i,j,k) for(int i = j;i >= k;i --)
#define mem(a,p) memset(a,p,sizeof(a))


const ll MOD = 1E9 + 7;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3f;
ll qmod(ll a,ll b,ll c) {ll res=1;a%=c; assert(b>=0); for(;b;b>>=1){if(b&1)res=res*a%c;a=a*a%c;}return res;}
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;}

template<typename T, typename S>
void upmax(T& a,S b){if(a<b) a=b;}
template<typename T, typename S>
void upmin(T& a,S b){if(a>b) a=b;}
template<typename T>
void W(T m) {cout << m <<endl;}
void gettle() {while(1);}
void getre() {int t=0;t/=t;}








int prime[100];
int tot;
int ans;
ll n;

bool isp(int n)
{
    for(int i = 2;i * i <= n;i ++) if(n % i == 0) return 0; return 1;
}

void dfs(int k, ll num, ll now, int limit)
{
    if(k >= tot) return;
    upmax(ans, num);
    rep (i,1,limit) {
        if(now <= n / prime[k]) dfs(k + 1, num * (i + 1), now *= prime[k], i);
        else break;
    }
}

int main()
{
    rep (i,2,60) {
        if(isp(i)) prime[tot++] = i;
    }
    int T;
    scanf("%d", &T);
    while(T --) {
        ans = 0;
        scanf("%lld", &n);
        dfs(0,1,1,30);
        printf("%lld\n", ans);
    }
    return 0;
}