1101: [POI2007]Zap

原创 2016年08月30日 21:54:50

1101: [POI2007]Zap

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 2172  Solved: 839
[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

Sample Output

3
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。

HINT

Source

[Submit][Status][Discuss]

莫比乌斯反演,此题不会爆int,切勿作死上long long
#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<bitset>
#include<algorithm>
#include<cstring>
#include<map>
#include<stack>
#include<set>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

typedef int LL;
const int maxn = 5E4 + 50;

int T,a,b,d,tot,mu[maxn],prime[maxn];
bool not_prime[maxn];

int getint()
{
	char ch = getchar();
	int ret = 0;
	while (ch < '0' || '9' < ch) ch = getchar();
	while ('0' <= ch && ch <= '9')
		ret = ret*10 + ch - '0',ch = getchar();
	return ret;
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	mu[1] = 1;
	for (int i = 2; i < maxn; i++) {
		if (!not_prime[i]) {
			prime[++tot] = i;
			mu[i] = -1;
		}
		for (int j = 1; j <= tot; j++) {
			int Nex = prime[j]*i;
			if (Nex >= maxn) break;
			not_prime[Nex] = 1;
			mu[Nex] = -mu[i];
			if (i % prime[j] == 0) {
				mu[Nex] = 0;
				break;
			}
		}
	}
	for (int i = 1; i < maxn; i++) mu[i] += mu[i-1];
	
	T = getint();
	while (T--) {
		a = getint(); b = getint(); d = getint();
		LL Ans = 0,A = a/d,B = b/d;
		if (A > B) swap(A,B);
		int last;
		for (int i = 1; i <= A; i = last + 1) {
			last = min(A/(A/i),B/(B/i));
			Ans += 1LL*(mu[last] - mu[i-1])*(A/i)*(B/i);
		}
		printf("%d\n",Ans);
	}
	return 0;
}

版权声明:这个人很懒什么都没有留下

相关文章推荐

莫比乌斯反演 bzoj1101【Poi2007】Zap

题目大意: 求对于给定的整数a,b和d,有多少正整数对x,y,满足x
  • Todobe
  • Todobe
  • 2017年03月06日 15:41
  • 389

BZOJ 1101: [POI2007]Zap

题目 题目描述 题解 题目解答 代码 题目代码
  • WNJXYK
  • WNJXYK
  • 2014年10月29日 20:25
  • 511

【POI2007】【BZOJ1101】Zap

DescriptionFGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x...

莫比乌斯函数(bzoj 1101: [POI2007]Zap)

莫比乌斯反演: 转载自http://blog.csdn.net/acdreamers/article/details/8542292

bzoj2045: 双亲数&bzoj1101: [POI2007]Zap

传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1101 思路:设a 求的就是f(a,b,d)=ΣΣ[gcd(i,j)==d](1 ...

[BZOJ 1101] POI 2007 Zap · 莫比乌斯 & 分块 超详细题解

初学莫比乌斯反演,翻了大量的题解才搞懂这题,所以决定自己写一个最详细的题解,虽然有些繁琐,但是每一步推导都十分详细。神犇就不要嘲讽我了2333 首先,我们定义 题目即要求 由于d是给定的,所以另 ...
  • ycdfhhc
  • ycdfhhc
  • 2016年02月05日 11:12
  • 858

【BZOJ】【P1101】【POI2007】【Zap】【题解】【莫比乌斯反演】

传送门:www.lydsy.com/JudgeOnline/problem.php?id=1101 shish

Zap[Poi2007][bzoj1101]

题目描述 DescriptionFGD FGD 正在破解一段密码,他需要回答很多类似的问题:对于给定的整数 a,b a,b 和 d d ,有多少正整数对 x,y x,y ,满足 x...

bzoj1101 [POI2007]Zap

Description FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x...

【POI2007】bzoj1101 Zap

狄利克雷卷积
  • sdfzyhx
  • sdfzyhx
  • 2017年01月13日 16:38
  • 141
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:1101: [POI2007]Zap
举报原因:
原因补充:

(最多只允许输入30个字)