BZOJ 2005 洛谷P4450 雙親數 容斥原理/莫比烏斯反演

容斥原理:

已知,輸入A,B,d,A內d的倍數有A/d個(設爲aa),B內d的倍數有B/d個(設爲bb)

gcd(a1,b1) = d 的有aa*bb,但是,這是需要去重的,如何去重呢,把上面的倍數次全減一遍.所以從min(a,b)到d,篩一遍,然後直接輸出d的.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 1000000
#define ll long long
using namespace std;
typedef long long LL;
LL f[N];
int n,m,d;
int main()
{
    cin>>n>>m>>d;
    int xiao = min(n,m);
    for(int i=xiao;i>=d;i--)
    {
        f[i] = 1ll*(n/i)*(m/i);
        for(int j=2;j*i<=xiao;j++)
            f[i] -= f[j*i];
    }
    cout<<f[d]<<'\n';
    return 0;
}

莫比烏斯的坑先挖着...