#include <chrono>

int const Size = 1024+1/*5*/;
int const Half = Size / 2;

int arr[Size][Size] =
{
     1,  2,  3,  4,  5,
     6,  7,  8,  9, 11,
    12, 14, 15, 16, 17,
    18, 19, 20, 21, 22,
    23, 24, 25, 26, 27
};

int summ_diamond_1()
{
    int res = 0;

    for (int i=0; i<Half; ++i)
    {
        for (int j=Half-i; j<=Half+i; ++j)
        {
            res += arr[i][j];
            res += arr[Size-i-1][j];
        }
    }

    for (int i=0; i<Size; ++i)
    {
        res += arr[Half][i];
    }

    return res;
}

int summ_diamond_2()
{
  constexpr size_t dim = Size;
  constexpr size_t m = dim / 2;

  int sum = 0;
  int
      *head = ((int*)arr) + m,
      *tail = ((int*)arr) + dim*dim - m - 1;
  size_t
      n = 1,
      c = n,
      d;

  while (head != tail)
  {
    if (c--)
    {
      sum += *head++;
      sum += *tail--;
    }
    else
    {
      d = dim - 1 - n;
      head += d;
      tail -= d;
      n += 2;
      c = n;
    }
  }
  sum += *head;

  return sum;
}

double test(int (&summ_func)())
{
    auto start = std::chrono::high_resolution_clock::now();
    for (size_t i=0; i<100000; ++i)
    {
        int const res = summ_func();
    }
    auto end = std::chrono::high_resolution_clock::now();
    std::chrono::duration<double> elapsed_seconds = end - start;
    double const res = elapsed_seconds.count();
    std::cout << "Elapsed time: " << res << "s\n";
    return res;
}

void calc_times()
{
    for (size_t i=0; i<Size; ++i)
    {
        for (size_t j=0; j<Size; ++j)
        {
            arr[i][j] = std::rand();
        }
    }

    double const d1_1 = test(summ_diamond_1);
    double const d2_1 = test(summ_diamond_2);
    std::cout << "Relation: " << (d2_1 / d1_1) << "\n";

    double const d1_2 = test(summ_diamond_2);
    double const d2_2 = test(summ_diamond_1);
    std::cout << "Relation: " << (d2_2 / d1_2) << "\n";

    double const d1_3 = test(summ_diamond_1);
    double const d2_3 = test(summ_diamond_2);
    std::cout << "Relation: " << (d2_3 / d1_3) << "\n";
}