acwing 215. 破译密码 (莫比乌斯反演+数论分块）

link

思路

我们需要了解一下莫比乌斯函数

下来就是 我们如何在更优 的复杂度内求 a/i b/i
我们可以想到 数论分块！！

#include <bits/stdc++.h>
using namespace  std;
#define  int long long
//typedef long long ll;
typedef pair<int,int> pii;
#define x first 
#define y second
#define pb  push_back
#define inf 1e18
#define endl '\n'
#define IOS   std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define  fer(i,a,b)  for(int i=a;i<=b;i++)
#define  der(i,a,b)  for(int i=a;i>=b;i--) 
const int maxn=1e5+10;
const int mod=1e9+7;
/*ll qsm(int a,int b)
{ll res=1;  while(b){  if(b&1) res=res*a%mod;  a=a*a%mod;  b>>=1; } return res;}
*/
const int N=5e4+100;
int dr[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
int n,k;
int m;
int mul[N];
int vis[N];
int p[N];
int sum[N];
int cnt;
void init(int n)
{
	mul[1]=1;
	fer(i,2,n)
	{
		if(!vis[i]){
			p[cnt++]=i;
			mul[i]=-1;
		}
		for(int j=0;p[j]<=n/i;j++)
		{
			vis[i*p[j]]=1;
			if(i%p[j]==0){
				mul[i*p[j]]=0;
				break;
			}
			mul[i*p[j]]=mul[i]*-1;
		}
	}
	fer(i,1,n)sum[i]=sum[i-1]+mul[i];
}
void solve()
{
	int a,b,d;
	cin>>a>>b>>d;
	a/=d;
	b/=d;
	int n=min(a,b);
	int res=0;
	for(int l=1,r;l<=n;l=r+1)
	{
		r=min(n,min(a/(a/l),b/(b/l)));
		res+=(sum[r]-sum[l-1])*(a/l)*(b/l);
	}
	cout<<res<<endl;
}
signed main()
{
	IOS;
	init(N-4);
	int _;
	cin>>_;
	while(_--) solve();
	return 0;
}