読者です 読者をやめる 読者になる 読者になる

AOJ 2016: プールの監視員

2016年になってしまい仕方なく解いた。
と思ったら拍子抜けするほど簡単だった……ええ何これは……。
まあ確かに点列挙して最短路みたいな幾何問題は無限に練習したしなあ……。

#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
#include<queue>
using namespace std;
const double EPS = 1e-10;
const double INF = 1e+10;
const double PI = acos(-1);
Pt p[20];
Pt h[20];
double ijk[1100];
int v[1100];
int main(){
	int a;
	while(scanf("%d",&a),a){
		double tg,tw;
		for(int i=0;i<a;i++){
			double X,Y;
			scanf("%lf%lf",&X,&Y);
			p[i]=Pt(X,Y);
		}
		p[a]=p[0];
		scanf("%lf%lf",&tg,&tw);
		Pt s,t;
		double X,Y;
		scanf("%lf%lf",&X,&Y);
		s=Pt(X,Y);
		scanf("%lf%lf",&X,&Y);
		t=Pt(X,Y);
		double th=asin(tg/tw);
		for(int i=0;i<a;i++){
			h[i]=hLP(p[i],p[i+1],t);
		}
		vector<Pt> hb;
		hb.push_back(s);
		hb.push_back(t);
		for(int i=0;i<a;i++){
			hb.push_back(p[i]);
		}
		for(int i=0;i<a;i++){
			for(int j=0;j<a;j++){
				if(i==j||i==j+1)continue;
				Pt tmp=hLP(p[j],p[j+1],p[i]);
				double dist=(tmp-p[i]).ABS();
				double len=dist*tan(th);
				Pt vec=p[j+1]-p[j];
				vec=vec/(vec.ABS())*len;
				Pt t1=tmp+vec;
				Pt t2=tmp-vec;
				if(iSP(p[j],t1,p[j+1])==2)hb.push_back(t1);
				if(iSP(p[j],t2,p[j+1])==2)hb.push_back(t2);
			}
		}
		for(int j=0;j<a;j++){
			Pt tmp=hLP(p[j],p[j+1],t);
			double dist=(tmp-t).ABS();
			double len=dist*tan(th);
			Pt vec=p[j+1]-p[j];
			vec=vec/(vec.ABS())*len;
			Pt t1=tmp+vec;
			Pt t2=tmp-vec;
			if(iSP(p[j],t1,p[j+1])==2)hb.push_back(t1);
			if(iSP(p[j],t2,p[j+1])==2)hb.push_back(t2);
		}
		for(int j=0;j<a;j++){
			Pt tmp=hLP(p[j],p[j+1],s);
			double dist=(tmp-s).ABS();
			double len=dist*tan(th);
			Pt vec=p[j+1]-p[j];
			vec=vec/(vec.ABS())*len;
			Pt t1=tmp+vec;
			Pt t2=tmp-vec;
			if(iSP(p[j],t1,p[j+1])==2)hb.push_back(t1);
			if(iSP(p[j],t2,p[j+1])==2)hb.push_back(t2);
		}
		int sz=hb.size();
		for(int i=0;i<sz;i++){
			v[i]=0;
			ijk[i]=999999999;
		}
		ijk[0]=0;
		for(int i=0;i<sz;i++){
			int at=0;
			double tmp=999999999;
			for(int j=0;j<sz;j++){
				if(!v[j]&&tmp>ijk[j]){
					at=j;tmp=ijk[j];
				}
			}
			v[at]=1;
			for(int j=0;j<sz;j++){
				if(v[j])continue;
				Pt M=(hb[at]+hb[j])/2;
				double toc=ijk[at];
				if(sGP(a,p,M)==0)toc+=(hb[at]-hb[j]).ABS()*tg;
				else toc+=(hb[at]-hb[j]).ABS()*tw;
				ijk[j]=min(ijk[j],toc);
			}
		}
		//for(int i=0;i<sz;i++)printf("(%f %f): %f\n",hb[i].x,hb[i].y,ijk[i]);
		printf("%.12f\n",ijk[1]);
	}
}