#include <iostream>
#include <cmath>

using namespace std;

class Point{
  double x, y, z;
  
public:
  Point(double _x, double _y, double _z){
    x=_x;
    y=_y;
    z=_z;
  }
  
  ~Point(){
    cout<<"Destructor"<<endl;
  }
  
  void PrintCoordinatesXYZ(){
    cout << "( " << x << ", " << y << ", " << z <<" )"<<endl;
  }
  
  double GetR(){
    return ( sqrt( x*x + y*y + z*z) );
  }
  
  double GetPhi(){
    if(x==0)
      return 0.0;
    return (atan(y/x));
  }
  
  double GetTheta(){
    if(GetR()==0)
      return 0.0;
    return (asin(z/GetR()));
  }
  
  void PrintCoordinatesRPhiTheta(){
    cout << "( " << GetR() << ", " << GetPhi() << ", " << GetTheta() <<" )"<<endl;
  }
  
  double GetDistance(Point &B){
    return ( sqrt( pow((z-B.z),2) + pow((y-B.y),2) + pow((x-B.x),2) ) ) ;
  }

};

int main(){
  Point A(1, 1, 1), B(0, 0, 0);
  
  A.PrintCoordinatesRPhiTheta();
  B.PrintCoordinatesRPhiTheta();
  A.PrintCoordinatesXYZ();
  
  cout<<A.GetDistance(B)<<endl;
  

  return 0;
}