HDU 1695

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

x是[1,b]，y是[1,d]，求GCD(x,y)=k的对数（x,y无序）

对x，y都除以k，则求GCD(x,y)=1

此时枚举x，问题转化为[1,d]区间内与x互素的数字个数，这个问题是hdu 4135

有一个特殊的地方是x，y无序，对于这点只要保证x始终小于y就可以了

特判k=0

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;

typedef __int64 ll;

vector <int> v;
int b,d,k;

int main(){
    int T;
    scanf("%d",&T);
    for(int cas=1;cas<=T;cas++){
        scanf("%*d%d%*d%d%d",&b,&d,&k);
        if(!k){
            printf("Case %d: 0
",cas);
            continue;
        }
        b/=k;d/=k;
        if(b>d)swap(b,d);
        int i;
        ll ans=0; 
        for(i=1;i<=b;i++){
            int temp=i;
            v.clear();
            for(int j=2;j*j<=i;j++){
                if(i%j==0){
                    v.push_back(j);
                    while(i%j==0)i/=j;
                }
            }
            if(i>1)v.push_back(i);
            int m=v.size();
            i=temp;
            ll res=0; 
            for(int j=1;j<(1<<m);j++){
                int cnt=0;
                ll val=1;
                for(int k=0;k<m;k++){
                    if(j&(1<<k)){
                        cnt++;
                        val*=v[k];
                    }
                }
                if(cnt&1)res+=d/val-(i-1)/val;
                else res-=d/val-(i-1)/val;
            }
            ans+=(d-i+1)-res;
        }
        printf("Case %d: %I64d
",cas,ans);
    }
    return 0;
}