bzoj1053 [HAOI2007]反素数ant（数论+搜索）

求一段区间内最大的反质数，其实就是求一段区间内约数个数最大的数，如果有多个取最小的。

我们知道
$x=\sum\limits_{i=1}^kp_i^{a_i}$
则约数个数为$\prod\limits_{i=1}^k(a_i+1)$
因此我们不难发现我们一定有$a_1>=a_2>=...>=a_k$
而且我们发现前10个质数的乘积就已经超过2e9了，因此我们爆搜前十个质数即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 13
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,ans=0,mx=0,prime[11]={0,2,3,5,7,11,13,17,19,23,29};
void dfs(int id,int res,int cnt,int k){
    if(id>10){
        if(cnt>mx) mx=cnt,ans=res;
        else if(cnt==mx) ans=min(ans,res);return;
    }ll tmp=res;
    for(int i=0;i<=k;++i){
        dfs(id+1,tmp,cnt*(i+1),i);
        tmp*=prime[id];if(tmp>n) break;
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();dfs(1,1,1,30);printf("%d\n",ans);
    return 0;
}