文章目录
二分法
二分法(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 块巧克力分给小朋友们。
切出的巧克力需要满足:
- 形状是正方形,边长是整数
- 大小相同
例如一块 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;
}