zoj2562（反素数）

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2562

 

这题很有意思，思路来源于

http://hi.baidu.com/speakless/blog/item/1be8d61b63711cd5ac6e751f.html

对于任何正整数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个素数构造:2,3,5,7,11,13,17,19,23,29

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

#include<iostream> using namespace std; const int prime[16]={1,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47}; __int64 maxsum,bestnum, n; void make(__int64 num, __int64 k,__int64 sum,int limit) {//num：当前枚举到的数，k：枚举到的第k大的质因子；sum：该数的约数个数；limit：质因子个数上限； int i; __int64 temp; if(sum > maxsum) { maxsum = sum; bestnum = num; //如果约数个数更多，将最优解更新为当前数； } if(sum==maxsum && bestnum > num) bestnum = num; //如果约数个数相同，将最优解更新为较小的数； if(k > 15) return; temp = num; for(i=1; i<=limit; i++) //开始枚举每个质因子的个数； { if(temp*prime[k] > n) break; temp = temp * prime[k]; //累乘到当前数； make(temp, k+1, sum*(i+1), i); //继续下一步搜索； } } int main() { while(scanf("%I64d",&n) != EOF) { maxsum == 0; bestnum = n; make(1,1,1,50); printf("%I64d/n",bestnum); } return 0; }