198. 反素数(dfs)

题目链接:https://www.acwing.com/problem/content/description/200/
知识共享许可协议
本作品采用知识共享署名-相同方式共享 4.0 国际许可协议进行许可。

对于任何正整数x,其约数的个数记作g(x),例如g(1)=1、g(6)=4。

如果某个正整数x满足:对于任意的小于x的正整数 i,都有g(x)>g(i) ,则称x为反素数。

例如,整数1,2,4,6等都是反素数。

现在给定一个数N,请求出不超过N的最大的反素数。
输入格式

一个正整数N。
输出格式

一个整数,表示不超过N的最大反素数。
数据范围

1≤N≤2∗109

输入样例:

1000

输出样例:

840

分析:

要不是有书,给我三十年,我也做不出这道题。 好菜呀。 关键有书我还能wa。

讲一下自己的理解吧。
明确几点概念:
正整数N被唯一分解为N=p1^c1* p2^2 * ~~pn ^cn
我们可得:
此正整数约数的个数为:
(c1 + 1) * (c2 +1) ~ ~(cn + 1);

第一反应求出1到N里面所有数的约数个数,并从左到右进行扫描找到最大的约数个数,如果最大数有多个,那应该是数小的那一个。到此方法必超时。
但可得引理一:
1到N中最大的反质数,就是1到N中约数个数最多的数中的最小一个数。
引理二:
看N的范围为1到2e9.那么1到N中任何数的不同质因子个数一定不会超过十个
eg:
2 * 3 * 5 * 7 * 11 * 13* 17 * 19 * 23 * 29 * 31 > 2e9.
那么同理,我们可的最高次幂不会超过31
eg:
2^31 > 2e9
引理三:
我们要求的那个数势必是下面形式:
2^c1 * 3^c2 * ~~* 31^cn。并且c1>= 2 >= cn
证明:不是非上升幂。那么我们通过交换质因子的方法,可以得到一个比他小的数。
最后在跑2e9的时候超时了。原因归咎是sum爆了longlong,成为了负数,所以特别加了一个判定是否为负数的情况。

#include"string.h"
#include"math.h"
#include"algorithm"
using namespace std;
typedef long long ll;

ll N,m,max_cnt;
ll prime[12] = {0,2,3,5,7,11,13,17,19,23,29,31};

void dfs(int id,int x,ll sum,ll cnt)
{
    if(sum > N || sum < 0) return ;
    //printf("sum = %lld m = %lld\n",sum,m);
    if(x == 0 || id > 11)
    {
        if(cnt > max_cnt)
        {
            m = sum; max_cnt = cnt;
        }
        else if(cnt == max_cnt)
            m = min(m,sum);
        return ;
    }

    dfs(id,x - 1,sum,cnt);
    dfs(id + 1,x,sum * (ll)pow(prime[id],x),cnt * (x + 1));
}

int main()
{
    scanf("%lld",&N);
    m = N + 1;
    max_cnt = 1;
    dfs(1,30,1LL,1LL);
    printf("%lld\n",m);
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值