码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ 1053】 [HAOI2007]反素数ant

时间:2015-08-27 14:47:26      阅读:144      评论:0      收藏:0      [点我收藏+]

标签:

1053: [HAOI2007]反素数ant

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1991  Solved: 1118
[Submit][Status][Discuss]

Description

 

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

Input

一个数N(1<=N<=2,000,000,000)。

Output

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

Sample Input

1000

Sample Output

840

HINT

 

呵呵,刚看到数据范围的时候我也是惊呆了。

但是到最后竟然是一道搜索题。

如果有点组合数学知识的人都知道,如果知道了一个数的素数分解式,那么就可以把每一个素数幂的指数加上1然后全部乘起来。这个我不多说了,想一想都会知道。

根据观察,我们得出一个结论:假如因数的数量相同的情况下,肯定数越小越小好。

也就是说,我们只需使用前几个素数就行了(动动脑筋~~~)

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

int fact[20];
int len;
long long n, last, ans;

bool is_prime(int m)
{
    for (int i = 2; i * i <= m; ++i)
        if (m % i == 0) return false;
    return true;
}

void init()    /* 把在范围内的素数选出来 */
{
    fact[len = 0] = 1;
    long long sum = 1;
    for (int i = 2;; ++i)
        if (is_prime(i))
        {
            if (sum * i > n) return;
            fact[len++] = i;
            sum *= i;
        }
}

void dfs(int s, long long sum, long long f)
{
    if (last < f) { ans = sum; last = f; }    /* 如果因数比目前的答案多就更新 */
    if (last == f) ans = min(ans, sum);    /* 相等的情况下肯定是越小越好 */
    if (s == len) return;
    long long nowp = 1;
    for (int i = 1;;++i)
    {
        dfs(s + 1, sum * nowp, f * i);
        if (sum * nowp * fact[s] > n) return;
        nowp *= fact[s];
    }
}

int main()
{
    cin >> n;
    init();
    ans = 0;
    last = 0;
    dfs(0, 1, 1);
    cout << ans << endl;
    return 0;
}

 

【BZOJ 1053】 [HAOI2007]反素数ant

标签:

原文地址:http://www.cnblogs.com/albert7xie/p/4763084.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!