分治算法

分治算法是将问题分解成几个子问题(一般是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;
}