机房测试3：太阳神 ra （莫比乌斯反演）

Posted 2022-12-22 mowanying

题目：

 

 

分析：
（终于在yyr大佬的援助下弄懂了这道题。。。）

首先lcm>n的限制太少，不好直接处理，转换成求补集，也就是lcm<=n，最后用n^n-ans即可。

考虑怎么求lcm<=n:

 

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define N 100005
#define ri register int
const ll mod=1e9+7;
int mu[N],lim,pri[N],su[N],cnt=0;
void init()

    mu[1]=1;
    for(ri i=2;i<=lim;++i)
        if(!pri[i]) mu[i]=-1,su[++cnt]=i;
        for(ri j=1;j<=cnt&&i*su[j]<=lim;++j)
            pri[i*su[j]]=1;
            if(i%su[j]==0)  mu[i*su[j]]=0; break; 
            else mu[i*su[j]]=-mu[i];
        
    

int main()

    freopen("ra.in","r",stdin);
    freopen("ra.out","w",stdout);
    ll n;    
    scanf("%lld",&n);
    lim=(int)(sqrt(n)+0.5);
    init();
    ll anss=0;
    for(ri d=1;d<=lim;++d)
        ll tmp=n/d/d,ans=0;
        for(ll a=1;a*a*a<=tmp;++a)//在极小的范围里面枚举a和b 就可以得到c的范围 
         for(ll b=a;b*b<=tmp/a;++b)
             ll c=tmp/a/b-b;//注意c是统计的范围 
             if(a==b) ans=( ans + 1 + c*3 ) %mod;
             else ans=( ans + 3 + c*6 ) %mod;
        
        anss=( anss + mu[d]*ans )%mod;
    
    n%=mod;
    printf("%lld\\n",(n*n%mod-anss+mod)%mod);

View Code