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

zoj2562:搜索+数论(反素数)

时间:2014-09-09 22:45:29      阅读:298      评论:0      收藏:0      [点我收藏+]

标签:style   blog   color   os   io   for   div   sp   代码   

题目大意:求n以内因子数量最多的数  n的范围为1e16

其实相当于求n以内最大的反素数。。。

由素数中的 算数基本原理

设d(a)为a的正因子的个数,则

              d(n)=(a1+1)(a2+1).....*(an+1);

又由反素数的性质2:

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

我们可以指定搜索策略和剪枝

详见代码和注释

#include <iostream>
#include<stdio.h>
#include<algorithm>
#include<math.h>
using namespace std;
long long n;
long long a[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};//这些素数全部乘积已经超过了1e16,所以不用再往下面乘了
long long ans,ans0;
void dfs(long long num,long long sum,long long limit,long long po)//limit 存储上一个素数的指数 相当于反素数性质中的t(i-1),由性质可知ti<=t(i-1)
{
    if(num>n)
        return;
    if(sum==ans0)
    {
        ans=min(ans,num);
    }
    if(sum>ans0)
    {
        ans0=sum;
        ans=num;
    }
    long long p=1;
    for(int i=1;i<=limit;i++)
    {
        p*=a[po];
        if(p*num>n)
            break;
        dfs(num*p,sum*(i+1),i,po+1);
    }
    return;
}
int main()
{
    while(scanf("%lld",&n)!=EOF)
    {
        ans0=0;
        ans=0;
        if(n==1)
        {
            puts("1");
            continue;
        }
        dfs(1,1,100,0);
        cout<<ans<<endl;
    }
    return 0;
}

 

zoj2562:搜索+数论(反素数)

标签:style   blog   color   os   io   for   div   sp   代码   

原文地址:http://www.cnblogs.com/lnever/p/3963371.html

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