BZOJ 2301 Problem b(莫比乌斯反演)
来源:互联网 发布:tiny core linux 中文 编辑:程序博客网 时间:2024/06/05 03:43
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
2
2 5 1 5 1
1 5 1 5 2
Sample Output
14
3
Solution
设f(x,y,k)为1<=i<=x,1<=j<=y且gcd(i,j)=k的(i,j)对数,则有
进而ans=f(b,d,k)-f(a-1,d,k)-f(b,c-1,k)+f(a-1,c-1,k)
Code
#include<cstdio>#include<iostream>#include<cstring>#include<algorithm>using namespace std;#define maxn 55555typedef long long ll;bool check[maxn];int prime[maxn],mu[maxn],sum[maxn];void Moblus(int n){ memset(check,0,sizeof(check)); mu[1]=sum[1]=1; int 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]]=1; if(i%prime[j]==0) { mu[i*prime[j]]=0; break; } else mu[i*prime[j]]=-mu[i]; } sum[i]=sum[i-1]+mu[i]; }}ll solve(int n,int m,int k){ n=n/k,m=m/k; if(n>m)swap(n,m); ll ans=0; for(int i=1,next=0;i<=n&&i<=m;i=next+1) { next=min(n/(n/i),m/(m/i)); ans+=1ll*(n/i)*(m/i)*(sum[next]-sum[i-1]); } return ans;}int T,a,b,c,d,k;int main(){ Moblus(maxn-10); scanf("%d",&T); while(T--) { scanf("%d%d%d%d%d",&a,&b,&c,&d,&k); a--,c--; ll ans=solve(b,d,k)-solve(a,d,k)-solve(b,c,k)+solve(a,c,k); printf("%lld\n",ans); } return 0;}
0 0
- BZOJ 2301 Problem b (莫比乌斯反演)
- BZOJ 2301 Problem b(莫比乌斯反演)
- BZOJ 2301-Problem b(莫比乌斯反演)
- 【bzoj 2301】Problem b(莫比乌斯反演)
- [BZOJ 2301][HAOI2011] Problem b 莫比乌斯反演
- bzoj 2301: [HAOI2011]Problem b 【莫比乌斯反演】
- BZOJ 2301: [HAOI2011]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 2301: [HAOI2011]Problem b (莫比乌斯反演)
- BZOJ 2301 Problem b(莫比乌斯反演+容斥)
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演 + 容斥原理 + 分块优化)
- BZOJ 2301: [HAOI2011]Problem b(莫比乌斯反演,分块,容斥)
- BZOJ 2301 [HAOI2011]Problem b (莫比乌斯反演)
- PL/SQL database character set(AL32UTF8) and Client character set(ZHS16GBK) are different
- 记录1
- 移动端前端UI库—Frozen UI、WeUI、SUI Mobile
- [NOI2016]游记
- ScrollView 内嵌套 ListVIew/GrideView
- BZOJ 2301 Problem b(莫比乌斯反演)
- java 连接Mysql 中文查询不到结果的问题解决方案,已经字符串拼接方法
- 第三方支付
- 关于OOP和AOP的释义
- 记录2
- 车联网概念——转载自互联网
- hdu 5791 dp简单
- Mac OS X 下如何配置才能通过 sz/rz 命令下载/上传文件
- HDU 1576 扩展欧几里得&&暴力