关闭

1101: [POI2007]Zap

222人阅读 评论(0) 收藏 举报
分类:

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;
}

0
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

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

莫比乌斯反演: 转载自http://blog.csdn.net/acdreamers/article/details/8542292
  • Jaihk662
  • Jaihk662
  • 2017-07-30 16:41
  • 259

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

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

【bzoj1101】[POI2007]Zap 莫比乌斯反演

DescriptionFGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd(x,y)=d。作为FGD的同学,FGD希望得...
  • LOI_DQS
  • LOI_DQS
  • 2016-01-19 19:38
  • 407

golang有用的库及工具 之 zap.Logger包 使用

zap.Logger  是go语言中相对日志库中性能最高的。那么如何开始使用?不多说直接上代码: import ( "encoding/json" "fmt" "log" ...
  • niyuelin1990
  • niyuelin1990
  • 2017-10-25 13:32
  • 657

OWASP ZAP上手体验

新领到一个任务,试下OWASP ZAP。 工具说明          ZAproxy是一个易于使用交互式的用于web应用程序漏洞挖掘的渗透测试工具,可以检测sql注入,跨站...
  • chinajobs
  • chinajobs
  • 2016-12-07 11:36
  • 2464

莫比乌斯反演 bzoj1101【Poi2007】Zap

题目大意: 求对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd(x,y)=d。题目分析:(莫比乌斯反演) 莫比乌斯反演公式: 常见的另一种形式(也是比较常用...
  • Todobe
  • Todobe
  • 2017-03-06 15:41
  • 427

BZOJ1101 [POI2007]Zap 【莫比乌斯反演】

1101: [POI2007]ZapTime Limit: 10 Sec Memory Limit: 162 MB Submit: 2813 Solved: 1213 [Submit][Sta...
  • qq_38678604
  • qq_38678604
  • 2017-12-18 14:02
  • 39

OWASP ZAP上手体验

新领到一个任务,试下OWASP ZAP。 工具说明          ZAproxy是一个易于使用交互式的用于web应用程序漏洞挖掘的渗透测试工具,可以检测sql注入,跨站脚本,跨站请求伪造以及路...
  • u013224189
  • u013224189
  • 2015-10-30 17:53
  • 11105

OWASP ZAP 2.7.0 版本全球发布|棉花哥的博客

OWASP Zed Attack Proxy (ZAP)工具是世界上最受欢迎的免费安全工具之一。ZAP可以帮助安全测试人员在开发和测试应用程序过程中,自动发现 Web应用程序中的安全漏洞。另外,它也是...
  • qq_18209847
  • qq_18209847
  • 2017-12-07 09:11
  • 116

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

传送门:www.lydsy.com/JudgeOnline/problem.php?id=1101 shish
  • u012732945
  • u012732945
  • 2014-05-31 07:36
  • 1004
    个人资料
    • 访问:182254次
    • 积分:8756
    • 等级:
    • 排名:第2551名
    • 原创:733篇
    • 转载:1篇
    • 译文:0篇
    • 评论:16条
    文章分类
    最新评论