马上开学了,加一个操作系统和数据库标签 不玩了,求1-n和1-m的lcm(i,j)和 首先想到把lcm(i,j)转化为i * j / gcd(i, j) 然后gcd,要素察觉,开始枚举d使得gcd(i,j) == d 之后莫比乌斯反演常规套路是miu(d),然后分母剩下d,分子剩下gcd(i,j) == 1,提出d我们得到d2 这时候柿子是µ(d) * d2 * ∑∑ij 其中i和j为除d,可以用等差数列求和公式O(1)求解 之后我们就好办了,枚举d 然后d2通过前缀和在线性筛里面处理,这一步可

题意:求\(\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)==1]\) \(assume\ n<m\) \(\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)==1]\) \(\Longrightarrow \sum_{x=1}^{n}\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x|gcd(i,j)}\mu(x)\) \(\Longrightarrow \sum_{x=1}^{n}\sum_{i=1}^{n}\sum_{

题目描述 求$\sum\limits_{i=1}^n\sum\limits_{j=1}^m|\mu(gcd(i,j))|lcm(i,j)$,即$gcd(i,j)$不存在平方因子的$lcm(i,j)$之和. 输入 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N.M 输出 T行 每行一个整数 表示第i组数据的结果 样例输入 4 2 4 3 3 6 5 8 3 样例输出 24 28 233 178 题解 莫比乌斯反演+线性筛 (为了方便,以下公式默认$n\le m$) $\ \ \ \

一个经典问题 求 \[ \sum_{k=1}^n\mathbb{lcm}(k,n) \] 一般的做法是使用\(\varphi(n)\)函数. 不经典的做法 \[ \begin{align*} \sum_{k=1}^n\mathbb{lcm}(k,n) &=\sum_{k=1}^n\frac{nk}{\gcd(n,k)}\tag{gcd与lcm的关系}\\ &=\sum_{d|n}\sum_{k=1}^{n/d}\frac{ndk}{d}[\gcd(n,dk)=d]\tag{枚举gcd(n,

题目大意: 表格中每一个位置(i,j)填的值是lcm(i,j) , 求n*m的表格值有多大 论文贾志鹏线性筛中过程讲的很好 最后的逆元我利用的是欧拉定理求解的 我这个最后线性扫了一遍,勉强过了,效率不是很高... /*bzoj 2154*/ #include <bits/stdc++.h> using namespace std; #define ll long long #define N 10000000 ; ] , prime[N+] , f[N+] , tot; ]; void ini

BZOJ 2154 crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如,LCM(6, 8) = 24.回到家后,Crash还在想着课上学的东西,为了研究最小公倍数,他画了一张N*M的表格.每个格子里写了一个数字,其中第i行第j列的那个格子里写着数为LCM(i, j).一个4*5的表格如下: 1 2 3 4 5 2 2 6 4

题意:求\(\sum\limits_{i=1}^n \sum\limits_{j=1}^m lcm(i,j)\ : gcd(i,j) 是sf 无平方因子数\) 无平方因子数?搞一个\(\mu(gcd(i,j))\)不就行了..不对不对有正负,是\(\mu^2\)才行 套路推♂倒 (ﾉ*･ω･)ﾉ \[ \begin{align*} \sum\limits_{i=1}^n \sum\limits_{j=1}^m \frac{ij}{gcd(i,j)} \mu(gcd(i,j))^2 &=\sum_

4833: [Lydsy1704月赛]最小公倍佩尔数 Time Limit: 8 Sec  Memory Limit: 128 MBSubmit: 240  Solved: 118[Submit][Status][Discuss] Description 令(1+sqrt(2))^n=e(n)+f(n)*sqrt(2),其中e(n),f(n)都是整数,显然有(1-sqrt(2))^n=e(n)-f(n)*sqrt(2).令g( n)表示f(1),f(2)…f(n)的最小公倍数,给定两个正整数n和

BZOJ2694 Lcm Description Input 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N.M Output T行 每行一个整数 表示第i组数据的结果 Sample Input 4 2 4 3 3 6 5 8 3 Sample Output 24 28 233 178 HINT T <= 10000 N, M<=4000000 文章链接:https://www.cnblogs.com/dream-maker-yk/p/9676383.html #include&

题意:给出n [1,3*1e6] 求 并模2^64. 思路:先手写出算式 观察发现可以化成 那么关键在于如何求得i为1~n的lcm(i,n)之和.可以知道lcm(a,b)为ab/gcd(a,b) 变换得(a/gcd) (b/gcd)gcd 由于GCD的性质,可以知道a/gcd 与 b/gcd是互质的两个质数.由此可以想到应用欧拉函数,并且由性质能够证明 nphi(n)/2为小于n所有与n互质数之和(证明:已知一个质数p那么显然n-p与它互质,那么phi(n)中有phi(n)/2对数,每对数和为n

原题 定义整数a,b,求所有满足条件的lcm(a,b)的和: 1<=a<=A 1<=b<=B ∀n>1,n2†gcd(a,b)(即任意n>1,\(n^2\)不是gcd(a,b)的约数) 输出答案对2^30取模. 要求gcd(a,b)不能含平方因子,所以gcd(a,b)一定是mu不等于0的数. 那么我们设所有满足条件的数为p 其余与bzoj 2693是一样的,推倒见这里! //敲公式累死了-- #include<cstdio> #include<algo

求 $\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)\mu(gcd(i,j))^2$   $\Rightarrow \sum_{d=1}^{n}\mu(d)^2\sum_{i=1}^{n}\sum_{j=1}^{m}\frac{ij}{d}[gcd(i,j)==d]$   $\Rightarrow \sum_{d=1}^{n}\mu(d)^2\frac{1}{d}\sum_{i=1}^{n}\sum_{j=1}^{m}ij[gcd(i,j)==d]$   $\Right

题目传送门 传送门I 传送门II 题目大意 设$F(n) = \sum_{i = 1}^{n}\sum_{j = 1}^{n}\left [ [i, j] + (i, j) \geqslant n \right ]$,求$\sum_{i = 1}^{n} F(i)$. 考虑设$f(n) = \sum_{i = 1}^{n}\sum_{i = 1}^{n}\left[ \left[i,j \right ] + (i, j) = n \right ]$,那么有$F(n) = F(n - 1) - f

2154: Crash的数字表格 Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 2924  Solved: 1091[Submit][Status][Discuss] Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如,LCM(6, 8) = 24.回到家后,Crash还在想着课上学的东西,为了研究

题意:求sigma{lcm(i,j)},1<=i<=n,1<=j<=m 不妨令n<=m 首先把lcm(i,j)转成i*j/gcd(i,j) 正解不会...总之最后化出来的莫比乌斯反演式子并没有除法- 本脑子有坑选手的做法:20101009是一个质数,而且n和m的范围小于20101009,这一定有其原因.经过仔细思考,我们发现这保证了每个1~n的数都有mod20101009意义下的乘法逆元.用inv[x]表示x的逆元,我们发现原先的式子等于sigma{inv[gcd(i,j)]

题目链接:http://www.spoj.com/problems/VLATTICE/ 题意:求gcd(x,y,z)=1,1<=x,y,z<=n,的个数. 开始做的时候枚举gcd(x,y),然后求z与gcd(x,y)互素的个数个数,O(n*sqrt(n))赌赌RP,然后TLE了... 后来才知道要用到莫比乌斯反演定理: 已知 f(n) = sigma(d|n, g(d)) 那么 g(n) = sigma(d|n, mu(d)*f(n/d)) 还有另一种形式更常用: 在某一范围内,已知 f(n)

传送门:GCD 题意:求[1,n],[1,m]gcd为k的对数. 分析:莫比乌斯入反演门题,gcd(x,y)==k等价于gcd(x/k,y/k)==1,求出[1,n][1,m]互质的对数,在减去[1,2][2,1]之类重复的个数即答案. 莫比乌斯反演资料: 贾志鹏线性筛 莫比乌斯反演入门 莫比乌斯反演介绍 莫比乌斯反演:46ms #pragma comment(linker,"/STACK:1024000000,1024000000") #include <cstdio>

[BZOJ2154]Crash的数字表格(莫比乌斯反演) 题面 BZOJ 简化题意: 给定\(n,m\) 求\[\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\] 题解 以下的一切都默认\(n<m\) 我们都知道\(lcm(i,j)=\frac{ij}{gcd(i,j)}\) 所以所求化简 \[\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}\] 看到\(gcd(i,j)\)很不爽,于是就再提出来 \[\sum_{d=1}^{n}\sum_

题链: http://www.lydsy.com/JudgeOnline/problem.php?id=2694 题解: 莫比乌斯反演 不难看出,造成贡献的(i,j)满足gcd(i,j)无平方因子. 其实也就是$\mu(gcd(i,j))!=0$ 先列出求ANS的式子 $\begin{align*}ANS&=\sum_{a=1}^{A}\sum_{b=1}^{B} lcm(a,b)\mu(gcd(a,b))^2\;(同样的,先枚举gcd的值g)\\&=\sum_{g=1}^{min(A,B

Description 求sigma lcm(x,y),x<=n,y<=m.n,m<=1e7. Solution lcm没有什么直接做的好方法,用lcm=x*y/gcd转成gcd来做 就是要求sigma d*f(x/d,y/d) f(x,y)为x和y以内gcd正好为1的对数 F为所有对数,于是有F(x,y)=x*(x+1)/2*y*(y+1)/2 f(x,y)=sigma (1<=i<=x) i*i*mu(i)*F(x/i,y/i) f用莫比乌斯反演解决,这两个式子都套上分块