题目:BZOJ1053.
题目大意:给定一个数N,求一个最大的不超过N的数,使得这个数的约数个数大于所有小于它的自然数的约数个数.
1
≤
n
≤
2
∗
1
0
9
1\leq n\leq 2*10^9
1≤n≤2∗109.
形式化一些,设 g ( x ) = ∑ d ∣ x 1 g(x)=\sum_{d|x}1 g(x)=∑d∣x1,求一个最大的数 x x x满足 x ≤ N x\leq N x≤N且对于 ∀ i < x , g ( x ) > g ( i ) \forall i<x,g(x)>g(i) ∀i<x,g(x)>g(i).
很明显这个数要求最大,那么这个数肯定是 [ 1 , N ] [1,N] [1,N]中约数个数最多的数之一,而且要求比它小的数中约数个数不能有比它大的,所以这个数肯定是所有 [ 1 , N ] [1,N] [1,N]中约数个数最多的数中最小的一个.
有了这个性质后,想到一个数 n n n的质因数个数不超过 log n \log n logn个,所以考虑暴力枚举质因数.又想到约数个数最多的数中最小的一个的每一个质因子的数目中,大的质因子的数量肯定不超过小的质因子的数量.
形式化的,若一个数n表示为唯一分解的形式时:
n
=
∏
i
=
1
k
p
i
c
i
n=\prod_{i=1}^{k}p_i^{c_i}
n=i=1∏kpici
那么若 p x < p y p_x<p_y px<py,则 c x ≥ c y c_x\geq c_y cx≥cy.
想到这里的时候再考虑一下一个数 x x x的质因子数量的上限?最多不会超过 log x \log x logx个.也就是说最后拆出来的数质因子分解后最多不超过 30 30 30个,而且不相同的质因子数量最多只有 10 10 10个.这表明我们可以考虑搜索,然后用上面的两条性质来优化了,事实证明这样子的时间复杂度可过.
至于具体怎么搜…其实只需要大力dfs每一个数出现的次数就可以了,然后每次搜的时候顺带上两个优化,在dfs的同时顺带求解一下约数个数,最后不断更新就好了.
代码如下:
#include<bits/stdc++.h>
using namespace std;
#define Abigail inline void
typedef long long LL;
const int N=10;
const int pr[N+9]={0,2,3,5,7,11,13,17,19,23,29};
int n,cnt[N+9],ans=1,num=1;
void dfs(int k,LL a,int c,int last){ //a有可能爆int,莫名80分就是因为这个
if (k>N){
if (a>ans&&c>num){ans=a,num=c;}
if (a<=ans&&c>=num){ans=a;num=c;}
return;
}
int p=1;
for (int i=0;i<=last;++i){
dfs(k+1,a*p,c*(i+1),i);
p=p*pr[k];
if (a*p>n) break;
}
}
Abigail into(){
scanf("%d",&n);
}
Abigail work(){
dfs(1,1,1,30);
}
Abigail outo(){
printf("%d\n",ans);
}
int main(){
into();
work();
outo();
return 0;
}