luogu P4619 [SDOI2018]旧试题

背景：

$\text{orz}$神仙题。

题意：

${\sum }_{i=1}^a{\sum }_{j=1}^b{\sum }_{k=1}^{c}d(ijk)$，其中$d(ijk)$表示$i\cdot j\cdot k$的因子的个数。

思路：

大力推式子。
$$\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^{c}d(ijk)$$

因为$d(ij)={\sum }_{u|i}{\sum }_{v|j}[\gcd (u,v)=1]$，所以有：$d(ijk)={\sum }_{u|i}{\sum }_{v|j}{\sum }_{w|k}[\gcd (u,v,k)=1]$，因此：

$$\begin{array}{rl}& =\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c\sum _{u|i}\sum _{v|j}\sum _{w|k}[\gcd (u,v,w)=1]\\ & =\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c\sum _{u|i}\sum _{v|j}\sum _{w|k}[\gcd (u,v)=1][\gcd (u,w)=1][\gcd (v,w)=1]\\ & =\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c[\gcd (i,j)=1][\gcd (i,k)=1][\gcd (j,k)=1]\lfloor \frac{a}{i}\rfloor \lfloor \frac{b}{j}\rfloor \lfloor \frac{c}{k}\rfloor \\ & =\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c\lfloor \frac{a}{i}\rfloor \lfloor \frac{b}{j}\rfloor \lfloor \frac{c}{k}\rfloor \sum _{d1|\gcd (i,j)}{\mu }_{d1}\sum _{d2|\gcd (i,k)}{\mu }_{d2}\sum _{d3|\gcd (j,k)}{\mu }_{d3}\\ & =\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c\lfloor \frac{a}{i}\rfloor \lfloor \frac{b}{j}\rfloor \lfloor \frac{c}{k}\rfloor \sum _{d1=1}^{\min (a,b)}{\mu }_{d1}[d1|\gcd (i,j)]\sum _{d2=1}^{\min (a,c)}{\mu }_{d2}[d2|\gcd (i,k)]\sum _{d3=1}^{\min (b,c)}{\mu }_{d3}[d3|\gcd (j,k)]\\ & =\sum _{d1=1}^{\min (a,b)}\sum _{d2=1}^{\min (a,c)}\sum _{d3=1}^{\min (b,c)}{\mu }_{d1}{\mu }_{d2}{\mu }_{d3}\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c\lfloor \frac{a}{i}\rfloor \lfloor \frac{b}{j}\rfloor \lfloor \frac{c}{k}\rfloor [d1|\gcd (i,j)][d2|\gcd (i,k)][d3|\gcd (j,k)]\\ & =\sum _{u=1}^{\min (a,b)}\sum _{v=1}^{\min (a,c)}\sum _{w=1}^{\min (b,c)}{\mu }_u{\mu }_v{\mu }_w\sum _{\text{lcm(u,v)}|i}\lfloor \frac{a}{i}\rfloor \sum _{\text{lcm(u,w)}|j}\lfloor \frac{b}{j}\rfloor \sum _{\text{lcm(v,w)}|k}\lfloor \frac{c}{k}\rfloor \end{array}$$

设$f_y(x)={\sum }_{x|i}\lfloor \frac{y}{i}\rfloor$，因此有：
$$=\sum _{u=1}^{\min (a,b)}\sum _{v=1}^{\min (a,c)}\sum _{w=1}^{\min (b,c)}{\mu }_u{\mu }_v{\mu }_w{f}_a(\text{lcm}(u,v))f_b(\text{lcm}(u,w))f_c(\text{lcm}(v,w))$$

我们可以预处理出$f_a$。
但是好像时间复杂度还是$3$次方级别。
考虑什么时候没有贡献：
$[1]$：${\mu }_u=0$或${\mu }_v=0$或${\mu }_w=0$；
$[2]$：$\text{lcm}(u,v)>a$或$\text{lcm}(u,w)>b$或$\text{lcm}(v,w)>c$。

怎么去除这两种废掉的情况呢？
设$u=kx,v=ky$，保证$u\perp v$。若${\mu }_u\ne 0,{\mu }_v\ne 0$，且$\text{lcm}(u,v)=kxy\le \max (a,b,c)$，则连边$u,v$。
最后我们三元环计数（据说建边用$\text{vector}$存储更快，因为访问连续的一段内存）即可。

代码：

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#define LL long long
#define mod 1000000007
using namespace std;
struct node{int u,v,lcm;} da[1000010];
vector <node> F[100010];
	int prime[100010],mu[100010],vis[100010];
	LL f[3][100010];
	bool bz[100010];
	int du[100010];
	int a,b,c,len=0;
	LL ans;
void init(int ma)
{
	bz[0]=bz[1]=true;
	mu[1]=1;
	int t=0;
	for(int i=2;i<=ma;i++)
	{
		if(!bz[i]) prime[++t]=i,mu[i]=-1;
		for(int j=1;j<=t&&i*prime[j]<=ma;j++)
		{
			bz[i*prime[j]]=true;
			if(!(i%prime[j]))
			{
				mu[i*prime[j]]=0;
				break;
			}
			mu[i*prime[j]]=-mu[i];
		}
	}
}
void clear(int limit)
{
	ans=0;
	len=0;
	for(int i=1;i<=limit;i++)
	{
		f[0][i]=f[1][i]=f[2][i]=0;
		du[i]=0;
		F[i].clear();
	}
}
void init_f(int op)
{
	int limit=(op==0?a:op==1?b:c);
	for(int i=1;i<=limit;i++)
		for(int j=i;j<=limit;j+=i)
			f[op][i]=(f[op][i]+limit/j)%mod;
}
void work_3(int limit)
{
	for(int i=1;i<=limit;i++)
		if(mu[i]) ans=(ans+(LL)mu[i]*mu[i]*mu[i]*f[0][i]*f[1][i]*f[2][i]%mod)%mod,ans=(ans+mod)%mod;
}
void work_2(int limit)
{
	for(int k=1;k<=limit;k++)
		for(int x=1;(LL)x*k<=limit;x++)
			if(mu[k*x])
				for(int y=x+1;(LL)x*y*k<=limit;y++)
					if(mu[k*y]&&__gcd(x,y)==1)
					{
						int u=x*k,v=y*k,lcm=x*y*k;
						ans=(ans+(LL)mu[u]*mu[u]*mu[v]*((f[0][u]*f[1][lcm]*f[2][lcm])%mod+(f[0][lcm]*f[1][u]*f[2][lcm])%mod+(f[0][lcm]*f[1][lcm]*f[2][u])%mod))%mod;
						ans=(ans+mod)%mod;
						ans=(ans+(LL)mu[v]*mu[v]*mu[u]*((f[0][v]*f[1][lcm]*f[2][lcm])%mod+(f[0][lcm]*f[1][v]*f[2][lcm])%mod+(f[0][lcm]*f[1][lcm]*f[2][v])%mod))%mod;
						ans=(ans+mod)%mod;
						du[u]++,du[v]++;
						da[++len]=(node){u,v,lcm};
					}
}
void work_0(int limit)
{
	for(int i=1;i<=len;i++)
	{
		int u=da[i].u,v=da[i].v,lcm=da[i].lcm;
		if((du[u]>du[v])||(du[u]==du[v]&&u<v)) F[u].push_back((node){u,v,lcm}); else F[v].push_back((node){v,u,lcm});
	}
	for(int u=1;u<=limit;u++)
	{
		if(!mu[u]) continue;
		for(int i=0;i<F[u].size();i++) vis[F[u][i].v]=F[u][i].lcm;
		
		for(int i=0;i<F[u].size();i++)
		{
			int v=F[u][i].v,uv=F[u][i].lcm;
			if(!mu[v]) continue;
			for(int j=0;j<F[v].size();j++)
			{
				int w=F[v][j].v,vw=F[v][j].lcm,uw=vis[w];
				if(!mu[w]||!vis[w]) continue;
				ans=(ans+mu[u]*mu[v]*mu[w]*((f[0][uv]*f[1][uw]*f[2][vw])%mod+(f[0][uv]*f[1][vw]*f[2][uw])%mod+
											(f[0][uw]*f[1][uv]*f[2][vw])%mod+(f[0][uw]*f[1][vw]*f[2][uv])%mod+
											(f[0][vw]*f[1][uw]*f[2][uv])%mod+(f[0][vw]*f[1][uv]*f[2][uw])%mod)%mod)%mod;
				ans=(ans+mod)%mod;
			}
		}
		for(int i=0;i<F[u].size();i++) vis[F[u][i].v]=0;
	}
}
int main()
{
	int T;
	init(100000);
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d %d %d",&a,&b,&c);
		int ma=max(max(a,b),c);
		clear(ma);
		init_f(0),init_f(1),init_f(2);
		work_3(ma);
		work_2(ma);
		work_0(ma);
		printf("%lld\n",ans);
	}
}