#include <lemon/time_measure.h>

 #include <vector>

 template<class T>

 class Matrix 

 {

 public:

   int m,n;

   std::vector<T> values;

   Matrix(int _m, int _n) : m(_m), n(_n), values(n*m) {}

   typename std::vector<T>::reference operator()(int i, int j)

   { return values[i*n+j];}

   class Row 

   {

   public:

     Matrix<T> &mat;

     int row;

     Row(Matrix<T> &_mat,int _row) : mat(_mat), row(_row) {}

     typename std::vector<T>::reference operator()(int j)

     { return mat.values[row*mat.n+j];}

   };

   class Column 

   {

   public:

     Matrix<T> &mat;

     int col;

     Column(Matrix<T> &_mat,int _col) : mat(_mat), col(_col) {}

     typename std::vector<T>::reference operator[](int i)

     { return mat.values[i*mat.n+col];}

   };

   Row operator[](int i) { return Row(*this,i); }

   Column operator()(int j) { return Column(*this,j); }

 };

 main() 

 {

   const int COUNT=1000000;

   const int SIZE=100;

   Matrix<int> A(SIZE,SIZE);

   for(int i=0;i<SIZE;i++) for(int j=0;j<SIZE;j++) A(i,j)=i+j;

   int s=0;

   lemon::Timer T;

   for(int c=0;c<COUNT;c++) {

     for(int i=0;i<SIZE;i++) for(int j=0;j<SIZE;j++) s+=A[i](j);

   }

   std::cout << "A[i](j): " << T << ' ' << s <<'\n';

   s=0;

   T.reset();

   for(int c=0;c<COUNT;c++) {

     for(int i=0;i<SIZE;i++) for(int j=0;j<SIZE;j++) s+=A(i,j);

   }

   std::cout << "A(i,j): " << T << ' ' << s << '\n';

   s=0;

   T.reset();

   for(int c=0;c<COUNT;c++) {

     for(int i=0;i<SIZE;i++) for(int j=0;j<SIZE;j++) s+=A(j)[i];

   }

   std::cout << "A(j)[i]: " << T << ' ' << s << '\n';

   s=0;

   T.reset();

   for(int c=0;c<COUNT;c++) {

     for(int i=0;i<SIZE;i++) for(int j=0;j<SIZE;j++) s+=A(i,j);

   }

   std::cout << "A(j,i): " << T << ' ' << s << "\n\n";

   return 0;

 }