aboutsummaryrefslogtreecommitdiff
path: root/impls/qsort.cpp
blob: 22efbfe5e7fa4f86d8df99489beb1d48553bf123 (plain)
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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
/* qsort.cpp
 * Date: 2024-06-22
 * Author: Aryadev Chavali
 * Commentary: The inplace O(nlog(n)) sortitng algorithm.  Originally designed
 * in C, hence the lack of iostream.
 */
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <stack>
#include <vector>

using u64 = uint64_t;

#define ARR_SIZE(XS) (sizeof(XS) / sizeof((XS)[0]))
#define SWAP(A, B) \
  do               \
  {                \
    (A) ^= (B);    \
    (B) ^= (A);    \
    (A) ^= (B);    \
  } while (0)

void print_arr(u64 *arr, u64 size)
{
  printf("{\n");
  for (u64 i = 0; i < size; ++i)
    printf("  %lu,\n", arr[i]);
  printf("}\n");
}

struct Partition
{
  u64 *arr;
  u64 arr_size;
};

void quicksort(u64 *arr, u64 arr_size)
{
  std::stack<Partition> partitions;
  partitions.push({arr, arr_size});

  while (!partitions.empty())
  {
    const auto partition = partitions.top();
    arr                  = partition.arr;
    arr_size             = partition.arr_size;
    partitions.pop();

    // TODO: Use a different sort for small enough sizes?
    if (arr_size > 1)
    {
      // TODO: Choose a better pivot maybe?
      u64 pivot       = arr[arr_size - 1]; // Choose the last element as pivot
      u64 pivot_index = -1;

      for (u64 j = 0; j < arr_size; ++j)
      {
        if (arr[j] < pivot)
        {
          ++pivot_index;
          std::swap(arr[pivot_index], arr[j]);
        }
      }

      std::swap(arr[pivot_index + 1], arr[arr_size - 1]);

      partitions.push({arr, pivot_index + 1});
      partitions.push({arr + pivot_index + 1, arr_size - pivot_index - 1});
    }
  }
}

int main(int argc, char *argv[])
{
  (void)argc;
  (void)argv;

#define TEST_SIZE 10
  // randomly generate some numbers
  // u64 *arr     = new u64[TEST_SIZE];
  // u64 arr_size = TEST_SIZE;

  // for (size_t i = 0; i < TEST_SIZE; ++i)
  //   arr[i] = rand() % (2L << 30);

  // Use nums.txt
  std::ifstream input(argc > 1 ? argv[1] : "nums.txt");
  std::vector<u64> num_vec;

  for (std::string line; std::getline(input, line);
       num_vec.push_back(std::stoi(line)))
    continue;
  u64 *arr     = &num_vec[0];
  u64 arr_size = num_vec.size();

  printf("Before: ");
  print_arr(arr, arr_size);
  printf("\n");

  quicksort(arr, arr_size);

  printf("After: ");
  print_arr(arr, arr_size);
  printf("\n");

  // delete[] arr;

  return 0;
}