【BZOJ】2820: YY的GCD

【题意】给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对。T<=10^4，N,M<=10^7。

【算法】数论（莫比乌斯反演）

【题解】公式推导见DQSSS。

推到ans= Σ_p是素数Σ_d≤minsμ(d) * (n/pd) * (m/pd)，mins=min(n/p,m/p)。

使用枚举取值的方法再枚举素数单次询问复杂度√n*(n/ln n)，显然不能满足要求。

问题在于枚举素数，令T=pd，则：

ans= Σ_T≤mins (n/T) * (m/T)*Σ_{p|T&&p是素数} μ(T/p)，mins=min(n,m)。

后面部分可以枚举素数的倍数预处理出μ前缀和，复杂度O((n/ln n)*ln n)即O(n)。

每次询问再枚举取值O(√n)解决。

总复杂度O(T*√n+n)。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=10000010,N=10000000;
int miu[maxn],mius[maxn],prime[maxn],tot;
ll s[maxn];
bool mark[maxn];
void pre(){miu[1]=1;for(int i=2;i<=N;i++){if(!mark[i])miu[prime[++tot]=i]=-1;for(int j=1;j<=tot&&i*prime[j]<=N;j++){mark[i*prime[j]]=1;if(i%prime[j]==0)break;miu[i*prime[j]]=-miu[i];}}for(int i=1;i<=tot;i++){for(int j=prime[i];j<=N;j+=prime[i])mius[j]+=miu[j/prime[i]];}for(int i=1;i<=N;i++)s[i]=s[i-1]+mius[i];
}
int main(){pre();int T,n,m;scanf("%d",&T);while(T--){scanf("%d%d",&n,&m);int pos=0,mins=min(n,m);ll ans=0;for(int i=1;i<=mins;i=pos+1){pos=min(n/(n/i),m/(m/i));ans+=(s[pos]-s[i-1])*(n/i)*(m/i);}printf("%lld\n",ans);}return 0;
}

View Code

将枚举两个数改为枚举乘积和其中一个，即T和p|T，后面的p|T可以O(n log n)贡献处理前缀和。

转载于:https://www.cnblogs.com/onioncyc/p/8274302.html

【BZOJ】2820: YY的GCD相关推荐

BZOJ 2820 YY的GCD 莫比乌斯反演
2820: YY的GCD Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y) ...
bzoj 2820: YY的GCD（莫比乌斯反演）
2820: YY的GCD Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2111 Solved: 1137 [Submit][Status][Di ...
BZOJ 2820: YY的GCD
2820: YY的GCD Time Limit: 10 Sec Memory Limit: 512 MB Submit: 1705 Solved: 908 [Submit][Status][Dis ...
bzoj 2820 YY的GCD - 莫比乌斯反演 - 线性筛
Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种傻×必 ...
[BZOJ]2820: YY的GCD
Time Limit: 10 Sec Memory Limit: 512 MB Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1< ...
BZOJ 2820 YY的GCD（莫比乌斯反演）
Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对,kAc这种傻×必 ...
BZOJ 2820 YY的GCD 莫比乌斯反演
题意:链接方法:莫比乌斯反演解析: 这题跟上一篇博客有一点差别,当然我们能够考虑枚举素数这个大暴力.只是当你A掉这道题后发现正解?都将近5s时.就放弃了这个念头. 相同的式子我们能够直接搬过来.p ...
[BZOJ 2820] YY的GCD
题意求下式的值: \[ \sum_{i=1}^n\sum_{j=1}^m \mathbb{P}(\gcd(i,j)) \] 其中 \(\mathbb{P}(x)\) 当 \(x\) 为质数时为 \( ...
【BZOJ 2820】 YY的GCD
2820: YY的GCD Time Limit: 10 Sec Memory Limit: 512 MB Submit: 807 Solved: 404 [ Submit][ Status] ...

【BZOJ】2820: YY的GCD

【BZOJ】2820: YY的GCD相关推荐

最新文章

热门文章