barvinok/NTL_QQ.h

   1 #ifndef NTL_QQ_H
   2 #define NTL_QQ_H
   3
   4 #include <NTL/ZZ.h>
   5 #include <NTL/vector.h>
   6
   7 #ifdef NTL_STD_CXX
   8 using namespace NTL;
   9 #endif
  10
  11 struct QQ {
  12     ZZ  n;
  13     ZZ  d;
  14
  15     QQ() {}
  16     QQ(int n, int d) {
  17         this->n = n;
  18         this->d = d;
  19     }
  20
  21     QQ& operator += (const QQ& a) {
  22         ZZ g = GCD(d, a.d);
  23         ZZ num = a.n * (d / g) + (a.d / g) * n;
  24         ZZ den = a.d / g * d;
  25         g = GCD(num, den);
  26         n = num/g;
  27         d = den/g;
  28         return *this;
  29     }
  30
  31     QQ& operator *= (const QQ& a) {
  32         n *= a.n;
  33         d *= a.d;
  34         return *this;
  35     }
  36
  37     QQ& operator *= (const ZZ& a) {
  38         n *= a;
  39         return *this;
  40     }
  41 };
  42
  43 NTL_vector_decl(QQ,vec_QQ);
  44
  45 vec_QQ& operator *= (vec_QQ& a, const ZZ& b);
  46
  47 std::ostream& operator<< (std::ostream& os, const QQ& q);
  48 std::istream& operator>> (std::istream& os, QQ& q);
  49
  50 NTL_io_vector_decl(QQ,vec_QQ);
  51
  52 #endif