HDU 1695 GCD

HDU 1695 GCD

题意：
给你a, b, c, d, k.
求在[a, b]中找到一个x,[c, d]中找到一个y,求x和y的最大公约数是k的对数。
思路：
题目中说明要找的是最大公约数是k的两个数，那么我们把区间都除以k,再在剩下的这些区间找出两个互质的数，不就可以了么。可以把原问题分成两个部分。
第一部分：[1, b / k]
由于a和c都假定成1，而且需要按照x <= y来找（为了不重复），那么在[1, b / k]，求出每个小于次数并且与此数互质的数的个数就可以了。而这个可以通过欧拉函数快速求出。
第二部分：(b / k, d / k]
对于此区间的每个数，找到小于次数且与此数互质的数的个数。
那么我们就可以应用容斥原理来计算。
此数一个质因数的个数 - 此数两个质因子的乘积的个数 + 此数三个质因子的乘积的个数 -….
Code:

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<bitset>
#include<queue>
#include<stack>
#include<list>
#include<map>
#include<set>

#define TEST

#define LL long long
#define Mt(f, x) memset(f, x, sizeof(f));
#define rep(i, s, e) for(int i = (s); i <= (e); ++i)
#ifdef TEST
    #define See(a) cout << #a << " = " << a << endl;
    #define See2(a, b) cout << #a << " = " << a << ' ' << #b << " = " << b << endl;
    #define debug(a, s, e) rep(_i, s, e) {cout << a[_i] << ' ';} cout << endl;
    #define debug2(a, s, e, ss, ee) rep(i_, s, e) {debug(a[i_], ss, ee)}
#else
    #define See(a)
    #define See2(a, b)
    #define debug(a, s, e)
    #define debug2(a, s, e, ss, ee)
#endif // TEST

const int MAX = 2e9;
const int MIN = -2e9;
const double eps = 1e-8;
const double PI = acos(-1.0);

using namespace std;

const int N = 1e5 + 5;

LL f[N];

vector<int> ou[N];

void init()//经典的初始化 求出了前N - 1项的欧拉函数，并分解了质因数。
{
    f[1] = 1;
    for(int i = 2; i < N; ++i)
    {
        if(!f[i])
        {
            for(int k = i; k < N; k += i)
            {
                if(!f[k]) f[k] = k;
                f[k] = f[k] * (i - 1) / i;
                ou[k].push_back(i);
            }
        }
        f[i] += f[i - 1];
    }
}

int dfs(int index, int b, int n)
{
    int ret = 0;
    for(int i = index; i < ou[n].size(); ++i)
    {
        int t = b / ou[n][i];
        ret += t - dfs(i + 1, t, n);
    }
    return ret;
}

int main()
{
    init();
    int T;
    cin >> T;
    for(int _ = 1; _ <= T; ++_)
    {
        int a, b, c, d, K;
        scanf("%d%d%d%d%d", &a, &b, &c, &d, &K);
        printf("Case %d: ", _);
        if(K == 0)
        {
            printf("0\n");
            continue;
        }
        if(b > d)
        {
            swap(b, d);
        }
        int st = b / K, ed = d / K;
        LL ans = f[st];
        for(int i = st + 1; i <= ed; ++i)
        {
            ans += st - dfs(0, st, i);
        }
        printf("%lld\n", ans);
    }
    return 0;
}