基于CLRS的合并排序C ++上的算法介绍(带反转计数)

时间:2019-07-10 04:40:54

标签: c++ arrays mergesort clrs inversion

我已经基于CLRS实现了对倒置进行计数的合并排序,但是答案不正确,没有对数组进行排序,也没有正确计算倒置。

我通过引用使用了相同的向量。

if(received_message_from_alice == decrypted_message)
  the signature proved the message is from alice

对于输入:{4,3,1,8,2}和5,正确的答案是6取反,而排序后的数组是{1,2,3,4,8}。它返回5个反演和{4,4,4,3,4}。

1 个答案:

答案 0 :(得分:1)

好几个月了,尽管我确实做了代码实现以返回排序后的数组的工作,但是在反转计数上仍然存在错误。今天,我正在研究并解决它,所以就在这里。

首先,在mergeInvCount方法的最后一个for中,使用基于1的索引访问arr,因此它不起作用,将其减去1即可访问基于0的索引。

第二,在比较两个要合并的辅助数组的条件下,必须计数一次反转的情况是不正确的。

当左辅助数组上的元素大于右辅助数组上的元素时,对于该数字,我们必须计算1个反转,之后的每个其他元素必须计算1个反转,但“无限” comodin除外。由于辅助数组是由于递归调用而排序的,因此是正确的。

它在左辅助数组从索引1开始时起作用,因为减法(mid-i)返回有序辅助数组中剩余的元素数量,而不考虑共生素。

但是当我们合并数组并且左边的数组不是从1开始时,减法无法在数组中的实际索引之后返回正确数量的元素。

因此,解决此问题的方法是使用n1,它在左辅助数组中存储元素数量。这样,(n1- i)返回正确的数字。

这是工作代码:

int mergeSortInvCount(vector<int> &arr, int izq, int der);
int mergeInvCount(vector<int> &arr, int izq, int mitad, int der);

void invCountRecursivo(vector<int> &arr, int n){

    int numInversiones = mergeSortInvCount(arr, 1, n);
    cout << "Num inversiones:" << numInversiones << endl;
    cout << "Ordered array, ascendant order" << endl;
    for(int i=0; i < n; i++){
        cout<<arr[i]<<endl;
    }
}

int mergeSortInvCount(vector<int> &arr, int izq, int der){

    int invCount = 0;

    if(izq < der){

        int mitad = (izq + der)/2;
        invCount = mergeSortInvCount(arr, izq, mitad);
        invCount += mergeSortInvCount(arr, mitad+1, der);
        invCount += mergeInvCount(arr, izq, mitad, der);
    }

    return invCount;
}

int infinito = numeric_limits<int>::max();

int mergeInvCount(vector<int> &arr, int izq, int mitad, int der){

    int n1 = mitad - izq + 1;
    int n2 = der - mitad;

    int invCount = 0;

    vector<int> vectorIzq;
    vector<int> vectorDer;

    for(int k=0;k<n1;k++){

        vectorIzq.push_back(arr[izq+k-1]);
    }

    vectorIzq.push_back(infinito);

    for(int k=0;k<n2;k++){

        vectorDer.push_back(arr[mitad+k]);
    }

    vectorDer.push_back(infinito);

    int i = 0;
    int j = 0;

    for(int k = izq; k <= der; k++){

        if(vectorIzq[i] <= vectorDer[j]){

            arr[k-1] = vectorIzq[i];
            i++;
        }
        else{

            arr[k-1] = vectorDer[j];
            j++;
            invCount += (n1 - i);
        }
    }

    return invCount;
}

int main(){

    vector<int> v = {4,3,1,8,2};
    invCountRecursivo(v, 5);
    // Returns 6, the correct # of inversions of A

    return 0;
}