洛谷 p1463 反素数 （数学方法）

题目链接

算法竞赛进阶指南P138

若n=p1^q1*p2^q2*...pn^qn  (pi是质数）,则n的因子个数为  （q1+1)*(q2+1)*...(qn+1),

1.由题意得，求n以内因子数最多的数且因子数相等时，取最小的数，

2.因为前11个素数的乘积>2e9,所以只用到前十个素数，因为2^31>2e9，所以指数和最多为30

3.由因子数的公式得，若公式一定，要让n最小，则按素数从小到大，他们的指数个数应该从大到小，且素数应为从2开始的连续的几个素数，

因此，可以用dfs搜索，从2开始确定他们的指数，且他们的指数非递增

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
#include <map>
using namespace std;
typedef long long ll;

int a[15]={0,2,3,5,7,11,13,17,19,23,29};
map<ll,ll>mp;
map<ll,ll>::iterator it;
void dfs(ll x,ll y,ll k,ll t){
    if(!y) return;
    if(x==11) return;
    if(k>2e9) return;
    ll ans=1;
    for(ll i=0;i<=y;i++){
        dfs(x+1,i,k*ans,t*(i+1));
        if(!mp[t*(i+1)]) mp[t*(i+1)]=k*ans;
        else mp[t*(i+1)]=min(k*ans,mp[t*(i+1)]);
        ans*=a[x];
        if(k*ans>2e9) break;
    }
}

int main()
{
    mp[1]=1;
    dfs(1,30,1,1);
    ll n;
    cin>>n;
    ll z=0,x=0;
    for(it=mp.begin();it!=mp.end();it++){
        if(it->second<=n){
            if(it->first>x || (it->first==x&&it->second<z)){
                x=it->first;
                z=it->second;
            }
        }
    }
    cout<<z<<endl;
    return 0;
}