分治法求逆序数

# include<iostream>
# include<stdlib.h>
# include<ctime>
# include<cmath>
using namespace std;

int Arr2[100] = { 0 };

/*
    在对序列进行二路归并排序的时候,要将序列拆分成若干子序列,先将子序列排序,再合并子序列构成最终排序后的序列。
    二路归并算法还有一个特点,在进行归并操作时候的两个子序列是有序序列,所以,我们可以利用这一点,
    在归并子序列的时候,其中的子序列内部的逆序数必然是0,这时候能产生逆序数的情况必然处于子序列之间,
    即:“位置靠后的子序列”中的元素小于“位置靠前的子序列”的元素。例如,有子序列x:2,4,5;子序列y:1,3,6,
    显然,子序列y中的元素1的逆序数为3,子序列y中的元素3的逆序数为2,其他元素的逆序数均为0。
    通过这样一种方法,我们可以在序列的二路归并排序的过程中将序列的逆序数计算出来。故其时间复杂度为O(nlogn)。

*/

// 归并排序,并统计两个有序子序列中的下逆序对数
int merge(int Arr2[], int p, int q, int r)
{
    int count=0;                    
    int len1, len2;
    int arrL[100] = { 0 }, arrR[100] = { 0 };
    int i, j;

    len1 = q - p + 1;
    len2 = r - q;

    for (i = 1; i <= len1; i++)
    {
        arrL[i] = Arr2[p + i - 1];
    }

    for (j = 1; j <= len2; j++)
    {
        arrR[j] = Arr2[q + j];
    }

    // 设置观察哨
    arrL[i] = INT_MAX;
    arrR[j] = INT_MAX;

    i = 1;
    j = 1;

    for (int k = 0; k < r - p + 1; k++)
    {
        if (arrL[i] > arrR[j])
        {
            // 计算逆序对数
            count = count + len1 - i + 1;
            // 排序 生成新的有序数组
            Arr2[p + k] = arrR[j];
            j++;
        }
        else
        {
            Arr2[p + k] = arrL[i];
            i++;
        }
    }

    return count;
}

// 分治法 计算每个子序列的逆序对数,再合并,求母序列的逆序对数
int sort_Count(int Arr2[], int p, int r)
{
    if (p == r)
    {
        return 0;
    }
    int q = (p + r) / 2;
    int r1, r2, r_total;
    r1 = sort_Count(Arr2, p, q);
    r2 = sort_Count(Arr2, q + 1, r);
    r_total = r1 + r2 + merge(Arr2, p, q, r);
    return r_total;
}

int main_4()
{
    //产生随机数组
    srand((unsigned int)(time(NULL)));

    for (int i = 0; i < 10; i++)
    {
        Arr2[i] = rand() % (10 - 0) + 0;
    }

    cout << "排序前" << endl;
    for (int i = 0; i < 10; i++)
    {
        cout << Arr2[i] << " ";
    }
    cout << endl;

    int count = sort_Count(Arr2, 0, 9);

    cout << "排序后" << endl;
    for (int i = 0; i < 10; i++)
    {
        cout << Arr2[i] << " ";
    }
    cout << endl;

    cout << "逆序对的个数为:" << count << endl;

    system("pause");
    return 0;
}

分治法逆序对个数c的实现如下: 1. 将数组a分成两个子数组a1和a2,分别对a1和a2递归地逆序对个数c1和c2,然后合并a1和a2并计算跨越a1和a2的逆序对个数c3。 2. 合并a1和a2时,设置两个指针i和j分别指向a1和a2的起始位置,然后将较小的元素放入新数组b中,并将指针向后移动,直到其中一个子数组遍历完毕。此时,将另一个子数组的剩余元素全部放入b中。 3. 计算跨越a1和a2的逆序对个数c3时,设置两个指针i和j分别指向a1和a2的末尾位置,并将它们向前移动,同时设置一个计数器count用于计算跨越a1和a2的逆序对个数。比较a1[i]和a2[j]的大小,如果a1[i]>a2[j],则说明a2[j]和a1[i+1]到a1[n-1]都构成了逆序对,因为a1[i+1]到a1[n-1]都比a2[j]大,此时将count加上i+1的值,并将a1[i]放入新数组b中,然后将指针i向前移动;否则,将a2[j]放入新数组b中,并将指针j向前移动。 4. 返回c1+c2+c3作为逆序对个数。 代码实现如下: ``` int mergeSort(vector<int>& a, int left, int right) { if (left >= right) return 0; int mid = (left + right) / 2; int c1 = mergeSort(a, left, mid); int c2 = mergeSort(a, mid + 1, right); vector<int> b(right - left + 1); int i = left, j = mid + 1, k = 0, c3 = 0; while (i <= mid && j <= right) { if (a[i] <= a[j]) b[k++] = a[i++]; else { b[k++] = a[j++]; c3 += mid - i + 1; } } while (i <= mid) b[k++] = a[i++]; while (j <= right) b[k++] = a[j++]; for (int i = left; i <= right; i++) a[i] = b[i - left]; return c1 + c2 + c3; } int countInversePairs(vector<int>& nums) { return mergeSort(nums, 0, nums.size() - 1); } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值