HDU 1695 GCD（欧拉函数+容斥原理）

传送门: http://acm.hdu.edu.cn/showproblem.php?pid=1695

求（1，b）区间和（1，d）区间里gcd（x，y）== k的个数。那么将b和d分别除以k求gcd==1的就好了。但题目要求对于（3，5）（5，3）这种只算一种，限制一下x<y就好了。

接下来就是枚举一个区间的x再确定另一个区间y的个数就好了。不能直接用欧拉函数，因为欧拉函数求的是从1开始的。区间里无非是于x互素的和不互素的。求出不互素的减掉得到的就是互素的个数了。

难点在于如何求与x不互素的个数。

若b<d，对于（1，b]的，直接用欧拉函数水出答案就好了。

对（b，d）我们先预处理一下得出每个x所有的素因子并存储起来，然后枚举（b，d）的所有x，显然，x/因子得到的就是不互素的个数，但很容易看出会有重复的。不过观察一下可以看到是有规律的，即奇数个因子的话就减掉，偶数个的话就加上，这样一来就很好用递归实现了。这个就是容斥原理。看一下紫书的话在推导欧拉函数过程中有一个公式就是这个意思。

代码：

// Header.
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <string>
#include <bitset>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <set>
#include <map>
using namespace std;

// Macro
typedef long long LL;
#define mem(a, n) memset(a, n, sizeof(a))
#define rep(i, n) for(int i = 0; i < (n); i ++)
#define REP(i, t, n) for(int i = (t); i < (n); i ++)
#define FOR(i, t, n) for(int i = (t); i <= (n); i ++)
#define ALL(v) v.begin(), v.end()
#define Min(a, b) a = min(a, b)
#define Max(a, b) a = max(a, b)
#define put(a) printf("%d\n", a)
#define ss(a) scanf("%s", a)
#define si(a) scanf("%d", &a)
#define sii(a, b) scanf("%d%d", &a, &b)
#define siii(a, b, c) scanf("%d%d%d", &a, &b, &c)
#define VI vector<int>
#define pb push_back
const int inf = 0x3f3f3f3f, N = 1e5 + 5, MOD = 1e9 + 7;
// Macro end

int T, cas = 0;
int a, b, c, d, k;
LL phi[N];
VI pr[N];
// Imp
int euler_phi(int n) {
    int ret = n, m = (int)sqrt(n + 0.5);
    FOR(i, 2, m) if(n % i == 0) {
        ret = ret / i * (i - 1);
        while(n % i == 0) n /= i;
    } if(n > 1) ret = ret / n * (n - 1);
    return ret;
}

void phi_table(int n) {
    mem(phi, 0);
    phi[1] = 1;
    FOR(i, 2, n) {
        if(!phi[i]) {
            for(int j = i; j <= n; j += i) {
                if(!phi[j]) phi[j] = j;
                phi[j] = phi[j] / i * (i - 1);
                pr[j].pb(i);
            }
        }
        phi[i] += phi[i-1];
    }
}

LL Inclusion_exclusion(int idx, int t, int n) {
    LL ret = 0, tmp;
    int sz = pr[n].size();
    REP(i, idx, sz) {
        tmp = t / pr[n][i];
        ret += tmp - cal(i + 1, tmp, n);
    }
    return ret;
}

int main(){
#ifdef LOCAL
    freopen("/Users/apple/input.txt", "r", stdin);
//  freopen("/Users/apple/out.txt", "w", stdout);
#endif
    
    phi_table(N - 1);
    si(T);
    while(T --) {
        sii(a, b), siii(c, d, k);
        printf("Case %d: ", ++ cas);
        if(b > d) swap(b, d);
        if(!k) {
            puts("0");
            continue;
        }
        LL ans = 0;
        b /= k, d /= k;
        ans = (LL)phi[b];
        FOR(i, b + 1, d) ans += (LL)b - Inclusion_exclusion(0, b, i);
        printf("%I64d\n", ans);
    }
    
    return 0;
}