[HAOI2011]Problem b（莫比乌斯反演）

最新推荐文章于 2024-05-13 09:55:25 发布

何况虚度光阴

最新推荐文章于 2024-05-13 09:55:25 发布

阅读量474

点赞数

分类专栏：数论文章标签： c++ 算法

本文链接：https://blog.csdn.net/weixin_63128867/article/details/130132233

版权

数论专栏收录该内容

41 篇文章 0 订阅

订阅专栏

[HAOI2011]Problem b

题目链接：https://www.luogu.com.cn/problem/P2522

题目描述

对于给出的 $n$ 个询问，每次求有多少个数对 $(x, y)$ ，满足 $\le x \le b$ ， $\le y \le d$ ，且 $\gcd(x,y) = k$ ， $\gcd(x,y)$ 函数为 $x$ 和 $y$ 的最大公约数。

输入格式

第一行一个整数 $n$ ，接下来 $n$ 行每行五个整数，分别表示 $a, b, c, d, k$ 。

输出格式

共 $n$ 行，每行一个整数表示满足要求的数对 $(x, y)$ 的个数。

样例 #1

样例输入 #1

2
2 5 1 5 1
1 5 1 5 2

样例输出 #1

14
3

提示

对于 $100\%$ 的数据满足： $\le n,k \le 5 \times 10^4$ ， $\le a \le b \le 5 \times 10^4$ ， $\le c \le d \le 5 \times 10^4$ 。

思路：
在这里插入图片描述

在这里插入代码片
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const long long N=60000;
long long s[N],m[N],sum[N],f[N];
long long cnt;
long long a[N];int k;
inline void read(int &x)
{
    x=0;
    static int p;p=1;
    static char c;c=getchar();
    while(!isdigit(c)){if(c=='-')p=-1;c=getchar();}
    while(isdigit(c)) {x=(x<<1)+(x<<3)+(c-48);c=getchar();}
    x*=p;	
}
void unit(long long n){
	m[1]=1;
	for(long long i=2;i<=n;i++){
		if(!a[i]){
			s[++cnt]=i;
			m[i]=-1;
		}
		for(long long j=1;s[j]*i<=n&&j<=cnt;j++){
//			int op=s[j]*i;
			a[s[j]*i]=1;
			if(i%s[j]==0){
				break;
			}
			m[s[j]*i]=m[i]*-1;
		 }
	}
	for(long long i=1;i<=n;i++){
		sum[i]=sum[i-1]+m[i];
	}
}
long long calc(int a,int b)
{
    static int max_rep;
    static long long ans;
    max_rep=min(a,b);ans=0;
    for(int l=1,r;l<=max_rep;l=r+1)
    {
        r=min(a/(a/l),b/(b/l));
        ans+=(1ll*a/(1ll*l*k))*(1ll*b/(1ll*l*k))*(sum[r]-sum[l-1]);
    }
    return ans;
}
int main(){
	unit(50000);
	int t;
	read(t);
	while(t--){
		int x,y,xx,yy;
		read(x);
		read(y);
		read(xx);
		read(yy);
		read(k);
		printf("%lld\n",calc(y,yy)-calc(y,xx-1)-calc(x-1,yy)+calc(x-1,xx-1));
	}
	return 0;
}