莫比乌斯反演、中国剩余定理题目总结

所有题目链接

B题

需要了解的公式：
μ(i)可通过莫比乌斯反演模板所求，需要注意若输入的k不等于1，即gcd(i,j) = k,等同于gcd(i/k,j/k) = 1。接下来问题就变得简单了，只需要注意(1,2)和(2,1)是等价的，所以需要排除。则如图 先定义一个ans算全部的对数，然后再用短边算ans2的对数/2，ans-ans2即所得答案

#include<iostream>
#include<algorithm>
#define ll long long
#define maxn 100010
using namespace std;
bool vis[maxn];
ll mu[maxn],prim[maxn];
void Mobius() { //打表算μ 
	mu[1] = 1;
	int cnt = 0;
	for(ll i=2; i<=maxn; i++) {
		if(!vis[i]){
			prim[cnt++] = i;
			mu[i] = -1;
		}
		for(ll j=0; j<cnt&&i*prim[j]<maxn; j++) {
			if(i*prim[j] > maxn)
				break;
			vis[i*prim[j]] = true;
			if(i%prim[j]==0) {
				mu[i*prim[j]] = 0;
				break;
			}
			mu[i*prim[j]] = -mu[i];
		}
	} 
}
int main() {
	Mobius();
	int T;
	cin >> T;
	int num = 0;
	while(T--) {
		int a1,a,b1,b;
		int k;  // 即f(k) 
		ll ans = 0,ans2= 0;
		cin >> a1 >> a >> b1 >> b >> k;
		cout<<"Case "<<++num<<": ";
		if(k==0) {
			cout << "0" << endl;
			continue;
		}
		else 
			a/=k,b/=k; //gcd(a,b)=c  -->  gcd(a/c,b/c)=1; 
		for(ll j=1;j<=min(a,b);j++) {
			ans += mu[j]*(a/j)*(b/j);  
		}
		for(ll j=1;j<=min(a,b);j++) {
        	ans2 += mu[j]*(min(a,b)/j)*(min(a,b)/j);  //画二维图进行观察 发现多出来的是边长最小的正方形个数的一半 
   		}
		ans2 /= 2;
		ans -= ans2;
		
		cout << ans << endl;
	} 
	return 0;
}