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]完全平方数相关推荐

  1. BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 3028  Solved: 1460 [Submit][S ...

  2. bzoj 2440: [中山市选2011]完全平方数(二分+莫比乌斯函数)

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 4018  Solved: 1952 [Submit][S ...

  3. bzoj 2440: [中山市选2011]完全平方数【莫比乌斯函数+二分】

    二分答案,然后用莫比乌斯函数作为容斥系数,计算当前枚举的mid内有几个满足要求的数 #include<iostream> #include<cstdio> #include&l ...

  4. bzoj2440: [中山市选2011]完全平方数

    自己写的第一个博客......... BZOJ 2440 [题意] 求第K个约数不含平方数的数 (1<=k<=10^9), 共有T组数据(T<=50). [题解] 首先题解并不是我独 ...

  5. BZOJ2440 [中山市选2011]完全平方数

    Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而这丝毫不影响他对其他数的热爱.  这天 ...

  6. 2017.3.11[bzoj2440][中山市选2011]完全平方数

    线性筛+莫比乌斯函数裸题 http://www.lydsy.com/JudgeOnline/problem.php?id=2440 考虑二分答案,对于每个数看它是第几个,即对于这个数看比他小的有多少个 ...

  7. [BZOJ2440][中山市选2011]完全平方数(莫比乌斯函数,二分)

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2440 题意:RT. 是https://wenku.baidu.com/view/fbec9 ...

  8. bzoj2440 [中山市选2011]完全平方数 容斥+莫比乌斯函数

    Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些 数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而 这丝毫不影响他对其他数的热爱. 这 ...

  9. [BZOJ 2438] [中山市选2011]杀人游戏 Tarjan缩点

    这个题很容易想到正解就是缩点找入度为零的点,那么我们考虑一种特殊情况就是,一个入度为零的点我们不访问他就知道他是不是凶手,那么这样的话就是:I. 他是一个真·孤立的点 II. 他在图里但是在他的强联通 ...

最新文章

  1. NSA释出逆向工程工具GHIDRA针对恶意软件逆向分析
  2. asp.net 6中的mini api和mvc api性能对比
  3. iphone NSTimer
  4. property属性学习
  5. 【原创工具 | NetSM】开源跨平台命令行网速监测(纯 Python 开发)
  6. Java UDP通信详解:单播、广播、组播
  7. Mac M1 安装JDK8
  8. PTC:能源互联网“双子星”,引领风电企业研发数字化变革
  9. php麻将,php实现麻将牌型的胡牌检测算法
  10. 致远oa系统unix 服务器,致远oa手机客户端服务器
  11. HCIA-GaussDB 华为认证数据库工程师
  12. 蔚来、小鹏、理想自动驾驶能力的纵向演进与横向比较
  13. NOIP2018初赛
  14. 机器学习系列(二)——回归模型
  15. 单目标优化、多目标优化
  16. ImageJ macro 批量vsi转tif
  17. 检测恶意软件分类模型中的概念漂移
  18. [开源名人访谈录] Philippe Gerum
  19. Image Processing for Embedded Devices 4
  20. era翻译_era是什么意思_era的翻译_音标_读音_用法_例句_爱词霸在线词典

热门文章

  1. if shell 条件语句简写_Shell if条件判断语句用法详解
  2. 【EOJ Monthly 2019.02 - E】中位数(二分 ,中位数 ,−1/1变换,dp求解DAG最长路)
  3. 笔试题 - 阶乘问题
  4. 【牛客 - 368D】动态连通块(并查集+bitset优化)
  5. Apollo进阶课程⑳丨Apollo感知之旅——机器学习与感知的未来
  6. mysql判断域为空_MySQL EXPLAIN 字段说明
  7. 买iphone不买android,为何宁可用4年前的苹果6s,也不买两三千的安卓呢?理由很真实...
  8. 我的家计算机教学反思,《我的家》教学反思
  9. Caused by: java.lang.IllegalStateException: StandardEngine[Tomcat].StandardHost[localhost].TomcatEmb
  10. 215. 数组中的第K个最大元素 BFPRT最牛解法