莫比乌斯反演

莫比乌斯函数

1、莫比乌斯函数定义:

对于x=p1^a1*p2^a2.......pn^an中pi为素数，ai>=1

对于u(x)有两种情况：

        情况1：存在a[i]>=2是u(x)=0;

        情况2：对于所有的a[i]=1,u(x)=(-1)^k其中k为p[i]的数目

2、

莫比乌斯反演定理

通过f[n]得到F(n)

通过F(n)得到f(n)常用

一般题目中F(n)容易得到，然后推导出f(n)

例题1：例1

给出两个点（a,c）,(b,d)两个点，用这两个点构成一个矩形，询问对于矩形里面的点（x,y）,其中gcd(x,y)=k的点一共有多少个；

思路：

AC代码：

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 50010;
int primes[N],cnt,mu[N],sum[N];
bool st[N];
void init()
{
    mu[1] = 1;
    for(int i = 2;i < N;i ++)
    {
        if(!st[i])primes[cnt ++ ] = i,mu[i] = -1;
        for(int j = 0;i*primes[j] < N;j ++)
        {
            st[primes[j]*i] = true;
            if(i%primes[j] == 0)break;
            mu[primes[j]*i] = -mu[i];
        }
    }
    for(int i = 1;i < N;i ++)sum[i] = sum[i-1] + mu[i];
}
int g(int k, int x)
{
    return k / (k / x);
}
ll f(int a, int b, int k)
{
    a = a/k,b = b/k;
    ll res  = 0;
    int n = min(a,b);
    for(int l = 1, r;l<=n;l=r+1)
    {
        r=min(n, min(g(a, l), g(b,l)));
        res += (ll)(sum[r] -  sum[l - 1]) * (a / l) * (b / l);
    }
    return res;
}
int main (){
    init();
    int T;
    scanf("%d", &T);
    while (T -- )
    {
        int a, b, c, d, k;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        printf("%lld\n",f(b,d,k) - f(a - 1,d,k) - f(b,c - 1,k) + f(a - 1,c - 1,k));
    }
    return 0;
}