BZOJ 2440: [中山市选2011]完全平方数
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 5160 Solved: 2515
[Submit][Status][Discuss]
Description
小 X 自幼就很喜欢数。但奇怪的是,他十分讨厌完全平方数。他觉得这些
数看起来很令人难受。由此,他也讨厌所有是完全平方数的正整数倍的数。然而
这丝毫不影响他对其他数的热爱。
这天是小X的生日,小 W 想送一个数给他作为生日礼物。当然他不能送一
个小X讨厌的数。他列出了所有小X不讨厌的数,然后选取了第 K个数送给了
小X。小X很开心地收下了。
然而现在小 W 却记不起送给小X的是哪个数了。你能帮他一下吗?
Input
包含多组测试数据。文件第一行有一个整数 T,表示测试
数据的组数。
第2 至第T+1 行每行有一个整数Ki,描述一组数据,含义如题目中所描述。
Output
含T 行,分别对每组数据作出回答。第 i 行输出相应的
第Ki 个不是完全平方数的正整数倍的数。
Sample Input
4
1
13
100
1234567
Sample Output
1
19
163
2030745
HINT
对于 100%的数据有 1 ≤ Ki ≤ 10^9
, T ≤ 50
题解
一看到完全平方数的倍数,自然而然就想到莫比乌斯函数,又因为数据范围较大
就想到二分来解决,所以先筛出莫比乌斯函数,之后根据容斥原理算出不是完全
平方数的倍数的个数。有些地方还是理解的不够透彻。
代码
#include<bits/stdc++.h>
#define LL long longusing namespace std;
const int MAXN = 1e5+5;int n,miu[MAXN],cnt,prime[MAXN],T;
bool vis[MAXN];inline LL check(LL x){int xx=(int)sqrt(x);LL ret=0;for(register int i=1;i<=xx;i++)ret+=miu[i]*(x/(i*i));return ret;
}signed main(){miu[1]=1;for(register int i=2;i<=MAXN;i++){if(!vis[i]) {vis[i]=1;prime[++cnt]=i;miu[i]=-1;}for(register int j=1;j<=cnt && (LL)prime[j]*i<=MAXN;j++){vis[(LL)i*prime[j]]=1;if(i%prime[j]==0){miu[(LL)i*prime[j]]=0;break;}else miu[(LL)i*prime[j]]=-miu[i];}}scanf("%d",&T);while(T--){scanf("%d",&n);LL l=0,r=n*2;while(l<r){LL mid=(l+r)>>1;if(check(mid)<n) l=mid+1;else r=mid;}printf("%lld\n",r);}return 0;
}
转载于:https://www.cnblogs.com/sdfzsyq/p/9677025.html
BZOJ 2440: [中山市选2011]完全平方数相关推荐
- BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 3028 Solved: 1460 [Submit][S ...
- bzoj 2440: [中山市选2011]完全平方数(二分+莫比乌斯函数)
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 4018 Solved: 1952 [Submit][S ...
- bzoj 2440: [中山市选2011]完全平方数【莫比乌斯函数+二分】
二分答案,然后用莫比乌斯函数作为容斥系数,计算当前枚举的mid内有几个满足要求的数 #include<iostream> #include<cstdio> #include&l ...
- bzoj2440: [中山市选2011]完全平方数
自己写的第一个博客......... BZOJ 2440 [题意] 求第K个约数不含平方数的数 (1<=k<=10^9), 共有T组数据(T<=50). [题解] 首先题解并不是我独 ...
- BZOJ2440 [中山市选2011]完全平方数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而这丝毫不影响他对其他数的热爱. 这天 ...
- 2017.3.11[bzoj2440][中山市选2011]完全平方数
线性筛+莫比乌斯函数裸题 http://www.lydsy.com/JudgeOnline/problem.php?id=2440 考虑二分答案,对于每个数看它是第几个,即对于这个数看比他小的有多少个 ...
- [BZOJ2440][中山市选2011]完全平方数(莫比乌斯函数,二分)
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2440 题意:RT. 是https://wenku.baidu.com/view/fbec9 ...
- bzoj2440 [中山市选2011]完全平方数 容斥+莫比乌斯函数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些 数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而 这丝毫不影响他对其他数的热爱. 这 ...
- [BZOJ 2438] [中山市选2011]杀人游戏 Tarjan缩点
这个题很容易想到正解就是缩点找入度为零的点,那么我们考虑一种特殊情况就是,一个入度为零的点我们不访问他就知道他是不是凶手,那么这样的话就是:I. 他是一个真·孤立的点 II. 他在图里但是在他的强联通 ...
最新文章
- NSA释出逆向工程工具GHIDRA针对恶意软件逆向分析
- asp.net 6中的mini api和mvc api性能对比
- iphone NSTimer
- property属性学习
- 【原创工具 | NetSM】开源跨平台命令行网速监测(纯 Python 开发)
- Java UDP通信详解:单播、广播、组播
- Mac M1 安装JDK8
- PTC:能源互联网“双子星”,引领风电企业研发数字化变革
- php麻将,php实现麻将牌型的胡牌检测算法
- 致远oa系统unix 服务器,致远oa手机客户端服务器
- HCIA-GaussDB 华为认证数据库工程师
- 蔚来、小鹏、理想自动驾驶能力的纵向演进与横向比较
- NOIP2018初赛
- 机器学习系列(二)——回归模型
- 单目标优化、多目标优化
- ImageJ macro 批量vsi转tif
- 检测恶意软件分类模型中的概念漂移
- [开源名人访谈录] Philippe Gerum
- Image Processing for Embedded Devices 4
- era翻译_era是什么意思_era的翻译_音标_读音_用法_例句_爱词霸在线词典
热门文章
- if shell 条件语句简写_Shell if条件判断语句用法详解
- 【EOJ Monthly 2019.02 - E】中位数(二分 ,中位数 ,−1/1变换,dp求解DAG最长路)
- 笔试题 - 阶乘问题
- 【牛客 - 368D】动态连通块(并查集+bitset优化)
- Apollo进阶课程⑳丨Apollo感知之旅——机器学习与感知的未来
- mysql判断域为空_MySQL EXPLAIN 字段说明
- 买iphone不买android,为何宁可用4年前的苹果6s,也不买两三千的安卓呢?理由很真实...
- 我的家计算机教学反思,《我的家》教学反思
- Caused by: java.lang.IllegalStateException: StandardEngine[Tomcat].StandardHost[localhost].TomcatEmb
- 215. 数组中的第K个最大元素 BFPRT最牛解法