HDOJ 1058 Humble Numbers

刚开始想到用类似素数筛法的方法打表，然后发现做不了。。。没忍住看了别人的题解，一眼就秒懂了。。。

果然我还是智硬啊。。。

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.
    }
}