寒假训练报告2.6（二分）

1.从有序数列中查找某个值 O(log n)
算法描述：将值与区间的中点进行比较，不断把解的范围缩小到原来的一半

int Binary_search(int x){
    int l = 0, high = n, mid;
    while(l <= h){
        mid = (l + h) >> 1;
        if(a[mid] == x) return mid;
        else if(a[mid] < x) l = mid + 1;
        else    h = mid - 1;
    }
}

也可直接用STL函数lower_bound， upper_bound

ForwardIterator lower_bound(ForwardIterator first,ForwardIterator last,
const T &val);  //从小到大
ForwardIterator upper_bound(ForwardIterator first,ForwardIterator last,
const T &val);  //从大到小

2.假定一个解并判断是否可行

void Binary_search(){
    for(int i = 0; i < 100; i++){
        double mid = (l + h) >> 1;
        if(judge(mid))  l = mid;
        else    h = mid;
    }
    return l;
}

e.g:POJ 1064

注意：输出小数问题，通常会指定误差范围。1次循环可以把区间的范围缩小一半，100次的循环可以达到10^(-30)的精度范围，或者可以把终止条件设为(h - l) > EPS这样，但如果EPS过小，可能因为浮点小数精度的原因导致死循环。

3.最大化最小值
对最小值进行二分搜索，不断判断，直到确定范围。

bool judge(int d){
    int last = 0;
    for(int i = 1; i < M; i++){
        int crt = last + 1;
        while(crt < n && a[crt] - a[last] < d)
            crt++;
        if(crt == n)    return false;
        last = cnt;
    }
    return true;
}
double Binary_search(){
    sort(a, a +n);
    int l = 0, h = INF;
    while(h - l > 1){
        int mid = (l + h) / 2;
        if(judge(mid))  l = mid;
        else    h = mid;
    }
    return l;
}

e.g:POJ 2456

4.最大化平均值
问题描述：有n个物品的重量和价值分别是wi, vi,从中选出k个物品使得单位重量的价值最大。
错误解法：将物品按照单位价值进行排序，从大到小贪心选取。错误原因在于选取以后权重不同。

bool judge(double x){
    for(int i = 0; i < n; i++)
        y[i] = a[i] -  x * b[i];
    sort(y, y + n);
    double sum = 0;
    for(int i = 0; i < n - k; i++)
        sum += y[n - i - 1];
    return sum >= 0;
}
double Binary_search(){
    double l = 0, h = INF;
    for(int i = 0; i < 100; i++){
        double mid= (l + h) / 2;
        if(judge(mid))  l = mid;
        else    h = mid;
    }
    return l;
}

e.g:POJ 2976

总结：由上述几种题型可知，二分的关键在于找好二分条件，即Judge()函数，然后进行查找即可。