BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]
2440: [中山市选2011]完全平方数
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 3028 Solved: 1460
[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
1
13
100
1234567
Sample Output
19
163
2030745
HINT
对于 100%的数据有 1 ≤ Ki ≤ 10^9, T ≤ 50
求第k个无平方因子数
二分这个数mid
小于sqrt(mid)的质数都可能成为平方因子,而一个数位平方因子数必定含有一个质数的组合(不一定是几个质数)的平方
根据容斥原理,[1,mid]中无平方因子数的个数为
- 0个质数乘积的平方的倍数的数的数量(1的倍数)
- -每个质数的平方的倍数的数的数量(9的倍数,25的倍数,...)
- +每2个质数乘积的平方的倍数的数的数量(36的倍数,100的倍数,...)-...
也就是容斥原理的变种“奇负偶正”
对于质因子的组合p,它的倍数的个数为mid/(p*p)
只有质因子的次数都是1才会用到,正好是莫比乌斯函数.....
#include <iostream> #include <cstdio> #include <algorithm> #include <cstring> #include <cmath> using namespace std; typedef long long ll; const int N=50000; inline int read(){char c=getchar();ll x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; } int n; bool notp[N]; int p[N],mu[N]; void sieve(){mu[1]=1;for(int i=2;i<=N-1;i++){if(!notp[i]) p[++p[0]]=i,mu[i]=-1;for(int j=1;j<=p[0]&&i*p[j]<=N-1;j++){int t=i*p[j];notp[t]=1;if(i%p[j]==0){mu[t]=0;break; }mu[t]=-mu[i];}} } int cal(int x){int ans=0,m=sqrt(x);for(int i=1;i<=m;i++) ans+=x/(i*i)*mu[i];return ans; } int sol(){int l=n,r=n<<1,ans=-1;while(l<=r){ll mid=l+((r-l)>>1),sum=cal(mid);//printf("hi %d %d\n",mid,sum);if(sum<n) l=mid+1;else ans=mid,r=mid-1;}return ans; } int main(){sieve();int T=read();while(T--){n=read();printf("%d\n",sol());} }
BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]相关推荐
- bzoj 2440: [中山市选2011]完全平方数【莫比乌斯函数+二分】
二分答案,然后用莫比乌斯函数作为容斥系数,计算当前枚举的mid内有几个满足要求的数 #include<iostream> #include<cstdio> #include&l ...
- bzoj 2440: [中山市选2011]完全平方数(二分+莫比乌斯函数)
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 4018 Solved: 1952 [Submit][S ...
- BZOJ 2440: [中山市选2011]完全平方数
Time Limit: 10 Sec Memory Limit: 128 MB Submit: 5160 Solved: 2515 [Submit][Status][Discuss] Descript ...
- [BZOJ2440][中山市选2011]完全平方数(莫比乌斯函数,二分)
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2440 题意:RT. 是https://wenku.baidu.com/view/fbec9 ...
- bzoj2440 [中山市选2011]完全平方数 容斥+莫比乌斯函数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些 数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而 这丝毫不影响他对其他数的热爱. 这 ...
- 2017.3.11[bzoj2440][中山市选2011]完全平方数
线性筛+莫比乌斯函数裸题 http://www.lydsy.com/JudgeOnline/problem.php?id=2440 考虑二分答案,对于每个数看它是第几个,即对于这个数看比他小的有多少个 ...
- bzoj2440: [中山市选2011]完全平方数
自己写的第一个博客......... BZOJ 2440 [题意] 求第K个约数不含平方数的数 (1<=k<=10^9), 共有T组数据(T<=50). [题解] 首先题解并不是我独 ...
- BZOJ2440 [中山市选2011]完全平方数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而这丝毫不影响他对其他数的热爱. 这天 ...
- bzoj2440:[中山市选2011]完全平方数
Pre 第一次提交的时候\(WA\)了. 百思(死)不得其解(求根公式:我还活着呢!). 最后发现交错题了(尴尬). 不过这题转换到莫比乌斯函数有一点神奇巧妙. 还有一些小细节需要注意. Soluti ...
最新文章
- Android自定义ViewGroup基本步骤
- 怎样用c语言写一个系统,用C语言写关于操作系统的一个问题。
- JAVA(利用jsp+javabean+servlet)实现简易计算器
- oracle 隐藏视图定义,【学习笔记】show hidden parameter 创建查看隐藏参数视图
- 数论--康托展开与逆康托展开模板
- 阅读《C陷阱与缺陷》的知识增量
- Linux下MySQL server和client安装
- POJ 1321-棋盘问题-简单搜索DFS
- mysql报1840_mysql 帮助手册 以及 warning: World-writable config file 以及 ERROR 1840 (HY000) at line 24:...
- 【辅助开发】游戏辅助开发全流程-golang
- ①万字《详解canvas api画图》小白前端入门教程(建议收藏)
- 解决:“redis.clients.jedis.exceptions.JedisClusterMaxAttemptsException: No more cluster attempts left”
- 开放、数字化、创新、合作……华为云发力ing!
- 热门智力题 过桥问题和倒水问题
- 织梦ajax登录界面,dede织梦后台登陆成功后又跳转到登陆页面
- ios给按钮图片根据不同的主题更改颜色
- 如何使用html、css制作一个期末作业网站【羽毛球体育运动主题html网页设计】
- 一周玩转示波器(六)
- html中用CSS设置背景
- MT360:工业级无线PDA(1维/2维条码/RFID/GPS/GSM/GPRS/WiFi/蓝牙)
热门文章
- 【Qt】ubuntu QtCreator的pro文件中使用pkg-config
- 如何使用matlab的siso,MIMO-SISO-MATLAB-program 和 的 仿真源代码,非常详细,可以参考学习 267万源代码下载- www.pudn.com...
- cocos creator 安卓原生平台环境_竞技对抗小游戏单挑篮球开发历程 | Cocos技术派第12期...
- 查看linux主机是否安装宋体码,Linux 安装宋体字体的简单办法
- centos java yum_CentOS7 使用yum命令安装Java SDK
- gff文件_GFF格式说明
- 16 导出pcb各网络的布线长度_设计PCB流程
- Java项目:嘟嘟校园一卡通系统(java+JSP+Servlet+html+css+JavaScript+JQuery+Ajax+mysql)
- 【Linux】Linux简单操作之vi与vim编辑器
- Jrebel 热部署插件的使用和破解