Zoj 2562 More Divisors

来源:互联网 发布:淘宝上衣服尺寸怎么看 编辑:程序博客网 时间:2024/06/05 10:53

题意:给定一个数n,求在n以内,一个数的约数最多,并且是最小的一个数。

分析:先把约数看成素因子了,导致算法和思路严重错误。

约数和素因子是不同的。比如12=2^2*3,它的素因子数为2,而约数为6,分别为:1,2,3,4,6,12.

不过,可以总结出约数的个数等于每个数的素因子的幂加1相乘。12=2^2*3.  m=(2+1)*(1+1)=3;

其实此题是一个关于反素数的问题,有关反素数:网上有这样的说明:

定义

对于任何正整数x,其约数的个数记做g(x).例如g(1)=1,g(6)=4.如果某个正整数x满足:对于任意i(0<i<x),都有g(i)<g(x),则称x为反素数.

性质

性质一:一个反素数的质因子必然是从2开始连续的质数.
性质二:p=2^t1*3^t2*5^t3*7^t4.....必然t1>=t2>=t3>=....
载自百度.....!
根据分析:题目只有10^16,故素因子越小的会成为最佳结果,故大楷只要15素数就可以了。
#include<cstdio>#include<cstring>typedef long long LL;using namespace std ;LL Max,best,n;int prime[15]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};void dfs(LL sum,LL cur,int k,int t){if(Max<sum){ Max=sum; best=cur;}if(sum==Max&&cur<best) best=cur;if(k>14) return;LL tmp=cur;for(int i=1;i<=t;i++){if(tmp*prime[k]>n) break;tmp*=prime[k];dfs(sum*(i+1),tmp,k+1,i);}}int main(){   while(scanf("%lld",&n)!=EOF){  Max=best=1;  dfs(1,1,0,50);  printf("%lld\n",best);   } return 0 ;}