[51nod 1222] 最小公倍数计数（莫比乌斯反演）

题意

• 给定$\mathrm{a},\mathrm{b}$，求有多少对$\mathrm{i}\le \mathrm{j}$满足$\mathrm{l}\mathrm{c}\mathrm{m}(\mathrm{i},\mathrm{j})\in [\mathrm{a},\mathrm{b}]$（$\mathrm{a}\le \mathrm{b}\le 10^{11})$。

参考了ww3113306的这篇博客。

其实也比较套路，枚举$\mathtt{g}\mathtt{c}\mathtt{d}$后反演就转化成了求有多少个三元组满足$\mathtt{a}\mathtt{b}\mathtt{c}\le \mathtt{n}$的问题，这个问题也很简单，强制给$\mathtt{a},\mathtt{b},\mathtt{c}$定下大小关系，暴力枚举可以做到$\mathtt{O}({\mathtt{n}}^{\frac{2}{3}})$的复杂度。

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

const int maxn=sqrt(1e11)+9;

int notp[maxn],prime[maxn],mui[maxn],cnt;

void Get_Prime(int n)
{
	notp[0]=notp[1]=1,mui[1]=1;
	for(int i=2;i<=n;++i)
	{
		if(!notp[i]) prime[++cnt]=i,mui[i]=-1;
		for(int j=1;j<=cnt&&prime[j]*i<=n;++j)
		{
			notp[prime[j]*i]=1;
			if(i%prime[j]==0) break;
			mui[prime[j]*i]=-mui[i];
		}
	}
}

ll h(ll n)
{
	ll res=0;
	for(ll i=1;i*i*i<=n;++i,++res)
		for(ll j=i+1;j<n/(i*j);++j)
			res+=(n/(i*j)-j)*6;
	for(ll i=1;i<n/(i*i);++i)
		res+=(n/(i*i)-i)*3;
	for(ll i=1;i*i*i<=n&&i<sqrt(n/i);++i)
		res+=((ll)sqrt(n/i)-i)*3;
	return res;
}

ll f(ll n)
{
	ll res=0;
	for(ll i=1;i*i<=n;++i)
		res+=mui[i]*h(n/(i*i));
	return (res+n)>>1;
}

int main()
{
	freopen("1222.in","r",stdin);
	freopen("1222.out","w",stdout);

	ll a,b;

	Get_Prime(maxn-1);

	cin>>a>>b;
	cout<<f(b)-f(a-1)<<endl;

	return 0;
}