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

[BZOJ 1053] [HAOI 2007]反素数ant

时间:2015-05-24 23:06:57      阅读:222      评论:0      收藏:0      [点我收藏+]

标签:

1053: [HAOI2007]反素数ant

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1857  Solved: 1034
[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
 
【题解】
我一个蒟蒻怎么能自己想出来这题呢……不过总算理解了……
我还有好多数论知识没有掌握啊TAT
===========================================
首先,可以用组合数学证得,如果一个数n=a1x1×a2x2×……×akxk
那么,设函数ex-phi(n)表示n的约数的个数
可推导出ex-phi(n)=(x1+1)(x2+1)...(xk+1) 
上面的结论称为定理1
===========================================
通过计算可以得出当N在2000000000以内时,最多只有10个素因子
证明:2*3*5*7*11*13*17*19*23*29≈60e>20e
上面的结论称为定理2
===========================================
小素因子多一定比大素因子多要优秀
小素因子多那么总因子多,ex-phi(n)肯定比大素因子多的大。
上面的结论称为定理3
===========================================
然而推出这些神奇的结论后,我们发现——
    我们维护素因子从小到大数量的单调递减性即可。
然后就AC了……
 
看了黄学长的博客果断明白了这题……然而我觉得我这种蒟蒻应该多写点来锻炼。
贴上代码
技术分享
 1 #include<bits/stdc++.h>
 2 #define ll long long
 3 using namespace std;
 4 int n;
 5 int prime[11]={0,2,3,5,7,11,13,17,19,23,29};
 6 int ans=-1,ysans;
 7 void tryy(int zs,ll sum,int dq,int sy,int res) {
 8     if (ysans == res*(dq+1) && sum<ans) ans=sum;
 9     if (res*(dq+1) > ysans) {ysans = res*(dq+1); ans=sum; }
10     if (dq+1<=sy && sum*prime[zs]<=n)    
11         tryy(zs,sum*prime[zs],dq+1,sy,res);
12     for (int i=zs+1;i<=10;++i) 
13         if(sum*prime[i]<=n) tryy(i,sum*prime[i],1,dq,res*(dq+1));
14 }
15 int main() {
16     scanf("%d",&n);
17     tryy(1,1,0,2147483000,1);
18     printf("%lld\n",ans);
19     return 0;
20 }
View Code

那个2147483000是我口胡的……因为你只要维护单调递减,而第一个数是什么并不影响。

可以证明,只需大于28即可……(有可能是错的,因为我看有人写20也A了。。)

 

[BZOJ 1053] [HAOI 2007]反素数ant

标签:

原文地址:http://www.cnblogs.com/TonyNeal/p/bzoj1053.html

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