蓝书（算法竞赛进阶指南）刷题记录——【HAOI2007】BZOJ1053反素数（数论+搜索）

题目：BZOJ1053.
题目大意：给定一个数N，求一个最大的不超过N的数，使得这个数的约数个数大于所有小于它的自然数的约数个数.
$1\le n\le 2\ast 10^9$.

形式化一些，设$g(x)={\sum }_{d|x}1$，求一个最大的数$x$满足$x\le N$且对于$\forall i<x,g(x)>g(i)$.

很明显这个数要求最大，那么这个数肯定是$[1,N]$中约数个数最多的数之一，而且要求比它小的数中约数个数不能有比它大的，所以这个数肯定是所有$[1,N]$中约数个数最多的数中最小的一个.

有了这个性质后，想到一个数$n$的质因数个数不超过$\log n$个，所以考虑暴力枚举质因数.又想到约数个数最多的数中最小的一个的每一个质因子的数目中，大的质因子的数量肯定不超过小的质因子的数量.

形式化的，若一个数n表示为唯一分解的形式时：
$$n=\prod _{i=1}^k{p}_i^{c_i}$$

那么若$p_x<p_y$，则$c_x\ge c_y$.

想到这里的时候再考虑一下一个数$x$的质因子数量的上限？最多不会超过$\log x$个.也就是说最后拆出来的数质因子分解后最多不超过$30$个，而且不相同的质因子数量最多只有$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;
}