蓝桥杯：C++二分算法

在基本算法中，二分法的应用非常广泛，它是一种思路简单、编程容易、效率极高的算法。蓝桥杯软件类大赛中需要应用二分法的题目很常见。

二分法有整数二分和实数二分两种应用场景

二分法的概念

二分法的概念很简单，每次把搜索范围缩小为上一次的1/2，直到找到答案为止。

二分法的效率很高，只需计算log(n)次。

下面介绍二分法的模板代码bin_search()函数：

我们用猜数字的例子，先给数组初始化，然后定义你要猜的数，用二分法效率高。

对于二分法的讲解非常细致，都在注释中。

#include<bits/stdc++.h>
using namespace std;
int a[1000];
int bin_search(int *a, int n, int x) {   //在数组a中查找数字x，返回位置
	int left = 0, right = n;	//left 通常初始化为 0，表示搜索范围的左边界是数组的第一个元素；right通常初始化为 n（数组的长度），表示搜索范围的右边界是数组的最后一个元素的下一个位置。
	while (left < right) {
		int mid = left+(right-left)/2;   //mid的标准写法，建议这样写，不能用(left+right)/2,有可能会整数溢出的。 
		if (a[mid] >= x) right = mid;   //x小，在左边，把右边的一半砍掉，这里就不用加1了，我们本身就是大于等于x。 
		else             left = mid+1;	//加1的原因是我们要跳过 a[mid] 这个元素，因为它小于 x，我们要的是等于x的元素 
		cout<<a[mid]<< " ";              //输出猜数的过程    如果你想省略过程，可以注释掉这一行的输出语句。 
	}
	return left;    //返回left所在的索引，不要牵扯到right，避免混淆，right一开始是索引的下一个位置。 
}
int main() {
	int n = 100;
	for(int i=0; i<n; i++) a[i]=i+1;    //赋值，数字1～100
	int test = 54;                      //猜54这个数
	int pos = bin_search(a,n,test);
	cout<<"\n"<<"test="<<a[pos];    
}

bin_search()有3个重要点：区间左端点left、区间右端点right、二分的中位数mid。每次把区间缩小一半，把left或right移动到mid；直到left = right为止，即找到答案所处的位置。

二分法的作用：

二分法可以把一个长度为n的有序序列上O(n)的查找时间优化到O(logn)。

注意应用二分法的前提：序列是有序的，按从小到大或从大到小排序。

无序的序列无法二分，如果是无序的序列，则应该先排序再对其进行二分，先排序再二分，排序的复杂度是O(nlog2(n))，二分的复杂度是O(log2(n))。排序加二分的总复杂度是O(nlog2(n))。如果使用暴力法，直接在无序的n个数里面查找，最多查找n次，复杂度是O(n)的，比先排序再二分快。如果不是查找一个数，而是查找m个数，那么先排序再做m次二分的计算复杂度是O(nlog2(n)+ mlog2(n))，而暴力法的复杂度是O(mn)，此时二分法远好于暴力法。

整数二分

在单调递增序列中查找x或者x的后继：

前面介绍的bin_search()函数就是“在单调递增序列中查找x或者x的后继”的模板代码。

二分函数都是一摸一样的，测试数据可以改一下，看看能不能查找后继：

int main() {
	int n = 100;
	for(int i=0; i<n; i++) a[i]=2*i+2;    //赋值，数字2～200，偶数
	int test = 55;                        //查找55或55的后继
	int pos = bin_search(a,n,test);
	cout<<"test="<<a[pos];//56   55没有，只能找56了。
}

在单调递增序列中查找x或者x的前驱：

#include<bits/stdc++.h>
using namespace std;
int a[1000];
int bin_search2(int *a, int n, int x) {
	int left = 0, right = n;
	while (left < right) {
		int mid = left + (right-left + 1)/2 ;   //+1是为了确保在 left 和 right 之间的元素数量是奇数时，mid 会指向中间元素；当元素数量是偶数时，mid 会指向中间两个元素的右侧那个元素。
		//这样做的原因是，我们希望在存在重复元素时，mid 尽可能向右偏移，从而找到最右侧的那个等于或小于 x 的元素。
		if (a[mid] <= x)  left = mid;
		else  right = mid - 1;
	}
	return left;
}
int main() {
	int n = 100;
	for(int i=0; i<n; i++) a[i]=2*i+2;   //赋值，数字2～200,偶数
	int test = 55;                       //查找55或55的前驱
	int pos = bin_search2(a,n,test);
	cout<<"test="<<a[pos]; //54
}

整数二分例题

例题1.分巧克力

2017年（第八届）省赛，lanqiaoOJ题号99

先试试暴力法：从边长为1开始到最大边长d，每个值都试一遍，一直试到刚好够分的最大边长为止。编程思路：边长初始值d = 1，然后d = 2、3、4……一个一个地试 。

代码：

#include<bits/stdc++.h>
using namespace std;
int h[100010],w[100010];//多申请10个空间 
int n,k;
bool check(int d) {             //检查够不够分
	int num=0;
	for(int i=0; i<n; i++)  num += (h[i]/d)*(w[i]/d);//假如，将6×5的巧克力的长边（6个单位）和宽边（5个单位）分别除以2×2的小正方形的边长（2个单位）。
	//这样可以得到长边可以切出3个2×2的巧克力，宽边可以切出2个2×2的巧克力。
	//接着，将长边和宽边切出的巧克力块数相乘，即3（长边切出的块数）× 2（宽边切出的块数）= 6。所以，一块6×5的巧克力可以切出6块2×2的巧克力。
	if(num>=k) return true;     //够分
	else       return false;    //不够分
}
int main() {
	cin >>n>>k;
	for(int i=0; i<n; i++)  cin>>h[i]>>w[i];  //长宽各自存在各自的数组中 
	int d=1;                    //正方形边长
	while(1) {
		if(check(d))  d++;      //边长从1开始，一个一个地试
		else          break;
	}
	cout << d-1;
	return 0;     //暴力求解只能过75的测试数据 ，最后两个测试数据错了，暂时不知道什么原因
}

整数二分法求解：

#include<bits/stdc++.h>
using namespace std;
int n,k;
const int N=100010;
int h[N],w[N];
bool check(int d) {
	int num=0;
	for(int i=0; i<n; i++)  num += (h[i]/d)*(w[i]/d);
	if(num>=k) return true;      //够分
	else       return false;     //不够分
}
int main() {
	cin >> n >> k;
	for(int i=0; i<n; i++)   cin>>h[i]>>w[i];
	int L=1, R=N;                //R的初值是100010
	//第一种写法：
	while(L<R) {
		int mid=(L+R+1) / 2;      //除以2，向右取整      不会整数溢出，直接L+R
		if(check(mid))  L=mid;   //新的搜索区间是右半部分，R不变，调整L=mid
		else            R=mid-1; //新的搜索区间是左半部分，L不变，调整R=mid–1
	}
	cout << L;
	//第二种写法：
	/*  while(L<R) {
	        int mid=(L+R) / 2;        //除以2，向左取整   不会整数溢出，直接L+R
	        if(check(mid)) L=mid+1;  //新的搜索区间是右半部分，R不变，更新L=mid+1
	        else           R=mid;    //新的搜索区间是左半部分，L不变，更新R=mid
	    }
	    cout << L-1;    */
	return 0;
}

实数二分

与整数二分相比，实数二分的编程就容易多了，不用考虑整数的取整问题。实数二分的模板代码如下。

const double eps = 1e-7;                    //精度。
while(right - left > eps) {                 
	double mid = left+(right-left)/2;
	if (check(mid)) right = mid;          //判定，然后继续二分，check(mid)为true执行此语句
	else            left  = mid;
}