bzoj 2820 YY的GCD 莫比乌斯反演

 题目大意:

给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对

 

这里就抄一下别人的推断过程了

后面这个g(x) 算的方法就是在线性筛的时候只考虑当前的数最小因子,如果进来的最小因子不存在,相当于在之前那个数的基础上的每个mu值都多加了一个质数,那么

这些mu值就要取反,如果已经包含了这个最小因子,我这里另外进行了跟之前类似的讨论方法,在代码中写着

因为这题目数据比较大,这里求解的时候不应该线性求,因为总是有一段区间的n/i*(m/i)值相同,将g[]数组求一个前缀和,记录一段区间得到的值,可以缩小到

sqrt(n)的复杂度

 

 1 /*bzoj2820 YY的GCD*/
 2 #include <bits/stdc++.h>
 3 
 4 using namespace std;
 5 #define ll long long
 6 #define N 10000000
 7 int mu[N+5] , prime[N+5] , tot , f[N+5] , sum[N+5];
 8 bool check[N+5];
 9 
10 void get_mu()
11 {
12     mu[1] = 1;
13     for(int i=2 ; i<=N ; i++){
14         if(!check[i]){
15             prime[tot++] = i;
16             mu[i] = -1;
17             f[i] = 1;
18         }
19         for(int j=0 ; j<tot ; j++){
20             if((ll)prime[j]*i>N) break;
21             check[prime[j]*i] = true;
22             if(i%prime[j]){
23                 mu[i*prime[j]] = -mu[i];
24                 f[i*prime[j]] = -f[i]+mu[i];
25             }else{
26                 /*
27                 本身i中已经含有素数prime[j]
28                 在if中表示含有至少3个因子prime[j],那么最后不管怎么样,当前数值除以一个因子都至少
29                 含有两个prime[j],那么必然为0
30                 else 只有2个prime[j]的情况,那么就是除了除以prime[j]其他情况得到的都是至少有两个因子
31                 的,那么mu[]必然为0,而因为从当前增加了一个因子,那么就是讲那个取到的取反就可以了
32                 */
33                 if((i/prime[j])%prime[j]==0) f[i*prime[j]] = mu[i];
34                 else f[i*prime[j]] = -mu[i/prime[j]];
35                 break;
36             }
37         }
38     }
39     for(int i=1 ; i<=N ; i++) sum[i] = sum[i-1]+f[i];
40 }
41 int n , m;
42 
43 ll solve()
44 {
45     int mn = min(n , m) , last ;
46     ll ret=0;
47     for(int i=1 ; i<=mn ; i=last+1){
48         last = min(n/(n/i) , m/(m/i));
49         ret += (ll)(sum[last]-sum[i-1])*(n/i)*(m/i);
50     }
51     return ret;
52 }
53 
54 int main()
55 {
56     freopen("in.txt" , "r" , stdin);
57     get_mu();
58     int T;
59     scanf("%d" , &T);
60     while(T--){
61         scanf("%d%d" , &n , &m);
62         printf("%lld\n" , solve());
63     }
64     return 0;
65 }

 

转载于:https://www.cnblogs.com/CSU3901130321/p/4787151.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/405097.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux open函数使用方法记录

1、作用 linux一切皆文件&#xff0c;那么如何打开文件&#xff0c;靠的就是open函数。打开和创建文件&#xff0c;对于文件来说存在文件则打开&#xff0c;不存在则创建&#xff1b;对于设备&#xff0c;可以打开相应设备&#xff08;串口&#xff0c;SPI&#xff0c;I2C等&a…

[react] 写个例子说明什么是JSX的内联条件渲染

[react] 写个例子说明什么是JSX的内联条件渲染 ifTrue ? <ComA /> : <ComB>ifTrue && <ComC/> 个人简介 我是歌谣&#xff0c;欢迎和大家一起交流前后端知识。放弃很容易&#xff0c; 但坚持一定很酷。欢迎大家一起讨论 主目录 与歌谣一起通关前…

tar.xz文件如何解压

From: http://blog.csdn.net/rheostat/article/details/7614451 XZ压缩最新压缩率之王 xz这个压缩可能很多都很陌生&#xff0c;不过您可知道xz是绝大数linux默认就带的一个压缩工具。 之前xz使用一直很少&#xff0c;所以几乎没有什么提起。 我是在下载phpmyadmin的时候…

JSON 之 SuperObject(11): TSuperTableString、TSuperAvlEntry

通过 ISuperObject.AsObject 可获取一个 TSuperTableString 对象.TSuperTableString 的常用属性: count、GetNames、GetValues varjo: ISuperObject;jts: TSuperTableString; beginjo : SO({A:1, B:2, C:3, D:{x:4, y:5, z:6}});jts : jo.AsObject;ShowMessage(IntToStr(jts.co…

linux重定向到文件permission denied

1. 执行重定向命令报权限错误&#xff0c;如下&#xff1a; sudo nohup ${obj} >${log_file} 2>&1 &./autoRun.sh: line 23: /home/netmeas/log-20200414-164657.txt: Permission denied2.分析原因 sudo 给了obj超级执行权限&#xff0c;但是重定向到log_file没…

[react] 在React中怎么阻止事件的默认行为?

[react] 在React中怎么阻止事件的默认行为&#xff1f; event.preventDefault();阻止浏览器默认行为&#xff0c; 例如标签不跳转 event.stopPropagation();阻止冒泡&#xff1b; 例如上级点击事件不生效 个人简介 我是歌谣&#xff0c;欢迎和大家一起交流前后端知识。放弃很…

构建第一个Flex的Mobile APP

Flash Builder 4.5已经支持直接创建Flex Mobile Project&#xff0c;写一个最简单的例子 1、建立工程 右击--》新建--》输入工程名“MyFirstMobileApp” 点击“Next”进入下一步 修改初始化的标题文本信息为“Home”&#xff08;默认为HomeView&#xff09;&#xff0c;勾选“G…

由于可能不会将凭据发送到远程计算机,因此将不会进行连接。若要获得协助,请与您的系统管理员联系。...

windows10系统下&#xff0c;在通过VPN连接到堡垒机的时候&#xff0c;出现“由于可能不会将凭据发送到远程计算机&#xff0c;因此将不会进行连接。若要获得协助&#xff0c;请与您的系统管理员联系。” 以下是解决方法&#xff1a; 方法一.组策略 请按照下列步骤操作&#xf…

Linux多线程开发-线程同步-互斥锁pthread_mutex_t

1、互斥锁 同一时刻只允许一个线程对临界区进行访问。POSIX库中用类型pthread_mutex_t来定义互斥锁&#xff0c;类型在pthreadtypes.h中定义。 2、如何声明一个互斥锁 #include <pthread.h>pthread_mutex_t mutex; 3、如何初始化一个互斥锁 //静态初始化互斥锁 pthr…

[react] React的事件和普通的HTML事件有什么不同

[react] React的事件和普通的HTML事件有什么不同 区别&#xff1a; 对于事件名称命名方式&#xff0c;原生事件为全小写&#xff0c;react 事件采用小驼峰对于事件函数处理语法&#xff0c;原生事件为字符串&#xff0c;react 事件为函数react 事件不能采用 return false 的方…

cmb

if (this.cmbMonitorLevelPollutionSource.Items.Count 0) { this.cmbMonitorLevelPollutionSource.Items.Add(new ComboBoxItem() { Content "国控", Tag "00" }); this.cmbMonitorLevelPollutionSource.Ite…

C++面试宝典2011版

From: http://blog.csdn.net/jingxuewang110/article/details/6789557/ 1.new、delete、malloc、free关系 delete会调用对象的析构函数,和new对应free只会释放内存&#xff0c;new调用构造函数。malloc与free是C/C语言的标准库函数&#xff0c;new/delete是C的运算符。它们都可…

hdu1466(dp)

Problem Description平面上有n条直线&#xff0c;且无三线共点&#xff0c;问这些直线能有多少种不同交点数。比如,如果n2,则可能的交点数量为0(平行)或者1(不平行)。Input输入数据包含多个测试实例,每个测试实例占一行,每行包含一个正整数n&#xff08;n<20&#xff09;,n表…

[react] 你最喜欢React的哪一个特性(说一个就好)

[react] 你最喜欢React的哪一个特性&#xff08;说一个就好&#xff09; virtual domHOC 个人简介 我是歌谣&#xff0c;欢迎和大家一起交流前后端知识。放弃很容易&#xff0c; 但坚持一定很酷。欢迎大家一起讨论 主目录 与歌谣一起通关前端面试题

Linux多线程开发-线程同步-读写锁pthread_rwlock_t

1、读写锁概念 对资源的访问抽象为两种类型&#xff0c;即独占和共享。独占资源具有排他性&#xff0c;例如写操作。共享资源可以同时由多个线程访问&#xff0c;不对该资源执行写操作。针对资源的读和写操作分别加锁&#xff0c;写操作加锁与互斥锁相同&#xff0c;但读操作加…

Wireshark抓包工具使用教程以及常用抓包规则

From: http://fangxin.blog.51cto.com/1125131/735178 Wireshark是一个非常好用的抓包工具&#xff0c;当我们遇到一些和网络相关的问题时&#xff0c;可以通过这个工具进行分析&#xff0c;不过要说明的是&#xff0c;这只是一个工具&#xff0c;用法是非常灵活的&#xff0c;…

SSH整合--1

简单说明&#xff1a;整个整合过程使用mysql数据库、Myeclipse 8.5&#xff0c;框架使用struts2.1.6,hibernate3.x,spring2.5 功能&#xff1a;实现简单的用户登录 0. 简单的流程 1. Resister.jsp <% page language"java" import"java.util.*" pageEnco…

[react] 在React中什么时候使用箭头函数更方便呢?

[react] 在React中什么时候使用箭头函数更方便呢&#xff1f; 合成事件当中使用箭头函数比较方便 注意点在于当事件作为props传递给子组件时&#xff0c;会造成额外的性能损耗 个人简介 我是歌谣&#xff0c;欢迎和大家一起交流前后端知识。放弃很容易&#xff0c; 但坚持一…

poj3750

简单题 View Code #include <iostream>#include <cstdio>#include <cstdlib>#include <cstring>usingnamespacestd;#definemaxn 100intn, w, s;charname[maxn][100];boolout[maxn];intmain(){//freopen("t.txt", "r", stdin);sca…

Linux多线程开发-线程同步-条件变量pthread_cond_t

1、条件变量的概念 一个线程A的执行需要另一个线程B来唤醒&#xff0c;否则A挂起等待。线程B可以产生线程A继续执行的信号。条件变量常用在共享数据状态变化的场景中&#xff0c;例如&#xff1a;生产则和消费者问题。POSIX线程库提供了条件变量这种同步机制。使用条件变量需要…