【51NOD 1238】最小公倍数之和 V3

最新推荐文章于 2022-02-25 20:59:32 发布

HOWARLI

最新推荐文章于 2022-02-25 20:59:32 发布

阅读量2.3k

点赞数

CC 4.0 BY-SA版权

分类专栏： # 莫比乌斯反演 # 杜教筛文章标签：莫比乌斯反演欧拉函数杜教筛

本文链接：https://blog.csdn.net/howarli/article/details/61926079

莫比乌斯反演同时被 2 个专栏收录

23 篇文章

订阅专栏

杜教筛

8 篇文章

订阅专栏

本文提出了一种高效的算法解决方案，用于计算特定数学问题。通过使用反演和杜教筛等技巧，该方案能够在较短时间内求解涉及大量数值运算的问题，复杂度达到了O(n^(2/3))。

Description

A n s = \sum i = 1 n \sum j = 1 n l c m (i, j)

$Ans=\sum_{i=1}^n\sum_{j=1}^n lcm(i,j)$

Solution1

设 $g_d$ 表示 $\gcd(i,j)$ 为d的倍数的 $i*j$ 之和，
反演过程省略（反演时先除一下，后面可以消掉）
设： $S_n=(\frac{1}{2}*(n+1)n)^2$

A n s = \sum i = 1 n μ (i) i 2 \sum d = 1 ⌊ n i ⌋ d * S ⌊ n i d ⌋

$Ans=\sum_{i=1}^n\mu(i)i^2\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}d*S_{\lfloor\frac{n}{id}\rfloor}$
用杜教筛做

∑xi=1μ(i)i2 $\sum_{i=1}^x\mu(i)i^2$

\sum d = 1 n d 2 \sum i | d μ (i) = 1

$\sum_{d=1}^nd^2\sum_{i|d}\mu(i)=1$
这个式子显然是等于1的，
这个还可以转化一下：

\sum i = 1 n i 2 \sum j = 1 ⌊ n i ⌋ μ (j) j 2 = 1

$\sum_{i=1}^ni^2\sum_{j=1}^{\lfloor\frac{n}{i}\rfloor}\mu(j)j^2=1$
于是

\sum i = 1 n μ (i) i 2 = 1 - \sum i = 2 n i 2 \sum j = 1 ⌊ n i ⌋ μ (j) j 2

$\sum_{i=1}^n\mu(i)i^2=1-\sum_{i=2}^ni^2\sum_{j=1}^{\lfloor\frac{n}{i}\rfloor}\mu(j)j^2$
求出这个以后整个用分块搞一波，

发现这个怎么卡都过不了啊QwQ。。。

Solution2

A n s = (2 * \sum i = 1 n \sum j = 1 i - 1 i j gcd ( i , j )) + n * (n + 1) / 2

$Ans=(2*\sum_{i=1}^n\sum_{j=1}^{i-1}\frac{ij}{\gcd(i,j)})+n*(n+1)/2$
设

ans′=∑ni=1∑i−1j=1ijgcd(i,j) $ans'=\sum_{i=1}^n\sum_{j=1}^{i-1}\frac{ij}{\gcd(i,j)}$ （注意j的范围）
设：

Sn=12∗(n+1)n $S_n=\frac{1}{2}*(n+1)n$
（先枚举两个互质的数，再算它们的倍数）（后面有phi的括号是互质的数的和）

a n s' = \sum i = 1 n i * S ⌊ n i ⌋ * (φ (i) * i * 1 2)

$ans'=\sum_{i=1}^ni*S_{\lfloor\frac{n}{i}\rfloor}*(\varphi(i)*i*\frac{1}{2})$

a n s' = \sum i = 1 n 1 2 * S ⌊ n i ⌋ * φ (i) * i 2

$ans'=\sum_{i=1}^n\frac{1}{2}*S_{\lfloor\frac{n}{i}\rfloor}*\varphi(i)*i^2$
用杜教筛求

∑xi=1φ(i)∗i2 $\sum_{i=1}^x\varphi(i)*i^2$ （下面的n和上面的不一样）

\sum d = 1 n d 2 \sum i | d φ (i)

$\sum_{d=1}^nd^2\sum_{i|d}\varphi(i)$
这个式子显然是立方和，
这个还可以转化一下：

\sum i = 1 n i 2 \sum j = 1 ⌊ n i ⌋ φ (j) j 2 = (S n) 2

$\sum_{i=1}^ni^2\sum_{j=1}^{\lfloor\frac{n}{i}\rfloor}\varphi(j)j^2=(S_n)^2$
于是

\sum i = 1 n φ (i) i 2 = (S n) 2 - \sum i = 2 n i 2 \sum j = 1 ⌊ n i ⌋ φ (j) j 2

$\sum_{i=1}^n\varphi(i)i^2=(S_n)^2-\sum_{i=2}^ni^2\sum_{j=1}^{\lfloor\frac{n}{i}\rfloor}\varphi(j)j^2$

回到刚才的式子，用分块搞一波即可

复杂度： $O(n^{\frac{2}{3}})$

Code

#include <cstdio>
#include <cstdlib>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define sqr(q) ((q)*(q))
using namespace std;
typedef long long LL;
const int N=20000500,mo=1e9+7,M=224737;
const LL eni=500000004,sni=166666668;
LL n,ans;
bool prz[N+10];
int pr[N/2];
int phi[N+10];
int Hx[M+1][2];
int co;
int HX(LL q)
{
    int i=q%M;
    while(Hx[i][0]&&Hx[i][0]!=q)i=(i+1)%M;
    return i;
}
LL SM(LL s,LL t)
{
    if((s+t)&1)return (s+t)%mo*(((t-s+1)>>1)%mo)%mo;
    return ((s+t)>>1)%mo*((t-s+1)%mo)%mo;
}
LL PF(LL q)
{
    q%=mo;
    return q*(q+1)%mo*(2*q+1)%mo*sni%mo;
}
LL Gphi(LL q)
{
    if(q<=N)return phi[q];
    int t=HX(q);
    if(Hx[t][0])return Hx[t][1];
    Hx[t][0]=q;
    LL ans=0;
    for(LL i=2,nx;i<=q;i=nx+1)
    {
        nx=q/(q/i);
        ans=(ans+(PF(nx)-PF(i-1))*Gphi(q/i))%mo;
    }
    q%=mo;
    return Hx[t][1]=(sqr(SM(1,q))-ans)%mo;
}
int main()
{
    phi[1]=1;
    fo(i,2,N)
    {
        if(!prz[i])pr[++pr[0]]=i,phi[i]=i-1;
        fo(j,1,pr[0])
        {
            int t=pr[j]*i;
            if(t>N)break;
            prz[t]=1;
            phi[t]=phi[i]*pr[j];
            if(i%pr[j]==0)break;
            phi[t]=phi[i]*(pr[j]-1);
        }
    }
    fo(i,2,N)phi[i]=((LL)phi[i]*(sqr((LL)i)%mo)%mo+phi[i-1])%mo;
    scanf("%lld",&n);
    for(LL i=2,nx;i<=n;i=nx+1)
    {
        nx=n/(n/i);
        ans=(ans+SM(1,n/i)*(Gphi(nx)-Gphi(i-1))%mo*eni%mo)%mo;
    }
    printf("%lld\n",((ans*2+SM(1,n))%mo+mo)%mo);
    return 0;
}