#ifndef POINT_H
#define POINT_H
#include <iostream> //for class ostream
#include <cmath>    //for the sqrt function
using namespace std;

template <typename T>
class point; //forward declaration for a template class

template <typename T>
ostream& operator<<(ostream& ost, const point<T>& p);

template <typename T>
bool operator==(const point<T>& p1, const point<T>& p2);

template <typename T>
class point {
	T x;
	T y;
public:
	point<T>(T init_x, T init_y): x {init_x}, y {init_y} {}

	T r() const;   //distance from origin

	friend ostream& operator<< <T>(ostream& ost, const point<T>& p);

	friend bool operator==(const point<T>& p1, const point<T>& p2) {
		return p1.x == p2.x && p1.y == p2.y;
	}
};

template <typename T>
inline bool operator!=(const point<T>& p1, const point<T>& p2) {
	return !(p1 == p2);
}

template <typename T>
T point<T>::r() const {
	return sqrt(x*x + y*y);
}

template <typename T>
ostream& operator<<(ostream& ost, const point<T>& p) {
	return ost << "(" << p.x << ", " << p.y << ")";
}
#endif