Problem: Spinning tops are one of the most popular and the most traditional toys. Not only spinning them, but also making one’s own is a popular enjoyment.

One of the easiest way to make a top is to cut out a certain shape from a cardboard and pierce an axis stick through its center of mass. Professionally made tops usually have three dimensional shapes, but in this problem we consider only two dimensional ones.

Usually, tops have rotationally symmetric shapes, such as a circle, a rectangle (with 2-fold rotational symmetry) or a regular triangle (with 3-fold symmetry). Although such symmetries are useful in determining their centers of mass, they are not definitely required; an asymmetric top also spins quite well if its axis is properly pierced at the center of mass.

When a shape of a top is given as a path to cut it out from a cardboard of uniform thickness, your task is to find its center of mass to make it spin well. Also, you have to determine whether the center of mass is on the part of the cardboard cut out. If not, you cannot pierce the axis stick, of course.

Hints

An important nature of mass centers is that, when an object O can be decomposed into parts O1 , . . . , On with masses M1 , . . . , Mn , the center of mass of O can be computed by:

<image>

where Gk is the vector pointing the center of mass of Ok.

A circular segment with its radius r and angle θ (in radian) has its arc length s = rθ and its chord length c = r√(2 - 2cosθ). Its area size is A = r2(θ - sinθ)/2 and its center of mass G is y = 2r3sin3(θ/2)/(3A) distant from the circle center.

<image>

Figure 2: Circular segment and its center of mass

<image>

Figure 3: The first sample top



Input

The input consists of multiple datasets, each of which describes a counterclockwise path on a cardboard to cut out a top. A path is indicated by a sequence of command lines, each of which specifies a line segment or an arc.

In the description of commands below, the current position is the position to start the next cut, if any. After executing the cut specified by a command, the current position is moved to the end position of the cut made.

The commands given are one of those listed below. The command name starts from the first column of a line and the command and its arguments are separated by a space. All the command arguments are integers.


start x y


Specifies the start position of a path. This command itself does not specify any cutting; it only sets the current position to be (x, y).


line x y


Specifies a linear cut along a straight line from the current position to the position (x, y), which is not identical to the current position.


arc x y r


Specifies a round cut along a circular arc. The arc starts from the current position and ends at (x, y), which is not identical to the current position. The arc has a radius of |r|. When r is negative, the center of the circle is to the left side of the direction of this round cut; when it is positive, it is to the right side (Figure 1). The absolute value of r is greater than the half distance of the two ends of the arc. Among two arcs connecting the start and the end positions with the specified radius, the arc specified is one with its central angle less than 180 degrees.


close


Closes a path by making a linear cut to the initial start position and terminates a dataset. If the current position is already at the start position, this command simply indicates the end of a dataset.

The figure below gives an example of a command sequence and its corresponding path. Note that, in this case, the given radius -r is negative and thus the center of the arc is to the left of the arc. The arc command should be interpreted as shown in this figure and, not the other way around on the same circle.

<image>

Figure 1: A partial command sequence and the path specified so far

A dataset starts with a start command and ends with a close command.

The end of the input is specified by a line with a command end.

There are at most 100 commands in a dataset and at most 100 datasets are in the input. Absolute values of all the coordinates and radii are less than or equal to 100.

You may assume that the path does not cross nor touch itself. You may also assume that paths will never expand beyond edges of the cardboard, or, in other words, the cardboard is virtually infinitely large.

Output

For each of the dataset, output a line containing x- and y-coordinates of the center of mass of the top cut out by the path specified, and then a character ‘+’ or ‘-’ indicating whether this center is on the top or not, respectively. Two coordinates should be in decimal fractions. There should be a space between two coordinates and between the y-coordinate and the character ‘+’ or ‘-’. No other characters should be output. The coordinates may have errors less than 10-3 . You may assume that the center of mass is at least 10-3 distant from the path.

Example

Input

start 0 0
arc 2 2 -2
line 2 5
arc 0 3 -2
close
start -1 1
line 2 1
line 2 2
line -2 2
arc -3 1 -1
line -3 -2
arc -2 -3 -1
line 2 -3
line 2 -2
line -1 -2
line -1 -1
arc -1 0 2
close
start 0 0
line 3 0
line 5 -1
arc 4 -2 -1
line 6 -2
line 6 1
line 7 3
arc 8 2 -1
line 8 4
line 5 4
line 3 5
arc 4 6 -1
line 2 6
line 2 3
line 1 1
arc 0 2 -1
close
end


Output

1.00000 2.50000 +
-1.01522 -0.50000 -
4.00000 2.00000 +

Solution (in CPP):

Solution:
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
using namespace std;
const double EPS = 1e-10;
const double INF = 1e+10;
const double PI = acos(-1);
int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
inline double ABS(double a){return max(a,-a);}
struct Pt {
	double x, y;
	Pt() {}
	Pt(double x, double y) : x(x), y(y) {}
	Pt operator+(const Pt &a) const { return Pt(x + a.x, y + a.y); }
	Pt operator-(const Pt &a) const { return Pt(x - a.x, y - a.y); }
	Pt operator*(const Pt &a) const { return Pt(x * a.x - y * a.y, x * a.y + y * a.x); }
	Pt operator-() const { return Pt(-x, -y); }
	Pt operator*(const double &k) const { return Pt(x * k, y * k); }
	Pt operator/(const double &k) const { return Pt(x / k, y / k); }
	double ABS() const { return sqrt(x * x + y * y); }
	double abs2() const { return x * x + y * y; }
	double arg() const { return atan2(y, x); }
	double dot(const Pt &a) const { return x * a.x + y * a.y; }
	double det(const Pt &a) const { return x * a.y - y * a.x; }
};
double tri(const Pt &a, const Pt &b, const Pt &c) { return (b - a).det(c - a); }
int iSP(Pt a, Pt b, Pt c) {
	int s = sig((b - a).det(c - a));
	if (s) return s;
	if (sig((b - a).dot(c - a)) < 0) return -2; // c-a-b
	if (sig((a - b).dot(c - b)) < 0) return +2; // a-b-c
	return 0;
}
int iLL(Pt a, Pt b, Pt c, Pt d) {
	if (sig((b - a).det(d - c))) return 1; // intersect
	if (sig((b - a).det(c - a))) return 0; // parallel
	return -1; // correspond
}
double dSP(Pt a, Pt b, Pt c) {
	if (sig((b - a).dot(c - a)) <= 0) return (c - a).ABS();
	if (sig((a - b).dot(c - b)) <= 0) return (c - b).ABS();
	return ABS(tri(a, b, c)) / (b - a).ABS();
}

bool iLS(Pt a, Pt b, Pt c, Pt d) {
	return (sig(tri(a, b, c)) * sig(tri(a, b, d)) <= 0);
}
bool iSS(Pt a, Pt b, Pt c, Pt d) {
	return (iSP(a, b, c) * iSP(a, b, d) <= 0 && iSP(c, d, a) * iSP(c, d, b) <= 0);
}
bool iSSstrict(Pt a, Pt b, Pt c, Pt d) {
	return (sig(tri(a, b, c)) * sig(tri(a, b, d)) < 0 && sig(tri(c, d, a)) * sig(tri(c, d, b)) < 0);
}
Pt pLL(Pt a, Pt b, Pt c, Pt d) {
	b = b - a; d = d - c; return a + b * (c - a).det(d) / b.det(d);
}
bool iCS(Pt a, double r, Pt b, Pt c) {
	return (sig(r - dSP(b, c, a)) >= 0 && sig(r - max((b - a).ABS(), (c - a).ABS())) <= 0);
}
pair<Pt,Pt> pCC(Pt a, double r, Pt b, double s) {
	double d = (b - a).ABS();
	double x = (d * d + r * r - s * s) / (d * 2);
	Pt e = (b - a) / d, w = e * Pt(0, 1) * sqrt(max(r * r - x * x, 0.0));
	return make_pair(a + e * x - w, a + e * x + w);
}
pair<Pt,Pt> pCL(Pt a, double r, Pt b, Pt c) {
	Pt h = b + (c - b) * (c - b).dot(a - b) / (c - b).abs2(); // perp(b, c, a)
	double d = (h - a).ABS();
	double y = sqrt(max(r * r - d * d, 0.0));
	Pt e = (c - b) / (c - b).ABS();
	return make_pair(h - e * y, h + e * y);
}

int sGP(int n, Pt p[], Pt a) {
	int side = -1, i;
	p[n] = p[0];
	for (i = 0; i < n; ++i) {
		Pt p0 = p[i] - a, p1 = p[i + 1] - a;
		if (sig(p0.det(p1)) == 0 && sig(p0.dot(p1)) <= 0) return 0;
		if (p0.y > p1.y) swap(p0, p1);
		if (sig(p0.y) <= 0 && 0 < sig(p1.y) && sig(p0.det(p1)) > 0) side = -side;
	}
	return side;
}
int sAP(Pt a, Pt b, Pt c) {
	return sig(a.det(c)) - sig(b.det(c)) - sig(a.det(b));
}
char in[10];
Pt p[210];
int isl[210];
int main(){
	while(1){
		int sz=0;
		scanf("%s",in);
		if(in[0]=='e')break;
		for(int i=0;i<210;i++)isl[i]=0;
		double X,Y;
		scanf("%lf%lf",&X,&Y);
		p[sz++]=Pt(X,Y);
		vector<pair<Pt,int> > left;
		vector<pair<Pt,int> > right;
		while(1){
			scanf("%s",in);
			if(in[0]=='c'){
				if((p[0]-p[sz-1]).ABS()<EPS){sz--;break;}
				isl[sz-1]=1;
				p[sz]=p[0];break;
			}
			if(in[0]=='l'){
				scanf("%lf%lf",&X,&Y);
				isl[sz-1]=1;
				p[sz++]=Pt(X,Y);
			}
			if(in[0]=='a'){
				double r;
				scanf("%lf%lf%lf",&X,&Y,&r);
				p[sz++]=Pt(X,Y);
				Pt M=(p[sz-1]+p[sz-2])/2;
				double dist=(p[sz-1]-p[sz-2]).ABS()/2;
				double h=sqrt(r*r-dist*dist);
				if(r<0){
					Pt v=(p[sz-1]-p[sz-2])*Pt(0,1);
					v=v/v.ABS()*h;
					left.push_back(make_pair(M+v,sz-2));
				}else{
					Pt v=(p[sz-1]-p[sz-2])*Pt(0,-1);
					v=v/v.ABS()*h;
					right.push_back(make_pair(M+v,sz-2));
				}
			}
		}
		double S=0;
		Pt G=Pt(0,0);
		for(int i=0;i<sz;i++){
			double det=p[i].det(p[i+1])/2;
			S+=det;
			G=G+(p[i]+p[i+1])/3*det;
		}
		for(int i=0;i<left.size();i++){
			Pt c=left[i].first;
			int at=left[i].second;
			double r=(p[at]-c).ABS();
			double th=asin(((p[at+1]-p[at]).ABS()/2)/r)*2;
			double A=r*r*(th-sin(th))/2;
			double y=2*r*r*r*sin(th/2)*sin(th/2)*sin(th/2)/A/3;
			Pt v=(p[at+1]+p[at])/2-c;
			v=v/v.ABS()*y;
			v=v+c;
			S+=A;
			G=G+v*A;
		}
		for(int i=0;i<right.size();i++){
			Pt c=right[i].first;
			int at=right[i].second;
			double r=(p[at]-c).ABS();
			double th=asin(((p[at+1]-p[at]).ABS()/2)/r)*2;
			double A=r*r*(th-sin(th))/2;
			double y=2*r*r*r*sin(th/2)*sin(th/2)*sin(th/2)/A/3;
			Pt v=(p[at+1]+p[at])/2-c;
			v=v/v.ABS()*y;
			v=v+c;
			S-=A;
			G=G-v*A;
		}
		G=G/S;
		//printf("%f ",S);
		printf("%f %f ",G.x,G.y);
		Pt os=Pt(114,514);
		bool ok=false;
		for(int i=0;i<sz;i++){
			if(isl[i]){
				if(iSS(p[i],p[i+1],os,G)){
				//	printf("line %d\n",i);
					ok=!ok;
				}
			}
		}
		for(int i=0;i<left.size();i++){
			Pt v1=p[left[i].second]-left[i].first;
			Pt v2=p[left[i].second+1]-left[i].first;
			if(v1.det(v2)>EPS)swap(v1,v2);
			double r=v1.ABS();
			if(iCS(left[i].first,r,os,G)){
				pair<Pt,Pt>hb=pCL(left[i].first,r,os,G);
				//printf("left %d\n",i);
				if(iSP(os,hb.first,G)==2&&sAP(v1,hb.first-left[i].first,v2)==1){ok=!ok;}
				if(iSP(os,hb.second,G)==2&&sAP(v1,hb.second-left[i].first,v2)==1){ok=!ok;}
			}
		}
		for(int i=0;i<right.size();i++){
			Pt v1=p[right[i].second]-right[i].first;
			Pt v2=p[right[i].second+1]-right[i].first;
			if(v1.det(v2)>EPS)swap(v1,v2);
			double r=v1.ABS();
			if(iCS(right[i].first,r,os,G)){
				pair<Pt,Pt>hb=pCL(right[i].first,r,os,G);
				if(iSP(os,hb.first,G)==2&&sAP(v1,hb.first-right[i].first,v2)==1){ok=!ok;}
				if(iSP(os,hb.second,G)==2&&sAP(v1,hb.second-right[i].first,v2)==1){ok=!ok;}
			}
		}
		if(ok)printf("+\n");
		else printf("-\n");
	}
}