题意 求\(\sum_{i=1}^{n} \sum_{j=1}^{m} lcm(i, j)^{gcd(i, j)}\)(\(n, m<=500000\)) 分析 很显然要死推莫比乌斯 题解 设\(n \le m\) \[ \begin{aligned} ans & = \sum_{i=1}^{n} \sum_{j=1}^{m} lcm(i, j)^{gcd(i, j)} \\ & = \sum_{i=1}^{n} \sum_{j=1}^{m} (\frac{ij}{gcd(i, j)…

3561: DZY Loves Math VI Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 205  Solved: 141 Description 给定正整数n,m.求 Input 一行两个整数n,m. Output 一个整数,为答案模1000000007后的值. Sample Input 5 4 Sample Output 424 HINT 数据规模: 1<=n,m<=500000,共有3组数据. Source By Jcvb [分析]…

3309: DZY Loves Math Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b). Input 第一行一个数T,表示询问数. 接下来T行,每行两个数a,b,表示一个询问. Output 对于每一个询问,输出一行一个非负整数作为回答. Sample In…

题意 \(T(T \le 10000)\)次询问,每次给出\(a, b(1 \le a, b \le 10^7)\),求 \[\sum_{i=1}^{a} \sum_{j=1}^{b} f((i, j))\] 其中\(f(n)\)表示\(n\)所含质因子的最大幂指数.\(f(1)=0\). 分析 以下默认\(a \le b\) $$\begin{align} & \sum_{i=1}^{a} \sum_{j=1}^{b} f((i, j)) \\= & \sum_{d=1}^{a} f(d…

BZOJ 3561 DZY Loves Math VI 求\(\sum_{i=1}^{n}\sum_{j=1}^{m}\text{lcm}(i,j)^{\gcd(i,j)}\),钦定\(n\leq m\) \(\sum_{i=1}^{n}\sum_{j=1}^{m}(\frac{ij}{{\gcd(i,j)}})^{\gcd(i,j)}\) 按套路,提出\(\gcd(i,j)\),枚举的\(i\),\(j\)都除\(g\) \(\sum_{g=1}^ng^g\sum_{i=1}^{n/g}\su…

题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3561 题解: 莫比乌斯反演 $$\begin{aligned}ANS&=\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)^{gcd(i,j)}\\&=\sum_{g=1}^{min(n,m)}\sum_{i=1}^{\frac{n}{g}}\sum_{j=1}^{\frac{m}{g}}g^gi^gj^g[gcd(i,j)==1]\\&=\sum_{g…

推到了一个推不下去的形式,然后就不会了 ~ 看题解后傻了:我推的是对的,推不下去是因为不需要再推了. 复杂度看似很大,但其实是均摊 $O(n)$ 的,看来分析复杂度也是一个能力啊 ~ code: #include <bits/stdc++.h> #define ll long long #define N 500006 #define mod 1000000007 #define setIO(s) freopen(s".in","r",stdin) us…

题意 \(m * m\)的网格,有\(n\)个点.\(t\)个询问:操作一:第\(x\)个点向四个方向移动了\(d\)个单位.操作二:询问同行同列其他点到这个点的曼哈顿距离和.强制在线.(\(n \le 10^5,m \le 10^{18}\)) 分析 没啥好分析的,就是推一下能推出每行每列的一个式子来,然后套两个区间维护的结构就行了. 题解 set + 线段树 #include <bits/stdc++.h> using namespace std; typedef long long ll…

[BZOJ3561]DZY Loves Math VI (数论) 题面 BZOJ 题解 \[\begin{aligned} ans&=\sum_{i=1}^n\sum_{j=1}^m\sum_{d=1}^n[gcd(i,j)=d](\frac{ij}{d})^d\\ &=\sum_{d=1}^nd^d\sum_{i=1}^{n/d}\sum_{j=1}^{m/d}[gcd(i,j)=1]i^dj^d\\ &=\sum_{d=1}^nd^d\sum_{i=1}^{n/d}\sum_{…

3560: DZY Loves Math V Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 241  Solved: 133 Description 给定n个正整数a1,a2,…,an,求 的值(答案模10^9+7). Input 第一行一个正整数n. 接下来n行,每行一个正整数,分别为a1,a2,…,an. Output 仅一行答案. Sample Input 3 6 10 15 Sample Output 1595 HINT 1<=n<=1…