mooc北大算法课第四周：二分算法

ps：其实不算完全写完，倒数几题搞死我了（说的就是矩形分割），先交为敬，不然真的会死。。
程序算法的时间复杂度

1. 一个程序或算法的时间效率，也称时间复杂度，有时简称复杂度。
2. 复杂度常用大的字母0和小写字母n表示，比如0(n),0(n^2)等。n代表问题的规模。
3. 时间复杂度是用算法运行过程中，某种时间固定的操作需要被执行的次数和m的关系来度量的。在无序数列中查找某个数，复杂度是0(n)。
4. 计算复杂度时，止痛剂执行次数最多的(n足够大时)那种固定操作的次数。比如某个算法需要执行加法n^2次， 除法n次，那么就计复杂度是0(n^2)的。
5. 复杂度有“平均复杂度”和“最坏复杂度”两种。两者可能一致，也可能不一致。
6. 如果复杂度是多个n的函数之和，则只关心随n的增长增长的最快的那个函数。
   O(n^3 + n^2) => O(n^3)
   O(n^3 + 2^n) => O(2^n)
   O(n! + 3^n) => O(n!)
   常数复杂度：O(1) // 时间（操作次数）和问题的规模无关
   对数复杂度：O(log(n))
   线性复杂度：O(n)
   多项式复杂度：O(n^k)
   指数复杂度：O(a^n)
   阶乘复杂度：O(n!)
   在无需数列中查找某个数（顺序查找）：O(n)
   平面上有n个点，要求出任意两点之间的距离：O(n^2)
   插入排序，选择排序，冒泡排序：O(n^2)
   快速排序：O(n*log(n))
   二分查找：O(log(n))

例题1：插入法排序（从小到大）

代码：

void InsertionSort(int a[],int size)
{
   
	for(int i = 1; i < size; ++i) {
   
		//a[i]是最左的无序元素，每次循环将a[i]放到合适位置
		for(int j = 0; j < i; ++j)
		  if(a[j] > a[i]) {
   
			  //要把a[i]放到位置j，原下标j到i-1的元素都往后移一位，直接覆盖a[i]数字
			    int tmp = a[i];
				for(int k = i; k > j; --k)
				  a[k] = a[k-1];
				a[j] = tmp;
				break;
		  }
	}
}// 复杂度O(n^2)

总结

1. 理解问题
   1）由于比较，i起点为1，j起点为0
   2）调整位置时，采用倒序方法，值不会被覆盖
   3）注意最终a[j]的值为tmp，而不是a[i]

二分查找

例题2：二分查找函数

描述：
写一个函数BinarySearch，在包含size个元素的，从小到大排序的int数组a里查找元素p，如果找到，则返回元素下标，如果找不到，则返回-1。要求复杂度O(log(n))

代码：

int BinarySearch(int a[], int size, int p)
{
   
	int L = 0; //查找区间的左端点
	int R = size - 1; //查找区间的右端点
	while(L <= R) {
    //如果查找区间不为空就继续查找
		int mid = L + (R - L) / 2; //去查找区间正中元素的小标
        if(p == a[mid])
		  return mid;
		else if(p > a[mid])
		  L = mid + 1; //设置新的查找区间的左端点
		else
		  R = mid - 1; //设置新的查找区间的右端点
	}
	return -1;
}// 复杂度O(log(n))

总结：

1. 理解问题
   1）返回元素下标

例题3：二分查找函数

描述：
写一个函数LowerBound，在包含size个元素的，从小到大排序的int数组a里查找比给定整数p小的，下标最大的元素。如果找到，则返回元素下标，如果找不到，则返回-1。要求复杂度O(log(n))

代码：

int LowerBound(int a[], int size, int p) //复杂度O(log(n))
{
   
	int L = 0; //查找区间的左端点
	int R = size - 1;//查找区间的右端点
	int lastPos = -1;//到目前为止找到的最优解
	while(l <= R) {
    //如果查找区间不为空就继续查找
		int mid = L + (R - L) / 2; //取查找区间正中元素的下标
		if(a[mid] >= p)
		  R = mid - 1;
		else {
   
			lastPos = mid;
			L = mid + 1;
		}
	}
	if(a[lastPos] >= p)
	  return -1;
	else 
	  return lastPos;
}

总结：

例题4：二分法求方程的根
求下面方程的一个根，f(x) = x^3 - 5x^2 + 10x - 80 = 0,若求出的根是a，则要求|f(a)| <= 10^(-6)。

解题思路：
对f(x)求导，得f’(x) = 3x^2 - 10x + 10。由一元二次方程求根公式知方程f’(x) = 0无解，因此f’(x) 恒大于0。故f(x)单调递增。易知f(0) < 0且f(100) > 0，所以区间[0,100]内必然有且只有一个根。由于f(x)在[0,100]内是单调的，所以可以用二分的办法在区间[0,100]中寻找根。

代码：

#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
double EPS = 1e-6;
double f(double x) 
{
   
	return x*x*x - 5*x*x + 10*x - 80;
}

int main()
{
   
	double root, x1 = 0, x2 = 100, y;
	root = x1