AC 数的范围 + 数的三次方根（二分）

给定一个按照升序排列的长度为 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]中包含答案点res

1.初始

显然[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 (向下取整) 一定小于原来的 r

if语句分支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;
}