注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

fudq's AC Road

何以解忧,唯有AC!

 
 
 

日志

 
 

hdu 3400 Line belt  

2013-09-12 12:17:00|  分类: ACM-hdu |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
有两条线段AB和CD,在AB上运行的速度为P,在CD上运行的速度为Q,在其它地方运行的速度为R,求A运动到D需要的最短时间。
运行轨迹可以归结为先在AB上运行一段,然后在其它地方运行一段,然后在CD上运行一段。
先三分在AB上的点,然后三分在CD上的点,这样三分嵌套三分解决。
*注:刚开始三分时是按照点在线段上运行的距离三分,会有精度问题,wa数次;改成按照在线段上运行的距离对线段总长的比值三分就过了。。。

/*
* pro.cpp
*
* Created on: 2013-09-12
* Author: fudq
*/
#include <functional>
#include <algorithm>
#include <iostream>
//#include <fstream>
#include <sstream>
#include <iomanip>
#include <numeric>
#include <cstring>
#include <cassert>
#include <cstdio>
#include <string>
#include <vector>
#include <bitset>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <list>
#include <set>
#include <map>
using namespace std;
//#pragma comment(linker,"/STACK:102400000,102400000")

#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define LL __int64

const int N=310;
const int M=2000000000;
const int MOD=20100501;
const int INF=0x7fffffff;
const int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
//const int dir[8][2]={{-1,0},{1,0},{0,-1},{0,1},{-1,1},{1,-1},{-1,-1},{1,1}};
const double eps=1e-8;
const double PI=acos(-1.0);

inline int sign(double x){return (x>eps)-(x<-eps);}
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T lcm(T a,T b,T d){return a/d*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
template<class T> inline T Min(T a,T b,T c){return min(min(a, b),c);}
template<class T> inline T Max(T a,T b,T c){return max(max(a, b),c);}
template<class T> inline T Min(T a,T b,T c,T d){return min(min(a, b),min(c,d));}
template<class T> inline T Max(T a,T b,T c,T d){return max(max(a, b),max(c,d));}
/*************************/

double ax,ay,bx,by,cx,cy,dx,dy,p,q,r;
double disab,discd,tax,tay,tbx,tby;

double dist(double ax,double ay,double bx,double by){
return sqrt((ax-bx)*(ax-bx)+(ay-by)*(ay-by));
}

double Calc2(double a)
{
tbx=dx+a*(cx-dx);
tby=dy+a*(cy-dy);
double t3=dist(dx,dy,tbx,tby)/q;
double t2=dist(tax,tay,tbx,tby)/r;
return t2+t3;
}

double solve2(double sta,double end)
{
double Left, Right;
double mid, midmid;
double mid_area, midmid_area;
Left = sta; Right = end;
while (Left + eps < Right)
{
mid = (Left + Right) / 2.0;
midmid = (mid + Right) / 2.0;
mid_area = Calc2(mid);
midmid_area = Calc2(midmid);
if (mid_area+eps < midmid_area) Right = midmid;
else Left = mid;
}
return mid_area;
}

double Calc(double a)
{
tax=ax+(bx-ax)*a;
tay=ay+(by-ay)*a;
return dist(ax,ay,tax,tay)/p+solve2(0,1.0);
}

void solve(double sta,double end)
{
double Left, Right;
double mid, midmid;
double mid_area, midmid_area;
Left = sta; Right = end;
while (Left + eps < Right)
{
mid = (Left + Right) / 2.0;
midmid = (mid + Right) / 2.0;
mid_area = Calc(mid);
midmid_area = Calc(midmid);
if (mid_area+eps < midmid_area) Right = midmid;
else Left = mid;
}
printf("%.2lf\n",mid_area);
}

int main()
{
#ifndef ONLINE_JUDGE
freopen("testin.txt","r",stdin);
// freopen("testout.txt","w",stdout);
#endif
int T;
scanf("%d",&T);
while(T--)
{
scanf("%lf%lf%lf%lf",&ax,&ay,&bx,&by);
scanf("%lf%lf%lf%lf",&cx,&cy,&dx,&dy);
scanf("%lf%lf%lf",&p,&q,&r);
disab=dist(ax,ay,bx,by);
discd=dist(cx,cy,dx,dy);
solve(0,1.0);
}
return 0;
}



  评论这张
 
阅读(56)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017