HYSBZ 2301 Problem b 莫比乌斯 分块优化
来源:互联网 发布:sql语句查询重复数据 编辑:程序博客网 时间:2024/06/06 04:03
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
22 5 1 5 11 5 1 5 2
Sample Output
143
Hint
100%的数据满足:1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000
#include<iostream>#include<cstring>#include<cstdio>#include<ostream>#include<istream>#include<algorithm>#include<queue>#include<string>#include<cmath>#include<set>#include<map>#include<stack>#include<vector>#define fi first#define se second#define ll long long#define pii pair<int,int>#define inf (1<<30)#define eps 1e-8#define pb push_backusing namespace std;const int maxn=100010;int sum[maxn];int mu[maxn],prime[maxn],tot;bool check[maxn];void getmu(){ int N=100000; memset(check,0,sizeof(check)); mu[1]=1,tot=0; for(int i=2;i<=N;i++){ if(!check[i]) prime[tot++]=i,mu[i]=-1; for(int j=0;j<tot;j++){ if(i*prime[j]>N) break; check[i*prime[j]]=true; if(i%prime[j]==0){ mu[i*prime[j]]=0; break; } else mu[i*prime[j]]=-mu[i]; } }}int main(){ getmu(); int n,m; while(~scanf("%d%d",&n,&m)) { if(n>m) swap(n,m); ll ans=0; for(int i=1;i<=n;i++) { ll u=0; for(int j=i;j<=n;j+=i) { u=u+(ll)mu[j/i]*(n/j)*(m/j); } ans+=u*(2*(i-1)+1); } printf("%lld\n",ans); } return 0;}
0 0
- HYSBZ 2301 Problem b 莫比乌斯 分块优化
- HYSBZ 2301 Promblem b 莫比乌斯反演+分块
- [HYSBZ/BZOJ2301]Problem b [莫比乌斯反演+分块] 【组合数学】
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演 + 容斥原理 + 分块优化)
- 【莫比乌斯反演】[HYSBZ/BZOJ2301]Problem b
- BZOJ 2301 [HAOI2011]Problem b (容斥+莫比乌斯反演+分块优化 详解)
- [BZOJ 2301] Problem b【莫比乌斯反演/容斥原理/分块】
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演,分块,容斥)
- Bzoj 2301: [HAOI2011]Problem b(莫比乌斯反演+除法分块)
- BZOJ 2301 Problem b (莫比乌斯反演+区间容斥+分块)
- 【BZOJ】2301 [HAOI2011]Problem b && 【BZOJ】1101 [POI2007]Zap 莫比乌斯函数+数论分块
- 【BZOJ2301】【HAOI2011】Problem b 莫比乌斯反演+分块+前缀和
- 【bzoj2301】[HAOI2011]Problem b 莫比乌斯反演+线性筛法+数论分块
- [莫比乌斯反演+容斥+分块求和] BZOJ2301: [HAOI2011]Problem b
- 【莫比乌斯函数+除法分块】BZOJ2301(HAOI2011)[Problem b]题解
- CSU 1325 A very hard problem (莫比乌斯反演+分块求和优化)
- CSU 1325 A very hard problem (莫比乌斯反演+分块求和优化)
- bzoj 2301 -莫比乌斯函数反演+分块优化
- bzoj4102【Usaco2015 Open】Bessie
- android studio制作9-patch图片
- LM2596的肖特基二极管选择的注意点
- Scala模式匹配(一)
- 从源码角度深入理解Toast
- HYSBZ 2301 Problem b 莫比乌斯 分块优化
- C#入门经典
- 为Python配置Vim编辑器(GUI/非GUI皆可)
- 九度OJ 1136:Number Steps(步数) (基础题)
- 支付宝 免签约即时到账原理 ecshop等 订单状态返回
- JSP中的编码问题
- UVA Josephus Problem 11089 (约瑟夫环) 数学好题
- 九度OJ 1137:浮点数加法 (大数运算)
- android中drawable的各种state