题目:洛谷P1463、BZOJ1053、Vijos P1172、codevs2912。
题目大意:对于任何正整数x,其约数的个数记作g(x)。例如g(1)=1、g(6)=4。
如果某个正整数x满足:g(x)>g(i) 0<i<x,则称x为反质数。例如,整数1,2,4,6等都是反质数。
现在给定一个数N(≤2,000,000,000),求出不超过N的最大的反质数。
解题思路:N那么大,第一反应是找规律。
首先得打出一个范围较小的表(三秒内能打完的),发现:范围大,但反素数很少。
于是果断选择了打表大法!
不过打表也得讲究技巧。首先发现所有反素数(除1以外)都是偶数,那么在循环判断时只需循环偶数即可。
其次发现从60开始,每个反素数都是20的倍数,那么此时每次加20即可。
但二十亿个数,还是很多。怎么办?我们发现,两个反素数直接总是相差很多,那我们每求出一个反素数,就跳过一部分,使得减少循环次数。
最终,成功打出表,见代码(最后那个$2^{32}-1$是占位用的)。
C++ Code:
#include<cstdio> const int a[]={1,2,4,6,12,24,36,48,60,120,180,240,360,720,840,1260,1680,2520, 5040,7560,10080,15120,20160,25200,27720,45360,50400,55440,83160,110880,166320, 221760,277200,332640,498960,554400,665280,720720,1081080,1441440,2162160,2882880, 3603600,4324320,6486480,7207200,8648640,10810800,14414400,18378360,21621600,32432400, 36756720,43243200,61261200,73513440,110270160,122522400,147026880,183783600,245044800, 328648320,410810400,551350800,698377680,735134400,1102701600,1396755360,2147483647}; int main(){ int n; scanf("%d",&n); for(int i=0;;++i) if(a[i]>n){ printf("%d ",a[i-1]); return 0; } }
下附打表代码,在我家里的渣机上跑了33分02秒279(O2优化下),不过总归是能打完的。
#include<cstdio> #define reg register #define ull unsigned long long int main(){ freopen("output.txt","w",stdout); printf("1,"); reg int t=2; reg ull max=1; for(reg ull i=2;i<=2000000200;i+=t){ if(i==60)t=20; reg ull p=1,k=i; for(reg int j=2;j*j<=k;++j) if(k%j==0){ reg ull f=1; for(;k%j==0;k/=j,++f); p*=f; } if(k>1)p<<=1; if(p>max){ max=p; printf("%llu,",i); if(i>=800000000)i+=100000000;else if(i>=410810400)i+=30000000;else if(i==328648320)i=410810380;else if(i>=200000000)i+=50000000;else if(i==147026880)i+=35000000;else if(i>=122522400)i+=20000000;else if(i>=43243200)i+=10000000;else if(i>=21621600)i+=4000000;else if(i>=10000000)i+=3000000;else if(i>=1400000)i+=600000;else if(i>=1000000)i+=300000;else if(i>=100000)i+=50000; } } return 0; }