筛选法

记录点滴。

/*
2015.6    HT
ACM Work_6

*/
#include <iostream>
#include <cstdio>
using namespace std;


/*
七夕节

*/
//int main()
//{
//    int T, a, i;
//    cin >> T;
//    while (T--)
//    {
//        cin >> a;
//        int sum = 1;
//        int k = (int)sqrt((double)a);
//        for (i = 2; i <= k; i++)
//        {
//            int t;
//            if (a%i == 0)
//            {
//                sum += i;
//                t = a / i;
//                if (t != i)
//                    sum += t;
//            }
//        }
//        cout << sum << endl;
//    }
//    return 0;
//}



/*
找新朋友
和N有大于1的公约数都是新朋友

用筛选法, 因为num如果能整除 i  ,i > 1, 那么对 i 的倍数, 肯定有大于1的公约数
*/
//int p[40000];
//int euler(int num)
//{
//    memset(p, 0, sizeof (p));
//    int cnt = 0;
//    for (int i = 2; i <= num / 2; ++i)
//    {
//        if (num % i == 0 && p[i] == 0)
//        {
//            for (int j = i; j < num; j += i)
//            {
//                if (p[j] == 0)
//                    cnt++;
//                p[j] = 1;
//            }
//        }
//    }
//    return num - cnt - 1;
//}
//
//int main()
//{
//    int T;
//    cin >> T;
//    while (T--)
//    {
//        int num;
//        cin >> num;
//        cout << euler(num) << endl;
//    }
//    return 0;
//}



/*
平方和与立方和

*/
//int main()
//{
//    int n, m, sum1, sum2, i, t;
//    while (scanf_s("%d%d", &n, &m) && n && m)
//    {
//        if (n>m)
//        {
//            t = n;
//            n = m;
//            m = t;
//        }
//        sum1 = 0;
//        sum2 = 0;
//        for (i = n; i <= m; i++)
//        {
//            if (i % 2 == 0)
//                sum1 += i*i;
//            else
//                sum2 += i*i*i;
//        }
//        printf_s("%d %d
", sum1, sum2);
//    }
//    return 0;
//}



/*
Hamming Problem
三个质数p1,p2,p3，因子只包含他们的数构成连续的一列数，问第K个数是什么?

*/
__int64 h[10000];
__int64 min(__int64 a, __int64 b)
{
    return a<b ? a : b;
}

int main()
{
    int a, b, c, i, k;
    // LONGLONG
    __int64 p, q, r;
    while (scanf_s("%I64d%I64d%I64d%d", &p, &q, &r, &k) != EOF)
    {
        a = b = c = 1;
        h[1] = 1;
        for (i = 2; i <= k + 1; i++)
        {
            h[i] = min(p * h[a], min(q * h[b], r * h[c]));
            if (h[i] == p * h[a])
                a++;
            if (h[i] == q * h[b]) 
                b++;
            if (h[i] == r * h[c])
                c++;
        }
        printf_s("%I64d
", h[k + 1]);
    }
    return 0;
}