Разработать и испытать рекурсивную функцию для вычисления общего делителя двух введенных натуральных чисел

#include <iostream>
#include <string>
using namespace std;
// ru.wikipedia.org/wiki/Бинарный_алгоритм_вычисления_НОД
/*
НОД(0, n) = n; НОД(m, 0) = m; НОД(m, m) = m;
НОД(1, n) = 1; НОД(m, 1) = 1;
Если m, n чётные, то НОД(m, n) = 2*НОД(m/2, n/2);
Если m чётное, n нечётное, то НОД(m, n) = НОД(m/2, n);
Если n чётное, m нечётное, то НОД(m, n) = НОД(m, n/2);
Если m, n нечётные и n > m, то НОД(m, n) = НОД(m, (n-m)/2) = НОД((n-m)/2, m);
Если m, n нечётные и n < m, то НОД(m, n) = НОД((m-n)/2, n);
*/
long long CSD (long long m, long long n) {
  if (m == 0) return n;
  if (n == 0 || m == n) return m;
  if (m == 1 || n == 1) return 1;
  long long mEven = ! (m & 1);
  long long nEven = ! (n & 1);
  if (mEven && nEven)
    return 2 * CSD (m >> 1, n >> 1);
  if (mEven)
    return CSD (m >> 1, n);
  if (nEven)
    return CSD (m, n >> 1);
  if (m < n)
    return  CSD (m, (n - m) >> 1);
  return CSD ((m - n) >> 1, n);
}
int main (int argc, char** argv) {
  if (argc != 3) {
    cout << "Usage: CSD m n";
    return 0;
  }
  try {
    long long m = stoll (argv[1]);
    long long n = stoll (argv[2]);
    cout << CSD (abs (m), abs (n));
  } catch (exception err) {
    cout << err.what();
  }
}