pattern.cpp

#include <cstdlib>
#include <cmath>
#include <iostream>
#include <iomanip>
#include <string>
#include "DGtal/arithmetic/LighterSternBrocot.h"
#include "DGtal/arithmetic/Pattern.h"


using namespace DGtal;


void usage( int, char** argv )
{
  std::cerr << "Usage: " << argv[ 0 ] << " <p> <q> [SUB]" << std::endl;
  std::cerr << "\t - computes the pattern of slope p / q and displays it." << std::endl;
  std::cerr << "\t - the optional [SUB] parameter displays the pattern with its recursive Berstel decomposition." << std::endl;
}

int main( int argc, char** argv )
{
  if ( argc < 3 )
    {
      usage( argc, argv );
      return 1;
    }

  typedef DGtal::int32_t Integer;
  typedef DGtal::int32_t Quotient;
  typedef LighterSternBrocot<Integer, Quotient, StdMapRebinder> SB; // the type of the Stern-Brocot tree
  typedef SB::Fraction Fraction; // the type for fractions
  typedef Pattern<Fraction> MyPattern; // the type for patterns

  DGtal::int32_t p = atoi( argv[ 1 ] );
  DGtal::int32_t q = atoi( argv[ 2 ] );
  MyPattern pattern( p, q );

  bool sub = ( argc > 3 ) && ( std::string( argv[ 3 ] ) == "SUB" );
  std::cout << ( ! sub ? pattern.rE() : pattern.rEs( "(|)" ) ) << std::endl;
  return 0;
}