JZYZOJ 1375 双亲数 莫比乌斯反演

http://172.20.6.3/Problem_Show.asp?id=1375

网上搜推理图。

有一段没有写莫比乌斯反演都快忘了。。数学能力--，定理完全不会推，但是这道题整体来说应该是比较好写的（虽然我没写出来）

#include<iostream>  
#include<cstdio>  
#include<cstring>  
#include<algorithm>  
#include<cmath>
using namespace std;
const int maxn=1000010;
long long a,b,d;
long long f[maxn]={},su[maxn]={};
long long w[maxn][2]={};
long long tot=0,cnt1=0,cnt2=0;
bool vis[maxn]={};
int main(){
    scanf("%I64d%I64d%I64d",&a,&b,&d);
    a/=d,b/=d;
    if(b<a)swap(a,b);
    f[1]=1;   int z;
    for(int i=2;i<=a;i++){
        if(!vis[i])su[++tot]=i,f[i]=-1;
        for(int j=1;j<=tot;j++){
            z=i*su[j];
            if(z>a)break;
            vis[z]=1;
            if(i%su[j]) f[z]=-f[i];
            else break;
        }
    }
    long long ans=0;
    for(int i=1;i<=a;i++){
        ans+=f[i]*(a/i)*(b/i);
    }
    printf("%I64d\n",ans);
    return 0;
}