【bzoj 2820】YY的GCD(莫比乌斯反演)
传送门biu~
参考博客:莫比乌斯反演
#include<bits/stdc++.h>
#define N 10000000
using namespace std;
int T,n,m,prime[N+5],mu[N+5],f[N+5];
bool mark[N+5];
void GetPrime(){mu[1]=1;for(int i=2;i<=N;i++){if(!mark[i]) prime[++prime[0]]=i,mu[i]=-1,f[i]=1;for(int j=1;j<=prime[0] && i*prime[j]<=N;j++){mark[i*prime[j]]=1;if(i%prime[j]==0){mu[i*prime[j]]=0;f[i*prime[j]]=mu[i];break;}mu[i*prime[j]]=-mu[i];f[i*prime[j]]=mu[i]-f[i];}}for(int i=1;i<=N;i++) f[i]+=f[i-1];
}
long long calc(int n,int m){if(n>m) swap(n,m);long long res=0;for(int r,i=1;i<=n;i=r+1){r=min(n/(n/i),m/(m/i));res+=(long long)(f[r]-f[i-1])*(n/i)*(m/i);}return res;
}
int main() {GetPrime();scanf("%d",&T);while(T--){scanf("%d%d",&n,&m);printf("%lld\n",calc(n,m));}return 0;
}
【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 莫比乌斯反演
题意:链接 方法:莫比乌斯反演 解析: 这题跟上一篇博客有一点差别,当然我们能够考虑枚举素数这个大暴力.只是当你A掉这道题后发现正解?都将近5s时.就放弃了这个念头. 相同的式子我们能够直接搬过来.p ...
- 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(莫比乌斯反演)
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 ...
- 【bzoj2820】YY的GCD 莫比乌斯反演
题目描述 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种 傻×必然不会了,于是 ...
- [洛谷P2257] YY的GCD (莫比乌斯反演)
YY的GCD 题目描述 神犇YY虐完数论后给傻×kAc出了一题 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻× ...
- BZOJ 2820 YY的GCD(莫比乌斯反演)
Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对,kAc这种傻×必 ...
- YY的GCD 莫比乌斯反演
---题面--- 题解: $ans = \sum_{x = 1}^{n}\sum_{y = 1}^{m}\sum_{i = 1}^{k}[gcd(x, y) == p_{i}]$其中k为质数个数 $$ ...
- [BZOJ]2820: YY的GCD
Time Limit: 10 Sec Memory Limit: 512 MB Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1< ...
最新文章
- 阿里云联合8家芯片商推“全平台通信模组”,加速物联网生态建设...
- 如何在Python中反转列表?
- Windows10 + Visual Studio 2017 + CMake +OpenCV编译、开发环境配置及测试
- mysql 权重搜索
- 微博遭鹿晗恋情暴击瘫痪,如何快速提高系统性能?
- ALAssetsLibrary获取相册列表
- java 擦出,Java擦除和转换实例分析
- MAC电脑8款常用设计工具推荐
- MySQL索引原理、失效情况
- mysql 事物隔离级别解读
- 【ES10(2019)】Object.fromEntries()
- 诗与远方:无题(九)
- linux 查找文件 mysql数据库_Linux下MySQL数据库目录多了好多文件
- Windows API ——GetLogicalDriveStrings——获取逻辑驱动器
- mysql中drop语法错误_MySQL DROP TABLE操作以及 DROP 大表时的注意事项
- Dns数据包内容分析
- 菜鸟教程android数据库,Flutter - Android 菜鸟教程打包
- linux泰语语言包,linux安装中文语言包(示例代码)
- 拼多多的“最初一公里”战事
- Android性能优化(一):APP启动优化