二分算法题解

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

#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;

double f(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    double l,r,m,x1,x2;
    int s=0,i;
    cin>>a>>b>>c>>d;
    for(i=-100;i<100;i++)
    {
        l=i;
        r=i+1;
        x1=f(l);
        x2=f(r);
        if(!x1)
        {
            cout<<fixed<<setprecision(2)<<l<<" ";
            s++;
        }//判断左端点，是零点直接输出。

        //不能判断右端点，会重复
        if(x1*x2<0)
        {
            while(r-l>=0.001)//二分控制精度
            {
                m=(l+r)/2;//中间值
                if(f(m)*f(r)<=0)
                    l=m;
                else
                    r=m;//计算中点处
            }
            cout<<fixed<<setprecision(2)<<r<<" ";
            s++;
        }
        if(s==3)
            break;
        //找到三个就退出
    }
    return 0;
}

P2678 [NOIP2015 提高组] 跳石头

最短距离尽可能长的问题
我们找最小距离就是二分的mid
那么显然mid的初始距离为(l+r)/2，然后的工作判断mid是不是最大的最短距离
二分循环结束的条件是移走的岩石数达到给定最大的数量
二分跳跃距离，以这个距离为标准移动石头
judge函数每个题有每个题的写法
等全部拿完再把拿走的数量和限制进行比对，如果超出限制，那么这就是一个非法解
反之就是一个合法解

#include <bits/stdc++.h>

using namespace std;

const int N=500010;
int d,n,m;
int a[N];
int l,r,mid,ans;
bool judge(int x)
{
    int tot=0;//tot代表计数器，记录以前答案需要移走的实际石头数
    int i=0;//代表下一块石头的编号
    int now =0;//now代表模拟挑石头的人当前在什么位置
    while(i<n+1) //千万注意不是n，n不是终点，n+1才是
    {
        i++;
        if(a[i]-a[now]<x) //判断距离，看两者之间差值
        {
            tot++;
        }//如果小把石头搬走
        else
        {
            now = i;//判断失败，这块石头不用搬走，直接跳过 考虑下一个
        }
    }
    if(tot>m)
            return false;
        else
            return true;
}

int main()
{
    cin>>d>>n>>m;
    for(int i=1; i<=n; i++)
    {
        cin>>a[i];
    }
    a[n+1]=d;//n不是终点
    l=1;
    r=d;
    while(l<=r)
    {
        mid=(l+r)/2;
        if(judge(mid))
        {
            ans=mid;
            l = mid + 1;//继续扩大
        }
        else
        {
            r=mid -1;
        }
    }
    cout<<ans<<endl;
    return 0;
}

P1102 A-B 数对

方法一：

已知n,c将a-b=c转化成a-c=b，首先将a数组每个元素出现的次数统计出来
用map映射，最后将数组a每次减去一个c，再将a数组扫描一遍，将所有的映射数加一起就是答案

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N=2e5+10;
ll a[N];
int n;
ll c,ans;
map<ll,ll> m;
//建立一个数字到出现次数的映射
//map<num,times>

int main()
{
    cin>>n>>c;
    for(int i=0;i<n;i++)
    {
        cin>>a[i];
        m[a[i]]++;
        //a[i]-=c;
    }
    for(int i=0;i<n;i++)
    {
        a[i]-=c;
        ans+=m[a[i]];
    }
    cout<<ans<<endl;
    return 0;
}

方法二：

哈希
实现：1.结构体存储x为这个位置的对应的数，y为这个数出现了几次
2.find(ll x) 找到哈希表h[]中x的位置
3.先找到此数在hash表中的位置，将其放入，并将这个位置对应的数量+1
4.check(ll x) 输出这个数在hash表中出现的次数；

#include <bits/stdc++.h>
#define ll long long
#define p 1000003//这个数越大越好最好是质数

using namespace std;
const int N=2e5+10;
int n;
ll m,ans;
ll a[p];
//哈希
struct node
{
    ll x;
    ll y;//x为这个位置的对应的数，y为这个数出现了几次
} h[p];
int find (ll x)
{
    //找到x位置
    int y=abs(x)%p;
    while(h[y].x&&h[y].x!=x)
    {
        y=(++y)%p;//不断地寻找y
    }
    return y;
}
void push(ll x)
{
    int y=find(x);
    h[y].y++;
    h[y].x=x;
    //先找到此数在hash表中的位置，并将这个位置对应的数量+1
}
int check(ll x)
{
    return h[find(x)].y;
    //输出这个数在hash表中出现的次数；
}
int main()
{
    cin>>n>>m;
    for(int i=0; i<n; i++)
    {
        cin>>a[i];
        push(a[i]);
        //输入并放入
    }
    for(int i=0; i<n; i++)
    {
        ans+=check(a[i]-m);//统计
    }
    cout<<ans<<endl;
    return 0;
}

方法三：

stl 二分函数
upper_bound(a+1,a+1+n,k)-lower_bound(a+1,a+1+n,k);

#include <bits/stdc++.h>
#define ll long long
//来个正宗的二分
//虽然我也不会

using namespace std;
const int N=2e5+10;
int n;
ll c,ans;
ll a[N];

int main()
{
    cin>>n>>c;
    for(int i=1; i<=n; i++)
    {
        cin>>a[i];
    }
    sort(a+1,a+1+n);//注意对这个数组进行排序
    for(int i=1; i<=n; i++)
    {
        int k=a[i]-c;
        ans+=upper_bound(a+1,a+1+n,k)-lower_bound(a+1,a+1+n,k);
        //大于这个元素的最小指针减去等于这个元素的最小值指针就是这个元素的个数
        //当然需要有序的情况下
        if(a[i]-c==a[i]) k--;
        //这里当注意的是数对的情况要求a-b=c；查找B出现的次数,当C=0的时候就会多算上a[i]的个数
        //所以要减去
    }
    cout<<ans<<endl;
    return 0;
}

方法四：典型二分模板

#include <bits/stdc++.h>
#define ll long long
//题解中说道是双指针的好题
//就是在说二分呗
//为了加深印象抄一下别人的解题思路
//我们考虑题目的要求求出所有a-b=c;的数对，首先将数组排序然后就会发现数组中的
//每一个a，对应的数都是一段连续的区间
//然后我们考虑如何去找这个区间
//我们显然是要找这个连续区间的左右端点。
//考虑到排序之后序列是有序的，我么枚举每个数，他们的左端点和右端点都是单调不降的，因此
//到这里我发现只是在模拟stl函数，并不是二分模板
//换一个
using namespace std;
const int N=2e5+10;
int n;
ll m,t;
ll l,r,mid,q,w;
ll a[N];
//来个正宗的二分查找
//略略略
//怎么找的不明白我服了这边界问题
int main()
{
    cin>>n>>m;
    for(int i=1; i<=n; i++)
    {
        cin>>a[i];
    }
    sort(a+1,a+1+n);//注意对这个数组进行排序
    for(int i=1; i<=n; i++)
    {
        l=i+1;//又是边界问题
        r=n;
        while(l<=r)//应该是闭区间
        {
            mid =(l+r)/2;
            if(a[i]+m<=a[mid])
                r=mid-1;
            else
                l=mid+1;
        }
        w=l;//找出后面连续数的第一位
        l=i+1;
        r=n;
        while(l<=r)
        {
            mid =(l+r)/2;
            if(a[i]+m<a[mid])
                r=mid -1;
            else
                l=mid +1;
        }
        q=l;//找出后面所制定连续数的末位
        t+=q-w;
    }
    cout<<t<<endl;
    return 0;
}

对于二分的边界问题，再理解一下整理下来