BZOJ2301:[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)的个数

Sample Input

2 5 1 5 1

1 5 1 5 2

Sample Output

HINT

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

Solution

和BZOJ1101一样……只不过简单容斥一下就好了……假设下界为1，答案为$ans_{b,d}-ans_{a-1,d}-ans_{b,c-1}+ans_{a-1,c-1}$

Code

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4 #define N (100000+1000)
 5 using namespace std;
 6
 7 int T,a,b,c,d,k,vis[N],prime[N],sum[N],mu[N],cnt;
 8
 9 void Get_mu()
10 {
11     mu[1]=1;
12     for (int i=2; i<=50000; ++i)
13     {
14         if (!vis[i]){prime[++cnt]=i,mu[i]=-1;}
15         for (int j=1; j<=cnt && prime[j]*i<=50000; ++j)
16         {
17             vis[prime[j]*i]=true;
18             if (i%prime[j]==0) break;
19             mu[prime[j]*i]=-mu[i];
20         }
21     }
22     for (int i=1; i<=50000; ++i) sum[i]=sum[i-1]+mu[i];
23 }
24
25 int Calc(int n,int m)
26 {
27     int ans=0; if (n>m) swap(n,m);
28     for (int l=1,r; l<=n; l=r+1)
29     {
30         r=min(n/(n/l),m/(m/l));
31         ans+=(sum[r]-sum[l-1])*(n/l)*(m/l);
32     }
33     return ans;
34 }
35
36 int main()
37 {
38     scanf("%d",&T);
39     Get_mu();
40     while (T--)
41     {
42         scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
43         printf("%d\n",Calc(b/k,d/k)-Calc((a-1)/k,d/k)-Calc(b/k,(c-1)/k)+Calc((a-1)/k,(c-1)/k));
44     }
45 }

转载于:https://www.cnblogs.com/refun/p/9554132.html

BZOJ2301:[HAOI2011]Problem b(莫比乌斯反演,容斥)相关推荐

BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演容斥原理]【学习笔记】
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MB Submit: 4032 Solved: 1817 [Submi ...
BZOJ2301: [HAOI2011]Problem b(莫比乌斯反演)
Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数 ...
莫比乌斯反演/容斥 +2020ICPC 江西省大学生程序设计竞赛 A Simple Math Problem
题目描述输入描述: 输出描述: 示例1 输入 3 输出 5 分析: 1.这个题其实考的是一个莫比乌斯反演题,但是由于我知识储备不够,没有看出来,题目给的范围可以瞎搞一下,所以下面容斥可以过. 2.转 ...
[POI2007]ZAP-Queries [HAOI2011]Problem b 莫比乌斯反演
1,[POI2007]ZAP-Queries ---题面--- 题解: 首先列出式子:$$ans = \sum_{i = 1}^{n}\sum_{j = 1}^{m}[gcd(i, j) == d]$ ...
BZOJ 2301 - Problem b（莫比乌斯反演+容斥）
题目链接 https://cn.vjudge.net/problem/HYSBZ-2301 [题意] 对于给出的 n 个询问,每次求有多少个数对(x,y)(x,y) ,满足 a≤x≤b,c≤y≤da≤ ...
P2522 HAOI2011 Problem b [莫比乌斯反演,数论分块]
P2522 HAOI2011 题意对于给出的n个询问,每次求有多少个数对(x,y)(x,y)(x,y),满足a≤x≤ba≤x≤ba≤x≤b,c≤y≤dc≤y≤dc≤y≤d,且gcd(x,y)=kgc ...
Comet OJ - Contest #8 E神奇函数（莫比乌斯函数容斥）
Comet OJ - Contest #8 E神奇函数(莫比乌斯函数容斥) 题目大意定义d(x)d(x)d(x)的值为x的最小素因子,定义 f(x)={1x=1d(x)f(xd2(x))x>1 ...
BZOJ2301: [HAOI2011]Problem b
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MB Submit: 6435 Solved: 2986 [Submi ...
[bzoj2301][HAOI2011]Problem b
题目大意:给你$a,b,c,d,k(1\leqslant a,b,c,d,k\leqslant 5\times10^4)$,求$\displaystyle\sum\limits_{x=a}^b\dis ...