Странная ошибка в слиянии

Я только что запрограммировал следующую простую функцию слияния в mergesort, которая точно соответствует той, что написана в книге CLRS.

#include<iostream>
#include<vector>
#include<limits>

using namespace std;

//Define the numerical positive infinity
const int INFP = numeric_limits<int>::max();

void MERGE(vector<int> A, int p, int q, int r){
//Create the size of left and right array
int n1 = q-p+1;
int n2 = r-q;
//Create the left array
vector<int> L(n1+1);
for(int i=0; i<n1;i++){
L[i] = A[i];
}
L[n1] = INFP; //Insert sentinel here!
//Create the right array
vector<int> R(n2+1);
for(int i=0; i<n2;i++){
R[i] = A[q+i+1];
}
L[n2] = INFP; //Insert sentinel here!
int i = 0;
int j = 0;
for(int k = 0; k <= r; k++){
if(L[i]<=R[j]){
A[k] = L[i];
i=i+1;
}
else{
A[k] = R[j];
j=j+1;
}
}
for(int m=0;m<4;m++){
cout<< A[m] << " ";
}
cout << endl;
}

int main(){
//test for merge function:
vector<int> A(4);
A[0]=1;
A[1]=3;
A[2]=2;
A[3]=4;
MERGE(A,0,1,3);
for(int m=0;m<4;m++){
cout<< A[m] << " ";
}
cout << endl;
return 0;
}

Тем не менее, он дал мне следующую распечатку, которая действительно смутила меня:

1 2 3 4
1 3 2 4

Я не знаю, является ли проблема функцией void, что я не могу использовать функцию void для векторов или что-то еще.

Очень надеюсь, что кто-нибудь может мне помочь. Спасибо!