反素数ant

最新推荐文章于 2019-08-14 21:07:41 发布

RevolIA

最新推荐文章于 2019-08-14 21:07:41 发布

阅读量212

点赞数

本文链接：https://blog.csdn.net/Du_Mingm/article/details/84797501

版权

问题 B: 反素数ant

时间限制: 1 Sec 内存限制: 128 MB

题目描述

对于任何正整数x，其约数的个数记作g(x)。例如g(1)=1、g(6)=4。如果某个正整数x满足：g(x)>g(i) 0<i<x，则称x为反质数。例如，整数1，2，4，6等都是反质数。现在给定一个数N，你能求出不超过N的最大的反质数么？

输入

一个数N（1<=N<=2,000,000,000）。

输出

不超过N的最大的反素数。

样例输入

样例输出

唯一分解得到，不同的素因子超不过12个，（前12个素数乘起来大于2e9）

一个数的约数的个数，可以这样算:

设 $x=a_{1}^{p_{1}}*a_{2}^{p_{2}}*a_{3}^{p_{3}}****a_{n}^{p_{n}}$

那么 $Fun(x)=\prod_{i=1}^{n}(p_{i}+1)$

然后暴力跑dfs，一个优化就是2*2*5和2*5*5，显然约数相等，但是明显的前面的小，

所以，优化就是后面的数（数值大）的幂要小于前面的数的幂

代码来源：我是超链接、

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll prim[] = {0,2,3,5,7,11,13,17,19,23,29,31,37};
ll Cot[100];
ll n,m,ans,Mx;
void dfs(int x,ll cnt,ll mult){///第x个素数，约数个数，值
    if(x>12)return;
    if(cnt>Mx || cnt == Mx && mult < ans){
        Mx = cnt;
        ans = mult;
    }
    Cot[x] = 0;
    while(mult*prim[x]<=n && Cot[x] < Cot[x-1]){
        Cot[x]++;
        mult *= prim[x];
        dfs(x+1,cnt*(Cot[x]+1),mult);
    }
}
int main(){
    cin>>n;
    Cot[0] = 1e9;
    dfs(1,1,1);
    cout<<ans<<endl;
    return 0;
}