推狮子的部分
\[ \sum_{i=1}^A\sum_{j=1}^B\sum_{k=1}^C\sigma(ijk) =\sum_{i=1}^A\sum_{j=1}^B\sum_{k=1}^C\sum_{x|i}\sum_{y|j}\sum_{z|k}\epsilon(\gcd(x,y))\epsilon(\gcd(y,z))\epsilon(\gcd(x,z))\\ =\sum_{i=1}^A\sum_{x|i}\sum_{j=1}^B\sum_{y|j}\sum_{k=1}^C\sum_{z|k}\epsilon(\gcd(x,y))\epsilon(\gcd(y,z))\epsilon(\gcd(x,z))\\ =\sum_{x=1}^A\lfloor\frac{A}{x}\rfloor\sum_{y=1}^B\lfloor\frac{B}{y}\rfloor\sum_{z=1}^C\lfloor\frac{C}{z}\rfloor\epsilon(\gcd(x,y))\epsilon(\gcd(y,z))\epsilon(\gcd(x,z))\\ =\sum_{x=1}^A\lfloor\frac{A}{x}\rfloor\sum_{y=1}^B\lfloor\frac{B}{y}\rfloor\sum_{z=1}^C\lfloor\frac{C}{z}\rfloor\sum_{d|x,d|y}\mu(d)\sum_{p|y,p|z}\mu(p)\sum_{q|x,q|z}\mu(q)\\ =\sum_{d=1}^{\min(A,B)}\mu(d)\sum_{p=1}^{\min(B,C)}\mu(p)\sum_{q=1}^{\min(A,C)}\mu(q)\sum_{d|x,q|x}^A\lfloor\frac{A}{x}\rfloor\sum_{d|y,p|y}^B\lfloor\frac{B}{y}\rfloor\sum_{p|z,q|z}^C\lfloor\frac{C}{z}\rfloor\\ =\sum_{d=1}^{\min(A,B)}\mu(d)\sum_{p=1}^{\min(B,C)}\mu(p)\sum_{q=1}^{\min(A,C)}\mu(q)\sum_{lcm(d,q)|x}^A\lfloor\frac{A}{x}\rfloor\sum_{lcm(d,p)|y}^B\lfloor\frac{B}{y}\rfloor\sum_{lcm(p,q)|z}^C\lfloor\frac{C}{z}\rfloor\\ \text{define } f(n,t)=\sum_{t|x}\lfloor\frac{n}{x}\rfloor ,N=\max(A,B,C)\\ \cdots=\sum_{d=1}^N\mu(d)\sum_{p=1}^N\mu(p)\sum_{q=1}^N\mu(q)f(A,lcm(d,q))f(B,lcm(d,p))f(C,lcm(p,q))\\ \]
计算答案
其中\(f(n,t)\)可以\(O(n\log n)\)预处理。
考虑对\(T\)个点连边建图，\(u\)、\(v\)之间有边当且仅当\(\mu(u)\not=0,\mu(v)\not=0,lcm(a,b)\not>T\)。那么图中的每个三元环都能算入答案，这里的三元环还包括只有俩点的和只有单点的。
对于包含三个点的\(<d,p,q>\)的贡献为
\[ \mu(d)\mu(p)\mu(q)\times\\ (f(A,lcm(d,q))f(B,lcm(d,p))f(C,lcm(p,q))+\\ f(A,lcm(d,q))f(B,lcm(p,q))f(C,lcm(d,p))+\\ f(A,lcm(d,p))f(B,lcm(d,q))f(C,lcm(p,q))+\\ f(A,lcm(d,p))f(B,lcm(p,q))f(C,lcm(d,q))+\\ f(A,lcm(p,q))f(B,lcm(d,p))f(C,lcm(d,q))+\\ f(A,lcm(p,q))f(B,lcm(d,q))f(C,lcm(d,p))) \]
对于包含两个、一个的环同理。统计三元环的方法参照不常用的黑科技——「三元环」。
参考实现
#include <bits/stdc++.h>
using namespace std;
const int N=1e5+10;
const int mod=1e9+7;
//[SDOI2018]旧试题
int pr[N],cnt;
int mu[N],bk[N],deg[N];
long long fa[N],fb[N],fc[N];
bool vis[N];
struct Node {int u,v,w;};
struct Edge {int ver,len;};
Node ech[N*300];
vector<Edge> e[N];
int solve(int a,int b,int c) {
long long ans=0;
int n=max(a,max(b,c)),m=0;
for(int i=1; i<=n; ++i) {
deg[i]=fa[i]=fb[i]=fc[i]=0;
e[i].clear();
}
for(int i=1; i<=n; ++i) {
for(int x=i; x<=n; x+=i) {
fa[i]+=a/x,fb[i]+=b/x,fc[i]+=c/x;
}
}
for(int i=1; i<=a&&i<=b&&i<=c; ++i) {
if(mu[i]) ans+=mu[i]*mu[i]*mu[i]*fa[i]*fb[i]*fc[i];
}
for(int g=1; g<=n; ++g) {
for(int i=1; i*g<=n; ++i) if(mu[i*g]) {
for(int j=i+1; 1LL*i*j*g<=n; ++j) if(mu[j*g]&&__gcd(i,j)==1) {
int u=i*g,v=j*g,w=i*j*g;
deg[u]++,deg[v]++,ech[++m]=(Node){u,v,w};
ans+=mu[u]*mu[u]*mu[v]*(fa[u]*fb[w]*fc[w]+fa[w]*fb[u]*fc[w]+fa[w]*fb[w]*fc[u]);
ans+=mu[u]*mu[v]*mu[v]*(fa[v]*fb[w]*fc[w]+fa[w]*fb[v]*fc[w]+fa[w]*fb[w]*fc[v]);
}
}
}
for(int i=1; i<=m; ++i) {
if(deg[ech[i].u]<deg[ech[i].v]||(deg[ech[i].u]==deg[ech[i].v]&&ech[i].u<ech[i].v))
swap(ech[i].u,ech[i].v);
e[ech[i].u].push_back((Edge){ech[i].v,ech[i].w});
}
#define veit vector<Edge>::iterator
for(int i=1; i<=n; ++i) if(mu[i]) {
for(veit j=e[i].begin(); j!=e[i].end(); ++j) bk[j->ver]=j->len;
for(veit j=e[i].begin(); j!=e[i].end(); ++j) {
for(veit k=e[j->ver].begin(); k!=e[j->ver].end(); ++k) {
if(!bk[k->ver]) continue;
ans+=mu[i]*mu[j->ver]*mu[k->ver]*(
fa[j->len]*fb[k->len]*fc[bk[k->ver]]+fa[j->len]*fb[bk[k->ver]]*fc[k->len]+fa[k->len]*fb[j->len]*fc[bk[k->ver]]+
fa[k->len]*fb[bk[k->ver]]*fc[j->len]+fa[bk[k->ver]]*fb[j->len]*fc[k->len]+fa[bk[k->ver]]*fb[k->len]*fc[j->len]
);
}
}
for(veit j=e[i].begin(); j!=e[i].end(); ++j) bk[j->ver]=0;
}
return (ans%mod+mod)%mod;
}
void sieve() {
mu[1]=1;
for(int i=2; i<N; ++i) {
if(!vis[i]) mu[pr[++cnt]=i]=-1;
for(int j=1; j<=cnt&&i*pr[j]<N; ++j) {
vis[i*pr[j]]=1;
if(i%pr[j]==0) break;
else mu[i*pr[j]]=-mu[i];
}
}
}
int main() {
sieve();
int T,a,b,c;
scanf("%d",&T);
while(T--) {
scanf("%d%d%d",&a,&b,&c);
printf("%d\n",solve(a,b,c));
}
return 0;
}