题目类型:莫比乌斯反演 传送门:>Here< 题意:求有多少对正整数对\((a,b)\),满足\(0<a<A\),\(0<b<B\),\(gcd(a,b)=d\) 解题思路 学了莫比乌斯反演,就以这道题来介绍一下莫比乌斯反演的题的应用(下文中,对数表示在规定范围内满足特定条件的数对数量,不是\(log\)的那个对数) 一般碰到有关\(gcd\)的题,一般地,设\(f(n)\)表示\(gcd=n\)的对数,\(F(n)\)表示\(n|gcd\)的对数 根据定义,满足\[F…
题意 有一张 \(n\times m\) 的数表,其第\(i\)行第\(j\)列的数值为能同时整除\(i\)和\(j\)的所有自然数之和. \(T\)组数据,询问对于给定的 \(n,m,a\) , 计算数表中\(\leq a\) 的数之和. \(T \leq 2\times 10^4,1 \leq n,m\leq 10^5\). 题解 令\(\sigma(x)\)表示\(x\)的约数和,容易写出答案的式子: \[\sum_{i=1}^n \sum_{j=1}^m\sigma(\gcd(i,j))…
题意 设\(d(x)\)为\(x\)的约数个数,求\(\sum_{i=1}^{n}\sum_{j=1}^{m}d(ij)\). 题解 首先证个公式: \[d(ij) = \sum_{x|i}\sum_{y|j} [gcd(x,y)=1]\] 可以这么考虑:利用唯一分解定理把\(i,j\)分解,即: $i=\prod_{k = 1}^{m} p_k^{c_k},j=\prod_{k=1}^m p_k^{d_k} $ 那等式左边显然为\(\prod(c_k+d_k+1)\), 然后考虑等式右边在干什…
一个结论:(从二维扩展来的,三维也是对的,证明可以考虑质因数分解) \[ d(ijk)=\sum_{i'|i}\sum_{j'|j}\sum_{k'|k}[\gcd(i',j')=1][\gcd(i', k')=1][\gcd(j', k')=1] \] \[ \sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^c\sum_{i'|i}\sum_{j'|j}\sum_{k'|k}[\gcd(i',j')=1][\gcd(i', k')=1][\gcd(j', k')=1] \]…
[Luogu3455][POI2007]ZAP-Queries(莫比乌斯反演) 题面 题目描述 FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd(x,y)=d.作为FGD的同学,FGD希望得到你的帮助. 输入输出格式 输入格式: The first line of the standard input contains one integer nn (1\le n\le 50 0001≤n≤50 000),de…
\(\mathcal{Preface}\)   单位根反演,顾名思义就是用单位根变换一类式子的形式.有关单位根的基本概念可见我的这篇博客. \(\mathcal{Formula}\)   单位根反演的公式很简单: \[[k|n]=\frac{1}k\sum_{i=0}^{k-1}\omega_k^{ni} \] \(\mathcal{Proof}\)   分类讨论: \(k|n\). 那么 \((\forall i)(\omega_k^{ni}=1)\),所以右侧为 \(\frac{1}k\su…
题目: http://www.lydsy.com/JudgeOnline/problem.php?id=1101 题解: http://www.cnblogs.com/mrha/p/8203612.html 数学公式太难打了,核心思想是化成gcd(i,j)==1,然后用莫比乌斯反演变成枚举约数d,然后再搞式子 #include<cstdio> #include<algorithm> #define N 50005 typedef long long ll; using namesp…
link ms是莫比乌斯反演里最水的题... 题意:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd(x,y)=d. 多组询问, T<=50000,d,a,b<=50000 稍微推下shizi \(\sum_{i=1}^a\sum_{j=1}^b[\gcd(i,j)=d]\) \(=\sum_{i=1}^{a/k}\sum_{j=1}^{b/k}[\gcd(i,j)=1]\) \(=\sum_{i=1}^{a/k}\sum_{j=1}^{b/k}\…
目录 Preface 数论函数 积性函数 Dirichlet 卷积 Dirichlet 卷积中的特殊函数 Mobius 函数 & Mobius 反演 Mobius 函数 Mobius 反演 基础应用 约数个数 欧拉函数 反演魔法 例一 例二 例三 魔法中的 tricks 线性筛 trick 筛 筛 筛 刷表 trick Conclusion   UPD:修改了 Euler 筛法代码框架.   若无特别说明,\(x,y\) 等形式变量均 \(\in\mathbb N_+\):\(p\) 为素数.…
[BZOJ2301][HAOI2011]Problem B(莫比乌斯反演) 题面 Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数n,接下来n行每行五个整数,分别表示a.b.c.d.k Output 共n行,每行一个整数表示满足要求的数对(x,y)的个数 Sample Input 2 2 5 1 5 1 1 5 1 5 2 Sample Outp…