hdu 1695 GCD(欧拉函数+容斥原理）

题目大意：给你 a , b , c , d , k 五个值 （题目说明了 你可以认为 a=c=1）  x 属于 [1,b] ，y属于[1,d]  让你求有多少对这样的 （x,y）满足gcd（x,y）==k。给你的时间是 3000 MS。   0 < a <= b <= 100,000, 0 < c <= d <= 100,000, 0 <= k <= 100,000

分析：1由于gcd（x，y）=k满足（x/k，y/k）=1，所有区间可以缩小为[1/k,b/k],[1/k,d/k],问题也转化成为取两区间的元素，使得gcd(x,y)=1。

2.由于gcd(x,y)=1和gcd(y,x)=1是相同的,假设x<y,于是可以先取小区间进行讨论,由于gcd(x,y)=1利用欧拉函数，所以小区间可以计算为ans += phi[1] + phi[2] + ...+ phi[b/k](用欧拉函数的递推形式保存起来即可)

3.对于[b/k+1,d/k]的讨论.设y为区间的一个元素，则可以对y进行素因子分解，于是得到集合{p1,p2,p3,...}，其中pi为素数.虽然是要求gcd(x,y)=1的组合，但反过来也可以求gcd(x,y)!=1的组合，于是可以用容斥原理进行统计能被这些素数整除的数的个数，最后相减求补数即可加到ans了。（对[1,b/k]考虑素数分解也是可以的，只不过麻烦一些，因为筛[b/k+1,d/k]比筛[1,b/k]复杂一点）

代码：

#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
const int maxn = 1e6 + 10;
bool vis[maxn];
int phi[maxn],prim[maxn/10],cnt;

void getphi(){
    phi[1] = 1;
    for(int i=2;i<=maxn;i++){
        if(!vis[i]) {
            phi[i] = i-1;
            prim[++cnt] = i;
        }
        for(int j=1;j<=cnt;j++){
            int x = prim[j];
            if(i*x>maxn) break;
            vis[i*x]=true;
            if(i%x==0) {
                phi[i*x]=x*phi[i];
                break;
            }else phi[i*x]=phi[i]*phi[x];
        }
    }
}

int t,a,b,c,d,k,tt;
int main(){
    getphi();
    cin>>t;
    int cas = 1;
    while(t--){
        cin>>a>>b>>c>>d>>k;
        if(k==0) {
            cout<<"Case "<<cas++<<": 0"<<endl;
            continue;
        }
        b/=k,d/=k;
        if(b>d) tt=b,b=d,d=tt;
        ll ans = 0;
        for(int i=1;i<=b;i++)
            ans += phi[i];
        for(int i=b+1;i<=d;i++){
            ans += b;
            int tmp=i;
            int ct=0,tp[20];
            for(int j=2;j*j<=tmp;j++){
                if(tmp%j==0) {
                    while(tmp%j==0) tmp/=j;
                    tp[ct++]=j;
                }
            }
            if(tmp>1) tp[ct++]=tmp;
            for(int k=1;k< (1<<ct);k++){
                int flag=0,s=1;
                for(int l=0;l<ct;l++){
                    if(k&(1<<l)){
                        flag++;
                        s*=tp[l];
                    }
                }
                if(flag&1) ans -= b/s;
                else ans += b/s;
            }
        }
        cout<<"Case "<<cas++<<": "<<ans<<endl;
    }
}