[BZOJ2440][中山市选2011]完全平方数(莫比乌斯函数,二分)
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2440
题意:RT。
是https://wenku.baidu.com/view/fbec9c63ba1aa8114431d9ac.html中的一道例题,拿来切了练手。
mark一个博客:http://blog.csdn.net/acdreamers/article/details/8542292
1 #include <bits/stdc++.h> 2 using namespace std; 3 4 typedef long long LL; 5 const LL maxn = 1001000; 6 LL vis[maxn], mu[maxn], prime[maxn], cnt; 7 LL n; 8 9 void init() { 10 memset(vis,0,sizeof(vis)); 11 mu[1] = 1; 12 cnt = 0; 13 for(LL i=2; i<maxn; i++) { 14 if(!vis[i]) { 15 prime[cnt++] = i; 16 mu[i] = -1; 17 } 18 for(LL j=0; j<cnt&&i*prime[j]<maxn; j++) { 19 vis[i*prime[j]] = 1; 20 if(i%prime[j]) mu[i*prime[j]] = -mu[i]; 21 else { 22 mu[i*prime[j]] = 0; 23 break; 24 } 25 } 26 } 27 } 28 29 bool check(LL x) { 30 LL xx = (LL)sqrt(x); 31 LL ret = 0; 32 for(LL i = 1; i <= xx; i++) { 33 ret += mu[i]*(LL)(x/(i*i)); 34 } 35 return ret < n; 36 } 37 38 signed main() { 39 // freopen("in", "r", stdin); 40 init(); 41 LL T; 42 scanf("%lld", &T); 43 while(T--) { 44 scanf("%lld", &n); 45 LL lo = 1, hi = (LL)1e10; 46 LL ret = 0; 47 while(lo <= hi) { 48 LL mid = (lo + hi) >> 1; 49 if(check(mid)) lo = mid + 1; 50 else hi = mid - 1; 51 } 52 printf("%lld\n", lo); 53 } 54 return 0; 55 }
转载于:https://www.cnblogs.com/kirai/p/7247964.html
[BZOJ2440][中山市选2011]完全平方数(莫比乌斯函数,二分)相关推荐
- bzoj2440 [中山市选2011]完全平方数 容斥+莫比乌斯函数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些 数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而 这丝毫不影响他对其他数的热爱. 这 ...
- BZOJ2440 [中山市选2011]完全平方数
Description 小 X 自幼就很喜欢数.但奇怪的是,他十分讨厌完全平方数.他觉得这些数看起来很令人难受.由此,他也讨厌所有是完全平方数的正整数倍的数.然而这丝毫不影响他对其他数的热爱. 这天 ...
- bzoj2440:[中山市选2011]完全平方数
Pre 第一次提交的时候\(WA\)了. 百思(死)不得其解(求根公式:我还活着呢!). 最后发现交错题了(尴尬). 不过这题转换到莫比乌斯函数有一点神奇巧妙. 还有一些小细节需要注意. Soluti ...
- 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). [题解] 首先题解并不是我独 ...
- bzoj 2440: [中山市选2011]完全平方数(二分+莫比乌斯函数)
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 4018 Solved: 1952 [Submit][S ...
- BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 3028 Solved: 1460 [Submit][S ...
- bzoj 2440: [中山市选2011]完全平方数【莫比乌斯函数+二分】
二分答案,然后用莫比乌斯函数作为容斥系数,计算当前枚举的mid内有几个满足要求的数 #include<iostream> #include<cstdio> #include&l ...
- BZOJ 2440: [中山市选2011]完全平方数
Time Limit: 10 Sec Memory Limit: 128 MB Submit: 5160 Solved: 2515 [Submit][Status][Discuss] Descript ...
最新文章
- 关系型数据库管理系统和SQL介绍
- python默认参数 可变对象_当心Python函数可变默认参数(list,set,dict…)的陷阱
- Linux web服务器初始化设置
- STATISTICS TIME ON
- SAP CRM BSP URL 302 redirect
- linux下缓存命中测试,linux为什么报表缓存缓存未命中?_linux_开发99编程知识库...
- Lightroom 如何安装lrplugin格式插件
- 向量范数证明例题_第八课:向量的范数
- linux目录文件变化,Linux下监测目录或文件的变化---inotify
- mysql查询id为5的数据类型_MySql学习05---数据类型
- laravel手册链接
- 2021年常规赛NBA球员数据分析
- 计算机专业实训图片,实训一图片的简单处理_计算机软件及应用_IT计算机_专业资料...
- 0045-量化第十天:QMT-以DMA指标为例调用系统指标
- 错误:.TemplateProcessingException: Exception evaluating SpringEL expression:
- 这是我见过的最垃圾的代码,没有之一
- APICloud 原生模块、H5模块、多端组件使用教程
- xshell报错-要继续使用此程序,您必须应用最新的更新或使用新版本
- 2021年中国驴行业存栏量延续下降走势,存栏基数降低使得驴肉供应减少,价格持续走强[图]
- 华为起诉最新进展,国内航司暂停运行有关客机,收买家庭不得继续抚养被解救儿童,脸书创始人后悔没早学微信,这就是今天的大新闻...
热门文章
- 第二章 HTML5存储 JSON格式
- Python入门--闭包,工程函数
- Python入门--面向过程,面向对象,类与对象
- 计算机组成原理—双口RAM和多模块存储器
- RayMarching3:组合与变幻
- Codecraft-18 and Codeforces Round #458: D. Bash and a Tough Math Puzzle(线段树)
- bzoj 4260: Codechef REBXOR(01字典树)
- bzoj 1121: [POI2008]激光发射器SZK
- bzoj 1691: [Usaco2007 Dec]挑剔的美食家(multiset贪心)
- [高光谱] GitHub开源项目Hyperspectral-Classification的解析