传送门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(莫比乌斯反演)相关推荐

  1. BZOJ 2820 YY的GCD 莫比乌斯反演

    2820: YY的GCD Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y) ...

  2. BZOJ 2820 YY的GCD 莫比乌斯反演

    题意:链接 方法:莫比乌斯反演 解析: 这题跟上一篇博客有一点差别,当然我们能够考虑枚举素数这个大暴力.只是当你A掉这道题后发现正解?都将近5s时.就放弃了这个念头. 相同的式子我们能够直接搬过来.p ...

  3. bzoj 2820 YY的GCD - 莫比乌斯反演 - 线性筛

    Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种 傻×必 ...

  4. bzoj 2820: YY的GCD(莫比乌斯反演)

    2820: YY的GCD Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 2111  Solved: 1137 [Submit][Status][Di ...

  5. BZOJ 2820: YY的GCD

    2820: YY的GCD Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 1705  Solved: 908 [Submit][Status][Dis ...

  6. 【bzoj2820】YY的GCD 莫比乌斯反演

    题目描述 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种 傻×必然不会了,于是 ...

  7. [洛谷P2257] YY的GCD (莫比乌斯反演)

    YY的GCD 题目描述 神犇YY虐完数论后给傻×kAc出了一题 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻× ...

  8. BZOJ 2820 YY的GCD(莫比乌斯反演)

    Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对,kAc这种傻×必 ...

  9. YY的GCD 莫比乌斯反演

    ---题面--- 题解: $ans = \sum_{x = 1}^{n}\sum_{y = 1}^{m}\sum_{i = 1}^{k}[gcd(x, y) == p_{i}]$其中k为质数个数 $$ ...

  10. [BZOJ]2820: YY的GCD

    Time Limit: 10 Sec  Memory Limit: 512 MB Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1< ...

最新文章

  1. 阿里云联合8家芯片商推“全平台通信模组”,加速物联网生态建设...
  2. 如何在Python中反转列表?
  3. Windows10 + Visual Studio 2017 + CMake +OpenCV编译、开发环境配置及测试
  4. mysql 权重搜索
  5. 微博遭鹿晗恋情暴击瘫痪,如何快速提高系统性能?
  6. ALAssetsLibrary获取相册列表
  7. java 擦出,Java擦除和转换实例分析
  8. MAC电脑8款常用设计工具推荐
  9. MySQL索引原理、失效情况
  10. mysql 事物隔离级别解读
  11. 【ES10(2019)】Object.fromEntries()
  12. 诗与远方:无题(九)
  13. linux 查找文件 mysql数据库_Linux下MySQL数据库目录多了好多文件
  14. Windows API ——GetLogicalDriveStrings——获取逻辑驱动器
  15. mysql中drop语法错误_MySQL DROP TABLE操作以及 DROP 大表时的注意事项
  16. Dns数据包内容分析
  17. 菜鸟教程android数据库,Flutter - Android 菜鸟教程打包
  18. linux泰语语言包,linux安装中文语言包(示例代码)
  19. 拼多多的“最初一公里”战事
  20. Android性能优化(一):APP启动优化

热门文章

  1. https://developer.apple.com 苹果开发者
  2. Scarborough Fair(斯卡布罗集市)歌词赏析
  3. 数据挖掘作业学习学习笔记-电商产品评论数据情感分析
  4. Kindle产品评论数据分析---基于excel的分析过程
  5. 跟小米、特斯拉分“蛋糕”的优必选要IPO
  6. 【WPF学习】第三十五章 资源字典
  7. Java设计模式—备忘录模式
  8. 大二下期英语学期总结
  9. Outsmart Your Own Biases 战胜自己的偏见
  10. 科技融合教育!重庆本土英语机构荣获雅思论坛奖项