一、简介
到底什么是二分呢?二分二分就是一分为二。简单来说二分就是在有序序列中,通过不断的二分,进而不断地缩小范围去寻找满足我们条件的解。这只是对二分一个狭义上的理解,广义二分其实是如果有一个临界值使得临界值一边的数据满足一种性质,另一边满足另一种性质,即使不是有序的但也可以利用二分去寻找这个临界值。
二分题目主要分为二分查找、二分答案,二分类型分为整数二分、浮点数二分。
二、二分答案
题目要求我们求出某种条件的最大值的最小可能情况或者最小值的最大情况 ,已知答案在一个固定的区间内,但难以通过搜索来找到符合要求的值, 但给定一个值你可以很快的判断它是不是符合要求,这时我们便可对答案进行二分。
1、模板
(1)最小值的最大情况
while(left<=right)
{
int middle=(right+left)/2;
if(check(middle))
{
ans=middle;
left=middle+1;
}
else
right=middle-1;
}
(2)最大值的最小情况
while(left<=right)
{
int middle=(right+left)/2;
if(check(middle))
{
ans=middle;
right=middle-1;
}
else
left=middle+1;
}
(3)浮点答案二分
while(right-left>=eps) //答案所需的精度eps
{
double middle=(l+r)/2;
if(check(middle)) left=middle;
else right=middle;
}
2、例题
(1)洛谷 P2678 [NOIP2015 提高组] 跳石头
题目背景
一年一度的“跳石头”比赛又要开始了!
题目描述
这项比赛将在一条笔直的河道中进行,河道中分布着一些巨大岩石。组委会已经选择好了两块岩石作为比赛起点和终点。在起点和终点之间,有 N N N
块岩石(不含起点和终点的岩石)。在比赛过程中,选手们将从起点出发,每一步跳向相邻的岩石,直至到达终点。为了提高比赛难度,组委会计划移走一些岩石,使得选手们在比赛过程中的最短跳跃距离尽可能长。由于预算限制,组委会至多从起点和终点之间移走 M M M
块岩石(不能移走起点和终点的岩石)。输入格式
第一行包含三个整数 L , N , M L,N,M L,N,M,分别表示起点到终点的距离,起点和终点之间的岩石数,以及组委会至多移走的岩石数。保证 L ≥ 1 L \geq 1 L≥1
且 N ≥ M ≥ 0 N \geq M \geq 0 N≥M≥0。接下来 N N N 行,每行一个整数,第 i i i 行的整数 D i ( 0 < D i < L ) D_i( 0 < D_i < L) Di(0<Di<L), 表示第 i i i
块岩石与起点的距离。这些岩石按与起点距离从小到大的顺序给出,且不会有两个岩石出现在同一个位置。输出格式
一个整数,即最短跳跃距离的最大值。
样例 #1
样例输入 #1
25 5 2
2
11
14
17
21
样例输出 #1
4
提示
输入输出样例 1 说明
将与起点距离为 2 2 2和 14 14 14 的两个岩石移走后,最短的跳跃距离为 4 4 4(从与起点距离 17 17 17 的岩石跳到距离 21 21 21
的岩石,或者从距离 21 21 21 的岩石跳到终点)。数据规模与约定
对于 20 % 20\% 20%的数据, 0 ≤ M ≤ N ≤ 10 0 \le M \le N \le 10 0≤M≤N≤10。 对于 50 % 50\% 50% 的数据, 0 ≤ M ≤ N ≤ 100 0 \le M \le N \le 100 0≤M≤N≤100。 对于 100 % 100\% 100%的数据, 0 ≤ M ≤ N ≤ 50000 , 1 ≤ L ≤ 1 0 9 0 \le M \le N \le 50000,1 \le L \le 10^9 0≤M≤N≤50000,1≤L≤109。
代码
#include <bits/stdc++.h>
using namespace std;
int L,M,N,ans;
int a[50005];
bool check(int middle)
{
int cnt=0;
int pre=0;
int index=0;
for(int i=N;i>=1;--i)
{
if(L-a[i]>=middle)
{
index=i;
break;
}
else
cnt++;
}
for(int i=1;i<=index;++i)
{
if(a[i]-pre>=middle)
pre=a[i];
else
cnt++;
if(cnt>M)
return false;
}
return true;
}
int main()
{
cin>>L>>N>>M;
for(int i=1;i<=N;++i)
cin>>a[i];
int left=0;
int right=L;
while(left<=right)
{
int middle=(right+left)/2;
if(check(middle))
{
ans=middle;
left=middle+1;
}
else
right=middle-1;
}
cout<<ans<<endl;
return 0;
}
(2)洛谷 P3382 【模板】三分法
题目描述
如题,给出一个 N N N 次函数,保证在范围 [ l , r ] [l, r] [l,r] 内存在一点 x x x,使得 [ l , x ] [l, x] [l,x] 上单调增, [ x , r ] [x, r] [x,r]
上单调减。试求出 x x x 的值。输入格式
第一行一次包含一个正整数 N N N 和两个实数 l , r l, r l,r,含义如题目描述所示。
第二行包含 N + 1 N + 1 N+1 个实数,从高到低依次表示该 N N N 次函数各项的系数。
输出格式
输出为一行,包含一个实数,即为 x x x 的值。若你的答案与标准答案的相对或绝对误差不超过 1 0 − 5 10^{-5} 10−5 则算正确。
样例 #1
样例输入 #1
3 -0.9981 0.5
1 -3 -3 1
样例输出 #1
-0.41421
提示
对于 100 % 100\% 100% 的数据, 6 ≤ N ≤ 13 6 \le N \le 13 6≤N≤13,函数系数均在 [ − 100 , 100 ] [-100,100] [−100,100] 内且至多 15 15 15
位小数, ∣ l ∣ , ∣ r ∣ ≤ 10 |l|,|r|\leq 10 ∣l∣,∣r∣≤10 且至多 15 15 15 位小数。 l ≤ r l\leq r l≤r。【样例解释】
如图所示,红色段即为该函数 f ( x ) = x 3 − 3 x 2 − 3 x + 1 f(x) = x^3 - 3 x^2 - 3x + 1 f(x)=x3−3x2−3x+1 在区间 [ − 0.9981 , 0.5 ] [-0.9981, 0.5] [−0.9981,0.5] 上的图像。
当 x = − 0.41421 x = -0.41421 x=−0.41421 时图像位于最高点,故此时函数在 [ l , x ] [l, x] [l,x] 上单调增, [ x , r ] [x, r] [x,r] 上单调减,故 x = − 0.41421 x = -0.41421 x=−0.41421,输出 − 0.41421 -0.41421 −0.41421。
代码
#include <bits/stdc++.h>
using namespace std;
int N;
double l,r;
double a[15];
const double eps=1e-6;
double jisuan(double x)
{
double sum=0;
for(int i=1;i<=N+1;++i)
sum=sum*x+a[i];
return sum;
}
int main()
{
cin >> N >> l >> r;
for (int i = 1; i <= N + 1; ++i)
cin >> a[i];
while(r-l>=eps)
{
double m=(l+r)/2;
if(jisuan(m+eps)>=jisuan(m-eps)) l=m;
else r=m;
}
printf("%.5lf",l);
return 0;
}
三、二分查找
对一个已知的有序数据集上进行二分的查找。二分查找也称折半查找,顾名思义,就是每次查找去掉不符合条件的一半区间,直到找到答案(整数二分)或者和答案十分接近(浮点二分)。
1、模板
(1)最简单二分查找,尽可能往左边寻找,找到数第一次出现的位置。
while(left<=right)
{
int middle=(left+right)>>1;
if(a[middle]>num) right=middle-1;
else if(a[middle]<num) left=middle+1;
else
{
ans=middle;
right=middle-1;
}
}
(2)最简单二分查找,尽可能往右边寻找,找到数最后一次出现的位置。
while(left<=right)
{
int middle=(left+right)>>1;
if(a[middle]>num) right=middle-1;
else if(a[middle]<num) left=middle+1;
else
{
ans=middle;
left=middle+1;
}
}
2、例题
题目描述
输入 n n n 个不超过 1 0 9 10^9 109 的单调不减的(就是后面的数字不小于前面的数字)非负整数
a 1 , a 2 , … , a n a_1,a_2,\dots,a_{n} a1,a2,…,an,然后进行 m m m 次询问。对于每次询问,给出一个整数
q q q,要求输出这个数字在序列中第一次出现的编号,如果没有找到的话输出 − 1 -1 −1 。输入格式
第一行 2 2 2 个整数 n n n 和 m m m,表示数字个数和询问次数。
第二行 n n n 个整数,表示这些待查询的数字。
第三行 m m m 个整数,表示询问这些数字的编号,从 1 1 1 开始编号。
输出格式
输出一行, m m m 个整数,以空格隔开,表示答案。
样例 #1
样例输入 #1
11 3
1 3 3 3 5 7 9 11 13 15 15
1 3 6
样例输出 #1
1 2 -1
提示
数据保证, 1 ≤ n ≤ 1 0 6 1 \leq n \leq 10^6 1≤n≤106, 0 ≤ a i , q ≤ 1 0 9 0 \leq a_i,q \leq 10^9 0≤ai,q≤109, 1 ≤ m ≤ 1 0 5 1 \leq m \leq 10^5 1≤m≤105
本题输入输出量较大,请使用较快的 IO 方式。
代码
#include <bits/stdc++.h>
using namespace std;
int n,m;
int a[1000005];
int main()
{
cin>>n>>m;
for(int i=1;i<=n;++i)
scanf("%d",&a[i]);
int left=1;
int right=a[n]-a[1];
while(m--)
{
int num;
scanf("%d",&num);
int left=1;
int right=n;
int ans=-1;
while(left<=right)
{
int middle=(left+right)>>1;
if(a[middle]>num) right=middle-1;
else if(a[middle]<num) left=middle+1;
else
{
ans=middle;
right=middle-1;
}
}
cout<<ans<<" ";
}
return 0;
}
3、low_bound()和upper_bound()
lower_bound( )和upper_bound( )都是利用二分查找的方法在一个排好序的数组中进行查找的,算法复杂度为logn。
在从小到大的排序数组中,
lower_bound( begin,end,num)是从数组的begin位置到end-1位置二分查找第一个大于或等于num的数字,找到返回该数字的地址(想获得数字记得在前面加),不存在则返回end。通过返回的地址减去起始地址begin,得到找到数字在数组中的下标。在函数参数列表中添加greater< type >()则查找第一个小于或等于num的数字*,也可以用自定义的cmp函数。
upper_bound( begin,end,num)则是二分查找第一个大于num的数字,找到返回该数字的地址,不存在则返回end。通过返回的地址减去起始地址begin,得到找到数字在数组中的下标。同理,upper_bound()也可使用greater< type >()和cmp。
而在从在从大到小的排序数组中,
lower_bound( begin,end,num)是从数组的begin位置到end-1位置二分查找第一个小于或等于num的数字,同理,lower_bound()和upper_bound()的其他用法变化而推出,笔者不再过多赘述。