算法笔记 二分法介绍与应用

算法笔记 二分法介绍与应用

二分法是解决以下一些问题的重要算法，如：判断函数与x轴交点的具体位置（也可以说是方程组的解）、查找、判断无理数的大小等。这一算法的思想其实很简单，就是每次都进行折半寻找，直到找到想要的解。下面用几个具体的例子来对这一方法的应用进行说明。
头文件均为
一、二分查找
这个理解比较简单，看代码就好。

  int binarySearch(int arr[],int left,int right,int findNum) {
    //此处的arr[]应该是排好序的，以从小到大排序。由大到小改写顺序的话也比较简单，此处不进行赘述。
    	while (left <= right) {
    		int mid = (left + right) / 2;
    		if (findNum == arr[mid]) return mid;
    		if (arr[mid]<findNum) {
    			left = mid + 1;
    		}
    		else {
    			right = mid-1;
    		}
    	}
    	return -1;
    }

二、计算无理数的值

#include<iostream>
const double eps = 1e-9;//精度
double fun(double x) {
	return x*x;
}
double calValue(double a,double b) {//这个数左边最接近数为a,求sqrt(b)
	double left = a, right = b,mid;
	while (right - left > eps) {
		mid = (left + right) / 2;
		if (fun(mid) > b) {
			right = mid;
		}
		else {
			left = mid;
		}
	}
	return mid;
}
int main() {
	printf("%.9lf",calValue(1, 2));
	system("pause");
}

三、函数求根

//一元函数求根（此处为y=sqrt(3)x-sqrt(2))  实验结果仅相差最后一位数字
/*
输出结果：0.816496581
	     0.816496582
*/
double fun2(double x) {
	return sqrt(3.0)*x - sqrt(2.0);
}
double calRoot(double left, double right) {
	double mid;
	while (fabs(right) - fabs(left) > eps) {
		mid = (left + right) / 2.0;
		if (fun2(mid)> 0.0) {
			right = mid;
		}
		else {
			left = mid;
		}

	}
	return mid;
}
int main() {
	printf("%.9lf\n", sqrt(6)/3.0);
	int i;
	for (i = -2;; i++) {//给一个范围，可以大一点，循环找就好
		if (i > (sqrt(2) / sqrt(3)))
			break;
	}
	double left = (double)(i - 1);
	double right = (double)i;
	printf("%.9lf\n", calRoot(left, right));
	system("pause");
}

四、木棒切割问题
给3根木棒，长度为10， 24， 15，要切割得到7根长度相等的木棒， 则7根木棒的长度最长为6，其组合为16+46+2*6。
我的解决办法是将长度先进行排序，再进行二分解决。

//木棒切割问题
int stick(int arr[],int len,int cuts) {//数组，数组长度，切分段数
	int minLen=0;//当前可切分的最小长度
	int maxLen=arr[0];//当前可分的最短木棒的长度(可视为能进行切割的最大长度）
	int mid = maxLen;
	while (maxLen-minLen>0) {
		int cur_cuts = 0;//当前可将木棒切分的段数
		for (int i = 0; i < len; i++) {
			cur_cuts += arr[i] / mid;
		}
		if (cur_cuts == cuts) {
			return mid;
		}
		else if(cur_cuts>cuts){
			minLen = mid;
		}
		else {
			maxLen = mid;
		}
		mid = (minLen + maxLen) / 2;
	}
	return -1;
}
int main() {
//在题目成立的条件下我自己找了一组数据，算了一下应该是满足条件的
	int arr[] = { 15,16,24,35,51 };
	printf("%d\n",stick(arr, 5,10));
	system("pause");
}
//如果以上几个代码有错误，欢迎大家指出来