莫比乌斯反演1003 SPOJ VLATTICE

题意:
有一个n*n*n的格子,这个格子是(0~n,0~n,0~n)的
从x点能看见y点,代表从x到y连线上没有其他点
问从(0,0,0)能看到多少个点
思路:
之前有个能量收集那道题跟这道题差不多…
考虑gcd(x,y,z)=1
但是这道题要考虑一些特殊情况,有0的存在
假设x,y,z都!=0,Ans1=∑(1<=T<=n,mu(T)(n/T)(n/T)*(n/T))
假设x,y,z有一个0,Ans2=3* ∑(1<=T<=n,mu(T)(n/T)(n/T))因为是x,y,z其中一个,有三种情况
假设x,y,z有两个0,那么只可能是(0,0,1),(0,1,0),(1,0,0)三种情况
所以Ans=Ans1+Ans2+3->∑(1<=T<=n,mu(T)(n/T)(n/T)*(n/T+3)) + 3

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
#define maxp 1000005
bool notprimes[maxp];
int primes[maxp];
int mu[maxp];
void get_mu()
{
    memset(notprimes,false,sizeof(notprimes));
    primes[0]=0;
    mu[1]=1;
    for(int i=2;i<maxp;++i)
    {
        if(!notprimes[i])
        {
            primes[++primes[0]]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=primes[0];++j)
        {
            if((LL)primes[j]*i>=maxp) break;
            notprimes[i*primes[j]]=true;
            if(i%primes[j])
                mu[i*primes[j]]=-mu[i];
            else
            {
                mu[i*primes[j]]=0;
                break;//代表i不是素数,mu[i*primes[j]]必然是0 
            }
        }
    }
}
LL Pre[maxp];
int main()
{
    get_mu();
    Pre[0]=0;
    for(int i=1;i<maxp;++i)
        Pre[i]=Pre[i-1]+mu[i];
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        LL Ans=0;
        int last;
        for(int i=1;i<=n;i=last+1)
        {
            last=n/(n/i);
            Ans+=(Pre[last]-Pre[i-1])*(n/i)*(n/i)*(n/i+3);
        }
        printf("%lld\n",Ans+3);
    }
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值