二分查找与二分答案

二分模板

由于二分不只是简单的查找一个数。它可以用于查找满足性质和不满足性质的分界点的位置。因此在记模板的时候，我们要用性质来解释。

二分步骤：
1.划分性质
2.描述check在符合性质处的代码语句
3.确定是否+1

举个例子：现在给定一串有序数字，询问n次，分别返回该数字的初始位置和结束位置

1 1 1 2 3 3 3 3 5 7 9

比如数字3的开始位置是4 和 7（下标从0开始）
要二分出来，就必须要找到对应的性质。

注：
1.划分性质的时候要有二段性，即分界点是空集，没有元素的。
2.check函数永远是在描述对象符合性质的条件

起始位置：

在查找起始位置这个问题来说：>=x是满足性质的，<x是不满足性质的。因此check函数是a[mid]>=x

结束位置：

最后考虑是否+1问题。

左缩（r = mid）不加1， 右缩（l = mid）加1

代码如下：

int find1(int 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) return l;
    else return -1;
}

int find2(int x)//最大下标
{
    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;
    }
    if(a[l] == x) return l;
    else return -1;
}

A-B数对

A-B数对
这道题有两种做法，一种是哈希枚举，一种是先枚举再二分。

---

哈希思路：用哈希初始化数组，得到每个数字出现的次数。由于C已知。因此A-B = C可以通过枚举A，得到B。B = A - C，然后通过哈希得到B在数组中出现的次数，加入ans中即可。

时间复杂度分析：哈希初始化为O(N),枚举A时间复杂度为O(N),查询哈希时间复杂度为O(1),因此总的时间复杂度为O(N)

---

二分思路：先枚举A，然后得到对应的B的值，用二分去数组里面查找是否存在B这样的值，并统计出现了多少次（因此要先把数组排序），加入ans即可。

时间复杂度分析：排序O(NlogN）+ 枚举O(N) +二分O(logN)。因此总时间复杂度为O(NlogN)

二分ac代码：

#include <iostream>
#include <algorithm>
using namespace std;

const int N = 200010;

int n, c;
int a[N];
long long res;

int find1(int 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) return l;
	else return -1;
}

int find2(int x)
{
	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;
	}
	if(a[l] == x) return l;
	else return -1;
}

int main()
{
	cin >> n >> c;
	for(int i = 0; i < n; i++) scanf("%d", &a[i]);
	sort(a, a+n);
	for(int i = 0; i < n; i++)
	{
		int x = a[i];
		int l = find1(x - c), r = find2(x - c);
		if(l != -1 && r != -1) res += (r - l + 1);
	}
	cout << res;
}

砍树

砍树
考点：二分。
注：题目要求砍下来的木材尽可能长，因此我们的锯子就要尽可能短。
因此划分的性质是

左边<= x是满足锯子尽可能短的性质的，右边是不满足锯子尽可能短的性质的.
check函数就看锯子是否尽可能短了。锯子尽可能短的体现是木材大于等最小值

---

ac代码

#include <iostream>

using namespace std;

const int N = 1e6 + 10;

int n, m;
int a[N];
int ans;

bool check(int x)
{
	long long sum = 0;
	for(int i = 0; i < n; i++) sum += (a[i] - x > 0 ? a[i] - x : 0);
	if(sum >= m) return true;
	else return false;
}

int find()
{
	int l = 0, r = 1e9;
	while(l < r)
	{
		int mid = l + r + 1 >> 1;
		if(check(mid)) l = mid;
		else r = mid - 1;
	}
	return l;
}

int main()
{
	cin >> n >> m;
	for(int i = 0; i < n; i++) scanf("%d", &a[i]);
	ans = find();
	cout << ans;
}

进击的奶牛

进击的奶牛
二分的难点其实不在于划分性质，难点在于check函数的编写。即如何验证性质的正确性。这一题就体验出来了。

考点:二分
1.题目问的是最大的最近距离，暗示了可以用二分来解决。
2.最大的最近距离就是求x的最大值。因此性质可以划分成<=x和>x
3.x是最近距离，x越小，棚子就能装下越多牛，x越大，棚子就有可能装不下所有牛。
4.因此我们可以用棚子装的牛的数量来判断x是否满足既能满足最大的最近距离，又能放下所有牛的性质
5.从最左边的棚子开始放牛，每隔大于等于x的距离就放牛，可以放下最多牛的数量（贪心思想），这就是check函数的思路

ac代码：

#include <iostream>
#include <algorithm>
using namespace std;

const int N = 100010;
int a[N];

int n, m;

bool check(int x)
{
	int k = 1, last = 1;
	
	for(int i = 1; i < n; i++)
		if(a[i] - last >= x)
		{
			k++;
			last = a[i];	
		}
	
	return k >= m;
}

int find()
{
	int l = a[0], r = a[n - 1];
	while(l < r)
	{
		int mid = l + r + 1 >> 1;
		if(check(mid)) l = mid;
		else r = mid - 1;
	}
	return l;
}

int main()
{
	cin >> n >> m;
	for(int i = 0; i < n; i++) scanf("%d", &a[i]);
	sort(a, a + n);
	int ans = find();
	cout << ans;
}

浮点数二分模板

bool check(double x) {/* ... */} // 检查x是否满足某种性质

double bsearch_3(double l, double r)
{
    const double eps = 1e-6;   // eps 表示精度，取决于题目对精度的要求。
    while (r - l > eps)//也可以写成for(int i = 0;i < 100;i++)这样相当于二分了100次，N/2的100次方
    {
        double mid = (l + r) / 2;
        if (check(mid)) r = mid;
        else l = mid;
    }
    return l;
}

注：由于double没有向下取整，因此不需要考虑是否+1了。

一元三次方程

一元三次方程
考点：浮点数二分
1.三次方程有三个解，如果只二分整一个区间，只能得到一个值。
2.技巧：将整个区间分成[i, i+1)这样的小区间 （题目说了每一个根之间的差值至少为1，因此这样分）。对每一个小区间进行二分
3.区间分为左开右闭的原因是防止出现重复的解
4.如何确定区间有根？零点存在定理
5.如何二分？画几个例子可以发现规律：

f(l) * f(mid) > 0,则l = mid

f(r) * f(mid) > 0,则r = mid 

注：二分只能找到区间内部的根，不能找到端点的根。因此得把端点是否为根也要判断一下。

ac代码：

#include <iostream>
#include <cmath>

using namespace std;

double esp = 1e-4;
double a, b, c, d;

double f(double x)
{
	return a * x * x * x + b * x * x + c * x + d;
}

int main()
{
	cin >> a >> b >> c >> d;
	for(int i = -100; i <= 100; i++)
	{
		double l = i, r = i + 1;
		if(abs(f(l)) <= esp) printf("%.2lf ", l);
		else if(abs(f(r)) <= esp) continue;
		else if(f(l) * f(r) < 0)
		{
			while(r - l > esp)
			{
				double mid = (l + r) / 2;
				if(f(mid) * f(l) > 0) l = mid;
				else r = mid;
			}
			printf("%.2lf ", l);
		}
	}
}