码迷,mamicode.com
首页 > 其他好文 > 详细

SPOJ 7001

时间:2014-10-03 12:58:24      阅读:279      评论:0      收藏:0      [点我收藏+]

标签:blog   io   os   使用   for   sp   div   c   on   

这次求的是三维的。如上篇一样,使用那个特殊的莫比乌斯反演来做。

需要注意的是,求得的只是从(1,1,1)到(n,n,n)的互质的个数,还要注意墙壁三面如(0,1,1)~(0,n,n)等。最后+3。因为是(0,0,1);

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 1000005

using namespace std;
typedef long long LL;
int mobi[N];
bool vis[N];
int n;

void initial(){
    int i,j;
    for(i=1;i<N;i++) mobi[i]=1,vis[i]=false; 
    for(i=2;i<N;i++) {
        if(vis[i]) continue;
        for(j=i;j<N;j+=i){
            vis[j]=true;
            if((j/i)%i==0){
                mobi[j]=0; continue;
            }
            mobi[j]=-mobi[j];
        }
    }
}

int main(){
	initial();
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		LL ans=0;
		for(int i=1;i<=n;i++){
			ans+=((LL)mobi[i]*(LL)(n/i)*(LL)(n/i)*(LL)(n/i)+3*(LL)mobi[i]*(LL)(n/i)*(LL)(n/i));
		}
		ans+=3;
		printf("%lld\n",ans);
	}
	return 0;
}

  

SPOJ 7001

标签:blog   io   os   使用   for   sp   div   c   on   

原文地址:http://www.cnblogs.com/jie-dcai/p/4004741.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!