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;
}

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