题意:
给出nnn,求出
max(d(1),d(2),...,d(n))
max(d(1),d(2),...,d(n))
max(d(1),d(2),...,d(n))
Solution:
由约数个数定理,只需要枚举质因数pip_{i}pi和他的幂次aia_{i}ai,使得乘积≤n\leq n≤n即可,此时约数个数即
(a1+1)(a2+1)...(an+1)
(a_{1}+1)(a_{2}+1)...(a_{n}+1)
(a1+1)(a2+1)...(an+1)
dfs枚举即可,不过需要剪枝
-
从小到大枚举素数,如果两个素数p,qp,qp,q满足p<qp<qp<q,我们分配给qqq有kkk次幂,不如让ppp多kkk次幂,这样可以使得原数更大,因为两数之积尽量相同积会大些,并且此时乘积也会小一些,能分配更多幂给其他的素数,所以一个优化是分配给当前素数的幂不超过上一个素数的幂
-
当上一次分配的幂是000是,接下来就不分配了,直接返回
只要有上述两个剪枝,就算筛出10710^7107次方之内的素数来枚举都不会改变复杂度
// #include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#include<bitset>
#include<map>
using namespace std;
using ll=long long;
const int N=1e7+5,inf=0x3fffffff;
const long long INF=0x3fffffffffffffff,mod=1e9+7;
ll cnt,prime[N];
bitset<N>nt;
void make_prime()
{
for(int i=2;i<N;i++)
{
if(!nt[i]) prime[++cnt]=i;
for(int j=1;j<=cnt&&i*prime[j]<N;j++)
{
nt[i*prime[j]]=true;
if(i%prime[j]==0) break;
}
}
}
ll dfs(int x,ll now,ll ans,ll last)
{
if(x==cnt+1||last==0) return ans;
ll tmp=now,tot=0,ret=-INF;
while(tmp&&tot<=last)
{
ret=max(ret,dfs(x+1,tmp,ans*(tot+1),tot));
tmp/=prime[x]; tot++;
}
return ret;
}
int main()
{
#ifdef stdjudge
freopen("in.txt","r",stdin);
#endif
make_prime();
int t; cin>>t;
while(t--)
{
ll n; cin>>n;
cout<<dfs(1,n,1,INF)<<'\n';
}
return 0;
}
598



被折叠的 条评论
为什么被折叠?



