莫比乌斯反演练习

最新推荐文章于 2021-08-15 23:30:41 发布
最新推荐文章于 2021-08-15 23:30:41 发布
阅读量152 收藏
点赞数
CC 4.0 BY-SA版权
原文链接：http://www.cnblogs.com/uid001/p/10563110.html
博客给出线性筛、整除分块等数论算法板子，还介绍了 $O(n)$ 求 $g$ 函数的方法。通过多个洛谷练习题，如 YY 的 GCD、杜教筛等，详细讲解了反演题目转化为卷积、交换求和次序，以及利用杜教筛求积性函数前缀和等解题思路和复杂度分析。

线性筛板子

const int N = 1e6+10;
int phi[N], mu[N], p[N], cnt, vis[N];
void init() {
	phi[1] = mu[1] = 1;
	for (int i=2; i<N; ++i) {
		if (!vis[i]) p[++cnt]=i,mu[i]=-1,phi[i]=i-1;
		for (int j=1,t; j<=cnt&&i*p[j]<N; ++j) {
			vis[t=i*p[j]] = 1;
			if (i%p[j]==0) {mu[t]=0,phi[t]=phi[i]*p[j];break;}
			mu[t]=-mu[i],phi[t]=phi[i]*phi[p[j]];
		}
	}
}

$O(nloglogn)$求$g=\mu * f$ 板子

void calc(int n, int *f, int *g) {
	REP(i,1,n) g[i] = f[i];
	REP(i,1,cnt) PER(j,1,n/p[i]) g[j*p[i]] -= g[j];
}

整除分块板子, 优化了一下除法次数

int mx = sqrt(n);
REP(i,1,mx) ans += (ll)(s[i]-s[i-1])*f(n/i);
for (int i=mx+1,j,k=n/i; i<=n; i=j+1,--k) {
    j = n/k;
    ans += (ll)(s[j]-s[i-1])*f(k);
}

 

 

练习1. luogu P2257 YY的GCD

大意: 求$\sum\limits_{i=1}^n \sum\limits_{j=1}^m[gcd(i,j)为素数]$

这种反演题目一般都是先转化为某个$g$函数的卷积, 再交换求和次序.

设$f(x)=[x为素数]$,  $g=\mu *f=\sum\limits_{p|n}\mu(\frac{n}{p})$, 则有

$\begin{align*} \sum\limits_{i=1}^n \sum\limits_{j=1}^mf(gcd(i,j)) &= \sum\limits_{i=1}^n \sum\limits_{j=1}^m \sum\limits_{d|gcd(i,j)} g(d) \\ &= \sum\limits_{d=1}^{min(n,m)}g(d)\sum\limits_{i=1}^n\sum\limits_{j=1}^m [d|gcd(i,j)] \\ &= \sum\limits_{d=1}^{min(n,m)}g(d) \lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor \end{align*}$

所以求出$g$以后整除分块即可, 考虑如何求$g$, $O(nloglogn)$当然很容易求出

这里思考一下$O(n)$的求法

对于$g(x)$, 假设$y$为$x$的最小素因子, 即$x=iy$

若$x$为素数, 显然有$g(x)=\mu(1)=1$

若$i%y==0$, $g(x)=\mu(i)$

若$i%y!=0$, $g(x)=-g(i)+mu(i)$

欧拉筛一次即可$O(n)-O(q\sqrt{n})$

const int N = 1e7+10;
int mu[N], p[N], cnt, vis[N], g[N];
void init() {
	mu[1] = 1;
	for (int i=2; i<N; ++i) {
		if (!vis[i]) p[++cnt]=i,mu[i]=-1,g[i]=1;
		for (int j=1,t; j<=cnt&&i*p[j]<N; ++j) {
			vis[t=i*p[j]] = 1;
			if (i%p[j]==0) {mu[t]=0,g[t]=mu[i];break;}
			mu[t]=-mu[i],g[t]=-g[i]+mu[i];
		}
		g[i]+=g[i-1];
	}
}

int main() {
	init();
	int t;
	scanf("%d", &t);
	REP(i,1,t) {
		int n, m;
		scanf("%d%d", &n, &m);
		ll ans = 0;
		if (n>m) swap(n,m);
		int lim = min(int(sqrt(n)*2),n);
		REP(i,1,lim) ans+=(ll)(g[i]-g[i-1])*(n/i)*(m/i);
		for (int i=lim+1,j; i<=n; i=j+1) {
			j = min(n/(n/i),m/(m/i));
			ans+=(ll)(g[j]-g[i-1])*(n/i)*(m/i);
		}
		printf("%lld\n", ans);
	}
}

练习2. Luogu P4213 杜教筛

大意: 求$\sum\limits_{i=1}^{n}\varphi(i),\sum\limits_{i=1}^{n}\mu(i)$
杜教筛是用来求一类积性函数前缀和的算法, 核心思想是利用公式
$$g(1)S(n)=\sum\limits_{i=1}^n(f*g)(i)-\sum\limits_{i=2}^n g(i)S(\lfloor \frac{n}{i}\rfloor)$$
其中$S(n)$是我们要求的$f(n)$的前缀和, $g(n)$是我们需要构造的一个函数, 需要满足$f*g$和$g$的前缀和可以$O(1)$求出.
可以证明上面这个递归式调用的不同的$S(n)$不超过$2\sqrt{n}$个, 即
$$S(1),S(2),...,S(\sqrt{n}),S(\frac{n}{1}),S(\frac{n}{2}),...,S(\frac{n}{\sqrt{n}})$$
对于$i\le n^{\frac{2}{3}}$的$S(i)$, 先暴力求出.
而对于$i>n^{\frac{2}{3}}$的$S(i)$, 即$S(\frac{n}{1}),S(\frac{n}{2}),...,S(\frac{n}{n^{\frac{1}{3}}})$, 只有$n^{\frac{1}{3}}$项, 可以每一项$\sqrt{n}$求出
所以总复杂度为$O(n^{\frac{2}{3}})$

对于本题来说, 对于$\sum\limits_{i=1}^n\mu(i)$取$g=\boldsymbol 1$, 就有$f*g=\varepsilon$, 对于$\sum\limits_{i=1}^{n}\varphi(i)$同理

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <math.h>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <string.h>
#include <bitset>
#define REP(i,a,n) for(int i=a;i<=n;++i)
#define PER(i,a,n) for(int i=n;i>=a;--i)
#define hr putchar(10)
#define pb push_back
#define lc (o<<1)
#define rc (lc|1)
#define mid ((l+r)>>1)
#define ls lc,l,mid
#define rs rc,mid+1,r
#define x first
#define y second
#define io std::ios::sync_with_stdio(false)
#define endl '\n'
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int P = 1e9+7, INF = 0x3f3f3f3f;
ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}
ll qpow(ll a,ll n) {ll r=1%P;for (a%=P;n;a=a*a%P,n>>=1)if(n&1)r=r*a%P;return r;}
ll inv(ll x){return x<=1?1:inv(P%x)*(P-P/x)%P;}
//head


const int N1 = ~0u>>1, N2 = 2e6, N3 = 1300;

int N;
ll phi[N2], phi2[N3];
int mu[N2], mu2[N3];
int vis[N3], v[N2];
int p[N2], cnt;

inline void init(int n) {
    phi[1] = mu[1] = 1;
    REP(i,2,n) {
        if (!v[i]) p[++cnt]=i,mu[i]=-1,phi[i]=i-1;
        for (int j=1,t; j<=cnt&&i*p[j]<=n; ++j) {
            v[t=i*p[j]] = 1;
            if (i%p[j]==0) {mu[t]=0,phi[t]=phi[i]*p[j];break;}
            mu[t]=-mu[i],phi[t]=phi[i]*phi[p[j]];
        }
    }
    REP(i,2,N2-1) mu[i]+=mu[i-1],phi[i]+=phi[i-1];
}

inline ll sum_phi(int n) {
    if (n<N2) return phi[n];
    int x = N/n;
    if (vis[x]) return phi2[x];
    vis[x] = 1;
    ll &ans = phi2[x];
    ans = (ll)n*(n+1)/2;
    int mx = sqrt(n);
    REP(i,2,mx) ans -= sum_phi(n/i);
    for (int i=mx+1,j,k=n/i; i<=n; i=j+1,--k) {
        j = n/k;
        ans -= (j-i+1)*sum_phi(k);
    }
    return ans;
}
inline int sum_mu(int n) {
    if (n<N2) return mu[n];
    int x = N/n;
    if (vis[x]) return mu2[x];
    vis[x] = 1;
    int &ans = mu2[x];
    ans = 1;
    int mx = sqrt(n);
    REP(i,2,mx) ans -= sum_mu(n/i);
    for (int i=mx+1,j,k=n/i; i<=n; i=j+1,--k) {
        j = n/k;
        ans -= (j-i+1)*sum_mu(k);
    }
    return ans;
}

int main() {
    init(N2-1);
    int t;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &N);
        if (!N) {puts("0\n0");continue;}
        memset(vis,0,sizeof vis);
        ll ans1 = sum_phi(N);
        memset(vis,0,sizeof vis);
        int ans2 = sum_mu(N);
        printf("%lld %d\n", ans1, ans2);
    }
}

练习3. Luogu P3768 简单数学题

大意: 求$\sum\limits_{i=1}^n\sum\limits_{j=1}^nijgcd(i,j)$

又是个杜教板子题, 简单推一下.

$\begin{align*} \sum\limits_{i=1}^n\sum\limits_{j=1}^nijgcd(i,j) &= \sum\limits_{i=1}^n\sum\limits_{j=1}^nij\sum\limits_{d|gcd(i,j)} \varphi(d) \\&= \sum\limits_{d=1}^n\varphi(d)\sum\limits_{i=1}^n\sum\limits_{j=1}^nij[d|gcd(i,j)] \\ &= \sum\limits_{d=1}^nd^2\varphi(d)\sum\limits_{i=1}^{\lfloor \frac{n}{d}\rfloor}i\sum\limits_{j=1}^{\lfloor \frac{n}{d}\rfloor}j\end{align*}$

对这个式子用整除分块, 转化为求$\sum\limits_{i=1}^ni^2\varphi(i)$

按照杜教的套路, 这里可以凑一个函数$g=Id^2$, 这样$f*g=Id^3$, 总的复杂度仍然是$O(n^{\frac{2}{3}})$

const int N2 = 5e6+10, N3 = 2500;
int P, inv6, inv2, cnt;
int phi[N2], vis[N2], p[N2];
int S1[N2], S2[N3], v[N3];
ll N;

ll inv(ll x){return x<=1?1:inv(P%x)*(P-P/x)%P;}
void init(int n) {
	phi[1] = 1;
	REP(i,2,n-1) {
		if (!vis[i]) p[++cnt]=i,phi[i]=i-1;
		for (int j=1,t;j<=cnt&&i*p[j]<=n; ++j) {
			vis[t=i*p[j]]=1;
			if (i%p[j]==0) {phi[t]=phi[i]*p[j];break;}
			phi[t]=phi[i]*phi[p[j]];
		}
	}
	REP(i,1,n) S1[i]=(S1[i-1]+(ll)i*i%P*phi[i]%P)%P;
}
ll sqr(ll n) {n%=P;return n*n%P;}
ll s_1(ll n) {n%=P;return n*(n+1)%P*inv2%P;}
ll s_2(ll n) {n%=P;return n*(n+1)%P*(2*n+1)%P*inv6%P;}
ll s_3(ll n) {n%=P;return sqr(s_1(n));}
int sum(ll n) {
	if (n<N2) return S1[n]; int x = N/n;
	if (v[x]) return S2[x]; v[x] = 1;
	int &ans = S2[x] = s_3(n), mx = sqrt(n);
	REP(i,2,mx) (ans-=sqr(i)*sum(n/i)%P)%=P;
	for (ll i=mx+1,j,k=n/i; i<=n; i=j+1,--k) {
		(ans-=(s_2(j=n/k)-s_2(i-1))*sum(k)%P)%=P;
	}
	return ans;
}

int main() {
	scanf("%d%lld", &P, &N);
	inv2 = inv(2), inv6 = inv(6);
	init(N2-1);
	int mx = sqrt(N), ans = 0;
	REP(i,1,mx) (ans+=(S1[i]-S1[i-1])*s_3(N/i)%P)%=P;
	for (ll i=mx+1,j,k=N/i; i<=N; i=j+1,--k) {
		(ans+=(sum(j=N/k)-sum(i-1))*s_3(k)%P)%=P;	
	}
	if (ans<0) ans+=P;
	printf("%d\n", ans);
}

练习4. luogu P3172 [CQOI2015]选数

大意: 求选择$n$个范围[L,R]的数, 且$gcd=k$的方案数

相当于求$\sum\limits_{L\le a_1 \le R}\cdots\sum\limits_{L\le a_n \le R}[gcd=k]$

等价于$\sum\limits_{l\le a_1 \le r}\cdots\sum\limits_{l\le a_n \le r}[gcd=1]$, 其中$l=\lfloor \frac{L-1}{K}\rfloor+1, r=\frac{R}{K}$

再反演一下就没了.....$\sum\limits_{d=1}^{r}\mu(d)(\lfloor \frac{r}{d}\rfloor-\lfloor \frac{l-1}{d} \rfloor)^n$

这题因为有$\lfloor \frac{r}{d}\rfloor$和$\lfloor \frac{l-1}{d} \rfloor$两部分, 所以为了方便直接用map记忆化了

这题刚开始想错了, 反演成$\varphi$了, 结果改了半天................

const int N2 = 2e6+10, N3 = 2000, P = 1e9+7;
int mu[N2], vis[N2], p[N2];
int n, k, l, r, cnt;
unordered_map<int,int> S2;
ll qpow(ll a,ll n) {ll r=1%P;for (a%=P;n;a=a*a%P,n>>=1)if(n&1)r=r*a%P;return r;} 
void init(int n) {
    mu[1] = 1;
    REP(i,2,n-1) {
        if (!vis[i]) p[++cnt]=i,mu[i]=-1;
        for (int j=1,t;j<=cnt&&i*p[j]<=n; ++j) {
            vis[t=i*p[j]]=1;
            if (i%p[j]==0) {mu[t]=0;break;}
            mu[t] = -mu[i];
        }
        mu[i] += mu[i-1];
    }
}
int sum(int n) {
    if (n<N2) return mu[n];
    if (S2.count(n)) return S2[n];
    int ans = 1, mx = sqrt(n);
    REP(i,2,mx) (ans-=sum(n/i))%=P;
    for (int i=mx+1,j,k=n/i; i<=n; i=j+1,--k) {
        (ans-=((j=n/k)-i+1)*sum(k)%P)%=P;
    }
    return S2[n]=ans;
}

int main() {
    scanf("%d%d%d%d", &n, &k, &l, &r);
    l = (l-1)/k, r /= k;
    init(N2-1);
    int ans = 0;
    for (int i=1,j; i<=r; i=j+1) {
        j = r/(r/i);
        if (l/i) j = min(j,l/(l/i));
        (ans+=(sum(j)-sum(i-1))*qpow(r/i-l/i,n)%P)%=P;
    }
    if (ans<0) ans+=P;
    printf("%d\n", ans);
}

练习5. Luogu P3327 [SDOI2015]约数个数和

大意: 求$\sum\limits_{i=1}^n\sum\limits_{j=1}^md(ij)$, $d(n)$为除数函数

对于除数函数可以这样处理$d(ij)=\sum\limits_{x|i}\sum\limits_{y|i}[x\perp y]=\sum\limits_{x|i}\sum\limits_{y|i}\sum\limits_{d|gcd(x,y)}\mu(d)$

所以就可以得到$\sum\limits_{i=1}^n\sum\limits_{j=1}^md(ij)=\sum\limits_{d=1}^{min(n,m)}\varphi(d)\sum\limits_{i=1}^{\lfloor \frac{n}{d}\rfloor}\lfloor \frac{n}{id} \rfloor \sum\limits_{j=1}^{\lfloor \frac{m}{d}\rfloor}\lfloor \frac{m}{jd} \rfloor$

然后就没了....主要这题的数据范围太小, 不需要杜教筛也行

 

转载于:https://www.cnblogs.com/uid001/p/10563110.html

《英雄编程体验课》第 15 课 | C语言中的数学库
英雄哪里出来
07-05 5万+
C语言中的数学库
莫比乌斯反演
长颜草的博客
08-03 815
莫比乌斯反演莫比乌斯反演是数论数学中很重要的内容，可以用于解决很多组合数学的问题。莫比乌斯函数莫比乌斯函数，数论函数，由德国数学家和天文学家莫比乌斯提出。梅滕斯首先使用μ(n)作为莫比乌斯函数的记号。莫比乌斯函数是指以下的函数：在这里，λ(n)是刘维尔函数莫比乌斯函数是一个数论函数，它同时也是一个积极函数（μ(ab) =μ(a)μ(b), a,b互质 ）当n不等于1时，n的所有因子的莫比乌斯函数值的
莫比乌斯反演习题总结
weixin_30732487的博客
02-24 483
一：常用卷积 \[1 * \mu = \epsilon\] \[Id * \mu = \varphi\] \[1 * \varphi = Id \] 最后一个是欧拉反演. 根据这些卷积之间的相互转化，在推式子时只要够敏感，就能巧妙地化简。（如P3768简单的数学题） 二： 推式子的技巧 1. 改变枚举对象、枚举顺序 比如 \[ \begin{aligned} 1. &\sum_{i=1}^...
一道莫比乌斯反演练习题
xiayuyang
01-07 298
题意 n<=1e9n<=1e9n<=1e9,答案对2322^{32}232取膜 部分分n<=1e6n<=1e6n<=1e6,n<=1e7n<=1e7n<=1e7 首先发现mu(i∗j)=mu(i)∗mu(j)∗[gcd(i,j)==1]mu(i*j)=mu(i)*mu(j)*[gcd(i,j)==1]mu(i∗j)=mu(i)∗mu(j)∗[...
【luogu】简单的数学题（莫比乌斯反演、杜教筛）
xyyxyyx的博客
07-18 209
题意 给出nnn和p,(n≤1010)p,(n \le 10^{10})p,(n≤1010)，求 ∑i=1n∑j=1ni⋅j⋅gcd(i,j)&ThickSpace;mod&ThickSpace;p\sum_{i=1}^{n}\sum_{j=1}^{n}i \cdot j \cdot gcd(i,j) \; mod \; pi=1∑n​j=1∑n​i⋅j⋅gcd(i,j)modp ...
2018蓝桥杯C语言b组国赛 6.矩阵求和 欧拉函数+莫比乌斯反演
Nicolas的博客
04-27 719
第六题 标题：矩阵求和 经过重重笔试面试的考验，小明成功进入 Macrohard 公司工作。 今天小明的任务是填满这么一张表： 表有 n 行 n 列，行和列的编号都从1算起。 其中第 i 行第 j 个元素的值是 gcd(i, j)的平方， gcd 表示最大公约数，以下是这个表的前四行的前四列： 1 1 1 1 1 4 1 4 1 1 9 1 1 4 1 16 小明突然冒出一个奇怪的想法，他...
【洛谷 P3768】简单的数学题（莫比乌斯反演+杜教筛）
Stargazer的博客
01-22 284
传送门 首先简单莫比乌斯反演可以得到 ∑i=1n∑j=1nijgcd(ij)\sum_{i=1}^{n}\sum_{j=1}^nijgcd(ij)i=1∑n​j=1∑n​ijgcd(ij) =∑d=1nd3∑k=1ndμ(k)k2S(nkd)2,S(n)=∑i=1ni=\sum_{d=1}^nd^3\sum_{k=1}^{\frac nd }\mu(k)k^2S(\frac n {kd})^2,S...
莫比乌斯反演一些式子的推导(附题目)
流水
04-10 1429
前置内容 1.a∣ba|ba∣b意为aaa整除bbb 2.[a=b]={1,a=b0,a≠b[a=b]=\begin{cases} 1,a=b\\0,a ≠b \end{cases}[a=b]={1,a=b0,a​=b​ 3.⌊a⌋\lfloor a \rfloor⌊a⌋意为对aaa向下取整 4.当n>k,n>dn>k,n>dn>k,n>d，有⌊⌊n/d⌋k⌋...
莫比乌斯反演公式与数论函数：数论中的隐藏力量
本文首先介绍了莫比乌斯反演公式的概念以及数论的基础知识，然后深入探讨了莫比乌斯函数的定义与性质，并阐述了莫比乌斯变换及其在数论函数中的应用。接着，本文拓展了莫比乌斯反演公式的理论，分析了其在组合数学中...
牛客练习赛76 F phi and phi（莫比乌斯反演）
qq_51682942的博客
08-15 386
牛客练习赛76 F phi and phi（莫比乌斯反演） 题解： 前置知识：前置知识：前置知识：∑i=1n∑j=1nϕ(ij)=∑i=1n∑j=1nϕ(i)ϕ(j)gcd(i,j)ϕ(gcd(i,j))\sum_{i=1}^{n}\sum_{j=1}^{n}\phi(ij)=\sum_{i=1}^{n}\sum_{j=1}^{n}\frac{\phi(i)\phi(j)gcd(i,j)}{\phi(gcd(i,j))}i=1∑n​j=1∑n​ϕ(ij)=i=1∑n​j=1∑n​ϕ(gcd(i,j))ϕ(
积性函数大全（欧拉函数、莫比乌斯反演、杜教筛……）
LZH的博客
09-11 685
fsdjgskg
莫比乌斯反演定理
海边拾贝，沧海一粟
06-02 496
[toc] // 先挖坑，有空填 莫比乌斯反演定理 百度百科 我感觉百度百科就讲的不错
【JAVA 课后习题 12.5】+ 自定义异常类
楚江枫
11-29 1607
package He;import java.util.Scanner;public class He { public static void main(String[] args) { Scanner input = new Scanner(System.in); System.out.print("Enter a hex number : ");
初涉莫比乌斯反演（附带例题）
热门推荐
litble的成(tui)长(fei)史
05-29 1万+
什么是莫比乌斯反演 关于莫比乌斯反演 莫比乌斯反演，又称懵逼钨丝繁衍，是一种看了就一脸懵逼的东西。 好吧好吧，严肃点。（不过因为本蒟蒻真的很懵逼，所以错误之处请大神指出） 莫比乌斯反演就是下面这个式子： 如果存在函数F(x)和f(x),满足 F(n)=∑d|nf(d)F(n)=∑d|nf(d)F(n)=\sum_{d|n} f(d) 那么就有： f(n)=∑d|nμ(d)F(n...
Excel表格通用模板：2000平方宾馆装修预算表.xls
09-10
Excel表格通用模板：2000平方宾馆装修预算表.xls
preludedb-tools-5.2.0-1.el8.tar.gz
09-10
# 适用操作系统：Centos8 #Step1、解压 tar -zxvf xxx.el8.tar.gz #Step2、进入解压后的目录，执行安装 sudo rpm -ivh *.rpm
计算机教学年度工作总结.doc
09-10
计算机教学年度工作总结.doc
网络推广方案：知己知彼-了解行业数据.doc
最新发布
09-10
网络推广方案：知己知彼-了解行业数据.doc
莫比乌斯反演详解与应用
“莫比乌斯反演-莫比乌斯反演经典讲稿” 莫比乌斯反演是一种在数论中用于解决计数问题的技巧，它由德国数学家莫比乌斯提出。在数论和组合数学中，它常常用于简化计算，特别是在处理关于数的乘积的问题时。讲稿来自...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值