UVA 11014Make a Crystal

题意：给定一个N*N*N的立方体，坐标范围为-N/2<=x,y<=N/2。求选出最多的点使得任意两点P,Q没有PQO共线，O为原点且POQ共线是可以的。。

分析：要求没有PQO共线，但是PQO是可以的，那么就是说可以跨原点共线，那么我们就可以把点分为3部分：(1)坐标轴上的点。(2)8个象限的点。(3)12个四分之一平面中的点。(1)中显然只有6个点。(2)是(3)的三维情况，我们先对(2)进行分析：令g(x)={(a,b,c),gcd(a,b,c)==x}的三元组的个数，令f(x)={(a,b,c),x|gcd(a,b,c)}的三元组的个数且f(x)=[n/x]*[n/x][n/x]([]为向下取整)，那么有f(x)=sigma(x|y)g(y)==>g(x)=sigma(x|y)*u(y/x)*f(y)。我们要求的就是g(1)啦，这就是莫比乌斯反演啦。(3)的情况只要令f(x)=[n/x]*[n/x]即可。O(t*sqrt(n))

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=100010;
const int MAX=151;
const int MOD1=1000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const ll INF=10000000010;
typedef unsigned long long ull;
int a[N],q[N],mu[N];
void deal() {
    int i,j,k=0,n=100000;
    memset(q,0,sizeof(q));
    q[1]=mu[1]=1;mu[0]=0;
    for (i=2;i<=n;i++) {
        if (!q[i]) { a[++k]=i;mu[i]=-1; }
        for (j=1;j<=k;j++) {
            if (a[j]*i>n) break ;
            q[a[j]*i]=1;mu[a[j]*i]=-mu[i];
            if (i%a[j]==0) { mu[a[j]*i]=0;break ; }
        }
    }
    for (i=2;i<=n;i++) mu[i]+=mu[i-1];
}
int main()
{
    int i,n,ca=0;
    ll ans;
    deal();
    while (scanf("%d", &n)&&n) {
        ans=6;n/=2;
        int i,t,last;
        for (i=1;i<=n;i=last+1) {
            last=n/(n/i);t=n/i;
            ans+=((ll)8*t*t*t+(ll)12*t*t)*(mu[last]-mu[i-1]);
        }
        printf("Crystal %d: %lld\n", ++ca, ans);
    }
    return 0;
}