莫比乌斯函数及其应用 - ZAP-Queries + Problem b + GCD

莫比乌斯函数及其应用

一、莫比乌斯函数

1、定义

$\mu 为莫比乌斯函数，定义为：$

$\mu (n)=\{\begin{array}{l}1n=1\\ 0n含有平方因子\\ (-1{)}^{k}k为n的本质不同质因子个数\end{array}$

$说明：$

$记n={\prod }_{i=1}^k{p}_i^{c_i}，其中p_i为质因子，c_i\ge 1。上述定义表示：$

$①、n=1时，\mu (n)=1;$

$②、当n\ne 1时：$

$a.若存在i\in [1,k]，使得c_i>1，\mu (n)=0，即若某个质因子次数大于0，\mu (n)就等于0。$

$$\begin{gathered} b.若任意i\in [1,k]，都有c_i=1，\mu (n)=(-1{)}^k，即若每个质因子都仅出现过一次， \\ \mu (n)等于-1的k次幂，k为不同质因子的个数。 \end{gathered}$$

2、应用

$莫比乌斯函数可用于解决如下问题：$

$计算[1,N]中与N互质的数的个数。我们记S_i表示[1,N]中，是i的倍数的数的个数。$

$则答案为：N-S_2-S_3-...+S_{2,3}+S_{2,5}+...-S_{2,3,5}-S_{2,3,7}-...+...$

$即：N-\lfloor \frac{N}{2}\rfloor -\lfloor \frac{N}{3}\rfloor -...+\lfloor \frac{N}{2\cdot 3}\rfloor +\lfloor \frac{N}{2\cdot 5}\rfloor +...-\lfloor \frac{N}{2\cdot 3\cdot 5}\rfloor -\lfloor \frac{N}{2\cdot 3\cdot 7}\rfloor -...+...$

$这里，S_i前面的系数就是莫比乌斯函数的第i项\mu (i)。$

3、莫比乌斯反演定理

定理： $F(n)和f(n)是定义在非负整数集合上的两个函数，并且满足条件：$
$$F(n)=\sum _{d|n}f(d)$$

$那么有：$

$$f(n)=\sum _{d|n}\mu (d)F(\lfloor \frac{n}{d}\rfloor )(1)$$

$或者：$

$$f(n)=\sum _{n|d}\mu (\lfloor \frac{d}{n}\rfloor )F(d)(2)$$

$利用性质：{\sum }_{d|n}\mu (d)=\{\begin{array}{l}1n=1\\ 0n=0\end{array}，容易证明。$

该定理的作用：

$对于一些函数f(n)，如果我们很难直接求出它的值，而容易求出$倍数和(式2)$或$约数和(式1)$F(n)，那么我们可以通过莫比乌斯反演来求得f(n)的值。$

$例：f(n)表示某一范围内(x,y)=n的数对的数量，F(n)表示某一范围内n|(x,y)的数对的数量。$

$当f(n)不好求时，我们可以通过对F(n)反演求f(n)。$

二、例题

1、P3455 [POI2007] ZAP-Queries

题意：

对于给定的整数a,b和d，有多少正整数对x,y，满足x≤a，y≤b，并且gcd(x,y)=d。

输入格式

第一行包含一个正整数n，表示一共有n组询问。

接下来n行，每行表示一个询问，每行三个正整数，分别为a,b,d。

输出格式

对于每组询问，输出一个正整数，表示满足条件的整数对数。

数据范围

1≤n≤50000,
1≤d≤a,b≤50000

输入样例：

2
4 5 2
6 4 3

输出样例：

3
2

---

分析：

$1\le x\le a，1\le y\le b，使得gcd(x,y)=d$

$等价于：$

$1\le x^{\prime }\le \lfloor \frac{a}{d}\rfloor ，1\le y^{\prime }\le \lfloor \frac{b}{d}\rfloor ，使得gcd(x^{\prime },y^{\prime })=1$

$我们的问题转化为：计算满足1\le x^{\prime }\le \lfloor \frac{a}{d}\rfloor ，1\le y^{\prime }\le \lfloor \frac{b}{d}\rfloor ，互质的数对的数量。$

$$\begin{gathered} 答案应当是：总方案数-不合法方案数 \\ =总方案数-有公因子2的方案数-有公因子3的方案数-... \\ +有公因子2和3的方案数+有公因子2和5的方案数+... \\ -有公因子2、3、5的方案数-有公因子2、3、7的方案总数-... \\ ... \end{gathered}$$

$方案总数：a^{\prime }\cdot b^{\prime }$

$$\begin{gathered} 只有有一个公因子p_i的方案总数:区间[1,a^{\prime }]中，p_i的倍数的个数\lfloor \frac{a^{\prime }}{p_i}\rfloor ， \\ 区间[1,b^{\prime }]中，p_i的倍数的个数\lfloor \frac{b^{\prime }}{p_i}\rfloor ， \\ 那么总的方案数为\lfloor \frac{a^{\prime }}{p_i}\rfloor \cdot \lfloor \frac{b^{\prime }}{p_i}\rfloor \end{gathered}$$

$有多个公因子的方案总数的计算同上。$

$答案即：$

$$\begin{gathered} a^{\prime }{b}^{\prime }-\lfloor \frac{a^{\prime }}{2}\rfloor \lfloor \frac{b^{\prime }}{2}\rfloor -\lfloor \frac{a^{\prime }}{3}\rfloor \lfloor \frac{b^{\prime }}{3}\rfloor -...+\lfloor \frac{a^{\prime }}{2\cdot 3}\rfloor \lfloor \frac{b^{\prime }}{2\cdot 3}\rfloor +\lfloor \frac{a^{\prime }}{2\cdot 5}\rfloor \lfloor \frac{b^{\prime }}{2\cdot 5}\rfloor +... \\ -\lfloor \frac{a^{\prime }}{2\cdot 3\cdot 5}\rfloor \lfloor \frac{b^{\prime }}{2\cdot 3\cdot 5}\rfloor -\lfloor \frac{a^{\prime }}{2\cdot 3\cdot 7}\rfloor \lfloor \frac{b^{\prime }}{2\cdot 3\cdot 7}\rfloor -...+... \end{gathered}$$

$能够发现，上式等价于:$

$$\sum _{i=1}^{min(a^{\prime },b^{\prime })}\lfloor \frac{a^{\prime }}{i}\rfloor \cdot \lfloor \frac{b^{\prime }}{i}\rfloor \cdot \mu (i)$$

$$\begin{gathered} 检验：当分母的k个质因子的次数都为1，系数为(-1{)}^k。 \\ 当分母的某个质因子次数大于1，系数为0。 \end{gathered}$$

$观察这个式子，前一半我们可以通过整除分块，将时间复杂度降到O(\sqrt{n})；$

$后一半借助前缀和来计算每一块的\sum \mu (i)。$

从莫比乌斯反演的角度来看待这个问题：

$$f(k)=\sum _{1\le x\le a,1\le y\le b}[gcd(x,y)==k]$$

$记F(i)为1\le x\le a，1\le y\le b且i|gcd(x,y)的数对(x,y)的个数$

$即F(i)是区间x\in [1,a]，y\in [1,b]上，满足gcd(x,y)是i的倍数的数对的个数：$

$$F(i)=\sum _{i|d}f(d)=\lfloor \frac{a}{i}\rfloor \cdot \lfloor \frac{b}{i}\rfloor$$

$利用反演定理，有：$

$$f(i)=\sum _{i|d}\mu (\frac{d}{i})F(d)=\sum _{i|d}\mu (\frac{d}{i})\lfloor \frac{a}{d}\rfloor \lfloor \frac{b}{d}\rfloor$$

$$则：f(k)=\sum _{k|d}\mu (\frac{d}{k})F(d)=\sum _{k|d}\mu (\frac{d}{k})\cdot \lfloor \frac{a}{d}\rfloor \cdot \lfloor \frac{b}{d}\rfloor =\sum _{j=1}^{min(\frac{a}{k},\frac{b}{k})}\mu (j)\cdot (\lfloor \frac{a}{k}\rfloor /j)\cdot (\lfloor \frac{b}{k}\rfloor /j)$$

$然后我们枚举k的所有倍数d，就能够计算得到f(k)$

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>

#define ll long long

using namespace std;

const int N = 50010;

int T, a, b, d;
int mu[N], primes[N], cnt, sum[N];
bool st[N];

void get_prime(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        {
            primes[cnt++]=i;
            mu[i]=-1;
        }
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            {
                mu[i*primes[j]]=0;
                break;
            }
            mu[i*primes[j]]=mu[i]*-1;
        }
    }
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+mu[i];
}

int main()
{
    get_prime(N-1);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&a,&b,&d);
        a/=d, b/=d;
        int n=min(a,b);
        
        ll res=0;
        for(int i=1, r=1;i<=n;i=r+1)
        {
            r=min(a/(a/i), b/(b/i));
            r=min(r,n);
            res+=(sum[r]-sum[i-1])*(ll)(a/i)*(b/i);
        }
        printf("%lld\n",res);
    }
    
    return 0;
}

---

2、P2522 [HAOI2011] Problem b

题目描述

对于给出的 n 个询问，每次求有多少个数对 (x,y)，满足 a≤x≤b，c≤y≤d，且 gcd(x,y)=k，gcd(x,y) 函数为 x 和 y 的最大公约数。

输入格式

第一行一个整数 n，接下来 n 行每行五个整数，分别表示 a,b,c,d,k。

输出格式

共 n 行，每行一个整数表示满足要求的数对 (x,y) 的个数。

输入输出样例
输入

2
2 5 1 5 1
1 5 1 5 2

输出

14
3

说明/提示
对于 100% 的数据满足：$1\le n,k\le 5\times 10^4$ ，$1\le a\le b\le 5\times 10^4$，$1\le c\le d\le 5\times 10^4$。

---

分析：

$记a^{\prime }=\frac{a}{k}，b^{\prime }=\frac{b}{k}，c^{\prime }=\frac{c}{k}，d^{\prime }=\frac{d}{k}$

$本题等价于计算：$

$区间[a^{\prime },b^{\prime }]和区间[c^{\prime },d^{\prime }]中，互质的数对的个数。$

我们可以借助前缀和和容斥的思想，来做一个巧妙的转化：

$分为4个部分：$

$先计算出区间[1,b^{\prime }]和区间[1,d^{\prime }]中互质的数对的个数C_1$

$再计算出区间[1,a^{\prime }-1]和区间[1,d^{\prime }]互质的数对的个数C_2，区间[1,b^{\prime }]和区间[1,c^{\prime }-1]互质的数对个数C_3$

$区间[1,a^{\prime }-1]和区间[1,c^{\prime }-1]互质的数对个数C_4$

$则答案为：$

$$C_1-C_2-C_3+C_4$$

$而上面每一部分，我们都可以借助例题1中的方法来解决。$

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>

#define ll long long

using namespace std;

const int N = 50010;

int T, a, b, c, d, k;
int mu[N], primes[N], cnt, sum[N];
bool st[N];

void get_prime(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        {
            primes[cnt++]=i;
            mu[i]=-1;
        }
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            {
                mu[i*primes[j]]=0;
                break;
            }
            mu[i*primes[j]]=mu[i]*-1;
        }
    }
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+mu[i];
}


ll cal(int a,int b)     //计算区间[1,a]和区间[1,b]中gcd(x,y)=k的数对数量
{
    a/=k, b/=k;
    int n=min(a,b);
    ll res=0;
    for(int i=1, r=1;i<=n;i=r+1)
    {
        r=min(a/(a/i), b/(b/i));
        r=min(r,n);
        res+=(ll)(sum[r]-sum[i-1])*(a/i)*(b/i);
    }
    return res;
}

int main()
{
    get_prime(N-1);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        printf("%lld\n",cal(b,d) - cal(a-1,d) - cal(b,c-1) + cal(a-1,c-1));
    }
    
    return 0;
}

---

3、GCD - HDU 1695

题意：

给出a,b,c,d,k,求出a<=x<=b, c<=y<=d 且gcd(x,y) == k 的（x,y）的对数，其中gcd(x,y)和gcd(y,x)被认为是同一种，应当舍去。

Input

样例个数T (T <= 3000) 每个样例输入a,b,c,d,k，保证所有的a和c都等于1. (a=1 , c=1 , 0 < b,d <= 100,000 , 0 <= k <= 100,000)

Output

对于每个样例，输出其对应的答案。

Sample Input

2
1 3 1 5 1
1 11014 1 14409 9

Sample Output

Case 1: 9
Case 2: 736427

---

分析：

$要考虑重复的情况，对于区间[1,b]和区间[1,d]，假设b<d，的情况。$

$重复的方案数应当是[1,b]与[1,b]的方案总数的一半。$

$因此我们计算出区间[1,b]和区间[1,d]的gcd(x,y)=k的对数cal(b,d)，再减去重复的方案数\frac{cal(b,b)}{2}$

本题需要注意： $k可能为0.$

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>

#define ll long long

using namespace std;

const int N = 100010;

int T, a, b, c, d, k;
int mu[N], primes[N], cnt, sum[N];
bool st[N];

void get_prime(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        {
            primes[cnt++]=i;
            mu[i]=-1;
        }
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            {
                mu[i*primes[j]]=0;
                break;
            }
            mu[i*primes[j]]=mu[i]*-1;
        }
    }
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+mu[i];
}


ll cal(int a,int b)     //计算区间[1,a]和区间[1,b]中gcd(x,y)=k的数对数量
{
    a/=k, b/=k;
    int n=min(a,b);
    ll res=0;
    for(int i=1, r=1;i<=n;i=r+1)
    {
        r=min(a/(a/i), b/(b/i));
        r=min(r,n);
        res+=(ll)(sum[r]-sum[i-1])*(a/i)*(b/i);
    }
    return res;
}

int main()
{
    get_prime(N-1);
    scanf("%d",&T);
    for(int i=1;i<=T;i++)
    {
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        printf("Case %d: ",i);
        if(k==0) puts("0");
        else printf("%lld\n",cal(b,d)-cal(min(b,d),min(b,d))/2);
    }
    
    return 0;
}