紫书第八章-----高效算法设计（二分查找、尺取法/滑动窗口法、三分查找）

二分查找

本文参考刘汝佳《算法竞赛入门经典》（第2版）

#include<iostream>
#include<cstdio>
using namespace std;

//二分查找（左闭右开区间）
int bin_search(int a[],int l,int r,int v){
    int mid;
    while(l<r){
        mid=l+(r-l)/2;
        if(v==a[mid]) return mid;
        else if(v>a[mid]) l=mid+1;
        else r=mid;
    }
    return -1;
}

//左闭右开区间，找最小的大于等于v的数在数组中的下标（越小越小）
int lower_bound(int a[],int l,int r,int v){
    int mid;
    while(l<r){
        mid=l+(r-l)/2;  //这样二分区间保证mid靠近区间起始l，不可能数组越界
        if(v<=a[mid]) r=mid;    //v小，一直把区间往左缩小，找到更小的且比v大的数
        else l=mid+1;   //v大，立马把区间右缩，要让数大于等于v，找到最小的大于等于v的数的数组下标
    }
    return l;
}
//左闭右开区间，找最小的大于v的数在数组中的下标（越小越小）
int upper_bound(int a[],int l,int r,int v){
    int mid;
    while(l<r){
        mid=l+(r-l)/2;//这样二分区间保证mid靠近区间起始l，不可能数组越界
        if(v>=a[mid]) l=mid+1;//v大，立马把区间右缩，要让数比v大，找到最小的大于v的数的数组下标
        else r=mid;//v小，一直把区间往左缩小，找到更小的且比v大的数
    }
    return l;
}


int main()
{
    int a[6]={1,2,4,5,6,7};
    cout<<"数组a: ";for(int i=0;i<6;i++) cout<<a[i]<<" ";cout<<endl;
    cout<<"4在数组a中的下标： ";cout<<bin_search(a,0,6,4)<<endl;
    cout<<"数组a中最小的大于等于4的数的下标： "; cout<<lower_bound(a,0,6,4)<<endl;
    cout<<"数组a中最小的大于4的数的下标是： ";cout<<upper_bound(a,0,6,4)<<endl;
    return 0;
}

尺取法/滑动窗口法

应用于有这么一类问题：需要在给的一组数据中找到不大于某一个上限的“最优连续子序列”

【例题 Subsequence POJ - 3061】

A sequence of N positive integers (10 < N < 100 000), each of them less than or equal 10000, and a positive integer S (S < 100 000 000) are given. Write a program to find the minimal length of the subsequence of consecutive elements of the sequence, the sum of which is greater than or equal to S.
Input
The first line is the number of test cases. For each test case the program has to read the numbers N and S, separated by an interval, from the first line. The numbers of the sequence are given in the second line of the test case, separated by intervals. The input will finish with the end of file.
Output
For each the case the program has to print the result on separate line of the output file.if no answer, print 0.
Sample Input
2
10 15
5 1 3 5 10 7 4 9 2 8
5 11
1 2 3 4 5
Sample Output
2
3

本题应用尺取法思路：

分析：首先，序列都是正数，如果一个区间其和大于等于S了，那么不需要在向后推进右端点了，因为其和也肯定大于等于S但长度更长，所以，当区间和小于S时右端点向右移动，和大于等于S时，左端点向右移动以进一步找到最短的区间，如果右端点移动到区间末尾其和还不大于等于S，结束区间的枚举。

针对样例1，尺取法的具体过程如下：

【代码】

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=100005;

long long a[maxn];


int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        long long sum=0,s;
        int ans=maxn;
        scanf("%d%I64d",&n,&s);
        for(int i=0;i<n;i++)
            scanf("%I64d",&a[i]);
        int st=0,en=0;
        while(1)
        {
            while(en<n && sum<s)sum+=a[en++];
            if(sum<s)break;
            ans=min(ans,en-st);
            sum-=a[st++];
        }
        if(ans==maxn)ans=0;
        printf("%d\n",ans);
    }
    return 0;
}

【例题 Unique Snowflakes UVA - 11572】

输入异构长度为n（n<=1000000）的序列，找到尽量长的连续子序列，使得该序列中没有相同的元素。

举个例子窗口滑动过程如下：

【代码】
参考《算法竞赛入门经典》（第2版） 刘汝佳著

#include<iostream>
#include<set>
#include<algorithm>
using namespace std;

const int maxn=1000000+5;
int a[maxn];

int main()
{
    int T,n;
    cin>>T;
    while(T--){
        cin>>n;
        for(int i=0;i<n;i++) cin>>a[i];

        set<int>s;
        int L=0,R=0,ans=0;
        while(R<n){
            while(R<n && !s.count(a[R])) s.insert(a[R++]);
            ans=max(ans,R-L);
            s.erase(a[L++]);
        }
        cout<<ans<<endl;
    }
    return 0;
}

三分查找法

三分法应用于单峰凸函数，下面参考如下博客地址：http://blog.csdn.net/pi9nc/article/details/9666627

如图所示，已知左右端点L、R，要求找到白点的位置。
思路：通过不断缩小 [L,R] 的范围，无限逼近白点。
做法：先取 [L,R] 的中点 mid，再取 [mid,R] 的中点 mmid，通过比较 f(mid) 与 f(mmid) 的大小来缩小范围。
当最后 L=R-1 时，再比较下这两个点的值，我们就找到了答案。
1、当 f(mid) > f(mmid) 的时候，我们可以断定 mmid 一定在白点的右边。
反证法：假设 mmid 在白点的左边，则 mid 也一定在白点的左边，又由 f(mid) > f(mmid) 可推出 mmid < mid，与已知矛盾，故假设不成立。
所以，此时可以将 R = mmid 来缩小范围。
2、当 f(mid) < f(mmid) 的时候，我们可以断定 mid 一定在白点的左边。
反证法：假设 mid 在白点的右边，则 mmid 也一定在白点的右边，又由 f(mid) < f(mmid) 可推出 mid > mmid，与已知矛盾，故假设不成立。
同理，此时可以将 L = mid 来缩小范围。

伪代码实现：

int SanFen(int l,int r) //找凸点
{
while(l < r-1)
{
int mid = (l+r)/2;
int mmid = (mid+r)/2;
if( f(mid) > f(mmid) )
r = mmid;
else
l = mid;
}
return f(l) > f(r) ? l : r;
}

【题目实战】

开始每个人都在一条数轴上的某个位置上，位置大于等于1，都是整数，每个人有个最大移动速度，在数轴上某个点集合，可以花费最短的时间，求出这个最短的时间。

代码：

#include<cstdio>
#include<algorithm>
#include<cmath>
#define INF 0x3f3f3f3f

using namespace std;
const double eps=0.000001;
const int maxn=60010;
double a[maxn],b[maxn];
int n;

double check(double x)
{
    double ans=0;
    for(int i=0;i<n;i++)
        ans=max(ans,fabs(x-a[i])/b[i]);
    return ans;
}

int main()
{
    double mmin=1.0*INF,mmax=-1.0*INF;
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%lf",&a[i]);
        mmin=min(mmin,a[i]);
        mmax=max(mmax,a[i]);
    }
    for(int i=0;i<n;i++)
        scanf("%lf",&b[i]);
    double l=mmin,r=mmax;
    while(l+eps<r)
    {
        double mid=(l+r)/2;
        double mmid=(mid+r)/2;
        if(check(mid)>check(mmid))l=mid;
        else r=mmid;
    }
    printf("%.12lf\n",check(r));
    return 0;
}