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