其他分享
首页 > 其他分享> > LJJ爱数数

LJJ爱数数

作者:互联网

LJJ爱数数

求\(\sum_{i=1}^n\sum_{j=1}^n\sum_{k=1}^n\epsilon(gcd(i,j,k))(\frac{1}{i}+\frac{1}{j}==\frac{1}{k})\)

显然无法用Mobius反演,问题在于\(\frac{1}{i}+\frac{1}{j}==\frac{1}{k}\),要将其转换为gcd条件。

法一:先约数拆分,再证明对应相等

分数我们无法处理,所以有

\[(i+j)k=ij\]

设\(g=gcd(i,j),I=i/g,J=j/g\),接着有

\[(I+J)k=IJg,gcd(I,J)=1\]

\[gcd(i,j,k)=gcd(g,k)=1\]

考虑分数反证,设

\[\frac{I+J}{g}=\frac{IJ}{k}=d\]

\(if(d>1)\)

\[d|IJ,d|(I+J)\]

\[gcd(I+J)=1\Rightarrow gcd(I+J,IJ)=1\]

所以与原命题矛盾,故\(d=1\),因此\(I+J=g,IJ=k,gcd(g,I)=1\),满足了等式,也满足了\(gcd(g,k)=1\),所以自然枚举g,再枚举I,以此k也就确定了,接下来问题在于范围,显然要回到\(i,j,k\)

\[g^2=i+j\Rightarrow g\leq\sqrt{2n}\Rightarrow k\leq n\]

现在考虑\(I\)

\[1\leq i=Ig,j=(g-I)g\leq n\Downarrow \]

\[I\in[max(1,\frac{g^2-n}{g}),min([n/g],g-1)]\]

于是总上,设\(L(g),R(g)\)为在g意义下\(I\)的范围,有

\[ans=\sum_{g=1}^{\sqrt{2n}}\sum_{I=L(g)}^{R(g)}(gcd(g,I)==1)\]

注意到后式取值范围受到g的影响,所以对后式单独维护

\[f(d)=\sum_{I=L(g)}^{R(g)}(gcd(g,I)==d)\]

\[F(d)=([R(g)/d]-[(L(g)-1)/d])(d|g)\]

由Mobius反演定理,我们有

\[f(d)=\sum_{d|x}(x|g)([R(g)/x]-[(L(g)-1)/x])\mu(x/d)\]

因此

\[f(1)=\sum_{x|g}([R(g)/x]-[(L(g)-1)/x])\mu(x)\]

所以

\[ans=\sum_{g=1}^{\sqrt{2n}}\sum_{x|g}([R(g)/x]-[(L(g)-1)/x])\mu(x)\]

\[=\sum_{x=1}^{\sqrt{2n}}\mu(x)\sum_{x|g}^{^{\sqrt{2n}}}([R(g)/x]-[(L(g)-1)/x])\]

显然后式可以暴力算,时间复杂度应为\(O(\sqrt{n}log(\sqrt{n}))\)。

法二:先证明对应相等再约数拆分

注意到\(ik+jk=ij\),类似因式分解形式,于是添项因式分解,

\[k^2+ij-ik-jk=k^2\]

\[k(k-i)+j(i-k)=k^2\]

\[(j-k)(i-k)=k^2\]

在进行约数拆分

\[j-k=ln^2,i-k=lm^2,k=lmn\]

\[gcd(i,j,k)=1\Rightarrow l=1,j-k=n^2,i-k=m^2,k=mn\]

\[j=mn+n^2,i=mn+m^2,k=mn\]

\[gcd(mn+n^2,mn+m^2,mn)=gcd(gcd(m,n)(m+n),mn)=1\]

\[gcd(m,n)=1\]

故充分条件为\(j-k=n^2,i-k=m^2,gcd(m,n)=1\),故考虑枚举n,m,为了与题目数据范围的n区分,把数据范围的n改为N,现在关键在于找数据范围,回到i,j,首先显然\(m\in[1,\sqrt{N}]\)

\[j=mn+n^2,i=mn+m^2,k=mn\in[1,N]\Downarrow\]

\[n\leq min(\frac{-m+\sqrt{m^2+4N}}{2},\frac{N-m^2}{m})\]

设其范围为\(r(m)\),于是有

\[ans=\sum_{i=1}^{\sqrt{n}}\sum_{j=1}^{r(i)}\epsilon(gcd(i,j))\]

同法一,后式是变的,维护后式,设

\[f(d)=\sum_{j=1}^{r(i)}(gcd(i,j)==d)\]

\[F(d)=(d|i)[r(i)/d]\]

由Mobius反演定理,我们有

\[f(d)=\sum_{d|x,x|i}[r(i)/x]\mu(x/d)\]

\[f(1)=\sum_{x|i}[r(i)/x]\mu(x)\]

代入原式即

\[ans=\sum_{i=1}^{\sqrt{n}}\sum_{x|i}[r(i)/x]\mu(x)=\sum_{x=1}^{\sqrt{n}}\mu(x)\sum_{x|i}[r(i)/x]\]

显然后式处理平均只要\(log(\sqrt(n))\),所以时间复杂度应为\(\sqrt{n}log(\sqrt{n})\)。

参考代码:

#include <iostream>
#include <cstdio>
#include <cmath>
#define il inline
#define ri register
#define ll long long
using namespace std;
bool check[1000001];
int prime[100001],pt,mu[1000001],
    li[1000001];
il void prepare(int);
template<class free>
il free Min(free,free);
int main(){
    ll n,ans(0),sn,i,j;
    scanf("%lld",&n),sn=sqrt(n),prepare(sn);
    for(i=1;i<=sn;++i)li[i]=Min((n-i*i)/i,(ll)(sqrt(i*i+4*n)-i)/2);
    for(i=1;i<=sn;++i)
        for(j=i;j<=sn;j+=i)
            ans+=mu[i]*(li[j]/i);
    printf("%lld",ans);
    return 0;
}
template<class free>
il free Min(free x,free y){
    return x<y?x:y;
}
il void prepare(int n){
    int i,j;check[1]|=mu[1]|=true;
    for(i=2;i<=n;++i){
        if(!check[i])prime[++pt]=i,mu[i]=-1;
        for(j=1;j<=pt&&i*prime[j]<=n;++j){
            check[i*prime[j]]|=true;
            if(!(i%prime[j]))break;
            mu[i*prime[j]]=~mu[i]+1;
        }
    }
}

标签:数数,LJJ,frac,gcd,sum,mn,sqrt,mu
来源: https://www.cnblogs.com/a1b3c7d9/p/10797996.html