GCD HDU - 1695

题意：

给出a,b,c,d,k,求出a<=x<=b, c<=y<=d 且gcd(x,y) == k 的（x,y）的对数。
求的是不同数量对的总数

题解：

和这个题一样P3455 [POI2007]ZAP-Queries，但是本题要求求不同数量对的总数，所以最后的结果要减去重复值
如果是莫比乌斯+分块做法，对于区间[1,b],[1,d],b<d,重复部分是[1,b]部分，所以减去solve(b,b)/2
如果是容斥做法，对于每个i∈[1,b],求[1,d]中互质的数量，那我们减去i与[1,i]中互质的数量
详细看代码

代码：

莫比乌斯+分块

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{#ifdef ONLINE_JUDGE
#elsestartTime = clock ();freopen("data.in", "r", stdin);
#endif
}
void Time_test()
{#ifdef ONLINE_JUDGE
#elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn=2e6+9;
int prime[maxn];
int mu[maxn];
int sum[maxn];
int vis[maxn];
int cnt=0;
void get_mu(int N){mu[1]=1;vis[1]=vis[0]=1;for(int i=2;i<=N;i++){if(!vis[i]){prime[++cnt]=i;mu[i]=-1;}for(int j=1;j<=cnt&&i*prime[j]<=N;j++){vis[i*prime[j]]=1;if(i%prime[j]==0)break;mu[i*prime[j]]=-mu[i];}}for(int i=1;i<=N;i++){sum[i]=sum[i-1]+mu[i];}
}
ll solve(int a,int b,int k){a/=k;b/=k;int minn=min(a,b);ll ans=0;for(int l=1,r;l<=minn;l=r+1){r=min(a/(a/l),b/(b/l));ans+=1ll*(sum[r]-sum[l-1])*(a/l)*(b/l);}return ans;
}
int main()
{//rd_test();get_mu(1000000);int t;read(t);int cas=0;while(t--){int a,b,c,d,k;read(a,b,c,d,k);if(k==0){printf("Case %d: 0\n",++cas);continue;}if(b>d)swap(b,d);
//      printf("%d %d %d %d\n",a,b,c,d); printf("Case %d: %lld\n",++cas,solve(b,d,k)-solve(b,b,k)/2);}//Time_test();
}

容斥

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{#ifdef ONLINE_JUDGE
#elsestartTime = clock ();freopen("data.in", "r", stdin);
#endif
}
void Time_test()
{#ifdef ONLINE_JUDGE
#elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn=3e5+9;
int prime[maxn];
int cnt=0;
void divide(int n){cnt=0;for(int i=2;i*i<=n;i++){if(n%i==0){prime[cnt++]=i;while(n%i==0)n/=i;}}if(n!=1)prime[cnt++]=n;
}
int solve(int S){int ans=0;for(int i=1;i<(1<<cnt);i++){int tmp=1;int num=0;for(int j=0;j<cnt;j++){if(i&(1<<j)){tmp*=prime[j];num++;}}if(num&1)ans+=S/tmp;else ans-=S/tmp;}return S-ans;
}
int main()
{//rd_test();int t;read(t);int cas=0;while(t--){int a,b,c,d,k;read(a,b,c,d,k);if(k==0){printf("Case %d: 0\n",++cas);continue;}if(b>d)swap(b,d);b/=k;d/=k;ll ans=0;for(int i=1;i<=b;i++){divide(i);ans+=solve(d)-solve(i-1);}printf("Case %d: %lld\n",++cas,ans);}return 0;//Time_test();
}

GCD HDU - 1695相关推荐

Gcd HDU - 6545 （基础数论）
wls 有一个整数 n,他想将 1 − n 这 n 个数字分成两组,每一组至少有一个数,并且使得两组数字的和的最大公约数最大,请输出最大的最大公约数. Input 输入一行一个整数 n. 2 ≤ n ...
HDU 1695（数论，筛选+素因子分解+容斥）
题意:已知给定k,x,y求 1<=a<=x 1<=b<=y 中满足 gcd(a,b)=k 的(a,b)对数.(注意数对是无序的). 1<=x,y<=1e5, 0&l ...
HDU 1695 BZOJ 2301 莫比乌斯反演
看了两个多小时,大概意思我觉得莫比乌斯是一个利用一个数的能被整除的一些数,通过容斥来计算出这个数的一些信息,很高级.这个题就当做是模板吧.然后把两种形式贴出来: 第一种: 第二种: 然后是容斥的函数: ...
hduoj 饭卡 java_饭卡 HDU - 2546
饭卡 HDU - 2546 原题链接考察:01背包这道题懵逼半天,看了大佬的代码终于懂了点思路: 这道题是和购买顺序有关的,我们先拿出最贵的菜,同时余额拿出5元,在剩余的余额尽可能买总和贵的菜. ...
UVA 11426 GCD-Extreme(II) ★ (欧拉函数)
题意求Σ{1<=i<N} Σ{i<j<=N} GCD(i, j) (N<=4000000) 分析原始思路暴力求明显是不行的,我们把式子简化形式一下发现它可以 ...
杭电OJ分类题目（1）
原题出处:HDOJ Problem Index by Type,http://acm.hdu.edu.cn/typeclass.php 杭电OJ分类题目(1) HDU Introduction HDU ...
莫比乌斯反演入门题目(详细)
目录 luoguP2568 GCD hdu1695 GCD luogu3455[POI2007]ZAP-Queries luogu2522 [HAOI2011]Problem b P4318 完全平方 ...
Mobius反演总结
刷了一些浅显的反演题目, 做一个总结知识点1 给定n个数(可以是1-n), 求有多少数与x互质 f(x)=∑d|xu[d]∗cnt[d] f(x) = \sum_{d|x}u[d]*cnt[d] 其 ...
求互质数对的两种方法:欧拉函数和容斥原理
1.HDU4135 首先,我们可以先求出1 ~ a - 1内与n互质的数,再求出1 ~ b内与n互质的数.最后利用前缀和算出来即可. #include <iostream> #includ ...

GCD HDU - 1695

题意：

题解：

代码：

莫比乌斯+分块

容斥

GCD HDU - 1695相关推荐

最新文章

热门文章