刚开始想到用类似素数筛法的方法打表,然后发现做不了。。。没忍住看了别人的题解,一眼就秒懂了。。。
果然我还是智硬啊。。。
humble number 是质因子只有2,3,5,7的数。
即有 hum[i] = pow(2,a) * pow(3,b) * pow(5,c) * pow(7,d) (a, b, c, d 属于 N)
任何一个 hum[i] 乘 2 或 3 或 5 或 7 得到的必然是另一个 hum[j] (j > i)
假设现在有4个序列,分别是所有 hum[i] 乘 2,3,5,7 的到的所有 humble number,那么我们要做的就是分别从这四个序列里找最小的加入我们的序列,然后将这个数从那四个序列里去掉。。。
记四个标志即可,不要忘记四个序列里有重复出现的数字,因为我们一直是找最小的,所以重复的话也一定是和现在所得序列里最小的那个重复。
总结:一定不要急着看题解。。。
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
long hum[5842];
void init()
{
int p2, p3, p5, p7;
hum[0] = 1;
p2 = p3 = p5 = p7 = 0;
for(int i = 1; i < 5842; i++) {
long t2 = hum[p2] * 2, t3 = hum[p3] * 3, t5 = hum[p5] * 5, t7 = hum[p7] * 7;
hum[i] = min(min(t2, t3), min(t5, t7));
if(hum[i] == t2) p2++;
else if(hum[i] == t3) p3++;
else if(hum[i] == t5) p5++;
else p7++;
if(hum[i] == hum[i-1]) i--;
}
}
int main()
{
int n;
init();
while(scanf("%d", &n) && n) {
printf("The %d", n);
if(n % 10 == 1 && (n / 10) % 10 != 1) printf("st");
else if(n % 10 == 2 && (n / 10) % 10 != 1) printf("nd");
else if(n % 10 == 3 && (n / 10) % 10 != 1) printf("rd");
else printf("th");
printf(" humble number is %ld.\n", hum[n-1]);
//The 5842nd humble number is 2000000000.
}
}