算法基础--二分(Binary Search)

算法基础–二分(Binary Search)

引入

高中数学中学习过求零点的二分法，算法二分与其大同小异，主要核心思想是一样的，其都是利用要查的目标左右的性质不同。

很多人对二分算法很头疼，可能是因为二分的边界问题很难掌握，也可能是判断条件很难写，这些问题在我刚接触到二分是也

很苦恼，所以我写下这篇文章，希望对你理解二分有所帮助。

介绍

在计算机科学中，二分算法，也称折半搜索，是一种在有序数组中查找某一特定元素的搜索算法。搜索过程从数组中间开始，如

果中间元素正好等于查找元素，则搜索过程结束(最理想状态)；如果要查找的元素大于或小于中间元素，则再次从大于或小于的半

边查找，递归循环直至找到或数组为空。这种搜索算法每一次将搜索范围缩小为一半。故其算法复杂度为O(log n),故也称对数搜索。

对于二分不仅仅适用于查找某个准确的数，更多的应用在于查找大于或小于某个值的最小值或最大值。

如果还不能理解的话，建议手动玩一下二分。

模板

用于查找>=x的最小值（向右查找）

while(l<r)
{
	mid=(l+r)>>1;
	//>>右移位操作，将左边操作数二进制形式右移一位，详参http://t.csdn.cn/o3Vp8
	if(check(mid))
		r=mid;
	else
		l=mid+1;//+1后面会提到
}
cout<<l<<endl;

用于查找<=x的最大值（向左查找）

while(l<r)
{
    mid=(l+r+1)>>1;
    if(check(mid))//check()判断函数，根据题目编写
        l=mid;
    else
        r=mid-1;//对于为什么-1，后面会提到
}
cout<<l<<endl;

这两个模板一定要牢牢记住哦！！！

用于浮点二分

whil(r-l<1e-5)//精度判断
{
    mid=(l+r)/2;
    if(check(mid))
        l=mid;
    else
        r=mid;//浮点数没有像整数一样的取整，故不用在加一减一
}

不要只死记模板，遇题还要灵活应用

STL

C++中的STL中有二分算法lower_bound(对应模板一)和upper_bound（模板二）

分类

二分查找

条件:查找区间是有序的。

光说不练纸上谈兵，先看一道例题

例1：查找

题目描述

输入 $n$ 个不超过 $10^9$ 的单调不减的（就是后面的数字不小于前面的数字）非负整数 $a_1,a_2,\dots ,a_n$，然后进行 $m$ 次询问。对于每次询问，给出一个整数 $q$，要求输出这个数字在序列中第一次出现的编号，如果没有找到的话输出 $-1$ 。

输入格式

第一行 $2$ 个整数 $n$ 和 $m$，表示数字个数和询问次数。

第二行 $n$ 个整数，表示这些待查询的数字。

第三行 $m$ 个整数，表示询问这些数字的编号，从 $1$ 开始编号。

输出格式

输出一行，$m$ 个整数，以空格隔开，表示答案。

样例 #1

样例输入 #1

11 3
1 3 3 3 5 7 9 11 13 15 15
1 3 6

样例输出 #1

1 2 -1

提示

数据保证，$1\le n\le 10^6$，$0\le a_i,q\le 10^9$，$1\le m\le 10^5$

本题输入输出量较大，请使用较快的 IO 方式。

code

#include<iostream>
using namespace std;
int n,m,a[1000010];
int main()
{
	int l,r,mid,q;
	cin>>n>>m;
	for(int i=0;i<n;i++){
		cin>>a[i];
	}
	while(m--)
	{
		cin>>q;
		l=0;r=n;
		while(l<r)
		{
			mid=l+r>>1;
			if(a[mid]>=q)//找第一次出现的位置，故运用模板一
				r=mid;
			else
				l=mid+1;//如果舍弃左半边，即a[mid]<q,故左半边从mid的下一位开始
		}
		if(a[l]==q)
			cout<<l+1<<" ";//找到了输出位置
		else
			cout<<-1<<" ";//未找到输出-1
	}
	return 0;
}

学会了自己也要练习

1.A-B 数对

2.烦恼的高考志愿

二分答案

二分答案就是已知答案区间，在答案区间上二分出正确答案。

**例2：**砍树

题目描述

伐木工人 Mirko 需要砍 $M$ 米长的木材。对 Mirko 来说这是很简单的工作，因为他有一个漂亮的新伐木机，可以如野火一般砍伐森林。不过，Mirko 只被允许砍伐一排树。

Mirko 的伐木机工作流程如下：Mirko 设置一个高度参数 $H$（米），伐木机升起一个巨大的锯片到高度 $H$，并锯掉所有树比 $H$ 高的部分（当然，树木不高于 $H$ 米的部分保持不变）。Mirko 就得到树木被锯下的部分。例如，如果一排树的高度分别为 $20,15,10$ 和 $17$，Mirko 把锯片升到 $15$ 米的高度，切割后树木剩下的高度将是 $15,15,10$ 和 $15$，而 Mirko 将从第 $1$ 棵树得到 $5$ 米，从第 $4$ 棵树得到 $2$ 米，共得到 $7$ 米木材。

Mirko 非常关注生态保护，所以他不会砍掉过多的木材。这也是他尽可能高地设定伐木机锯片的原因。请帮助 Mirko 找到伐木机锯片的最大的整数高度 $H$，使得他能得到的木材至少为 $M$ 米。换句话说，如果再升高 $1$ 米，他将得不到 $M$ 米木材。

输入格式

第 $1$ 行 $2$ 个整数 $N$ 和 $M$，$N$ 表示树木的数量，$M$ 表示需要的木材总长度。

第 $2$ 行 $N$ 个整数表示每棵树的高度。

输出格式

$1$ 个整数，表示锯片的最高高度。

样例 #1

样例输入 #1

4 7
20 15 10 17

样例输出 #1

15

样例 #2

样例输入 #2

5 20
4 42 40 26 46

样例输出 #2

36

提示

对于 $100\%$ 的测试数据，$1\le N\le 10^6$，$1\le M\le 2\times 10^9$，树的高度 $<10^9$，所有树的高度总和 $>M$。

分析：

锯片高度一定在最高树之下[0,max],即答案在此区间内，其要求满足M米木材的H最大值，故其是向左寻找（与你想的不同），因为高度与得到树木长度成反比。教你一个判断方向的秘诀，不考虑求最值，看那边可以满足条件，就向那边查找。这题锯片高度为0时可以满足得到的树木长度大于M，故左边满足条件，故向左查找。

code

#include<iostream>
#include<algorithm>
using namespace std;
int n,h,a[1000010];
long long m;
bool check(int aa)//判断函数
{
	int l=1,r=n,mid;
	long long temp=0;
	while(l<r)//对第一个高于锯齿高度树木进行查找，向左查找模板
	{
		mid=l+r>>1;
		if(a[mid]>=aa)
			r=mid;
		else
			l=mid+1;
	}
	for(int i=l;i<=n;i++){
		temp+=a[i]-aa;
	}
	if(temp>=m)
		return true;
	else
		return false;
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=n;++i)
		cin>>a[i];
	sort(a+1,a+1+n);//对树木高度进行排序
	int l=0,r=a[n],mid=0;
	while(l<r)//向右查找模板
	{
		mid=l+r+1>>1;
		if(check(mid))
			l=mid;
		else
			r=mid-1;
	}
	cout<<l;
}

练习：

1.木材加工

2.路标设置

3.数段分段

4.跳石头

浮点二分

显而易见，浮点二分的不同就是操作数为实数

例3：银行贷款

题目描述

当一个人从银行贷款后，在一段时间内他（她）将不得不每月偿还固定的分期付款。这个问题要求计算出贷款者向银行支付的利率。假设利率按月累计。

输入格式

三个用空格隔开的正整数。

第一个整数表示贷款的原值，第二个整数表示每月支付的分期付款金额，第三个整数表示分期付款还清贷款所需的总月数。

输出格式

一个实数，表示该贷款的月利率（用百分数表示），四舍五入精确到 $0.1\%$。

样例 #1

样例输入 #1

1000 100 12

样例输出 #1

2.9

提示

数据保证，$1\le$ 贷款的原值，分期付款金额 $\le 2^{31}-1$，$1\le$ 月数 $\le 3000$。

自己先思考下在看AC代码

code

#include<iostream>
#include<iomanip>
using namespace std;
long long  m,n,num;
bool check(double mid)
{
	double temp=m;
	for(int i=0;i<num;++i)
	{
		temp=temp*(mid+1)-n;
	}
	if(temp<=0)
		return true;
	else
		return false;
}
int main()
{
	cin>>m>>n>>num;
	double l=0,r=5,mid;
	while(r-l>1e-5)//精度控制
	{
		mid=(l+r)/2;
		if(check(mid))
			l=mid;
		else
			r=mid;
	}
	cout<<fixed<<setprecision(1)<<l*100;//控制小数位数
}

练习：

1.一元三次方程求解

总结：

相信看到这里你对二分算法有了一定的掌握，本文练习AC代码，私信发

第一次写博客，还请大佬多多指导，有什么问题提出来

以后还会更新其他算法基础内容，请多多支持！！！