❤️二分算法❤️

二分法

二分法（Bisection method），即假设目标值位于闭区间[l，r]中，每次将区间长度缩小一半，当l = r时，我们就找到了目标值。

对于区间[a，b]上连续不断且f（a）·f（b）<0的函数y=f（x），通过不断地把函数f（x）的零点所在的区间一分为二，使区间的两个端点逐步逼近零点，进而得到零点近似值的方法叫二分法。

如图，零点x0 即为我们所求答案，我们通过判断特定条件不断更新端点a，b的值，将端点不断向零点趋近，当得到近似值误差值小于一定范围内，我可认为此时所求值即为目标答案。

使用二分法的关键：二分后，能否判断出答案所在的区间，而不是给出的数据是否有序。

下面我们来看下Y总给的模板：

整数二分

模板一
当我们将区间[l, r]划分成[l, mid]和[mid + 1, r]时，其更新操作是r = mid或者l = mid + 1;，计算mid时不需要加1。

C++ 代码模板

// 区间[l, r]被划分成[l, mid]和[mid + 1, r]时使用：
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;
}

模板二
当我们将区间[l, r]划分成[l, mid - 1]和[mid, r]时，其更新操作是r = mid - 1或者l = mid;，此时为了防止死循环，计算mid时需要加1。

C++ 代码模板

// 区间[l, r]被划分成[l, mid - 1]和[mid, r]时使用：
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;
}

浮点数二分

判断的精度取决于题目要求，若要求保留k位小数，一般精确到1e^-(k+2)即可

C++ 代码模板

double bsearch_3(double l,double r)
{
	const double eps = 1e-6; // eps表示精度，取决于题目对精度的要求
	while(r-l > eps)
	{
		double mid = l + r >> 1;
		if(check(mid)) 
			r = mid;
		else
			l = mid;
	}
	return l;
}

二分法时间复杂度分析

二分法的操作是将每次区间[l,r]长度缩小为原来的二分之一。当区间内元素只有一个或者元素精度达到既定精度时停止。
第一次二分，区间长度：n/2；
第二次二分，区间长度：n/4；
…
第k次二分，区间长度：n/2^k；
最坏情况下，就剩一个元素，此时n/2^k=1；
k = logn；
所以时间复杂度为O(logn)。

---

acwing 680.剪绳子

有N根绳子，第i根绳子长度为Li，现在需要M根等长的绳子，你可以对N根绳子进行任意裁剪（不能拼接），请你帮忙计算出这M根绳子最长的长度是多少。

输入格式

第一行包含2个正整数N、M，表示原始绳子的数量和需求绳子的数量。

第二行包含N个整数，其中第 i 个整数LiLi表示第 i 根绳子的长度。

输出格式

输出一个数字，表示裁剪后最长的长度，保留两位小数。

数据范围

1≤N,M≤100000,
0<Li<109

输入样例：

3 4
3 5 4

输出样例：

2.50

样例解释

第一根和第三根分别裁剪出一根2.50长度的绳子，第二根剪成2根2.50长度的绳子，刚好4根。

题目分析

我们拿到这道题后发现，直接做的话没有办法下手，不知道如何最优裁剪绳子长度，才能达到需求绳子数量，因此我们可以换一种角度，假设我们知道我们要裁剪绳子长度为x，求最多能裁剪多少根绳子n，解决这个问题我们发现很简单，只需要将所有数据遍历一遍，求 ∑ [Li/x] 即可。
将这个问题返回到这道题，若ans(∑ [Li/x]) 从右端无限趋近于m(需求绳子数量) 即可，我们可以通过二分区间[l,r]判断mid = l + r >> 1;的是否符合判断条件（ans >= m）若符合，答案在[mid,r]区间内即更新左端点l = mid；若不符合，答案在[l,mid]区间内即更新右端点r = mid。

C++代码如下：

#include<iostream>

using namespace std;

const int N = 1e5+5;
int a[N];
int n,m;

bool check(double x)
{
    int cnt = 0;
    for(int i = 0; i < n; i ++)
        cnt += a[i]/x;
    return cnt >= m;
}

int main()
{
    cin >> n >> m;
    for(int i = 0; i < n; i++)
        cin >> a[i];
    double l = 0, r = 1e9;
    while(r - l > 1e-4)  //题中保留两位小数，即精确到1e-(k+2)即可
    {
        double mid = (l+r)/2;
        if(check(mid)) l = mid;
        else r = mid;
    }
    printf("%.2lf\n",r);
    return 0;
}

---

---

acwing 1227.分巧克力

儿童节那天有 K 位小朋友到小明家做客。

小明拿出了珍藏的巧克力招待小朋友们。

小明一共有 N 块巧克力，其中第 ii 块是 Hi×Wi 的方格组成的长方形。

为了公平起见，小明需要从这 N 块巧克力中切出 K 块巧克力分给小朋友们。

切出的巧克力需要满足：

1. 形状是正方形，边长是整数
2. 大小相同

例如一块 6×5 的巧克力可以切出 6 块 2×2 的巧克力或者 2 块 3×3 的巧克力。

当然小朋友们都希望得到的巧克力尽可能大，你能帮小明计算出最大的边长是多少么？

输入格式

第一行包含两个整数 N 和 K。

以下 N 行每行包含两个整数 Hi 和 Wi。

输入保证每位小朋友至少能获得一块 1×1 的巧克力。

输出格式

输出切出的正方形巧克力最大可能的边长。

数据范围

1≤N,K≤105,
1≤Hi,Wi≤105

输入样例：

2 10
6 5
5 6

输出样例：

2

题目分析

题意:就是让我们在1~100000中找到一个最大的ans,让其满足，我们能在大巧克力中划分出的小巧克力块数 >=m ;
思路:
假设n块巧克力的长和宽分别保存在两个一维数组a[100000]和b[100000]的n个元素里。
由于一块大小为 xy（长为x，宽为y）的巧克力，可以切割成的边长为mid的正方形巧克力的块数为：(x/mid)(y/mid)， 所以n块巧克力所切割成的总块数idxt可以按照下面的方式计算出来：

bool check(int x)
{
    int res = 0;
    for(int i = 0; i <= n; i ++)
    {
        res += (w[i]/x)*(h[i]/x);
    }
    return res >= k;
}

当idx大于等于m时，就可以切割出能够分给小朋友的巧克力，否则就切割不出。

C++代码

#include<iostream>
#include<algorithm>

using namespace std;

const int N = 1e5+5;

int n,k;
int h[N],w[N];

bool check(int x)
{
    int res = 0;
    for(int i = 0; i <= n; i ++)
    {
        res += (w[i]/x)*(h[i]/x);
    }
    return res >= k;
}

int main()
{
    cin >> n >> k;
    for(int i = 0; i <= n; i ++ )cin>> h[i] >> w[i];

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