poj1556 && nyoj227

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <stack>
#include <queue>
#include <vector>
#include <cmath>
#include <map>
#include <set>
#define ll long long
#define llu unsigned long long
using namespace std;
const double eps = 1e-6;
struct Point {
double x,y;
Point (double x = 0,double y = 0) : x(x),y(y) {}
};

Point p[150];

struct Line{
Point a,b;
};

Line line[100];
int indl,indp;
double Map[150][150];
bool dcmp(double x) {
if(fabs(x)-0.0 <= eps) return true;
return false;
}
void init(){
for(int i=1;i<=140; i++) {
for(int j=1; j<=140; j++) Map[i][j] = 100000000.0;
//            Map[i][j] = (i!=j) ? 100000000.0 : 0;
}
}
double dis(Point a,Point b) {
return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}
double cross(Point a,Point b,Point c) {
if(dcmp((b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y))) return 0;
return ((b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y));
}
bool f(Point a,Point b,Point c,Point d) {
if(min(a.x,b.x)<=max(c.x,d.x) && min(c.x,d.x)<=max(a.x,b.x)
&& min(a.y,b.y)<=max(c.y,d.y) && min(c.y,d.y)<=max(a.y,b.y)){
if(cross(a,b,c)*cross(a,b,d) < 0 && cross(c,d,a)*cross(c,d,b) < 0) return true;
}
return false;
}
bool pd(Point x,Point y){
for(int i=1; i<indl; i++) {
if(f(x,y,line[i].a,line[i].b)) {
return false;
}
}
return true;
}
void jian(){
for(int i=1;i<=indp; i++) {
if(i == 1 || i == indp) {
for(int j=1; j<=indp; j++) {
if(i!=j && pd(p[i],p[j])) {
Map[i][j] = min(Map[i][j],dis(p[i],p[j]));
}
}
continue;
}
for(int j=1; j<=indp; j++) {
if(i != j && (((i-2)/4+1) != ((j-2)/4+1)) && pd(p[i],p[j])) {
Map[i][j] = min(Map[i][j],dis(p[i],p[j]));
}
}
}
}
int main(){
int n;
while(scanf("%d",&n)!=EOF) {
if(n == -1) break;
init();
Point A,B;
indl = 1,indp = 1;
p[indp].x = 0, p[indp].y = 5;
indp++;
double a,b,c,d,e;
for(int i=1;i<=n;i++) {
scanf("%lf%lf%lf%lf%lf",&a,&b,&c,&d,&e);
p[indp].x = a; p[indp++].y = b;
p[indp].x = a; p[indp++].y = c;
p[indp].x = a; p[indp++].y = d;
p[indp].x = a; p[indp++].y = e;
A.x = a; A.y = 0;
B.x = a; B.y = b;
line[indl].a = A; line[indl++].b = B;
A.x = a; A.y = c;
B.x = a; B.y = d;
line[indl].a = A; line[indl++].b = B;
A.x = a; A.y = e;
B.x = a; B.y = 10;
line[indl].a = A; line[indl++].b = B;
}
p[indp].x = 10; p[indp].y = 5;
jian();  //建图
// 最短路
for(int k=1; k<=indp; k++) {
for(int i=1; i<=indp; i++) {
for(int j=1; j<=indp; j++) {
if(Map[i][k]!=100000000.0 && Map[k][j] != 100000000.0) {
Map[i][j] = min(Map[i][j],Map[i][k]+Map[k][j]);
}
}
}
}
printf("%.2lf\n",Map[1][indp]);
}
return 0;
}