[BZOJ]1101 莫比乌斯+分块优化
来源:互联网 发布:中国蝉的种类知乎 编辑:程序博客网 时间:2024/05/20 06:22
1101: [POI2007]Zap
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2731 Solved: 1163
[Submit][Status][Discuss]
Description
FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a
,y<=b,并且gcd(x,y)=d。作为FGD的同学,FGD希望得到你的帮助。
Input
第一行包含一个正整数n,表示一共有n组询问。(1<=n<= 50000)接下来n行,每行表示一个询问,每行三个
正整数,分别为a,b,d。(1<=d<=a,b<=50000)
Output
对于每组询问,输出到输出文件zap.out一个正整数,表示满足条件的整数对数。
Sample Input
2
4 5 2
6 4 3
4 5 2
6 4 3
Sample Output
3
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。
HINT
Source
HOME Back
公式推导过程:
令
用莫比乌斯函数的性质把求和的式子换掉,
其中,更换求和指标,
容易知道单调不上升,且最多有种不同的取值。所以按取值分成个段分别处理,一个连续段内的和可以用预处理出的莫比乌斯函数前缀和求出.
分块技能get.
#include<stdio.h>#include<algorithm>using namespace std;const int maxn = 50005;bool mark[maxn];int n, m, d, tot, T;int pr[maxn], mu[maxn], smu[maxn];inline void Linear_sieve(){smu[1] = mu[1] = 1;for(int i = 2; i <= maxn; ++i){if(!mark[i]) pr[++tot] = i, mu[i] = -1;for(int j = 1; j <= tot && pr[j] * i <= maxn; ++j){mark[i * pr[j]] = true;if(!(i % pr[j])){mu[i * pr[j]] = 0;break;}mu[i * pr[j]] = -mu[i];}smu[i] = smu[i - 1] + mu[i];}}int main(){Linear_sieve();scanf("%d", &T);while(T--){scanf("%d%d%d", &n, &m, &d);n /= d, m /= d;int lim = min(n, m), dd = 0, ans = 0;for(int i = 1; i <= lim; i = dd + 1){dd = min(n / (n / i), m / (m / i));ans += (smu[dd] - smu[i - 1]) * (n / i) * (m / i);}printf("%d\n", ans);}}
阅读全文
0 0
- [BZOJ]1101 莫比乌斯+分块优化
- bzoj 2301 -莫比乌斯函数反演+分块优化
- BZOJ 1101 莫比乌斯函数+分块
- bzoj 2820-莫比乌斯函数反演+分块优化 +换元优化
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演 + 容斥原理 + 分块优化)
- 【BZOJ】2301 [HAOI2011]Problem b && 【BZOJ】1101 [POI2007]Zap 莫比乌斯函数+数论分块
- BZOJ 2154(莫比乌斯函数+分块)
- BZOJ 2301 莫比乌斯函数+分块
- BZOJ 2301 [HAOI2011]Problem b (容斥+莫比乌斯反演+分块优化 详解)
- [BZOJ 1101] POI 2007 Zap · 莫比乌斯 & 分块 超详细题解
- BZOJ 1101 Luogu P3455 POI 2007 Zap (莫比乌斯反演+分块)
- HYSBZ 2301 Problem b 莫比乌斯 分块优化
- hdu 4746 Mophues 莫比乌斯 分块优化
- [BZOJ 2301] Problem b【莫比乌斯反演/容斥原理/分块】
- BZOJ 2301 - 莫比乌斯反演 + 前缀和 + 分块计算
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演,分块,容斥)
- Bzoj 2820: YY的GCD(莫比乌斯反演+除法分块)
- Bzoj 2301: [HAOI2011]Problem b(莫比乌斯反演+除法分块)
- NOI 2002 机器人M号 欧拉函数
- opencv-python(1):读入、显示、保存图像
- 2015ACM/ICPC亚洲区沈阳站-重现赛(感谢东北大学) Meeting(SPAF+拆点)
- 表格
- 最大的矩形
- [BZOJ]1101 莫比乌斯+分块优化
- 这是一篇没有教科书那样晦涩语言的晶体管原理的介绍
- CF 868 A. Bark to Unlock 【简单匹配】
- POJ3093 Margaritas on the River Wal ( 背包DP)
- java中的类型赋值常见错误
- CF 868 B. Race Against Time【模拟】
- MySql学习笔记(一)
- 上上次考试总结(2017.9.23考试)
- 寻找数组中的最大值和最小值