POJ-2773 欧几里得 + 二分 + 容斥

题意: 求与n互质的第k个数


开始看到这题很蒙 试了好几种做法都T了 后来才知道

由于GCD(a,b) =  GCD(b, a%b) = GCD(a%b,b)

所以 GCD(a,b) = GCD(a+b , b) = GCD( (a+b) %b, b) = GCD( a+b*t,b)


本题可以把小于n于n互质的数都求出来  后面的于n互质的是都是前面第一周期(也就是小于n于n互质的数)+ n的倍数了

所以先求出k是第一周期的多少倍 用倍数× n再加上 第k对n取余个(第一周期的互质数)

注意结果有可能爆int

注意取余结果有可能是0 需要返回p-1倍+第一周期最有一个数


code:

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
typedef long long ll;
const int lim = 1000005;
int z[lim];
ll gcd(ll a,ll b){return a%b==0?b:gcd(b,a%b);
}
int main()
{int n,k;while(~scanf("%d%d",&n,&k)){memset(z,0,sizeof(z));z[1]=1;int sze=1;for(int i=2;i<n;i++)if(gcd(i,n)==1)z[++sze]=i;ll p = k/sze;//注意int y = k%sze;if(y==0)p--,y=sze;//注意printf("%lld\n",p*n+z[y]);}return 0;
}


这道题还有一个做法 就是 二分那个第k个与n互质的数

利用容斥原理 去做二分查找的判定 那么这个与n互质的个数一定是随上界的上升而增多的

我们计算与n互质的个数 可以利用n-与n不互质的个数 = 与n互质的个数

找出n的所有因子 用 n - ( n/p1 + n/p2 +...n/pn - n/(p1*p2) -n/(p1*p3) - n/(p1*p4)...+n/(p1*p2*p3)+.....)

即可算出所有与n互质的数的个数

 

#include<iostream>
#include<cstdio>
#include<cctype>
#include<limits.h>
#include<vector>
using namespace std;
typedef unsigned long long ll;
const ll inf = LONG_MAX;
ll n,k;
ll sum;
vector<int>p;void devide(ll a)// 不必真的去那质数表分解 枚举分解法 因子也是质数 
{for(int i=2;i*i<=n;i++){if(a%i==0){p.push_back(i);do a/=i;while(a%i==0);}}if(a>1)p.push_back(a);
}
void dfs(ll pos,ll now,ll cnt,ll x){//容斥计算if(pos==p.size()){if(now==1)return;if(cnt&1)sum+=x/now;else sum-=x/now;return;}dfs(pos+1,now*p[pos],cnt+1,x);dfs(pos+1,now,cnt,x);
}
bool test(ll x)
{sum=0;dfs(0,1,0,x);return (x-sum)>=k; 
}
int main()
{while(~scanf("%d%d",&n,&k)){p.clear();devide(n);ll mid,l = 1,r = inf;if(k==1){printf("1\n");continue;}ll ans;while(l<r){		mid = (l+r)>>1; if(test(mid))r=mid;else l = mid+1; }printf("%lld\n",r);}	return 0;
} 

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

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

相关文章

html:(35):缩进和行高

段落排版--缩进 中文文字中的段前习惯空两个文字的空白&#xff0c;这个特殊的样式可以用下面代码来实现&#xff1a; p{text-indent:2em;} <p>1922年的春天&#xff0c;一个想要成名名叫尼克卡拉威&#xff08;托比?马奎尔Tobey Maguire 饰&#xff09;的作家&#x…

new/delete与malloc/free

C语言中使用malloc/calloc/realloc用来在堆上分配空间&#xff0c;free将申请的空间释放掉。 malloc&#xff1a; 原型&#xff1a;extern void *malloc(unsigned int num_bytes)。 功能&#xff1a;分配长度为num_bytes字节的内存块。 1 int *p(int*)malloc(sizeof(int));…

Java学习笔记5-1——多线程

目录前言核心概念线程创建继承Thread类实现Runnable接口上述两个方法小结实现Callable接口并发问题简介静态代理模式线程状态线程停止&#xff08;stop&#xff09;线程休眠&#xff08;sleep&#xff09;线程礼让&#xff08;yield&#xff09;线程强制执行&#xff08;join&a…

[Leedcode][JAVA][第128题][最长连续序列][Hash]

【问题描述】[困难] 给定一个未排序的整数数组&#xff0c;找出最长连续序列的长度。要求算法的时间复杂度为 O(n)。示例:输入: [100, 4, 200, 1, 3, 2] 输出: 4 解释: 最长连续序列是 [1, 2, 3, 4]。它的长度为 4。【解答思路】 1. SortCompare 先排序&#xff0c;题意要求连…

【应用】Lambda表达式

目录简介用法前提和语法作用例子推导lambda表达式&#xff08;内部类的简化过程&#xff09;带参数的情况以及lambda表达式的简化&#xff1a;简介 Lambda 表达式&#xff0c;也可称为闭包&#xff0c;它是推动 Java 8 发布的最重要新特性。Lambda 允许把函数作为一个方法的参…

html:(36):间距和对齐

段落排版--中文字间距、字母间距 中文字间隔、字母间隔设置&#xff1a; 如果想在网页排版中设置文字间隔或者字母间隔就可以使用 letter-spacing 来实现&#xff0c;如下面代码&#xff1a; h1{letter-spacing:50px; } ... <h1>了不起的盖茨比</h1>注意&…

二维凸包 Graham's Scan

凸包是啥&#xff1f;&#xff1f; 凸包就是&#xff1a; 给定二维平面上的点集&#xff0c;凸包就是将最外层的点连接起来构成的凸多边型&#xff0c;它能包含点集中所有的点。 外层的红线 就是凸包 一般题目会让你求利用最外层凸包的一些性质如下面这道题 来看一道模板题…

Spring boot项目集成Sharding Jdbc

环境    jdk&#xff1a;1.8    framework&#xff1a; spring boot&#xff0c; sharding jdbc    database&#xff1a; MySQL    搭建步骤    在pom 中加入sharding 依赖    <dependency>    <groupId>io.shardingjdbc</groupId…

[Leedcode][JAVA][第300题][最长上上子序列][动态规划][压缩空间]

【问题描述】[中等] 给定一个无序的整数数组&#xff0c;找到其中最长上升子序列的长度。示例:输入: [10,9,2,5,3,7,101,18] 输出: 4 解释: 最长的上升子序列是 [2,3,7,101]&#xff0c;它的长度是 4。 说明:可能会有多种最长上升子序列的组合&#xff0c;你只需要输出对应的…

Java学习笔记5-2——多线程

目录线程同步三大不安全案例一、不安全的买票过程二、不安全的取钱过程三、线程不安全的集合synchronized解决三大不安全案例一、解决不安全的买票过程二、解决不安全的取钱过程三、解决集合线程不安全死锁可重入锁&#xff08;ReentrantLock&#xff09;线程协作&#xff08;线…

html:(38):元素分类和块级元素

元素分类 在讲解CSS布局之前&#xff0c;我们需要提前知道一些知识&#xff0c;在CSS中&#xff0c;html中的标签元素大体被分为三种不同的类型&#xff1a;块状元素、内联元素(又叫行内元素)和内联块状元素。 常用的块状元素有&#xff1a; <div>、<p>、<h1…

Spyder kernel died 错误

Keras 2.2.4版本和 tensorflow1.2.1 版本不兼容导致的错误。降低Keras 为2.1.2版本 import keras 出现&#xff1a; Using TensorFlow backend. Kernel died, restarting conda uninstall kerasconda install keras2.1.2 转载于:https://www.cnblogs.com/maxiaodoubao/p/10421…

html:(39):块级元素和内联块级元素

元素分类--内联元素 在html中&#xff0c;<span>、<a>、<label>、 <strong> 和<em>就是典型的内联元素&#xff08;行内元素&#xff09;&#xff08;inline&#xff09;元素。当然块状元素也可以通过代码display:inline将元素设置为内联元素。…

HDU-6103

。 Kirinriki Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1204 Accepted Submission(s): 475 Problem Description We define the distance of two strings A and B with same length n is disA,B∑i0n−1|A…

Java学习笔记6——网络编程

目录基础知识查询IP地址端口通信协议TCP实现聊天TCP实现文件上传UDPUDP实现聊天UDP多线程在线咨询URL下载网络资源基础知识 请自行翻阅《计算机网络》 查询IP地址 import java.net.InetAddress; import java.net.UnknownHostException;public class TestInetAddress {public…

[剑指offer][JAVA]面试题第[15]题[二进制中1的个数][位运算]

【问题描述】[简单] 请实现一个函数&#xff0c;输入一个整数&#xff08;无符号数&#xff09;&#xff0c;输出该数二进制表示中 1 的个数。例如&#xff0c;把 9 表示成二进制是 1001&#xff0c;有 2 位是 1。因此&#xff0c;如果输入 9&#xff0c;则该函数输出 2。示例…

课外阅读(XHTML和XML简介)

课外阅读 一、XHTML简介 HTML从出现发展到今天&#xff0c;仍有些缺陷和不足。HTML的三个主要缺点如下。 &#xff08;1&#xff09;太简单。不能适应现在越来越多的网络设备和应用的需要&#xff0c;比如手机、PDA、信息家电都不能直接显示HTML&#xff1b; &#xff08;2…

Java学习笔记7-1——注解与反射

目录注解注解是啥内置注解元注解TargetRetentionDocumentedInherited自定义注解反射静态语言 和 动态语言Java 反射机制概述Java 反射机制提供的功能Java 反射优点和缺点反射相关的主要 API注解 注解是啥 注解&#xff08;Annotation&#xff09;是从 JDK 1.5 开始引入的新技…

C++11并发之std::thread

C11并发之std::thread知识链接&#xff1a;C11 并发之std::mutexC11 并发之std::atomic本文概要&#xff1a;1、成员类型和成员函数。2、std::thread 构造函数。3、异步。4、多线程传递参数。5、join、detach。6、获取CPU核心个数。7、CPP原子变量与线程安全。8、lambda与多线程…

极角排序的几种方法

转载自http://www.cnblogs.com/aiguona/p/7248311.html&#xff01;&#xff01; 关于极角排序&#xff1a; 在平面内取一个定点O&#xff0c;叫极点&#xff0c;引一条射线Ox&#xff0c;叫做极轴&#xff0c;再选定一个长度单位和角度的正方向&#xff08;通常取逆时针方向&a…