《算法笔记》第4章 入门篇（2）---算法初步 4.5二分

4.5.1二分查找

二分查找概念：

二分查找实例：

注意二分查找的过程与序列的下标从0开始还是从1开始无关，参考代码如下：

#include<iostream>
using namespace std;
int erfen(int a[],int left,int right,int x)
{
    int mid;
    while(left<=right)
    {
        mid=(left+right)/2;
        if(a[mid]>x)
            right=mid-1;
        else if(a[mid]<x)
            left=mid+1;
        else if(a[mid]==x)
            return mid;
    }
    return -1;
}
int main()
{
    int a[10]={1,3,4,6,7,8,10,11,12,15};
    cout << erfen(a,0,9,9);
    return 0;
}

注意如果二分上界超过int型数范围的一半，那么当欲查询的元素在序列较靠后的位置时，应该：

如果递增序列中元素有可能重复应该：

如果：求序列中第一个大于等于x的元素的位置

int erfen(int a[],int left,int right,int x)
{
    int mid;
    while(left<right)  //注意这里是left<right，不用深究
    {
        mid=(left+right)/2;
        if(a[mid]>=x)
            right=mid;
        else if(a[mid]<x)
            left=mid+1;
    }
    return left;
}

如果：求序列中第一个大于x的元素的位置

int erfen(int a[],int left,int right,int x)
{
    int mid;
    while(left<right)  //注意这里是left<right，不用深究
    {
        mid=(left+right)/2;
        if(a[mid]>x)
            right=mid;
        else if(a[mid]<=x)
            left=mid+1;
    }
    return left;
}

---

4.5.2 二分法拓展

#include<iostream>
using namespace std;
const double esp=1e-5;
int main()
{
    double left=1.0,right=2.00,mid;
    while(right-left>esp)
    {
        mid=(left+right)/2;
        if(mid*mid>2)
            right=mid;
        else
            left=mid;
    }
    printf("%f",mid);
}

上面的问题相当于求下面关于x方程的问题，如果f(x)递增或者递减，对应代码是有区别的

装水问题：

#include<iostream>
#include<cmath>
using namespace std;
const double PI=acos(-1.0);
const double eps=1e-5;
double fun(double R,double H)
{
    double alpha=2*acos((R-H)/R);
    double L=2*sqrt(R*R-(R-H)*(R-H));
    double S2=R*R*PI/2;
    double S1=alpha*R*R/2-(R-H)*L/2;
    return S1/S2;
}
double solve(double R, double r)
{
    double left=0,right=R,mid;
    while(right-left>eps)
    {
        mid=(left+right)/2;
        if(fun(right,mid)>r)
        {
            right=mid;
        }
        else
            left=mid;
    }
    return mid;
}
int main()
{
    double R,r;
    cin >> R >> r;
    printf("%.4f",solve(R,r));
    return 0;
}

木棒切割问题：

木棒切割问题存疑:现在想不出来是对木棒长度进行二分还是对k进行二分，代码暂时写不出来

4.5.3 快速幂

用快速幂的方法：

typedef int LL;
LL wsf(LL a, LL b, LL m)
{
    if(b==0)
        return 1;
    if(b%2==1)   //这里完全可以用b&1，当b为奇数时，返回1；否则返回0
        return a*wsf(a,b-1,m)%m;
    else
    {
        LL kk=wsf(a,b/2,m);
        return kk*kk%m;
    }
}

对于上面的代码需要注意的问题：

用快速幂的迭代：

typedef long long LL;
LL wsf(LL a, LL b, LL m)
{
    LL ans=1;
    while(b>0)
    {
        if(b&1)
            ans=ans*a%m;
        a=a*a%m;
        b>>=1;
    }
    return ans;
}