# include <cstdlib>
# include <iomanip>
# include <iostream>

using namespace std;

# include "polyomino_index.hpp"

int main ( );

//****************************************************************************80

int main ( )

//****************************************************************************80
//
//  Purpose:
//
//    POLYOMINO_INDEX_TEST tests POLYOMINO_INDEX.
//
//  Licensing:
//
//    This code is distributed under the GNU LGPL license.
//
//  Modified:
//
//    30 April 2018
//
//  Author:
//
//    John Burkardt
//
{
  int i;
  int j;
  int m = 3;
  int n = 4;
//
//  P is listed in column-major order;
//
  int p[3*4] = { 
    1, 1, 0,
    0, 1, 1,
    1, 1, 1,
    1, 0, 0 };
  int *pin;

  timestamp ( );
  cout << "\n";
  cout << "POLYOMINO_INDEX_TEST\n";
  cout << "  C++ version\n";
  cout << "  POLYOMINO_INDEX assigns an index to each nonzero entry\n";
  cout << "  of a polyomino.\n";

  polyomino_print ( m, n, p, "  The polyomino P:" );

  pin = polyomino_index ( m, n, p );

  cout << "\n";
  cout << "  PIN: Index vector for P:\n";
  cout << "\n";
  for ( i = 0; i < m; i++ )
  {
    for ( j = 0; j < n; j++ )
    {
      cout << " " << pin[i+j*m];
    }
    cout << "\n";
  }
//
//  Free memory.
//
  delete [] pin;
//
//  Terminate.
//
  cout << "\n";
  cout << "POLYOMINO_INDEX_TEST\n";
  cout << "  Normal end of execution.\n";
  cout << "\n";
  timestamp ( );

  return 0;
}