51nod1220 约数之和

本文针对一个具体的数学算法问题,详细推导了解题过程,并利用杜教筛法进行优化,最终给出了时间复杂度为O(n^2/3)的高效求解方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题解

  首先列出题目要求的式子
  

ans=i=1nj=1nd(ij)

  一看就不会做,百度一下我就知道有一个结论,于是套上结论。
  
ans=i=1nj=1na|ib|jajb[(a,b)=1]

  你问我为啥?我也不知道
  
ans=i=1nj=1na|ib|jajbd|(a,b)μ(d)=d=1nμ(d)d|and|bna|inb|jnajb=d=1nμ(d)d|and|bna|inb|jnajb=d=1nμ(d)d|adnd|bdnad|inbd|jnajb=d=1nμ(d)a=1ndb=1ndi=1nadj=1nbdajbdb=d=1nμ(d)a=1ndb=1ndi=1nadj=1nbdajd=d=1ndμ(d)a=1ndandab=1ndj=1nbdj

  最后面那两个可以考虑用算贡献的方法化简
  
b=1ndj=1ndbj=j=1ndjndj

  令D(n)表示约数和函数d()的前n项和。
  
D(n)=i=1nini

ans=d=1nμ(d)D(nd)2

  dμ(d)可以杜教筛,D则可以小数据线性筛+大数据暴力。
  时间复杂度O(n23)

代码

//杜教筛 
#include <cstdio>
#include <algorithm>
#include <cmath>
#define maxn 1000005
#define mod 1000000007ll
using namespace std;
typedef long long ll;
ll N, f[maxn+10], d[maxn+10];
int prime[maxn/10+10], mu[maxn+10], pr[maxn+10], prd[maxn+10];
bool mark[maxn+10];
void shai()
{
    int i, j;
    d[1]=mu[1]=1;
    for(i=2;i<maxn;i++)
    {
        if(!mark[i])prime[++*prime]=i, mu[i]=-1, d[i]=i+1, pr[i]=i, prd[i]=i;
        else
            if(prd[i]==i)d[i]=((ll)prd[i]*pr[i]-1)/(pr[i]-1);
            else d[i]=d[prd[i]]*d[i/prd[i]];
        for(j=1;i*prime[j]<maxn;j++)
        {
            mark[i*prime[j]]=1;
            pr[i*prime[j]]=prime[j];
            if(i%prime[j]==0){prd[i*prime[j]]=prd[i]*prime[j];break;}
            mu[i*prime[j]]=-mu[i];
            prd[i*prime[j]]=prime[j];
        }
    }
    for(i=1;i<maxn;i++)mu[i]=(mu[i-1]+(ll)i*mu[i])%mod, d[i]=(d[i-1]+d[i])%mod;
}
ll s1(ll x){return (x*(1+x)>>1)%mod;}
ll getf(ll n){return n<maxn?mu[n]:f[N/n];}
void djs(ll n)
{
    if(n<maxn or getf(n))return;
    ll i, last, t=N/n;
    f[t]=1;
    for(i=2;i<=n;i=last+1)
    {
        last=n/(n/i);
        djs(n/i);
        f[t]=(f[t]-(s1(last)-s1(i-1))*getf(n/i))%mod;
    }
}
ll sqr(ll x){return x*x%mod;}
ll getd(ll n)
{
    if(n<maxn)return d[n];
    ll ans=0, i, last;
    for(i=1;i<=n;i=last+1)
    {
        last=n/(n/i);
        ans=(ans+(s1(last)-s1(i-1))*(n/i))%mod;
    }
    return ans;
}
int main()
{
    ll i, last, ans=0, a, b, c;
    shai();
    scanf("%lld",&N);
    djs(N);
    for(i=1;i<=N;i=last+1)
    {
        last=N/(N/i);
        a=getf(last)-getf(i-1);
        b=sqr(getd(N/i));
        ans=(ans+a*b)%mod;
    }
    printf("%lld",(ans+mod)%mod);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值