分治算法的基本思想是将一个规模为N的问题分解为K个规模较小的子问题,这些子问题相互独立且与原问题性质相同。求出子问题的解,就可得到原问题的解。即一种分目标完成程序算法,简单问题可用二分法完成。
分治算法的求解步骤 (1) 划分:把规模为 n 的原问题划分为 k 个规模较小的子问题,并尽量使这 k 个子问题的规模大致相同。 (2) 求解子问题:各子问题的解法与原问题的解法通常是相同的,可以用递归的方法求解各个子问题,有时递归处理也可以用循环来实现。 (3) 合并:把各个子问题的解合并起来,合并的代价因情况不同有很大差异,分治算法的有效性很大程度上依赖于合并的实现。 分治算法的适用条件 (1) 问题的规模缩小到一定的程度就可以容易地解决。为问题的计算复杂性一般是随着问题规模的增加而增加,因此大部分问题满足该特征。 (2) 问题可以分解为若干个规模较小的相同问题,即该问题具有最优子结构性质。这是应用分治法的前提,反映了递归思想的应用。 (3) 问题分解出的子问题的解可以合并为该问题的解。能否利用分治法完全取决于问题是否具有这条特征,如果具备了前两条特征,而不具备这条特征,则可以考虑贪心算法或动态规划。 (4) 问题所分解出的各个子问题是相互独立的,即子问题之间不包含公共的子问题。这条特征涉及到分治法的效率,如果各子问题不独立,则分治法要做许多重复工作,重复地求解公共子问题,此时虽然也可用分治法,但一般用动态规划较好。接下来以一个实验例子来体现分治算法(用分治法查找数组元素的最大值和最小值)
实验内容 用分治法查找数组元素的最大值和最小值 (1)问题描述 给定任意几组数据,利用分治法的思想,找出数组中的最大值和最小值并输出。 (2)实验步骤 ① 先解决小规模的问题,如数组中只有 1 个元素或者只有两个元素时候的情况。 ② 将问题分解,如果数组的元素大于等于 3 个,将数组分为两个小的数组。 ③ 递归的解各子问题,将 中分解的两个小的数组再进行以上两个步骤最后都化为小规模问题。 ④ 将各子问题的解进行比较最终得到原问题的解。
实验所用语言为python3,参考代码如下:
import time import numpy def find_bestvalue(array): if len(array) <= 2: Max = max(array) Min = min(array) return [Max,Min] L = find_bestvalue(array[: len(array)//2]) R = find_bestvalue(array[len(array)//2 :]) return merge(L,R) def merge(L,R): Max = (L[0] if L[0]>R[0] else R[0]) Min = (L[1] if L[1]<R[1] else R[1]) return [Max,Min] if __name__ == '__main__': print("请输入随机数组的范围和长度(格式:下界 上界 长度):") cin = input().split(' ') info = [] info.append(int(cin[0])) info.append(int(cin[1])) info.append(int(cin[2])) array = numpy.random.randint(info[0],info[1],info[2]) print('产生的随机数组为:'+str(array)) true_max = max(array) true_min = min(array) print('真实的最大最小值分别为:'+str(true_max)+' '+str(true_min)) start = time.clock() res = find_bestvalue(array) end = time.clock() print('分治法得到的最大最小值分别为:'+str(res[0])+' '+str(res[1])) print('分治算法的运行时间为:'+str(end-start)+'s')C++代码实现:
#include <iostream> #include <windows.h> #include <cmath> #include <bits/stdc++.h>//后面计算时间用得上 using namespace std; void func(int *a , int left , int right , int &Max , int &Min){ int max1=-1; int max2=-1; int min1=1; int min2=1; if(left==right){ Max=Min=a[left]; }else if(right-left==1){ Max=max(a[left],a[right]); Min=min(a[left],a[right]); }else{ int m=(left+right)/2; func(a,left,m,max1,min1); func(a,m+1,right,max2,min2); if(max1>max2){ Max=max1; }else{ Max=max2; } if(min1<min2){ Min=min1; }else{ Min=min2; } } } int main(int argc, char** argv){ int n; int Max=-1; int Min=1; cout<<"input the length of array:"<<endl; cin>>n; int *a=new int[n]; for(int i=0;i<n;i++){ cin>>a[i]; } clock_t start=clock(); func(a,0,n-1,Max,Min); clock_t end=clock(); cout<<"the max value of the array is :"<<Max<<endl; cout<<"the min value of the array is :"<<Min<<endl; cout<<"the time of the fuction has runned:"<<(double)(end-start)*1000000/CLOCKS_PER_SEC <<"μs"<<endl; delete []a; system("pause"); return 0; }