Получить все делители числа q, взаимно простые к p

// факторизация
#include <vector>
#include <cmath>
 
std::vector
factorize(int num)
{
    std::vector factors;
    if (!num) throw 1; // что-нибудь другое можно сюда впихнуть
    num = num > 0 ? num : -num;
    int boundary = std::pow(num 1.0 / 3.0); // ходят слухи, что проверять нужно лишь до третьего корня
    for (int i = 2; i < boundary; ++i) { // тупо в лоб перебираем
        if (num % i == 0) {
            factors.push_back(i);
        }
    }
    return factors;
}
 
unsigned
GCD(unsigned u, unsigned v) { // НОД
    while ( v != 0) {
        unsigned r = u % v;
        u = v;
        v = r;
    }
    return u;
}
 
// в main сравниваешь каждое число из вектора с p в GCD, если выхлопнет 1, значит оно

#include <iostream>
#include <vector>
#include <algorithm>
#include <iterator>
 
std::vector<unsigned> factorize_unique(unsigned n)
{ // Факторизует n; возвращает только уникальные факторы
  std::vector<unsigned> factors;
 
  for (unsigned factor = 2; factor * factor <= n; )
    if (n % factor == 0)
    {
      if (factors.empty() || factors.back() != factor)
        factors.push_back(factor);
      n /= factor;
    }
    else
      ++factor;
 
  if (n > 1 && (factors.empty() || factors.back() != n))
    factors.push_back(n);
 
  return factors;
}
 
int main()
{
  unsigned q = 210;
  unsigned p = 12;
 
  std::vector<unsigned> q_factors = factorize_unique(q);
  std::vector<unsigned> p_factors = factorize_unique(p);
 
  std::vector<unsigned> q_only;
  std::set_difference(
    q_factors.begin(), q_factors.end(), 
    p_factors.begin(), p_factors.end(),
    std::back_inserter(q_only));
 
  if (q_only.empty())
  {
    std::cout << "Нет делителей" << std::endl;
    return 0;
  }
 
  for (unsigned pattern = 1, end_pattern = 1 << q_only.size();
       pattern != end_pattern;
       ++pattern)
  {
    unsigned divisor = 1;
 
    std::vector<unsigned>::iterator it;
    unsigned bit;
    for (it = q_only.begin(), bit = 1; it != q_only.end(); ++it, bit <<= 1)
      if ((pattern & bit) != 0)
        divisor *= *it;
 
    std::cout << divisor << " ";
  }
 
  std::cout << std::endl;
}