洛谷3455 【POI2007】ZAP-Queries（莫比乌斯反演）

传送门

【题目分析】

Problem b这道题的弱化版，证明看这篇博客吧。

直接粘代码。

【代码~】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=5e4+10;

int T;
int n,m,k,up,ans;
int prime[MAXN],vis[MAXN],tot;
int mu[MAXN];

int Read(){
	int i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void pre(){
	mu[1]=1;
	for(int i=2;i<=5e4;++i){
		if(!vis[i])
		  prime[++tot]=i,mu[i]=-1;
		for(int j=1;j<=tot&&i*prime[j]<=5e4;++j){
			vis[prime[j]*i]=1;
			if(i%prime[j]==0)
			  break;
			mu[prime[j]*i]=-mu[i];
		}
		mu[i]+=mu[i-1];
	}
}

int main(){
	T=Read();
	pre();
	while(T--){
		n=Read(),m=Read(),k=Read();
		n/=k,m/=k;
		ans=0;
		for(int i=1;i<=min(n,m);i=up+1){
			up=min(n/(n/i),m/(m/i));
			ans+=(mu[up]-mu[i-1])*(n/i)*(m/i);
		}
		cout<<ans<<'\n';
	}
	return 0;
}