P2257 YY的GCD

题意：

求 1≤x≤N,1≤y≤M1 \leq x \leq N,1 \leq y \leq M1≤x≤N,1≤y≤M 且gcd(x, y) 为质数的 (x,y) 有多少对。

题解：

莫比乌斯反演

代码：

#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 LOCALstartTime= clock();freopen("in.txt", "r", stdin);
#endif
}
void Time_test()
{#ifdef LOCALendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn= 1e7 + 9;
ll mu[maxn];
int vis[maxn], prim[maxn], cnt, f[maxn];
int sum[maxn];
void get_mu(int n)
{mu[1]= 1;for (int i= 2; i <= n; i++) {if (!vis[i]) {prim[++cnt]= i;mu[i]= -1;}for (int j= 1; j <= cnt && prim[j] * i <= n; j++) {vis[prim[j] * i]= 1;if (i % prim[j] == 0)break;elsemu[i * prim[j]]= -mu[i];}}for (int i= 1; i <= cnt; i++) {for (int j= 1; prim[i] * j <= n; j++) {f[j * prim[i]]+= mu[j];}}for (int i= 1; i <= n; i++)sum[i]= sum[i - 1] + f[i];
}
ll solve(int a, int b)
{ll ans= 0;for (int l= 1, r= 0; l <= a; l= r + 1) {r= min(a / (a / l), b / (b / l));ans+= 1ll * (sum[r] - sum[l - 1]) * 1ll * (a / l) * (b / l);}return ans;
}
int main()
{//rd_test();get_mu(10000000);int t;read(t);while (t--) {int n, m;read(n, m);if (n > m)swap(n, m);printf("%lld\n", solve(n, m));}//Time_test();
}

P2257 YY的GCD相关推荐

洛谷P2257 YY的GCD 莫比乌斯函数反演+线性筛
洛谷P2257 YY的GCD 标签莫比乌斯反演线性筛前言这题貌似和莫反没多大关系,就是用到了一个莫比乌斯函数的性质了,其他就是推公式,优化和式. 我的第一道懵逼反演-真的好难好难-而且套路特别 ...
[Luogu P2257] YY的GCD (莫比乌斯函数)
题面传送门:洛咕 Solution 推到自闭,我好菜啊显然,这题让我们求: \(\large \sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)\in prime]\) 根 ...
洛谷 P2257 YY的GCD
YY的GCD 题目描述神犇YY虐完数论后给傻×kAc出了一题给定\(N\), \(M\) ,求\(1 \le x \le N,1 \le y \le M\)且\(gcd(x, y)\)为质数的\( ...
[洛谷P2257] YY的GCD (莫比乌斯反演)
YY的GCD 题目描述神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻× ...
P2257 YY的GCD （莫比乌斯反演）
[题目链接] https://www.luogu.org/problemnew/show/P2257 // luogu-judger-enable-o2 /* -------------------- ...
洛谷P2257 YY的GCD
题目描述神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻×必然不会了,于 ...
并不对劲的bzoj2820:p2257:YY的GCD
题目大意 \(t\)(\(t\leq10^4\))组数据,给定\(n,m\)(\(n,m\leq10^6\))求 \[\sum_{x=1}^{n}\sum_{y=1}^{m}[gcd(x,y)=1]\ ...
acm-(数论、莫比乌斯反演)LuoGu P2257 YY的GCD
传送门本题就是让求这样一个式子: ∑ i = 1 n ⁣ ∑ j = 1 m [ g c d ( i , j ) i s p r i m e ] \sum_{i=1}^n\!\sum_{j=1}^m ...
luogu P2257 YY的GCD
嘟嘟嘟感觉这几道数论题都差不多,但这到明显是前几道的升级版. 推了一大顿只能得60分,不得不看题解. 各位看这老哥的题解吧我就是推到他用\(T\)换掉\(kd\)之前,然后枚举\(T\)的.这个转 ...

P2257 YY的GCD

题意：

题解：

代码：

P2257 YY的GCD相关推荐

最新文章

热门文章