洛谷题单算法1-6二分查找与二分答案

Cloud_Ivy37

已于 2024-03-24 11:39:35 修改

阅读量665

点赞数 16

分类专栏：洛谷题单文章标签：算法 c++

于 2024-03-24 09:00:00 首次发布

本文链接：https://blog.csdn.net/2301_76144982/article/details/136965149

版权

洛谷题单专栏收录该内容

6 篇文章 0 订阅

订阅专栏

本文介绍了多个C++编程实例，涉及二分查找算法、数组排序、一元三次方程求解、区间划分、设备充电问题等，展示了在实际编程中处理数学问题的方法和技术

摘要由CSDN通过智能技术生成

代码仅供参考，不足之处望理解。

P2249 【深基13.例1】查找

#include<iostream>
using namespace std;
const int N=1e6+5;
int n,m,a[N];
int main(){
    ios::sync_with_stdio(false);
    //输入
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    for(int i=1;i<=m;i++){
        int t=0;
        cin>>t;
        //整数的二分
        int left=1,right=n;
        while(left<right){
            int mid=(left+right)/2;
            if(a[mid]>=t) right=mid;
            else left=mid+1;
        }
        //输出
        if(a[left]==t) cout<<left<<" ";
        else cout<<"-1 ";
    }
    return 0;
}

P1102 A-B 数对

#include<iostream>
#include<algorithm>
using namespace std;
const int N=2e5+5;
int n,c,a[N];
long long ans;
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>c;
    //输入后sort排序
    for(int i=1;i<=n;i++) cin>>a[i]; sort(a+1,a+1+n);
    //双指针思想:i为快指针,j,k均为慢指针
    for(int i=1,j=1,k=1;i<=n;i++){
        while(j<=i&&a[i]-a[j]>=c){
            j++;
        }
        while(k<=i&&a[i]-a[k]>c){
            k++;
        }
        ans+=j-k;
    }
    cout<<ans<<endl;
    return 0;
}

P1024 [NOIP2001 提高组] 一元三次方程求解

#include<iostream>
#include<cmath>
using namespace std;
const double eps=1e-8;
double a,b,c,d;
//开口向上
void Dsearch1(double l,double r){
    while(r-l>eps){
        double mid=(l+r)/2;
        if(a*pow(mid,3)+b*pow(mid,2)+c*mid+d>=0) r=mid;
        else l=mid;
    }
    printf("%.2lf ",l);
}
//开口向下
void Dsearch2(double l,double r){
    while(r-l>eps){
        double mid=(l+r)/2;
        if(a*pow(mid,3)+b*pow(mid,2)+c*mid+d<=0) r=mid;
        else l=mid;
    }
    printf("%.2lf ",l);
}
int main(){
    ios::sync_with_stdio(false);
    cin>>a>>b>>c>>d;
    //将三元方程求导后,令二元方程等于0,求二元方程的解:t1,t2
    double temp=4.0*b*b-12.0*a*c;
    double t=sqrt(temp);
    double t1=(-2.0*b-t)/6.0/a;
    double t2=(-2.0*b+t)/6.0/a;
    //始终保持 t1 < t2
    if(t1>t2) {double tt=t1;t1=t2;t2=tt;}
    double flag=100.0;
    //根据二元函数开口方向即a的正负,让三个区间分别进行浮点数的二分
    if(a>0){Dsearch1(-flag,t1);Dsearch2(t1,t2);Dsearch1(t2,flag);}
    else {Dsearch2(-flag,t1);Dsearch1(t1,t2);Dsearch2(t2,flag);}
    return 0;
}

P1678 烦恼的高考志愿

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e6+5;
int n,m,a[N];
long long ans;
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>m;
    //输入完直接sort排序
    for(int i=1;i<=n;i++) cin>>a[i]; sort(a+1,a+1+n);
    //依次输入分数
    while(m--){
        int score=0;
        cin>>score;
        //整数的二分
        int l=1,r=n;
        while(l<r){
            int mid=(l+r)/2;
            if(a[mid]<score) l=mid+1;
            else r=mid;
        }
        //二分结束后应该是这样的数据状态: a[l-1] < score < a[l]
        int x=abs(score-a[l]);
        int y=abs(score-a[l-1]);
        //注意l==1的情况,此时a[l-1]应该省掉
        if(l==1) y=1000000; //赋值一个较大的数即可
        //累加ans
        ans+=(long long)min(x,y);
    }
    //输出
    cout<<ans<<endl;
    return 0;
}

P1163 银行贷款

#include<iostream>
#include<cmath>
using namespace std;
const double eps=1e-8;
int w0,w,m;
bool check(double x){
    double sum=0.0;
    //求和
    for(int i=1;i<=m;i++) sum+=w/pow((1.0+x/100),i);
    //判断
    if(sum<=w0) return true;
    else return false;
}
//浮点数的二分
void Dsearch(double l,double r){
    while(r-l>eps){
        double mid=(l+r)/2;
        if(check(mid)) r=mid;
        else l=mid;
    }
    //保留一位小数
    printf("%.1lf",l);
}
int main(){
    ios::sync_with_stdio(false);
    cin>>w0>>w>>m;
    //left有可能为0.0
    double left=0.0,right=300.0;
    Dsearch(left,right);
    return 0;
}

P1873 [COCI 2011/2012 #5] EKO / 砍树

#include<iostream>
using namespace std;
const int N=1e6+5;
int n,m,a[N];
bool check(int x){
    long long temp=0;
    for(int i=1;i<=n;i++){
        if(a[i]>x) temp+=a[i]-x;
        //再for循环里面就可以判断,节省时间
        if(temp>=m) return true;
    }
    return false;
}
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    int l=1,r=1e6;
    while(l<r){
        int mid=(r+l+1)>>1;
        if(check(mid)) l=mid;
        else r=mid-1;
    }
    cout<<l<<endl;
    return 0;
}

P2440 木材加工

#include<iostream>
using namespace std;
const int N=1e5+5;
int n,k,a[N];
bool check(int x){
    int temp=0;
    for(int i=1;i<=n;i++){
        temp+=a[i]/x;
        if(temp>=k) return true;
    }
    return false;
}
int main(){
    ios::sync_with_stdio(false);
    //输入
    cin>>n>>k;
    for(int i=1;i<=n;i++) cin>>a[i];
    //整数的二分
    int l=0,r=1e8;
    while(l<r){
        int mid=(l+r+1)/2;
        if(check(mid)) l=mid;
        else r=mid-1;
    }
    //输出
    cout<<l<<endl;
    return 0;
}

P2678 [NOIP2015 提高组] 跳石头

#include<iostream>
#include<algorithm>
using namespace std;
const int N=5e4+5;
int l,n,m,a[N];
bool check(int x){
    int last=0,cnt=0;
    for(int i=1;i<=n+1;i++){
        if(a[i]-last<x) cnt++; //凡是距离起点的值小于x,计数自增
        else last=a[i];        //否则更新起点值
    }
    //性质的划分应该是cnt与m的大小关系,而cnt<=m时,表示当前解可行,应该去右边查找更大的值
    if(cnt<=m) return true;
    else return false;
}
int main(){
    ios::sync_with_stdio(false);
    //输入
    cin>>l>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    a[n+1]=l;            //加上终点
    sort(a+1,a+1+n);     //排序
    //整数的二分
    int left=0,right=1e8;
    while(left<right){
        int mid=(left+right+1)/2;
        if(check(mid)) left=mid;
        else right=mid-1;
    }
    //输出
    cout<<left<<endl;
    return 0;
}

P3853 [TJOI2007] 路标设置

#include<iostream>
using namespace std;
const int N=1e5+5;
int l,n,k,a[N];
bool check(int x){
    //cnt负责计数
    long long cnt=0;
    for(int i=1;i<n;i++){
        //相邻路标中每x长度放一个路标,计算它们之间的路标个数
        long long temp=a[i+1]-a[i];
        if(temp%x==0) cnt+=temp/x-1;
        else cnt+=temp/x;
    }
    //性质的划分应该是cnt与k的大小关系,而cnt<=k时,表示当前解可行,应该去左边查找更小的值
    if(cnt<=k) return true;
    else return false;
}
int main(){
    ios::sync_with_stdio(false);
    //输入
    cin>>l>>n>>k;
    for(int i=1;i<=n;i++) cin>>a[i];
    //整数的二分
    int left=1,right=l;
    while(left<right){
        int mid=(left+right)/2;
        if(check(mid)) right=mid;
        else left=mid+1;
    }
    //输出
    cout<<left<<endl;
    return 0;
}

P1182 数列分段 Section II

#include<iostream>
using namespace std;
const int N=1e5+5;
int n,m,a[N];
bool check(int x){
    int cnt=0,temp=0;
    for(int i=1;i<=n;i++){
        //如果划分区域和大于x,cnt自增,更新temp
        if(temp+a[i]>x) cnt++,temp=a[i];
        //否则temp继续加值
        else temp+=a[i];
    }
    //划分次数会比组数少一,因此判断时cnt加一
    if(cnt+1<=m) return true;
    else return false;
}
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>m;
    int l=0,r=1e9;
    //左边界应该是数组中最大值,右边界应该是整个数组的和
    for(int i=1;i<=n;i++){
        cin>>a[i];
        l=max(l,a[i]);
        r+=a[i];
    }
    //整数的二分
    while(l<r){
        int mid=(l+r)>>1;
        if(check(mid)) r=mid;
        else l=mid+1;
    }
    //输出
    cout<<l<<endl;
    return 0;
}

P3743 小鸟的设备

#include<iostream>
using namespace std;
const int N=1e5+5;
int n;
double p,a[N],b[N];
bool check(double x){
    double maxx=p*x;//充电器最多提供的电量
    double sum=0.0; //需要的电量总和
    for(int i=1;i<=n;i++){
        //若已有的电量大于使用时间需要的电量则忽略
        if(a[i]*x<=b[i]) continue;
        //否则充电，并加上需要的电量。
        else sum+=(a[i]*x-b[i]);
    }
    //性质的划分应该是sum与maxx的大小关系,当sum<=maxx时,应该去右边查找更大的值
    if(sum<=maxx) return true;
    else return false;
}
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>p;
    //temp记录所有设备的耗电速度
    double temp=0.0;
    for(int i=1;i<=n;i++){
        cin>>a[i]>>b[i];
        temp+=a[i];
    }
    //若所有设备的耗点速度和小于充电速度，输出-1
    if(temp<=p){cout<<"-1"<<endl; return 0;}
    //浮点数的二分
    double l=0.0,r=5e9;
    while(r-l>1e-6){
        double mid=(l+r)/2;
        if(check(mid)) l=mid;
        else r=mid;
    }
    //输出
    printf("%.10lf",l);
    return 0;
}