zoj 2562 More Divisors(反素数)

原创 2013年12月02日 20:50:18
More Divisors

Time Limit: 2 Seconds      Memory Limit: 65536 KB

Everybody knows that we use decimal notation, i.e. the base of our notation is 10. Historians say that it is so because men have ten fingers. Maybe they are right. However, this is often not very convenient, ten has only four divisors -- 1, 2, 5 and 10. Thus, fractions like 1/3, 1/4 or 1/6 have inconvenient decimal representation. In this sense the notation with base 12, 24, or even 60 would be much more convenient.

The main reason for it is that the number of divisors of these numbers is much greater -- 6, 8 and 12 respectively. A good quiestion is: what is the number not exceeding n that has the greatest possible number of divisors? This is the question you have to answer.

Input:

The input consists of several test cases, each test case contains a integer n (1 <= n <= 1016).

Output:

For each test case, output positive integer number that does not exceed n and has the greatest possible number of divisors in a line. If there are several such numbers, output the smallest one.

Sample Input:
10
20
100
Sample Output:
6
12
60


题意:给出一个数n,求n以内的一个因子个数最多的且大小最小的数。

思路:

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

如果某个正整数x满足:对于任意i(0<i<x),都有g(i)<g(x),则称x为反素数.

题目转化为给一个N,求出不超过N的最大的反素数.

比如:输入1000 输出 840


求[1..N]中约数在大的反素数-->求约数最多的数

如果求约数的个数 756=2^2*3^3*7^1

(2+1)*(3+1)*(1+1)=24

基于上述结论,给出算法:按照质因数大小递增顺序搜索每一个质因子,枚举每一个质因子

为了剪枝:

性质一:一个反素数的质因子必然是从2开始连续的质数.

10^16以内的数最多只需要前16个素数构造

性质二:p=2^t1*3^t2*5^t3*7^t4.....必然t1>=t2>=t3>=....


AC代码:

#include <iostream>
#include <cmath>
#include <cstring>
#include <cstdio>
#include <queue>
#define ll long long

using namespace std;

const int prime[16] = {1,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
ll n, minnum, maxsum;

void dfs(ll num, ll sum, ll k, int cnt)
{
    if(sum > maxsum){
        minnum = num;
        maxsum = sum;
    }
    else if(sum == maxsum && num < minnum) minnum = num;
    if(k > 15) return;
    ll tmp = num;
    for(int i = 1; i <= cnt; i++)
    {
        tmp *= prime[k];
        if(tmp > n) break;
        dfs(tmp, sum * (i + 1), k + 1, i);
    }
}
int main()
{
    while(~scanf("%lld", &n))
    {
        minnum = n;
        maxsum = 0;
        dfs(1, 1, 1, 50);
        printf("%lld\n", minnum);
    }
    return 0;
}



版权声明:本文为博主原创文章,未经博主允许不得转载。

ZOJ 2562 More Divisors 解题报告(反素数)

More Divisors Time Limit: 2 Seconds      Memory Limit: 65536 KB Everybody knows that we use de...
  • kbdwo
  • kbdwo
  • 2014年04月21日 20:46
  • 449

反素数深度分析

今天要我要讲的是反素数,在ACM中也算是常见的考点,那么对于搞ACM的同学来说,很有必要搞清楚它,所以接下 来我会很详细地讲解。   在讲解反素数之前,我们先来看反素数的概念。   反素数的定义:对于...
  • ACdreamers
  • ACdreamers
  • 2014年05月05日 20:06
  • 10569

ACM里的反素数问题

定义 对于正整数 x ,其约数的个数记做 g(x) 。例如 g(1) = 1,g(6) = 4. 如果某个正整数x满足: 对于任意 i(0...
  • qq_15714857
  • qq_15714857
  • 2015年10月05日 22:48
  • 394

【HDU4392】【反素数强大的模版 java或者C++】

传送门:http://acm.hdu.edu.cn/showproblem.php?pid=4392 JAVA版本: import java.io.BufferedInputStream; im...
  • guhaiteng
  • guhaiteng
  • 2016年10月29日 21:19
  • 262

[bzoj 1053] [HAOI2007]反素数ant:数论,DAG上最短路

题意:设g(i) = i的约数个数,若对于某个正整数x,所有0
  • ruoruo_cheng
  • ruoruo_cheng
  • 2016年12月04日 14:52
  • 243

zoj 1562 反素数 附上个人对反素数性质的证明

反素数的定义:对于任何正整数,其约数个数记为,例如,如果某个正整数满足:对任意的正整             数,都有,那么称为反素数。 从反素数的定义中可以看出两个性质: (1)一个反素数的所...
  • u011026968
  • u011026968
  • 2014年08月24日 00:46
  • 1425

bzoj3085: 反质数加强版SAPGAP

传送门 没看过1053的走这儿 剪枝大法好。 1.前一个质数个数必定大于后一个质数个数。 2.前一个质数的幂必定小于后一个质数的幂。 在加上一些奇奇怪怪的乱搞就可以了。 他比我讲得好#in...
  • zhouyuyang233
  • zhouyuyang233
  • 2017年03月09日 15:50
  • 293

素数系列——反素数

Emirp Time Limit: 5000 MS Memory Limit: 100000 K Total Submit: 147(60...
  • sexgeek
  • sexgeek
  • 2017年01月12日 16:08
  • 135

提高篇第一讲项目6.4——编制函数输出1000以内的所有素数,回文数,回文素数,可逆素数

(4)编制main函数,调用上面定义的3个函数,完成 输出1000以内的所有素数。 输出1000以内的所有回文数。 输出1000以内的所有回文素数。 若一个素数的反序数仍为素数,则称它为可逆素数。求1...
  • xiaochenxiaoren
  • xiaochenxiaoren
  • 2017年03月08日 23:41
  • 581

ZOJ - 3496 Assignment 上下界最大流+二分

题目链接点这里 毒瘤题,,wa了无数发,,,才发现当m==0的时候,,会输出一些奇怪的东西。。。 以后,,不管什么题都要记得特判,,恩,,,不然会死的很惨。。 #include #inc...
  • qq_30927651
  • qq_30927651
  • 2017年03月31日 23:04
  • 193
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:zoj 2562 More Divisors(反素数)
举报原因:
原因补充:

(最多只允许输入30个字)