2019.10.3 noip模拟赛

T1 天空龙

巨水，不讲。

T3 太阳神

【问题描述】
太阳神拉很喜欢最小公倍数，有一天他想到了一个关于最小公倍 数的题目。
求满足如下条件的数对(a,b)对数:
a,b 均为正整数且 a,b<=n 而 lcm(a,b)>n。
其中的 lcm 当然表示最小公倍数。答案对 1,000,000,007 取模

【输入格式】
第一行一个正整数 n。

【输出格式】
一行一个整数表示答案，对 1,000,000,007 取模。

解析

问题先转化为求：
$$\sum _{i=1}^n\sum _{j=1}^{n}1[lcm(i,j)>n]$$

因为$i$和$j$总共的搭配有$n^2$种，所以问题可以转化为求：
$$n^2-\sum _{i=1}^n\sum _{j=1}^{n}1[lcm(i,j)\le n]$$
明显，${\sum }_{i=1}^n{\sum }_{j=1}^{n}1[lcm(i,j)\le n]$要更好求一些。下面我们探讨如何求该式。

首先，由$lcm(x,y)=\frac{x\cdot y}{gcd(x,y)}$，原式转化为：
$$\sum _{i=1}^n\sum _{j=1}^{n}1[\frac{i\cdot j}{gcd(i,j)}\le n]$$
设$gcd(i,j)=d$，则$\frac{i\cdot j}{gcd(i,j)}\le n$可转化为$\frac{i}{d}\cdot \frac{j}{d}\le \frac{n}{d}$，上式转化为：
$$\sum _{i=1}^n\sum _{j=1}^{n}1[\frac{i}{d}\cdot \frac{j}{d}\le \frac{n}{d}]$$

枚举$gcd(i,j)=d$，设$x=\frac{i}{d}$,$y=\frac{j}{d}$，由$gcd$的定义可知$gcd(x,y)=1$，原式转化为：
$$\sum _{d=1}^n\sum _{x=1}^n\sum _{y=1}^n[gcd(x,y)=1][xy\le \frac{n}{d}]$$
对$gcd(x,y)=1$进行莫比乌斯反演，得：
$$\sum _{d=1}^n\sum _{x=1}^n\sum _{y=1}^n\sum _{k|gcd(x,y)}\mu (k)[xy\le \frac{n}{d}]$$
将$k$提前，得：
$$\sum _k\mu (k)\sum _{d=1}^n\sum _{k|x}\sum _{k|y}[xy\le \frac{n}{d}]$$
因为$k$既是$x$的因数，又是$y$的因数，所以可以在中括号里同除两个$k$，从而抵消掉${\sum }_{k|x}{\sum }_{k|y}$，得到：
$$\sum _{d=1}^n\mu (k)[dxy\le \frac{n}{k^2}]$$

所以这道题就变成了枚举$k$，求$d,x,y$的可能取值组数。因为$k|gcd(x,y)$，而$gcd(x,y)\le \sqrt{n}$，所以$k$只有$\sqrt{n}$种取值。

因为不知道$x,y,d$的大小关系，所以不妨假设$x\le y\le d$，那么$x^3$一定小于$\frac{n}{k^2}$，$x^{2}y$也一定小于$\frac{n}{k^2}$。所以外层枚举$k$，内层枚举$x,y$注意$y$要从$x$开始枚举。然后讨论$d$的取值种类。$d$的取值种类有$\frac{n}{k^2\cdot x\cdot y}-y$种。由排列组合的知识可知，若$x==y$,总情况要乘3,；若$x!=y$，总情况要乘6.最后乘$\mu (k)$，计入总答案。复杂度约为$O(n^{\frac{2}{3}})$（我不会微积分）。

代码

#include<cstdio>
#include<cmath>
#define ll long long
#define maxn 100005
#define mod 1000000007
using namespace std;

ll n,m,tot,ans;
ll mu[maxn],prime[maxn];
bool mark[maxn];

void getmu()
{
	mu[1]=1;
	for(int i=2;i<=100000;++i)
	{
		if(!mark[i])
		{
			prime[++tot]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=tot;++j)
		{
			if(i*prime[j]>100000)break;
			mark[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				mu[i*prime[j]]=0;
				break;
			}
			else
			mu[i*prime[j]]=-mu[i];
		}
	}
}

int main()
{
	getmu();
	scanf("%lld",&n);
	m=sqrt(n+0.5);
	for(ll i=1;i<=m;++i)
	{
		ll nw=n/i/i,sum=0;
		for(ll a=1;a*a*a<=nw;++a)//枚举a
		{
			for(ll b=a;b*b<=nw/a;++b)//枚举b
			{
				ll c=nw/a/b-b;
				if(a==b)sum=(sum+c*3+1)%mod;
				else sum=(sum+c*6+3)%mod;
			}
		}
		sum*=mu[i];
		ans=(ans+sum)%mod;
	}
	ans=(((n%mod)*(n%mod)-ans)%mod+mod)%mod;
	printf("%lld",ans);
	return 0;
}

总结

对于莫比乌斯反演的题目，要先化公式（显然），尽量化出带等号的式子，然后莫比乌斯反演一下。不过我这辈子应该没机会再做这种题了吧。