Problem:
Problem E: Kuru-Kuru Robot

ロボット工学の研究者であるアシモフ博士は、新しいロボットの開発に成功した。このロボットは、床に張り巡らせた特殊なワイヤーに沿って自由に動くことができる。ロボットは常に現在いるワイヤーに平行な方向を向いて前進する。

すばらしいことに、このロボットはワイヤー上を移動するためにその距離にかかわらずエネルギーを消費しない。しかし、ワイヤーの端点やワイヤーの交点で向きを変えるためにはその点を中心に回転する必要があり、その回転角度だけエネルギーを消費する。

このロボットをスタート地点からゴール地点まで動かしたい。例えば、下図においてロボットを start に東向き（下図 x軸の正の向き）で配置し goal まで動かしてみよう。start から交点 A 交点 B を経由して goal にたどり着く経路と、交点 A 交点 C を経由する経路があるが、前者の方が移動距離が長いものの回転角度が少ないためエネルギー消費量は少ない。

あなたの仕事は、床にあるワイヤーの位置、スタート地点、ゴール地点を読み込み、スタート地点からゴール地点へ移動するために必要な回転角度の合計の最小値を報告するプログラムを作成することである。

ロボットはスタート地点においてどのような向きに配置してもよく、ゴール地点ではどの方向を向いていてもよい。

Input

入力は複数のデータセットからなる。各データセットの形式は以下のとおり：

n
x11 y11 x21 y21
x12 y12 x22 y22
.
.
.
x1n y1n x2n y2n
sx sy gx gy

n はワイヤーの数を示す整数である。x1i y1i は i 番目のワイヤーの一方の端点の座標、x2i y2i は i 番目のワイヤーのもう一方の端点の座標を示す。

sx sy、gx gy はそれぞれスタート地点とゴール地点の座標を示す。

与えられる座標はすべて整数値であり、x軸 y軸の値は -1000　以上 1000 以下である。また、スタート地点とゴール地点は与えられるワイヤーの端点上にあると仮定してよい。

n ≤ 50 であると仮定してよい。

n が 0 のとき入力の終了を示す。

Output

各データセットに対して、回転角度の合計の最小値を１行に出力せよ。ロボットがゴール地点にたどり着けない場合は "-1" と出力せよ。出力は 0.00001 以下の誤差を含んでもよい。

Sample Input

8
1 3 11 3
5 3 13 11
10 10 17 10
11 2 11 6
10 5 14 5
13 6 13 2
17 10 17 3
13 3 19 3
1 3 19 3
6
1 3 11 3
5 3 13 11
10 10 17 10
11 2 11 6
10 5 14 5
13 3 19 3
1 3 19 3
2
0 0 7 0
3 0 7 3
0 0 7 3
0

Output for the Sample Input

270.0000
-1
36.86989765

Solution in JAVA (judge verdict: Wrong Answer)

import java.util.ArrayList;
import java.util.Arrays;
import java.util.Collections;
import java.util.PriorityQueue;
import java.util.Scanner;

class Main {

	static int n;
	static double[] x1, y1, x2, y2;
	static double sx, sy, gx, gy;
	static Point[][] seg;
	static ArrayList<Integer>[] g;

	@SuppressWarnings("unchecked")
	public static void main(String[] args) {
		Scanner sc = new Scanner(System.in);
		ArrayList<Point> ps = new ArrayList<>();
		while (true) {
			System.out.println("In preparation.");
			n = sc.nextInt();
			if (n == 0)
				break;
			x1 = new double[n];
			y1 = new double[n];
			x2 = new double[n];
			y2 = new double[n];
			seg = new Point[n][2];
			for (int i = 0; i < n; ++i) {
				x1[i] = sc.nextDouble();
				y1[i] = sc.nextDouble();
				x2[i] = sc.nextDouble();
				y2[i] = sc.nextDouble();
				Point p1 = new Point(x1[i], y1[i]);
				Point p2 = new Point(x2[i], y2[i]);
				ps.add(p1);
				ps.add(p2);
				seg[i] = new Point[] { p1, p2 };
			}
			sx = sc.nextDouble();
			sy = sc.nextDouble();
			gx = sc.nextDouble();
			gy = sc.nextDouble();
			for (int i = 0; i < n; ++i) {
				for (int j = i + 1; j < n; ++j) {
					if (linesIntersect(seg[i][0], seg[i][1], seg[j][0], seg[j][1])) {
						Point p = crossPoint(seg[i][0], seg[i][1], seg[j][0], seg[j][1]);
						ps.add(p);
					}
				}
			}

			Collections.sort(ps);
			for (int i = 0; i < ps.size(); ++i) {
				while (i + 1 < ps.size() && Math.abs(ps.get(i).x - ps.get(i + 1).x) < 1e-6
						&& Math.abs(ps.get(i).y - ps.get(i + 1).y) < 1e-6) {
					ps.remove(i + 1);
				}
			}

			int N = ps.size();
			g = new ArrayList[N];
			for (int i = 0; i < N; ++i)
				g[i] = new ArrayList<>();

			for (int i = 0; i < N; ++i)
				ps.get(i).id = i;
			for (int i = 0; i < n; ++i) {
				ArrayList<Point> tmp = new ArrayList<>();
				for (Point p : ps) {
					if (ccw(seg[i][0], seg[i][1], p) == 0) {
						tmp.add(p);
					}
				}
				Collections.sort(tmp);
				for (int j = 0; j + 1 < tmp.size(); ++j) {
					g[tmp.get(j).id].add(tmp.get(j + 1).id);
					g[tmp.get(j + 1).id].add(tmp.get(j).id);
				}
			}

			int srcId = -1;
			int dstId = -1;
			for (int i = 0; i < N; ++i) {
				if ((Math.abs(ps.get(i).x - sx) + Math.abs(ps.get(i).y - sy)) < 1e-6) {
					srcId = i;
				}
				if ((Math.abs(ps.get(i).x - gx) + Math.abs(ps.get(i).y - gy)) < 1e-6) {
					dstId = i;
				}
			}

			double[][] dist = new double[N][N];
			for (int i = 0; i < N; ++i) {
				for (int j = 0; j < N; ++j) {
					dist[i][j] = Double.MAX_VALUE / 10;
				}
			}

			PriorityQueue<State> pq = new PriorityQueue<>();
			for (int dst : g[srcId]) {
				pq.add(new State(srcId, dst, 0));
				dist[srcId][dst] = 0;
			}
			while (!pq.isEmpty()) {
				State s = pq.poll();
				for (int dst : g[s.dst]) {
					double angle = angle(ps.get(s.src), ps.get(s.dst), ps.get(s.dst), ps.get(dst));
					if (dist[s.dst][dst] > s.cost + angle) {
						dist[s.dst][dst] = s.cost + angle;
						pq.add(new State(s.dst, dst, s.cost + angle));
					}
				}
			}

			double ans = Double.MAX_VALUE / 10;
			for (int i = 0; i < N; ++i) {
				ans = Math.min(ans, dist[i][dstId]);
			}
			if (ans >= Double.MAX_VALUE / 10) {
				System.out.println(-1);
			} else {
				ans /= Math.PI;
				ans *= 180;
				System.out.println(ans);
			}
		}
	}

	static class Point implements Comparable<Point> {
		double x, y;
		int id;

		public Point(double x1, double y1) {
			this.x = x1;
			this.y = y1;
		}

		double norm(Point p) {
			double dx = x - p.x;
			double dy = y - p.y;
			return Math.sqrt(dx * dx + dy * dy);
		}

		@Override
		public int compareTo(Point o) {
			if (x != o.x)
				return Double.compare(x, o.x);
			else
				return Double.compare(y, o.y);
		}
	}

	static class State implements Comparable<State> {
		int src, dst;
		double cost;

		public State(int src, int dst, double cost) {
			this.src = src;
			this.dst = dst;
			this.cost = cost;
		}

		public int compareTo(State o) {
			return Double.compare(cost, o.cost);
		};
	}

	static Point crossPoint(Point s1, Point d1, Point s2, Point d2) {
		double[][] mx = { { d1.x - s1.x, -(d2.x - s2.x) }, { d1.y - s1.y, -(d2.y - s2.y) } };
		double[] v = { -s1.x + s2.x, -s1.y + s2.y };
		double[][] rev = { { mx[1][1], -mx[0][1] }, { -mx[1][0], mx[0][0] } };
		double determinent = mx[0][0] * mx[1][1] - mx[0][1] * mx[1][0];
		rev[0][0] /= determinent;
		rev[0][1] /= determinent;
		rev[1][0] /= determinent;
		rev[1][1] /= determinent;
		double[] coe = mulVec(rev, v);
		Point p = new Point((1 - coe[0]) * s1.x + coe[0] * d1.x, (1 - coe[0]) * s1.y + coe[0] * d1.y);
		return p;
	}

	static double[] mulVec(double[][] mx, double[] v) {
		double[] retVec = new double[2];
		for (int i = 0; i < 2; ++i) {
			for (int j = 0; j < 2; ++j) {
				retVec[i] += mx[i][j] * v[j];
			}
		}
		return retVec;
	}

	static double[][] mul(double[][] m1, double[][] m2) {
		double[][] ret = new double[2][2];
		for (int i = 0; i < 2; ++i) {
			for (int j = 0; j < 2; ++j) {
				for (int k = 0; k < 2; ++k) {
					ret[i][j] += m1[i][k] * m2[k][j];
				}
			}
		}
		return ret;
	}

	static boolean linesIntersect(Point s1, Point d1, Point s2, Point d2) {
		long prd1 = ccw(s1, d1, d2) * ccw(s1, d1, s2);
		long prd2 = ccw(s2, d2, s1) * ccw(s2, d2, d1);
		return prd1 <= 0 && prd2 <= 0;
	}

	static long ccw(Point s, Point d, Point p) {
		double cp = crossProduct(s, d, s, p);
		if (Math.abs(cp) > 1e-6)// cp!=0
			return Double.compare(cp, 0);
		else {
			if (dotProduct(s, d, s, p) < 0)
				return -2;
			else if (dotProduct(d, s, d, p) < 0)
				return 2;
			else
				return 0;
		}
	}

	static double angle(Point s1, Point d1, Point s2, Point d2) {
		double dotprd = dotProduct(s1, d1, s2, d2);
		dotprd /= s1.norm(d1);
		dotprd /= s2.norm(d2);
		return Math.acos(dotprd);
	}

	static double crossProduct(Point s1, Point d1, Point s2, Point d2) {
		return (d1.x - s1.x) * (d2.y - s2.y) - (d1.y - s1.y) * (d2.x - s2.x);
	}

	static double dotProduct(Point s1, Point d1, Point s2, Point d2) {
		return (d1.x - s1.x) * (d2.x - s2.x) + (d1.y - s1.y) * (d2.y - s2.y);
	}

	static void tr(Object... objects) {
		System.out.println(Arrays.deepToString(objects));
	}
}