分治算法是将问题分解成几个子问题(一般是2个),之后对这些子问题进行操作,或者部分问题操作,然后再处理完成后的这一个或几个 部分的结果,实现整个任务的完成。
生活中的实例
16硬币,有可能有1枚假币,假币比真币轻。有一架天平,用最少称量次数确定有没有假币,若有的话,假币是哪一枚。
我们分部处理8和8。之后选取较轻的那一堆分4和4。之后同理最后分成1和1,确定最终的结果。
分治的典型应用:归并排序
递归将数组分为2半,之后对两边依次排序,之后两两结合返回到原来的数组中完成排序。
归并排序代码
#include<iostream>
using namespace std;
void merge1(int *a,int lo, int mid ,int hi, int *b)
{ int pb=0;
int p1=lo;
int p2=mid+1;
while(p1<=mid&&p2<=hi){
if(a[p1]>a[p2])
b[pb++]=a[p2++];
else
b[pb++]=a[p1++];
}
while(p1<=mid)
b[pb++]=a[p1++];
while(p2<=hi)
b[pb++]=a[p2++];
for(int i=0;i<hi-lo+1;i++)
a[lo+i]=b[i];
}
void mergesort(int *a,int lo,int hi,int *b)
{
if(lo<hi){
int m=lo+(hi-lo)/2;
mergesort(a,lo,m,b);
mergesort(a,m+1,hi,b);
merge1(a,lo,m,hi,b);
}
}
int main()
{
int a[10]={13,27,19,2,8,12,2,8,30,89};
int b[10];
mergesort(a,0,9,b);
for(int i=0;i<10;i++)
printf("%d ",a[i]);
cout<<endl;
return 0;
}
时间复杂度
分治的典型应用:快速排序
设k=a[0], 将k挪到适当位置,使得比k小的元素都在k左边,比k大的元素都在k右边,和k相等的,不关心在k左右出现均可 (O(n)时间完成)。之后把k左边的部分快速排序,再把k右边的部分快速排序。
快速排序代码
#include<iostream>
using namespace std;
void quicksort(int *a,int lo,int hi)
{
if(lo>=hi)
return ;
int i=lo;
int j=hi;
int k=a[lo];
while(i!=j){
while(i<j&&k<=a[j])
j--;
swap(a[i],a[j]);
while(i<j&&k>=a[i])
i++;
swap(a[i],a[j]);
}
quicksort(a,lo,i);
quicksort(a,i+1,hi);
}
int main()
{
int a[]={93,27,30,2,8,12,2,8,30,89};
quicksort(a,0,9);
for(int i=0;i<10;i++)
printf("%d ",a[i]);
cout<<endl;
}
时间复杂度
例题一 输出前m大的数
描述 给定一个数组包含n个元素,统计前m大的数并且把这m个数从大到小输出。
输入 第一行包含一个整数n,表示数组的大小。n < 100000。第二行包含n个整数,表示数组的元素,整数之间以一个空格分开,。每个整数的绝对值不超过100000000。 第三行包含一个整数m。m < n。
方法一 排序后再输出,复杂度 O(nlogn)。
方法二 用分治处理:复杂度 O(n+mlogm)
方法二 设key=a[0],进行操作是比key大的都在key右边,比key小的都在左边。之后在选择数组的前部或后部进行操作。
代码
#include<iostream>
using namespace std;
int s[100000];
int n,m;
void ac(int *s,int lo,int hi, int key)
{
int k=s[lo];
int i=lo,j=hi;
while(i!=j){
while(i<j&&s[j]>=k)
j--;
swap(s[i],s[j]);
while(i<j&&s[i]<=k)
i++;
swap(s[i],s[j]);
}
if(hi-i+1==key)
return ;
else if(hi-i+1>key)
ac(s,i+1,hi,key);
else
ac(s,lo,i-1,key-(hi-i+1));
}
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
scanf("%d",&s[i]);
ac(s,0,n-1,m);
for(int i=n-m;i<n;i++)
printf("%d ",s[i]);
cout<<endl;
return 0;
}
例题二 求排列的逆序数
考虑1,2,…,n (n <= 100000)的排列i1,i2,…,in,如果其中存在j,k,满足j < k 且 ij > ik, 那么就称(ij,ik)是这个排列的一个逆序。一个排列含有逆序的个数称为这个排列的逆序数。例如排列 263451 含有8个逆序(2,1),(6,3),(6,4),(6,5),(6,1),(3,1),(4,1),(5,1),因此该排列的逆序数就是8。现给定1,2,…,n的一个排列,求它的逆序数。
笨办法
分治
先将数组分成两半,分别求出左半边的逆序数和右半边的逆序数,再算有多少逆序是由左半边取一个数和右半边取一个数构成。
就相当于于归并算法,不过在每一步加上了统计逆序的算法,而且这个归并是由大到小排序的。
代码
#include<iostream>
#include<cstdio>
using namespace std;
int sum=0;
void merge1(int *a,int lo,int mid,int hi,int *b)
{
int i=lo,j=mid+1;
int k=0;
while(i<=mid&&j<=hi){
if(a[i]<=a[j])
b[k++]=a[j++];
else{
b[k++]=a[i++];
sum+=hi-j+1; // 逆序数统计
}
}
while(i<=mid)
b[k++]=a[i++];
while(j<=hi)
b[k++]=a[j++];
for(int i=0;i<hi-lo+1;i++)
a[lo+i]=b[i];
}
void mergesort(int *a,int lo,int hi,int *b)
{
if(lo<hi){
int mid=lo+(hi-lo)/2;
mergesort(a,lo,mid,b);
mergesort(a,mid+1,hi,b);
merge1(a,lo,mid,hi,b);
}
}
int main()
{
int a[100005],b[100005];
int n;
cin>>n;
for(int i=0;i<n;i++)
scanf("%d",&a[i]);
mergesort(a,0,n-1,b);
cout<<sum<<endl;
return 0;
}