luogu2522[HAOI2011]Problem b

最新推荐文章于 2023-01-10 19:36:46 发布

zz_ylolita

最新推荐文章于 2023-01-10 19:36:46 发布

阅读量148

点赞数

CC 4.0 BY-SA版权

分类专栏：数论莫比乌斯反演模板

本文链接：https://blog.youkuaiyun.com/zz_ylolita/article/details/82823055

数论同时被 3 个专栏收录

16 篇文章

订阅专栏

莫比乌斯反演

5 篇文章

订阅专栏

模板

5 篇文章

订阅专栏

本文介绍了一种使用容斥原理求解特定矩阵范围内问题的方法，通过预处理最小倍数函数并应用线性筛法优化质数检测，实现快速求解a<=x<=b,c<=y<=d范围内问题的算法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用一个简单的容斥就可以求出a<=x <=b, c<=y<=d范围的答案了

注意1LL呀qwq

#include <bits/stdc++.h>
using namespace std;
#define N 50005
bool vis[N];
int miu[N],prime[N];
typedef long long LL;
LL sum[N],ans;
int a,b,c,d,k,T,cnt;
void Miu()
{
	miu[1] = 1;
	cnt = 0;
	for (int i=2;i<=N-5;i++)
	{
		if (!vis[i]) 
		{
			prime[++cnt] = i;
			miu[i] = -1;
		}
		for (int j = 1;j<=cnt && prime[j] * i <= N-5;j++)
		{
			vis[prime[j] * i] = 1;
			if (i % prime[j] == 0) break;
			else miu[i * prime[j]] = -miu[i];
		}
	}
	for (int i=1;i<=N-5;i++)
		sum[i] = sum[i-1] + (LL)miu[i];
}
LL get_ans(int n, int m)
{
	LL ans = 0;
	if (n > m) swap(n,m);
	for (int l = 1,r;l <=n;l = r+1)
	{
		r = min(n / (n / l), m / (m / l));
		ans += 1LL * (n / l) * (m / l) * (sum[r / k] - sum[(l-1)/k]);//注意long long
	}
	return ans;
}
int main()
{
	Miu();
	scanf("%d", &T);
	while (T--)
	{
		scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);
		ans = get_ans(b,d) - get_ans(a-1,d) - get_ans(b,c-1) + get_ans(a-1,c-1);
		printf("%lld\n", ans);
	}
}