HDU - 1695 (GCD)

题意：给定 b , d , k ，求   

 

分析：转换一下就是求   这类题目可以用莫比乌斯反演和容斥做，前者会快一些，这里用容斥做一下，注意（x,y）和 (y,x) 是算一个 ；

 

代码：

#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

typedef long long ll;
const int N = 1E5+5;
struct node{
	int k;
	vector<int>fen;
}p[N];
bool is_p[N];

void init(){
	for(int i=0;i<N;i++) p[i].k=0;
	is_p[1]=1;
	for(int i=2;i<N;i++){
		if(!is_p[i]){
			for(int j=i;j<N;j+=i){
				is_p[j]=1;
				p[j].fen.push_back(i);
				p[j].k++;
			}
		}
	}
}

int r,m,n;
ll tot;

void dfs(int pos,int k,int cnt,ll sum){
	if(pos==p[k].k) return;
	ll res;
	for(int i=cnt;i<p[k].k;i++){
		res=sum*p[k].fen[i];
		if(pos&1) tot+=n/res;
	    else tot-=n/res;
	    dfs(pos+1,k,i+1,res);
	}
}

int main()
{
	init();
	int q,k;
	scanf("%d",&q);
	for(int cas=1;cas<=q;cas++){
		scanf("%d%d%d%d%d",&k,&r,&k,&m,&k);
		if(k==0){
			printf("Case %d: 0\n",cas);
            continue;
		}
		r/=k,m/=k;
		if(r>m) swap(r,m);
		ll ans=0;
		for(int i=1;i<=m;i++){
			n=min(i,r);
			tot=n;
			dfs(0,i,0,1);
			ans+=tot;
		}
		printf("Case %d: %lld\n",cas,ans);
	}
}