src/geometry/Z2/centroid/gravity.hpp

• View this file on GitHub
• Last update: 2023-05-23 10:28:02+09:00
• Include: #include "src/geometry/Z2/centroid/gravity.hpp"

Depends on

• 格子点 (整数幾何) (src/geometry/Z2/class/point.hpp)
• ベクトル (整数幾何) (src/geometry/Z2/class/vector.hpp)

Verified with

• test/aoj/2442.test.cpp

Code

#pragma once

#include "src/geometry/Z2/class/point.hpp"

#include <numeric>

namespace luz::Z2 {

  template < typename Z >
  Point< Z > ntimes_gravity(const Points< Z > &pts) {
    return std::accumulate(pts.begin(), pts.end(), Point< Z >());
  }

} // namespace luz::Z2

#line 2 "src/geometry/Z2/centroid/gravity.hpp"

#line 2 "src/geometry/Z2/class/point.hpp"

#line 2 "src/geometry/Z2/class/vector.hpp"

#include <vector>

namespace luz::Z2 {

  template < typename Z >
  class Vector {

    Z x_, y_;

   public:
    Vector(): x_(0), y_(0) {}
    Vector(Z x, Z y): x_(x), y_(y) {}

    Z x() const {
      return x_;
    }

    Z y() const {
      return y_;
    }

    bool operator==(const Vector &v) const {
      return x_ == v.x_ and y_ == v.y_;
    }

    bool operator!=(const Vector &v) const {
      return x_ != v.x_ or y_ != v.y_;
    }

    Vector &operator+=(const Vector &v) {
      x_ += v.x_;
      y_ += v.y_;
      return *this;
    }
    Vector &operator-=(const Vector &v) {
      x_ -= v.x_;
      y_ -= v.y_;
      return *this;
    }

    Vector operator+(const Vector &v) const {
      return Vector(*this) += v;
    }
    Vector operator-(const Vector &v) const {
      return Vector(*this) -= v;
    }

    Vector operator+() const {
      return *this;
    }
    Vector operator-() const {
      return Vector() - *this;
    }
  };

  template < typename Z >
  using Vectors = std::vector< Vector< Z > >;

} // namespace luz::Z2
#line 4 "src/geometry/Z2/class/point.hpp"

#line 6 "src/geometry/Z2/class/point.hpp"

namespace luz::Z2 {

  template < typename Z >
  using Point = Vector< Z >;

  template < typename Z >
  using Points = std::vector< Point< Z > >;

} // namespace luz::Z2
#line 4 "src/geometry/Z2/centroid/gravity.hpp"

#include <numeric>

namespace luz::Z2 {

  template < typename Z >
  Point< Z > ntimes_gravity(const Points< Z > &pts) {
    return std::accumulate(pts.begin(), pts.end(), Point< Z >());
  }

} // namespace luz::Z2

Back to top page