BZOJ1101 POI2007 Zap
来源:互联网 发布:微易企业号源码 编辑:程序博客网 时间:2024/06/07 04:55
BZOJ1101 POI2007 Zap
Time Limit: 10 Sec Memory Limit: 162 MB
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
Sample Output
3
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。
#include<iostream>#include<cstring>#include<cstdio>#include<cmath>using namespace std;#define LL long long#define MAXN 50010LL prime[MAXN],mu[MAXN];int tot=0,not_prime[MAXN],n,m,k,T;void Liner_Shaker(){//线性筛miu函数 mu[1]=1; for(int i=2;i<=50000;i++){ if(!not_prime[i]){ mu[i]=-1; prime[++tot]=i; } for(int j=1;j<=tot&&i*prime[j]<=50000;j++){ not_prime[i*prime[j]]=true; if(i%prime[j]==0){ mu[i*prime[j]]=0; break; } mu[i*prime[j]]=-mu[i]; } } for(int i=1;i<=50000;i++) mu[i]=mu[i]+mu[i-1];}LL Calc(int n,int m,int k){ LL ret=0; int last; n/=k;m/=k; for(int i=1;i<=m&&i<=n;i=last+1){ last=min(n/(n/i),m/(m/i)); ret+=(mu[last]-mu[i-1])*(m/i)*(n/i); } return ret;}int main(){ scanf("%d",&T); Liner_Shaker(); while(T--){ int n,m,k; scanf("%d%d%d",&n,&m,&k); printf("%lld\n",Calc(n,m,k)); } return 0;}
阅读全文
0 0
- 【bzoj1101】[POI2007]Zap
- 【POI2007】【BZOJ1101】Zap
- 【BZOJ1101】【POI2007】Zap
- Zap[Poi2007][bzoj1101]
- 【POI2007】bzoj1101 Zap
- 【bzoj1101】[POI2007]Zap
- bzoj1101: [POI2007]Zap
- BZOJ1101: [POI2007]Zap
- bzoj1101 [POI2007]Zap
- BZOJ1101 POI2007 Zap
- BZOJ1101 [POI2007]Zap
- 【BZOJ1101】[POI2007]Zap
- 【BZOJ1101】[POI2007]Zap【Mobius函数】
- BZOJ1101——[POI2007]Zap
- bzoj2045: 双亲数&bzoj1101: [POI2007]Zap
- bzoj1101/洛谷P3455 [POI2007]ZAP-Queries
- bzoj1101: [POI2007]Zap 莫比乌斯反演
- BZOJ1101: [POI2007]Zap 莫比乌斯反演
- 七参数对不同坐标系统的转换
- 【腾讯TMQ】Tensorflow 的 word2vec详细解释——basic篇
- 常用命令与优化思路
- sdut 字符统计2
- Activity设置成Dialog样式
- BZOJ1101 POI2007 Zap
- Android 沉浸式 完全沉浸
- Android开发中的常见异常
- unity多线程方案
- nginx-rtmp源码概述
- 利用进程快照枚举当前Windows运行进程的信息
- 网站添加微信登录功能
- -> 与 .的区别,&与*
- 使用grade查看某个库的版本,以及这个库被谁依赖