给定一个按照升序排列的长度为 n 的整数数组,以及 q 个查询。
对于每个查询,返回一个元素 k 的起始位置和终止位置(位置从 0 开始计数)。
如果数组中不存在该元素,则返回 -1 -1
。
输入格式
第一行包含整数 n 和 q,表示数组长度和询问个数。
第二行包含 n 个整数(均在 1∼10000 范围内),表示完整数组。
接下来 q 行,每行包含一个整数 k,表示一个询问元素。
输出格式
共 q 行,每行包含两个整数,表示所求元素的起始位置和终止位置。
如果数组中不存在该元素,则返回 -1 -1
。
数据范围
1≤n≤100000
1≤q≤10000
1≤k≤10000
输入样例:
6 3
1 2 2 3 3 4
3
4
5
输出样例:
3 4
5 5
-1 -1
#include<bits/stdc++.h>
using namespace std;
int n, k;
const int N = 1e5 + 10;
int a[N];
int main(){
cin >> n >> k;
for(int i = 0; i < n; i ++ ){
scanf("%d", &a[i]);
}
while(k -- ){
int x;
cin >> x;
int l = 0, r = n - 1;
while(l < r){
int mid = l + r >> 1;
if(a[mid] >= x)
r = mid;
else
l = mid + 1;
}
if(a[l] != x)
cout << "-1 -1" <<endl;
else{
cout << l << " " ;
int l = 0, r = n - 1;
while(l < r){
int mid = l + r + 1 >> 1;
if(a[mid] <= x)
l = mid;
else
r = mid - 1;
}
cout << l << endl;
}
}
return 0;
}
----------------------------------------------------------------------------------------------
y总二分模板
原文链接:https://www.acwing.com/blog/content/31/
一个题目,如果一个区间具有单调性质,那么一定可以二分,但是如果说这道题目没有单调性质,而是具有某种区间性质的话,我们同样可以使用二分. ——y总
二分模板一共有两个,分别适用于不同情况。
算法思路:假设目标值在闭区间[l, r]中, 每次将区间长度缩小一半,当l = r时,我们就找到了目标值。版本1
当我们将区间[l, r]划分成[l, mid]和[mid + 1, r]时,其更新操作是r = mid或者l = mid + 1;,计算mid时不需要加1。C++ 代码模板:
int bsearch_1(int l, int r) { while (l < r) { int mid = l + r >> 1; if (check(mid)) r = mid; else l = mid + 1; } return l; }
版本2
当我们将区间[l, r]划分成[l, mid - 1]和[mid, r]时,其更新操作是r = mid - 1或者l = mid;,此时为了防止死循环,计算mid时需要加1。C++ 代码模板:
int bsearch_2(int l, int r) { while (l < r) { int mid = l + r + 1 >> 1; if (check(mid)) l = mid; else r = mid - 1; } return l; }
不加1若l=r-1时,r=l+1 ,那么就是(l+l+1)/2向下取整就是l了l=l死循环
归结上面的两种二分方法,步骤为:
先写一个check函数
判定在check的情况下(true和false的情况下),如何更新区间。
在check(m)==true的分支下是:
l=mid的情况,中间点的更新方式是m=(l+r+1)/2
r=mid的情况,中间点的更新方式是m=(l+r)/2
这种方法保证了:
1. 最后的l==r
2. 搜索到达的答案是闭区间的,即a[l]是满足check()条件的。
---------------------------------------------------------------------------------------------------------------------------------
分治算法
原文链接:https://www.acwing.com/solution/content/16798/
二分算法属于分治算法,有三个步骤
整数二分算法
右分界点
此模板适用于寻找右边部分的分界点,即右分界点举例来说,对于排好序数组 {1 2 3 4 5 6}
对于二分的性质 ">=3",数组分成 {1 2} 和 {3 4 5 6} 两部分
左边部分不满足性质 ">=3",右边部分满足性质 ">=3"
算法寻找的点就是数组元素 3 的下标,即所谓的右分界点
如果寻找左分界点(例:对于性质"<3",左分界点就是数组元素2的下标),该算法就要做相应的修改
模板
int bsearch_1(int l, int r) { //第二步:递归处理子问题,用while循环来实现 while (l < r) { //第一步:分解成子问题,这是二分的核心--范围减半 int mid = l + r >> 1; if (check(mid)) r = mid; //向左边找 //if判断mid是否满足性质,注意该性质会划分数组的右边部分 else l = mid + 1; //向右边找 } //第三步:合并子问题.对二分算法来说,不需要这一步 return l; //l就是寻找的右分界点,如果数组中没有要找的点,l的值就是r,但这是一个错误答案 }
待证问题:循环结束后的l就是要找的点(默认包含答案)
证明
循环不变式:[l..r]中包含答案点res1.初始
显然[l..r]包含答案点res
2.保持
假设某轮循环开始之前,[l..r]包含答案点res
执行循环体
int mid = l + r >> 1 mid是向下取整得到的
if语句分支1:
如果mid满足性质,那么说明res在[l..mid]间(包括mid本身),令 r = mid
if语句分支2:
如果 mid 不满足性质,说明mid在左边部分,res在[mid+1..r]间,令l = mid + 1
∴ l和r更新之后,下一轮循环开始之前,循环不变式依然成立
3.终止
循环终止时, l >= r
易知 l不可能比 r 大 , 故 l = r
∴ 根据循环不变式,l 就是答案点 res
边界分析
问题:为什么 mid 是向下取整得到的,即 mid = l + r >> 1. 而不是向上取整,即 mid = l + r + 1 >> 1
答:mid向下取整是为了缩小范围,避免造成无限划分
证明:
if语句分支1: r = mid = l + r >> 1 (向下取整) 一定小于原来的 rif语句分支2: l = mid + 1 一定大于原来的 l
所以,mid向下取整的话,就不会造成无限划分
注:对于二分的另一种情况,即寻找左分界点, mid就需要向上取整了
寻找左分界点
int bsearch_2(int l, int r) { while (l < r) { int mid = l + r + 1 >> 1; // mid 向上取整 if (check(mid)) l = mid; else r = mid - 1; } return l; }
浮点数二分
浮点数二分算法和整数的两种算法一摸一样,只不过不需要考虑向上取整还是向下取整,以及边界处理
浮点数二分的本质也是边界, 唯一区别是浮点数没有整除, 区间长度可以严格的缩小一半
当区间长度足够小时, 便可以认为是一个数
浮点数寻找右分界点
double bsearch_3(double l, double r) { while (r - l > eps) eps 取决于题目对精度的要求, 一般比所求精度高 2 { double mid = (l + r) / 2; //这里不需要考虑取整,因为是浮点数 if (check(mid)) r = mid; //这里需不需要考虑 +1 -1 之类的,因为是浮点数 else l = mid; } return l; }
浮点数寻找左分界点
double bsearch_4(double l, double r) { while (r - l > eps) { double mid = (l + r) / 2; if (check(mid)) l = mid; //只有这两个地方与上面的模板相反,其余都一样 else r = mid; } return l; }
---------------------------------------------------------------------------------------------------------------------------------
给定一个浮点数 n,求它的三次方根。
输入格式
共一行,包含一个浮点数 n。
输出格式
共一行,包含一个浮点数,表示问题的解。
注意,结果保留 6 位小数。
数据范围
−10000≤n≤10000
输入样例:
1000.00
输出样例:
10.000000
#include<bits/stdc++.h>
using namespace std;
double n;
int main(){
cin >>n;
double l = - 10000, r = 10000;
while(r - l > 1e-8){
double mid = (l + r) / 2;
if(mid * mid * mid >= n)
r = mid;
else
l = mid;
}
printf("%.6lf\n", l);
return 0;
}
话说看到有很不讲武德的版本ORZ,那函数见都没见过QAQ
cbrt函数:cube root 立方根 cube 立方体,root 根
#include<cstdio>
#include<cmath>
int main(){
double a;
scanf("%lf",&a);
printf("%.6lf",cbrt(a));
return 0;
}