[一本通学习笔记] 二分与三分

最新推荐文章于 2022-04-26 20:48:32 发布
最新推荐文章于 2022-04-26 20:48:32 发布
阅读量464 收藏
点赞数
原文链接：http://www.cnblogs.com/mollnn/p/11583991.html
版权

摘要

本节主要关于二分与三分的问题，二分与三分时我们通常会区分整数和浮点而做不同的处理。

整数二分

最经典的一类，我们通常这样写

int L=... , R = ...;
while(R>L) {
    ...
    if(...) L=mid+1;
    else R=mid;
}
cout<<L (-1) <<endl;

需要注意的是在这里最终答案是否-1取决于L=mid+1这一切换进行时，mid是否为合法解。如果是，那么-1. 如果不是，就输出L。

当然还存在一些别的写法，比如将判断条件改为R>=L，切换语句写作L=mid+1,r=mid-1，同时记录最后一次合法的Ans。弄清楚一种即可。

 

浮点数二分

浮点类二分非常容易，只需要注意控制精度误差即可。

浮点数三分


三分法通常用于找凸函数极值（虽然理论上对导函数二分也挺好），这时候需要用到两个mid值，即m1和m2，分别对应区间[l,r]的两个三等分点，计算出两个三等分点，然后根据需要将左端点缩到左三等分点，或者右端点缩到右三等分点，这样来进行转移即可。在浮点数情况下实现是很容易的。

 

整数三分

有时候我们需要对整数三分，就需要注意一下取整误差的问题，提供参考写法如下。

int l=..., r=...;
while(r>l){
    int m1=(2*l+r)/3, m2=(2*r+l+2)/3;
    if(...) r=m2-1;
    else l=m1+1;
}
cout<<l<<endl;

 

二分答案+验证

如果问题的答案与决定性参变量的函数关系是非严格单调的，我们就可以用二分来将它转化为若干判定性问题，这在处理很多问题时是非常有效的。

 

二分套二分/三分套三分

当决定性参变量有多个，并且函数对每个参变量的偏导数都是恒非正或恒非负的时候，我们既可以用二分套二分来解决问题。即对f(x,y)，我们先二分一个x，将它转化为f(x,y)=g(y)，然后求出g(y)的答案，就得到了f(x,y)在这个x上的答案。

 

 

LOJ10011. 「一本通 1.2 例 1」愤怒的牛

#include <bits/stdc++.h>
using namespace std;
int n,m,a[1000005];
int main() {
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    sort(a+1,a+n+1);
    int l=0,r=1e+9;
    while(r>l) {
        int mid = (l+r)/2;
        int las = -1e+9, cnt = 0;
        for(int i=1;i<=n;i++)
            if(a[i]-las >= mid)
                ++cnt, las = a[i];
        if(cnt>=m) l=mid+1;
        else r=mid;
    }
    cout<<l-1<<endl;
}

LOJ10012. 「一本通 1.2 例 2」Best Cow Fences

这里涉及到二分答案的一个常用技巧，即对分数形式，比如这里的ans = (sum(...))/(sum(1))，我们二分一个ans，那么不等式ans > (sum(...))/(sum(1))就可以转化为sum(a[i]-ans)>0，通过后缀最大值处理一下就可以了。

浮点数问题一定要注意补偿舍入误差。

#include <bits/stdc++.h>
using namespace std;

int n,l,a[1000005];
double s[1000005];

int main() {
    cin>>n>>l;
    for(int i=1;i<=n;i++) cin>>a[i];
    double L=0, R = 1e+9;
    while(R-L > 1e-6) {
        double mid = (L+R) / 2.0;
        for(int i=1;i<=n;i++)
            s[i] = s[i-1] + a[i] - mid;
        int flag = 0;
        s[n+1] = -1e+19;
        for(int i=n-l;i>=0;--i)
            s[i+l] = max(s[i+l], s[i+l+1]),
            flag |= (s[i+l]-s[i]>-1e-6);
        if(flag) L=mid;
        else R=mid;
    }
    cout<<(int)(L*1000.0+0.005)<<endl;
}

 

#10013. 「一本通 1.2 例 3」曲线

#include <bits/stdc++.h>
using namespace std;

int T,n;
double a[10005],b[10005],c[10005],t1,t2,t3;

double f(double x) {
    double ret = -1e+18;
    for(int i=1;i<=n;i++)
        ret = max(ret, a[i]*x*x+b[i]*x+c[i]);
    return ret;
}

int main() {
    ios::sync_with_stdio(false);
    cin>>T;
    while(T--) {
        cin>>n;
        double ans = -1e+18;
        for(int i=1;i<=n;i++)
            cin>>a[i]>>b[i]>>c[i];
        double l=0,r=1000;
        while(r-l > 1e-11) {
            double m1 = l + (r-l)/3.0, m2 = r - (r-l)/3.0;
            if(f(m1)>f(m2)) l=m1;
            else r=m2;
        }

        cout<<setiosflags(ios::fixed)<<setprecision(4)<<f(l)<<endl;
    }
}

 

#10014. 「一本通 1.2 练习 1」数列分段 II

#include <bits/stdc++.h>
using namespace std;

int n,m,a[1000005];

int main() {
    cin>>n>>m;
    int mx=0;
    for(int i=1;i<=n;i++) cin>>a[i], mx=max(mx,a[i]);
    int l=mx,r=1e+9;
    while(r-l>0) {
        int mid = (l+r)/2;
        int sum = 0, cnt = 1;
        for(int i=1;i<=n;i++) {
            if(sum+a[i]>mid) ++cnt,sum=0;
            sum+=a[i];
        }
        if(cnt<=m) r=mid;
        else l=mid+1;
    }
    cout<<l<<endl;
}

#10015. 「一本通 1.2 练习 2」扩散

二分答案后考虑如何检验，暴力枚举点对看是否直接相连，并查集维护连通性即可。

#include <bits/stdc++.h>
using namespace std;

int n,x[100005],y[100005],f[105];

int find(int x) {
    return (f[x]==x)?x:(f[x]=find(f[x]));
}

int dis(int x1,int y1,int x2,int y2) {
    return abs(x1-x2) + abs(y1-y2);
}

int main() {
    cin>>n;
    for(int i=1;i<=n;i++) cin>>x[i]>>y[i];
    int l=0,r=1e+9;
    while(r-l) {
        int mid = (l+r)/2, cnt = n;
        for(int i=1;i<=n;i++) f[i]=i;
        for(int i=1;i<=n;i++)
            for(int j=1;j<i;j++)
                if(dis(x[i],y[i],x[j],y[j])<=2*mid)
                    if(find(i)!=find(j)) f[find(i)]=find(j), --cnt;
        if(cnt==1) r=mid;
        else l=mid+1;
    }
    cout<<l<<endl;
}

#10016. 「一本通 1.2 练习 3」灯泡

几何推公式

#include <bits/stdc++.h>
using namespace std;

int T; double H,h,D;

double f(double x) {
    double ret = h/(H-h)*x;
    if(ret+x > D)
        ret = D - x + (ret+x - D) * (H-h) / x;
    return ret;
}

int main() {
    cin>>T;
    while(T--) {
        cin>>H>>h>>D;
        double L=0,R=D;
        while(R-L > 1e-8) {
            double m1 = L + (R-L)/3, m2 = R - (R-L)/3;
            if(f(m1)>f(m2)) R=m2;
            else L=m1;
        }
        cout<<setiosflags(ios::fixed)<<setprecision(3)<<f(L)<<endl;
    }
}

#10017. 「一本通 1.2 练习 4」传送带

三分套三分，设两个比例k,u。

#include <bits/stdc++.h>
using namespace std;

double ax,ay,bx,by,cx,cy,dx,dy,p,q,r;

double f(double x1,double y1,double x2,double y2) {
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}

double calc(double k,double u) {
    double x1=k*ax+(1-k)*bx, y1=k*ay+(1-k)*by;
    double x2=u*cx+(1-u)*dx, y2=u*cy+(1-u)*dy;
    return f(ax,ay,x1,y1)/p+f(x1,y1,x2,y2)/r+f(x2,y2,dx,dy)/q;
}

double inner(double m) {
    double L=0,R=1;
    while(R-L > 1e-8) {
        double m1=L+(R-L)/3, m2=R-(R-L)/3;
        if(calc(m,m1)<calc(m,m2)) R=m2;
        else L=m1;
    }
    return calc(m,L);
}

int main() {
    cin>>ax>>ay>>bx>>by>>cx>>cy>>dx>>dy>>p>>q>>r;
    double L=0,R=1;
    while(R-L > 1e-8) {
        double m1=L+(R-L)/3, m2=R-(R-L)/3;
        if(inner(m1)<inner(m2)) R=m2;
        else L=m1;
    }
    cout<<setiosflags(ios::fixed)<<setprecision(2)<<inner(L)<<endl;
}

 

#include <bits/stdc++.h>using namespace std;
double ax,ay,bx,by,cx,cy,dx,dy,p,q,r;
double f(double x1,double y1,double x2,double y2) {    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));}
double calc(double k,double u) {    double x1=k*ax+(1-k)*bx, y1=k*ay+(1-k)*by;    double x2=u*cx+(1-u)*dx, y2=u*cy+(1-u)*dy;    return f(ax,ay,x1,y1)/p+f(x1,y1,x2,y2)/r+f(x2,y2,dx,dy)/q;}
double inner(double m) {    double L=0,R=1;    while(R-L > 1e-8) {        double m1=L+(R-L)/3, m2=R-(R-L)/3;        if(calc(m,m1)<calc(m,m2)) R=m2;        else L=m1;    }    return calc(m,L);}
int main() {    cin>>ax>>ay>>bx>>by>>cx>>cy>>dx>>dy>>p>>q>>r;    double L=0,R=1;    while(R-L > 1e-8) {        double m1=L+(R-L)/3, m2=R-(R-L)/3;        if(inner(m1)<inner(m2)) R=m2;        else L=m1;    }    cout<<setiosflags(ios::fixed)<<setprecision(2)<<inner(L)<<endl;}

 

转载于:https://www.cnblogs.com/mollnn/p/11583991.html

aodanchi3822
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
【51nod1281 山峰和棋子】 二分套二分 / 三分
我是谁 --- 卡比兽
07-15 188
51nod1281 一开始看到这题 想了下这不是可以直接三分么 因为随着k增大 能放的棋子先增大后减小 是一个单峰函数 如果是三分 就可以写二分套二分 这题主要是发现随着k增加 能放的棋子数先增加后减少 但是能放的棋子数 - k 一定是丛 0 逐渐变小 原来cal(x)-x 这个函数是单调递减的 那么我们就可以愉快的二分了 求答案里面有个二分 所以是二分套二分 滑稽 二分套二分 版本 /* ...
二分、三分入门
SZTU_ZJB的博客
01-15 430
二分、三分入门 二分、三分入门 小故事 1. 二分法简介 1.1 时间复杂度 1.2 工作原理 1.3 空间复杂度 1.3.1 递归版本 1.3.2 非递归版本 1.3.3 STL 的二分查找 2. 常见题型 2.1 二分查找 2.2 二分查找 2.2.1 [NOIP2015 提高组] 跳石头 3. 三分简介 小故事 一天DOGGOD_Q到图书馆借了 N 本书，出图书馆的时候，警报响了，于是保安把DOGGOD_Q拦下，要检查一下哪本书没有登记出借。DOGGOD_Q正准备把
[二分与三分]「一本通 1.2 练习 3」灯泡 「一本通 1.2 例 3」曲线]
wly_sh的博客
06-15 221
[一本通 1.2 练习 3」灯泡 简洁题意： 已知灯的高度 H，人的高度h，灯到墙距D，求人的地上的+墙上的影子L最长 题解： 在人远离墙的过程中，有 1.地上的影子先变长后变短，转折点为灯、人、墙角三点共线的时候 2.墙上的影子变短。灯、人、墙角三点共线之后为0 注意精读问题x #include<cstdio> double H,h,D,dd=1e-11; double check(double m) { double x=(D*h-m*H)/(H-h); double y=x*H/(D
一本通题解——1438：灯泡
努力中的老周的专栏
02-17 1123
题目链接 一本通OJ：http://ybt.ssoier.cn:8088/problem_show.php?pid=1438。 我的OJ：http://47.110.135.197/problem.php?id=4473。 题目 题目描述 相比 Wildleopard 的家，他的弟弟 Mildleopard 比较穷。他的房子是狭窄的，而且在他的房间里仅有一个灯泡。每天晚上，他徘徊在自己狭...
「一本通 1.2 练习 3」灯泡 --三分+三角函数
小元勋的博客
08-23 259
loj 10016 题目分析: 注意lll初值的设置 Code: #include <bits/stdc++.h> using namespace std; #define eps 1e-9 int T; double H,h,D; double F_(double x) { return D - x + H - D * ( ( H - h ) / x ); } v...
《算法图解》学习笔记一：二分查找
12-22
假设要在电话簿中找一个名字以K打头的人，（现在谁还用电话簿！）可以从头开始翻页，直到进入以K打头的部分。但你很可能不这样做，而是从中间开始，因为你知道以K打头的名字在电话簿中间。...二分查找是一种算
数据结构笔记：二分查找
12-22
对比序列的中间值，小于该值从左序列中二分查找，大于该值从右序列中二分查找，直到查找到目标值或子序列只有一个元素截止 实现 递归实现 def binarysearch(alist, item): '''二分查找_递归实现''' n = len(alist)...
运放与差分学习笔记.pdf
07-19
个人总结的一下运放和差分传输的知识！也有配合的仿真电路图，可以供大家参考！
计算机三级嵌入式学习笔记汇总.pdf
08-31
"计算机三级嵌入式学习笔记汇总.pdf" 本资源是一个关于计算机三级嵌入式学习笔记的汇总，涵盖了嵌入式系统概论、嵌入式系统开发的基础知识、数字媒体基础、网络通信技术等方面的知识点。 嵌入式系统概论 嵌入式...
单元测试学习笔记之三
01-31
我们就将学习两种单元测试的框架.第一种,就是鼎鼎大名的xUnit测试框架家族.xUnit测试框架有助于我们更加结构化地书写测试代码, 更方便地运行单元测试并检查运行结果.第二种框架称为隔离框架(isolationframewor
10个炫酷特效的网页写法（附源码），拿去就能用，nice
王春燕的博客
04-26 6万+
这是借鉴其他博主的，觉得挺实用的，给网页加个背景，给鼠标加个特效，这些“使用简单”，“效果爆炸”的页面，分享给大家，如果觉得有帮助可以点赞收藏支持一下，如果能关注一下就再好不过了。 目录 1、鼠标点击弹出爱心 2、鼠标点击弹出文字 3、鼠标点击弹出烟花波纹 4、鼠标小星星拖尾跟随 5、鼠标粒子随心拖尾跟随 6、鼠标笑脸跟随+仙女棒+泡泡+雪花+点击烟花效果（自由组合） 7、樱花特效 8、蜘蛛网特效 9、看板娘（左下角的小人） 10、烟花 1、鼠标点击弹出爱心 代码 &l
二分——数列分段 II
vina的博客
07-27 2830
问题 D: 数列分段 II 时间限制: 1 Sec  内存限制: 128 MB 提交: 11  解决: 7 [提交][状态][讨论版][命题人:add_cy] 问题描述：  对于给定的一个长度为N的正整数数列A[i]，现要将其分成M（M≤N）段，并要求每段连续，且每段和的最大值最小。 关于最大值最小： 例如一数列4 2 4 5 1要分成3段 将其如下分段： [4 2][4 5][1...
「一本通 1.2 例 3」曲线 --三分
小元勋的博客
08-22 262
Luogu 10013 题目分析： 本题是求最小值，三分时，令lmid=l+(r−l)/3,rmid=r−(r−l)/3,如果f(rmid)&gt;=f(lmid),就舍弃lmid右边那一截，令r=rmidl_{mid}=l+(r-l)/3,r_{mid}=r-(r-l)/3,如果f(r_{mid})&gt;=f(l_{mid}),就舍弃l_{mid}右边那一截，令r= r_{...
[贪心算法]「一本通 1.1 练习 2」数列分段 / [二分与三分]「一本通 1.2 练习 1」数列分段 II
wly_sh的博客
06-11 253
P1181 数列分段Section I luogu [一本通 1.1 练习 2」数列分段 简洁题意： 给n个数分组，每组和不超过m，求最小组数 题解： 从头开始枚，能塞进一组就塞 #include<cstdio> /*从前开始尽可能让每组分满 假设第1组有a，a+b<=m则化为一组，a+b>m则不化 .若a+b<=m , b+c<=m 但是 a+b+c>m ， 那当然选a,b一组,因为c+d<b+c+d */ int main() { int n,m
RGB、YUY2、YUYV、YVYU、UYVY、AYUV格式详解
zxh821112的专栏
01-29 1万+
以下摘自http://bbs.chinavideo.org/viewthread.php?tid=4143 　　I420是YUV格式的一种，而YUV有packed format和planar format两种，而I420属于planar format的一种。 　　同时I420表示了YUV的采样比例4:2:0。4:2:0的YUV并不是说没有V分量，而是指对于每一个行，只有一个U或者V分
#10015. 「一本通 1.2 练习 2」扩散
TJ
12-10 1491
目录    思路【floyd+dp】：  【代码实现：floyd+dp】    思路2:并查集+二分  【代码实现：并查集+二分】               最后来个小小的总结： 【题目描述】 一个点每过一个单位时间就会向 4 个方向扩散一个距离，如图所示：两个点 a 、b 连通，记作 e(a,b) 当且仅当 a 、b 的扩散区域有公共部分。连通块的定义是块内的任意两个点 u、v ...
给定平面上任意三个点的坐标(x​1​​,y​1​​)、(x​2​​,y​2​​)、(x​3​​,y​3​​)，检验它们能否构成三角形
qq_42475914的博客
11-12 2万+
1.题目来源  作者: 杨起帆 单位: 浙江大学 时间限制: 200 ms 内存限制: 64 MB 代码长度限制: 16 KB 2.题目简述 给定平面上任意三个点的坐标(x​1​​,y​1​​)、(x​2​​,y​2​​)、(x​3​​,y​3​​)，检验它们能否构成三角形。 输入格式: 输入在一行中顺序给出六个[−100,100]范围内的数字，即三个点的坐标x​1​​、y​1​...
第1部分 基础算法(提高篇)--第2章 二分与三分-1433：【例题1】愤怒的牛
zqhf123的博客
06-07 482
1433：【例题1】愤怒的牛 时间限制: 1000 ms 内存限制: 65536 KB 提交数: 2091 通过数: 1201 【题目描述】 农夫 John 建造了一座很长的畜栏，它包括N(2≤N≤100,000)个隔间，这些小隔间依次编号为x1,…,xN(0≤xi≤1,000,000,000). 但是，John的C(2≤C≤N)头牛们并不喜欢这种布局，而且几头牛放在一个隔间里，他们就要发生争斗。为了不让牛互相伤害。John决定自己给牛分配隔间，使任意两头牛之间的最小距离尽可能的大，
二分和三分总结与误区分析
litmxs的博客
10-08 692
今天学一个新算法的时候, 突然发现我以前的二分一直有点问题, 虽然以前的写的二分都能AC, 但还是有些问题, 如果题目给出的区间内没有解, 我的二分可能就会出问题. 所以干脆把所有二分三分思路整理一下. 加上更久之前发现自己三分写法的问题, 感觉自己能有些题目能AC真是运气好整数二分二分的坑点主要在整数的二分, 如果没有写好会出现死循环的问题.二分的原理是利用区间内值有序的特点, 不断让可行区间减半

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值