src/integer-geometry/ntimes_centroid.hpp

View this file on GitHub
Last update: 2022-05-14 03:51:17+09:00
Include: #include "src/integer-geometry/ntimes_centroid.hpp"

Depends on

Required by

src/integer-geometry/unit-test/unit-test-ntimes_centroid.cpp

Code

#pragma once

#include <numeric>

#include "point.hpp"
#include "points.hpp"

namespace intgeometry2d {

  template< typename Z >
  lattice_point< Z > ntimes_centroid(const lattice_points<Z> &pts) {
    return std::accumulate(pts.begin(), pts.end(), lattice_point<Z>(0, 0));
  }

} // intgeometry2d

#line 2 "src/integer-geometry/ntimes_centroid.hpp"

#include <numeric>

#line 2 "src/integer-geometry/point.hpp"

#line 2 "src/integer-geometry/alias.hpp"

#include <cstddef>
#include <cstdint>

namespace intgeometry2d {
  using isize = std::ptrdiff_t;
  using usize = std::size_t;

  using i32 = std::int_fast32_t;
  using i64 = std::int_fast64_t;
  using u32 = std::uint_fast32_t;
  using u64 = std::uint_fast64_t;
} // intgeometry2d
#line 4 "src/integer-geometry/point.hpp"

namespace intgeometry2d {

  template< typename Z >
  class lattice_point {
    Z x_, y_;

  public:
    lattice_point() {}
    lattice_point(Z x_, Z y_) : x_(x_), y_(y_) {}
    
    Z x() const { return x_; }
    Z y() const { return y_; }

    bool operator==(const lattice_point &p) const { return x_ == p.x_ and y_ == p.y_; }
    bool operator!=(const lattice_point &p) const { return x_ != p.x_ or  y_ != p.y_; }

    lattice_point operator+(lattice_point p) { return lattice_point(x_ + p.x_, y_ + p.y_); }
    lattice_point operator-(lattice_point p) { return lattice_point(x_ - p.x_, y_ - p.y_); }

    Z norm() const { return x_ * x_ + y_ * y_; }
  };

} // intgeometry2d
#line 2 "src/integer-geometry/points.hpp"

#include <vector>

#line 6 "src/integer-geometry/points.hpp"

namespace intgeometry2d {

  template< typename Z >
  using lattice_points = std::vector< lattice_point<Z> >;

} // intgeometry2d
#line 7 "src/integer-geometry/ntimes_centroid.hpp"

namespace intgeometry2d {

  template< typename Z >
  lattice_point< Z > ntimes_centroid(const lattice_points<Z> &pts) {
    return std::accumulate(pts.begin(), pts.end(), lattice_point<Z>(0, 0));
  }

} // intgeometry2d