题目传送门:http://acm.hdu.edu.cn/showproblem.php?pid=5878
题意:输入一个n,要求找出比n大的,因子只有2,3,5,7的最小的数。
题解:首先说一下什么是丑数,因子只有2,3,5,7的数就是丑数。。。。然后,该怎么做呢?肯定不能从n开始,一个一个暴力找,从题目的样例来看,当n是10的9次方时,要找两万多遍,会超时。如果做过丑数相关的题的话,就会知道其实10的9次方个数里面,丑数其实不多,可以直接打表,然后通过二分法把离n最近的丑数找出来。怎么打表就是关键了。
首先,第一个丑数为“1”,后面的每一个丑数都是有前一个丑数乘2、3、5或7而来,那么后一个丑数就是前一个乘这四个数得到的最小值,for example:第一个:1,第二个:1*2、1*3、1*5或1*7,显然为2,第三个:2*2,1*3,1*5或1*7,显然是3,第四个:2*2,,2*3,1*5,1*7为4,第五个:3*2,2*3,1*5,1*7……
以下代码枚举了5842个丑数,从小到大存在a数组里面,输出一下第5842个丑数,发现已经是2*10的9次方了
#include <iostream>
#include <cstdio>
#define min(a,b) ((a)<(b)?(a):(b))
#define min4(a,b,c,d) min(min(a,b),min(c,d))
using namespace std;
int a[5850];
int main()
{
int n=1;
int p2,p3,p5,p7;
p2=p3=p5=p7=1;
a[1]=1;
while(n<5843){
a[++n]=min4(2*a[p2],3*a[p3],5*a[p5],7*a[p7]);
if(a[n]==2*a[p2])p2++;
if(a[n]==3*a[p3])p3++;
if(a[n]==5*a[p5])p5++;
if(a[n]==7*a[p7])p7++;
}
return 0;
}
#include <iostream>
#include <cstdio>
#define min(a,b) ((a)<(b)?(a):(b))
#define min4(a,b,c,d) min(min(a,b),min(c,d))
using namespace std;
int a[5850];
int bSearch(int begin, int end, int e)
{
int mid, left = begin, right = end;
while(left <= right){
mid = (left + right) >> 1;
if(a[mid] >= e) right = mid - 1;
else left = mid + 1;
}
return left;
}
int main()
{
freopen("input.txt","r",stdin);
int n=1, T;
int p2,p3,p5,p7;
p2=p3=p5=p7=1;
a[1]=1;
while(n<5843){
a[++n]=min4(2*a[p2],3*a[p3],5*a[p5],7*a[p7]);
if(a[n]==2*a[p2])p2++;
if(a[n]==3*a[p3])p3++;
if(a[n]==5*a[p5])p5++;
if(a[n]==7*a[p7])p7++;
}
cin>>T;
while(T--){
scanf("%d",&n);
printf("%d\n",a[bSearch(1,5842,n)]);
}
return 0;
}