[bzoj2301][HAOI2011]Problem b

题目大意：给你$a,b,c,d,k(1\leqslant a,b,c,d,k\leqslant 5\times10^4)$，求$\displaystyle\sum\limits_{x=a}^b\displaystyle\sum\limits_{y=c}^d[(x,y)==k]$

题解：下文中令$n\leqslant m$先考虑求$\sum\limits_{i=1}^n\sum\limits_{j=1}^m[(i,j)==k]$
$$
\begin{align*}
\def\dsum{\displaystyle\sum\limits}
令f(p)&=\dsum_{i=1}^n\dsum_{j=1}^m[(i,j)==p]\\
令F(p)&=\dsum_{p|k}f(k)\\
   &=\dsum_{p|k}\dsum_{i=1}^n\dsum_{j=1}^m[(i,j)==k]\\
   &=\dsum_{i=1}^n\dsum_{j=1}^m[p|(i,j)]\\
   &=\left\lfloor\dfrac{n}{p}\right\rfloor\cdot \left\lfloor\dfrac{m}{p}\right\rfloor\\
莫比乌&斯反演得：\\
\therefore f(p)&=\dsum_{p|k}\mu\Big(\dfrac{k}{p}\Big)F(k)\\
           &=\dsum_{i=1}^n\mu(i)\left\lfloor\dfrac{n}{ip}\right\rfloor\cdot\left\lfloor\dfrac{m}{ip}\right\rfloor\\
\end{align*}\\
令g(p)=\dsum_{i=1}^p\mu(i)\\
然后容斥一下就好了\\
$$

卡点：无

C++ Code：

#include <cstdio>
#define maxn 50010
using namespace std;
int miu[maxn], plist[maxn], ptot;
bool isp[maxn];
void sieve(int n) {miu[1] = 1;for (int i = 2; i <= n; i++) {if (!isp[i]) {miu[i] = -1;plist[ptot++] = i;}for (int j = 0; j < ptot, i * plist[j] <= n; j++) {int tmp = i * plist[j];isp[tmp] = true;if (i % plist[j] == 0) {miu[tmp] = 0;break;}miu[tmp] = -miu[i];}}for (int i = 2; i <= n; i++) miu[i] += miu[i - 1];
}
inline int min(int a, int b) {return a < b ? a : b;}
int solve(int n, int m, int k) {n /= k, m /= k;int tmp = min(n, m);int ans = 0, l, r;for (l = 1; l <= tmp; l = r + 1) {r = min(n / (n / l), m / (m / l));ans += (miu[r] - miu[l - 1]) * (n / l) * (m / l);}return ans;
}
int Tim, a, b, c, d, k;
int main() {sieve(50000);scanf("%d", &Tim);while (Tim --> 0) {scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);printf("%d\n", solve(b, d, k) - solve(a - 1, d, k) - solve(b, c - 1, k) + solve(a - 1, c - 1, k));}return 0;
}

转载于:https://www.cnblogs.com/Memory-of-winter/p/9523943.html

[bzoj2301][HAOI2011]Problem b相关推荐

BZOJ2301: [HAOI2011]Problem b
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MB Submit: 6435 Solved: 2986 [Submi ...
BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演容斥原理]【学习笔记】
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MB Submit: 4032 Solved: 1817 [Submi ...
BZOJ2301:[HAOI2011]Problem b(莫比乌斯反演,容斥)
Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数 ...
bzoj2301: [HAOI2011]Problem b懵逼乌斯反演
属于结果的和好求但是结果不好求的题 (轻易能得到以k的倍数为最大公约数的对数,但是不好直接求k) 所以一波反演结束其实反演的时候完全没有反演的感觉,就是不停地恒等变形算是懵逼乌斯反演最简单的例题 ...
BZOJ2301: [HAOI2011]Problem b(莫比乌斯反演)
Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数 ...
[bzoj2301] [HAOI2011]Problem b
Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数 ...
BZOJ 2301 [HAOI2011]Problem b
2301: [HAOI2011]Problem b Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函 ...
[POI2007]ZAP-Queries [HAOI2011]Problem b 莫比乌斯反演
1,[POI2007]ZAP-Queries ---题面--- 题解: 首先列出式子:$$ans = \sum_{i = 1}^{n}\sum_{j = 1}^{m}[gcd(i, j) == d]$ ...
[HAOI2011]Problem c
链接 P2523 [HAOI2011]Problem c 想法还是很巧妙的. 其实只是问一个先后顺序,因为编号相同的话,那么$id$小的就在前面,$id$大的就在后面. 所以我们考虑的是到底有 ...

[bzoj2301][HAOI2011]Problem b

[bzoj2301][HAOI2011]Problem b相关推荐

最新文章

热门文章