HDU 6428 - Calculate （线性筛+莫比乌斯反演）

最新推荐文章于 2021-03-12 16:12:16 发布

啊诡是校草

最新推荐文章于 2021-03-12 16:12:16 发布

阅读量432

点赞数 1

分类专栏：莫比乌斯反演文章标签：莫比乌斯反演线性筛积性函数

本文链接：https://blog.csdn.net/RingoFF/article/details/82024911

版权

莫比乌斯反演专栏收录该内容

1 篇文章 0 订阅

订阅专栏

题解有空补……

#include<bits/stdc++.h>
#define N 10000005
using namespace std;
typedef long long ll;
int const MOD = 1 << 30;

bool vis[N];
int p[N],mu[N],cnt,phi[N];
int b[N],c[N],d[N];
ll sum[N];

void Getmu(){
    mu[1]=1;
    phi[1]=1;
    sum[1]=1;
    for(int i=1;i<N;i++)b[i]=c[i]=i;
    for(int i=2;i<N;i++){
        if(!vis[i]){
            p[++cnt]=b[i]=c[i]=i;
            mu[i]=-1;
            phi[i]=i-1;
            d[i]=1;
        }
        for(int j=1;j<=cnt;j++){
            int t=i*p[j];
            if(t>=N)break;
            vis[t]=1;
            if(i%p[j]){
                mu[t]=-mu[i];
                phi[t]=phi[i]*(p[j]-1);
                d[t]=1;
                b[t]=b[i]*p[j];
                c[t]=c[i]*p[j];
            }
            else{
                phi[t]=phi[i]*p[j];
                mu[t]=0;
                d[t]=d[i]+1;
                if(d[t]%2==1)b[t]=b[i]*p[j];
                else b[t]=b[i];
                if(d[t]%3==1)c[t]=c[i]*p[j];
                else c[t]=c[i];
                break;
            }
        }
    }
    for(int i=2;i<N;i++){
        if(!vis[i]){
            sum[i]=mu[1]*phi[i]+mu[i]*phi[1];
        }
        for(int j=1;j<=cnt;j++){
            int t=i*p[j];
            if(t>=N)break;
            if(i%p[j]){
                sum[t]=sum[p[j]]*sum[i]%MOD;
            }
            else{
                int x=i,f=p[j];
                while(x%p[j]==0)f*=p[j],x/=p[j];
                if(!sum[f]){
                    for(int l=1;;l*=p[j]){
                        sum[f]+=phi[l]*mu[f/l];
                        if(p[j]>f/l)break;
                    }
                    sum[f]%=MOD;
                    if(sum[f]<0)sum[f]+=MOD;
                }
                sum[t]=sum[f]*sum[x]%MOD;
                break;
            }
        }
    }
}

int main()
{
    Getmu();
    int t,x,y,z;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&x,&y,&z);
        ll ans=0;
        for(int f=1;f<=x;f++){
            ans+=sum[f]*(x/f)%MOD*(y/b[f])%MOD*(z/c[f])%MOD;
        }
        printf("%lld\n",(ans%MOD+MOD)%MOD);
    }
    return 0;
}