Main.cpp

#include "PolynomialDivision.cpp"


int main()
{
  // Define ordering on monomials
  const char* ordering = "lex";
  const int n = 3;

  // Define Grobner basis
  std::vector<Polynomial<n>> G = { Polynomial<n>("x + y + z - 3", ordering), Polynomial<n>("x^2 + y^2 + z^2 - 5", ordering), Polynomial<n>("x^3 + y^3 + z^3 - 7", ordering) };
  ConvertToGrobnerBasis(G);

  std::cout << "Grobner Basis:" << std::endl;
  for (unsigned int i = 0; i < G.size(); ++i)
    G[i].printP();
  std::cout << std::endl;

  Polynomial<n> f = Polynomial<n>("x^4 + y^4 + z^4 - 9", ordering);
  MultivariatePolynomialDivision(f, G, false).printP();
  std::cout << std::endl;

  Polynomial<n> g1 = Polynomial<n>("x^5 + y^5 + z^5", ordering);
  MultivariatePolynomialDivision(g1, G, false).printP();

  Polynomial<n> g2 = Polynomial<n>("x^6 + y^6 + z^6", ordering);
  MultivariatePolynomialDivision(g2, G, false).printP();
}