Problem:
Lifeguard in the Pool

プールの監視員

English text is not available in this practice contest.

Horton Moore はプールの監視員として働いている．彼が見回りのためにプールの縁を歩いていたところ，プールの中で一人の少女がおぼれていることに気づいた．もちろん，彼は直ちに救助に向かわなければならない．しかも，少女の身に何かあっては大変であるから，少しでも早く少女のもとにたどり着きたい．

あなたの仕事は，プールの形状（頂点数が 3〜10 の凸多角形），地上および水中における監視員の単位距離あたりの移動時間，そして監視員と少女の初期位置が与えられたときに，監視員が少女のところに到着するまでにかかる最短の時間を求めるプログラムを書くことである．

Input

入力は複数のデータセットの並びからなる．入力の終わりは 1 つの 0 だけを含む行によって示される．

各データセットは次の形式になっている．

n

x1 y1 x2 y2 ... xn yn

tg

tw

xs ys

xt yt

それぞれの記号の意味は次のとおりである．

n は凸多角形をしたプールの頂点数を示す．これは 3 以上 10 以下の整数である．

(xi, yi) はプールの i 番目の頂点の座標を示す．それぞれの座標値は絶対値が 100 以下の整数である．頂点は反時計回りの順番で与えられる．

tg は監視員が地上において移動するときにかかる単位距離あたりの時間を表す．tw は監視員が水中において移動するときにかかる単位距離あたりの時間を表す．これらはいずれも整数であり，さらに 1 ≦ tg < tw ≦ 100 を満たす．

(xs, ys) は監視員の初期位置の座標を表す．この座標はプールのちょうど辺上にある．

(xt, yt) は少女の初期位置の座標を表す．この座標はプールの内側にある．

同一の行にある数値と数値の間は 1 個の空白で区切られている．

この問題において，監視員および少女は点であるとみなす．また，監視員がプールの辺に沿って移動するときは地上を移動しているとみなす．監視員は，地上から水中に一瞬で入ることが，また水中から地上に一瞬で出ることができると仮定して構わない．監視員が地上から水中に，あるいは水中から地上に移るとき，監視員は同じ座標にとどまると考えること．したがって，たとえばプールの縁から離れたところに飛び込むことによって水中での移動距離を減らすことはできない．

Output

各データセットに対して，監視員が少女の元に到着するまでにかかる最短の時間を 1 行に出力しなさい．解答の誤差は 0.00000001 (10−8) を超えてはならない．精度に関する条件を満たしていれば，小数点以下は何桁数字を出力しても構わない．

Sample Input

4
0 0 10 0 10 10 0 10
10
12
0 5
9 5
4
0 0 10 0 10 10 0 10
10
12
0 0
9 1
4
0 0 10 0 10 10 0 10
10
12
0 1
9 1
8
2 0 4 0 6 2 6 4 4 6 2 6 0 4 0 2
10
12
3 0
3 5
0

Output for the Sample Input

108.0
96.63324958071081
103.2664991614216
60.0

Solution in C++ (judge verdict: Accepted)

#include<bits/stdc++.h>
#define MAX 100
#define inf 1<<29
#define linf 1e16
#define eps (1e-8)
#define mod 1000000007
#define pi acos(-1)
#define f first
#define s second
#define mp make_pair
#define pb push_back
#define all(a) (a).begin(),(a).end()
#define pd(a) printf("%.10f\n",(double)(a))
#define FOR(i,a,b) for(int i=(a);i<(b);i++)
#define RFOR(i,a,b) for(int i=(a)-1;(b)<=i;i--)
#define equals(a,b) (fabs((a)-(b))<eps)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<int,double> pid;
typedef pair<double,int> pdi;
typedef vector<int> vi;
typedef vector<pii> vpi;
const int dx[8]={1,0,-1,0,1,1,-1,-1};
const int dy[8]={0,1,0,-1,1,-1,1,-1};

  
class Point{
public:
  double x,y;
  Point(double x=0,double y=0):x(x),y(y){}

  Point operator+(Point p){ return Point(x+p.x,y+p.y);}
  Point operator-(Point p){ return Point(x-p.x,y-p.y);}
  Point operator*(double k){ return Point(x*k,y*k);}
  Point operator/(double k){ return Point(x/k,y/k);}
  bool operator<(Point p)const{ return (x!=p.x ? x<p.x : y<p.y);}
  bool operator==(Point p)const{ return fabs(x-p.x)<eps && fabs(y-p.y)<eps;}

  double abs(){ return sqrt(norm());}
  double norm(){ return (x*x+y*y);}
};
typedef Point Vector;
typedef vector<Point> Polygon;

class Segment{
public:
  Point p1,p2;
  Segment(Point p1=Point(),Point p2=Point()):p1(p1),p2(p2){}
};
typedef Segment Line;

double norm(Vector a){ return (a.x*a.x+a.y*a.y);}
double abs(Vector a){ return sqrt(norm(a));}
double dot(Vector a,Vector b){ return (a.x*b.x+a.y*b.y);}
double cross(Vector a,Vector b){ return (a.x*b.y-a.y*b.x);}

int ccw(Point p0,Point p1,Point p2){
  Vector a=p1-p0;
  Vector b=p2-p0;
  if(cross(a,b)>eps)return 1;
  if(cross(a,b)<-eps)return -1;
  if(dot(a,b)<-eps)return 2;
  if(a.norm()<b.norm())return -2;
  return 0;
}

vector<Point> unique(vector<Point> vp){
  vector<Point> res;
  if(vp.empty())return res;
  res.pb(vp[0]);
  FOR(i,0,vp.size())if(!(res.back()==vp[i]))res.pb(vp[i]);
  return res;
}

int n;
double tg,tw;
Polygon p;
vector<Point> vp;
vector<pid> g[MAX];
Point s,t;
int S,T;

void init(){
  p.clear();
  vp.clear();
  FOR(i,0,MAX)g[i].clear();
}

Point getPoint(Point a,Point b,Point c){
  Vector v=b-a;
  v=v/abs(v);
  double L=0.0,R=abs(b-a);
  FOR(i,0,100){
    double m1=(L+L+R)/3.0;
    double m2=(L+R+R)/3.0;
    double res1=(abs(v*m1))*tg+abs(a+v*m1-c)*tw;
    double res2=(abs(v*m2))*tg+abs(a+v*m2-c)*tw;
    if(res1<res2)R=m2;
    else L=m1;
  }
  return a+v*L;
}

double dijkstra(){
  double d[MAX];
  priority_queue<pdi,vector<pdi>,greater<pdi> > pq;
  fill(d,d+MAX,inf);
  d[S]=0;
  pq.push(mp(0,S));

  while(pq.size()){
    pdi u=pq.top();
    pq.pop();
    if(d[u.s]-u.f<-eps)continue;
    if(u.s==T)return d[u.s];

    FOR(i,0,g[u.s].size()){
      int next=g[u.s][i].f;
      double cost=d[u.s]+g[u.s][i].s;
      if(cost-d[next]<-eps){
        d[next]=cost;
        pq.push(mp(cost,next));
      }
    }
  }
  return inf;
}

double solve(){
  FOR(i,0,n){
    vp.pb(getPoint(p[i],p[(i+1)%n],s));
    vp.pb(getPoint(p[i],p[(i+1)%n],t));
    vp.pb(getPoint(p[(i+1)%n],p[i],s));
    vp.pb(getPoint(p[(i+1)%n],p[i],t));
  }
  vp.pb(s);
  vp=unique(vp);
  FOR(i,0,n){
    Segment s(p[i],p[(i+1)%n]);
    vector<pdi> list;
    FOR(j,0,vp.size()){
      if(ccw(s.p1,s.p2,vp[j])==0)list.pb(mp(abs(s.p1-vp[j]),j));
    }
    sort(all(list));
    FOR(j,0,list.size()-1){
      int a=list[j].s,b=list[j+1].s;
      g[a].pb(mp(b,abs(vp[a]-vp[b])*tg));
      g[b].pb(mp(a,abs(vp[a]-vp[b])*tg));
    }
  }
  vp.pb(t);
  FOR(i,0,vp.size()){
    if(vp[i]==s)S=i;
    if(vp[i]==t)T=i;
  }
  FOR(i,0,vp.size())if(i!=T)g[i].pb(mp(T,abs(vp[i]-t)*tw));
  FOR(i,0,vp.size())if(i!=S)g[S].pb(mp(i,abs(vp[i]-s)*tw));

  return dijkstra();
}

int main()
{
  while(cin>>n && n){
    init();
    FOR(i,0,n){
      int x,y;
      cin>>x>>y;
      p.pb(Point(x,y));
      vp.pb(Point(x,y));
    }
    cin>>tg>>tw;
    cin>>s.x>>s.y>>t.x>>t.y;
    pd(solve());
  }
  return 0;
}