bzoj2301【HAOI2011】Problem b

2301: [HAOI2011]Problem b

Time Limit: 50 Sec  Memory Limit: 256 MB

Submit: 2951  Solved: 1318

[Submit][Status][Discuss]

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

HINT

100%的数据满足:1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000

容斥原理+莫比乌斯反演

首先根据容斥原理,可以将一个询问转化成四个询问(a-1,c-1) (a-1,d) (b,c-1) (b,d),每次询问有多少个数对(x,y)满足1≤x≤n,1≤y≤m且gcd(x,y)=k。

这个问题等价于有多少个数对(x,y)满足1≤x小于等于(n/k),1≤y≤(m/k)且x与y互质。

这时候我们就可以考虑莫比乌斯反演了。

令f(i)表示1≤x≤n,1≤y≤m且gcd(x,y)=i的数对(x,y)的个数;g(i)表示1≤x≤n,1≤y≤m且i|gcd(x,y)的数对(x,y)的个数。

那么显然有g(i)=(n/i)*(m/i)。

根据莫比乌斯定理有f(i)=∑(i|d) mu(d/i)*g(d)。

这样枚举每一个k的倍数,就可以将每次询问做到O(n)。

但这还是有点慢,于是我们考虑进一步优化。我们维护莫比乌斯函数的前缀和,然后对于(n/i)和(m/i)相同的部分分块处理,最多有4*sqrt(n)块。所以单次询问复杂度就降到O(sqrt(n)),总复杂度为O(n*sqrt(n))。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 50005
using namespace std;
int t,a,b,c,d,k,tot;
int mu[maxn],sum[maxn],pri[maxn];
bool mark[maxn];
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
inline void getmu()
{
	mu[1]=1;
	F(i,2,50000)
	{
		if (!mark[i]){mu[i]=-1;pri[++tot]=i;}
		for(int j=1;j<=tot&&pri[j]*i<=50000;j++)
		{
			mark[i*pri[j]]=true;
			if (i%pri[j]==0){mu[i*pri[j]]=0;break;}
			else mu[i*pri[j]]=-mu[i];
		}
	}
	F(i,1,50000) sum[i]=sum[i-1]+mu[i];
}
inline int calc(int n,int m)
{
	if (n>m) swap(n,m);
	int ans=0,pos;
	for(int i=1;i<=n;i=pos+1)
	{
		pos=min(n/(n/i),m/(m/i));
		ans+=(sum[pos]-sum[i-1])*(n/i)*(m/i);
	}
	return ans;
}
int main()
{
	getmu();
	t=read();
	while (t--)
	{
		a=read();b=read();c=read();d=read();k=read();
		a--;c--;
		a/=k;b/=k;c/=k;d/=k;
		int ans=calc(a,c)+calc(b,d)-calc(a,d)-calc(b,c);
		printf("%d\n",ans);
	}
	return 0;
}
时间: 2024-10-29 16:39:57

bzoj2301【HAOI2011】Problem b的相关文章

【BZOJ2301】【HAOI2011】Problem b [莫比乌斯反演]

Problem b Time Limit: 50 Sec  Memory Limit: 256 MB[Submit][Status][Discuss] 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 Inp

BZOJ 2301 【HAOI2011】 Problem b

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)的个数 HINT 100%的数据满足:1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000 过了这么久终于写出了莫比乌斯反演的入门题TAT-- 这道

【BZOJ2298】【HAOI2011】problem a 动态规划

链接: #include <stdio.h> int main() { puts("转载请注明出处[vmurder]谢谢"); puts("网址:blog.csdn.net/vmurder/article/details/44979323"); } 题解: 一句话 (a,b) 可以理解成一个线段 (a,n?b] . 然后排个序去下重,最后一个线段的权值 x 就是表示 这 x 人互不冲突,一起算. 然后动态规划求若干条不相交线段的权值最大值,最后用总人数减

【HAOI2011】problem a

我的思维能力真的上来了好感动555 原题: 随便点的一道DP题,本来这题以前无数次想不出来,题解好像也看不懂,想的时候都打算放弃了 但是想着一定要死磕思维能力,最后居然还真的自己做出来了 果真是以前放弃得太轻易hhh 首先要确定一下题意,“有ai个人”我感觉也可以理解为>=ai个人 但是其实后来想一下,如果这样说得话就没必要强调可能有相同分数了 因为如果是>=ai的话,两个同分的和两个相邻分数的对于合法性其实都一样 所以认为题意为恰好有ai个人(当然考场上还是大力问监考) 那么研究一下这个约束

【HAOI2011】向量

题目描述: 给你一对数a,b,你可以任意使用(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), (-b,a), (-b,-a)这些向量,问你能不能拼出另一个向量(x,y). 说明:这里的拼就是使得你选出的向量之和为(x,y) 输入格式: 第一行数组组数t,(t<=50000) 接下来t行每行四个整数a,b,x,y (-2*10^9<=a,b,x,y<=2*10^9) 输出格式: t行每行为Y或者为N,分别表示可以拼出来,不能拼出来 分析: 我们看

【BZOJ2998】Problem A(动态规划)

[BZOJ2998]Problem A(动态规划) 题面 BZOJ 题解 一个人的成绩范围可以确定为一个区间 这样就变成了 选择若干区间,不重合, 每个区间有个权值,求最大权值和 这样就可直接\(dp\)了 #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #includ

【bzoj2298】【HAOI2011】【problem a】

2298: [HAOI2011]problem a http://www.lydsy.com/JudgeOnline/problem.php?id=2298 Time Limit: 10 Sec Memory Limit: 256 MB Submit: 696 Solved: 307 [Submit][Status][Discuss] Description 一次考试共有n个人参加,第i个人说:"有ai个人分数比我高,bi个人分数比我低."问最少有几个人没有说真话(可能有相同的分数)

【HAOI2011】【BZOJ2301】ProblemB

2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MB Submit: 1756 Solved: 755 [Submit][Status][Discuss] 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 Outpu

【bzoj2298】【HAOI2011】【problem a】【dp】

Description 一次考试共有n个人参加,第i个人说:"有ai个人分数比我高,bi个人分数比我低."问最少有几个人没有说真话(可能有相同的分数) Input 第一行一个整数n,接下来n行每行两个整数,第i+1行的两个整数分别代表ai.bi Output 一个整数,表示最少有几个人说谎 Sample Input 3 2 0 0 2 2 2 Sample Output 1 HINT 100%的数据满足: 1≤n≤100000 0≤ai.bi≤n 思路:其实我们可以很巧妙的把这道题转化