-
Notifications
You must be signed in to change notification settings - Fork 0
/
shell_sort.cpp
52 lines (44 loc) · 869 Bytes
/
shell_sort.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
#include <Rcpp.h>
using namespace Rcpp;
class Special{
public:
int nassign;
int ncomp;
bool Less(int a, int b){
ncomp++;
return a <= b;
}
void Assignment(){
nassign++;
}
Special(){
nassign = 0;
ncomp = 0;
}
};
// [[Rcpp::export]]
NumericVector ShellSort(NumericVector arr, NumericVector Gaps) {
int gap, j, i;
int n = arr.size();
Special special;
for(int t = Gaps.size()-1; t >= 0; t--) {
gap = Gaps[t];
for(j = gap; j<n; j++) {
int key = arr[j];
for(i = j-gap; i>=0; i -= gap) {
if(special.Less(arr[i], key))
break;
else{
arr[i+gap] = arr[i];
special.Assignment();
}
}
arr[i+gap] = key;
}
}
NumericVector res(2);
res[0] = special.nassign;
res[1] = special.ncomp;
//Rcout << arr << "\n";
return res;
}