Jacobi.cpp

#include "Jacobi.hpp"

void Jacobi::fac2(uintmax_t n, uintmax_t &a, uintmax_t &b)
{
    uintmax_t ret_a = 0;
    uintmax_t ret_b = n;
    while((ret_b&1) == 0)
    {
        ret_b = ret_b >> 1;
        ret_a++;
    }
    a = ret_a;
    b = ret_b;
}

Jacobi.hpp

#include <iostream>

class Jacobi
{
public:
    static void fac2(uintmax_t n, uintmax_t &a, uintmax_t &b);
};

main.cpp

#include <iostream>
#include "Jacobi.hpp"

using namespace std;

int main(int argc, char **argv)
{
    uintmax_t ret_a;
    uintmax_t ret_b;
    uintmax_t n = 1872;
    Jacobi::fac2(n, ret_a, ret_b);
    cout << "2^" << ret_a << " * " << ret_b << endl;
    return(0);
}