现在的位置: 首页 > 综合 > 正文

hdu3400 Line belt

2018年04月23日 ⁄ 综合 ⁄ 共 1891字 ⁄ 字号 评论关闭

如图所示,必然在AB,CD上存在两点e,f使A到D的时间最小,由此可以假设f点已知,由此确定下e点,所以采用嵌套三分

code

/*
ID: yueqiq
PROG: numtri
LANG: C++
*/
#include <set>
#include <map>
#include <ctime>
#include <queue>
#include <cmath>
#include <stack>
#include <limits>
#include <vector>
#include <bitset>
#include <string>
#include <cstdio>
#include <cstring>
#include <fstream>
#include <string.h>
#include <iostream>
#include <algorithm>
#define Si set<int>
#define LL long long
#define pb push_back
#define PS printf(" ")
#define Vi vector<int>
#define LN printf("\n")
#define lson l,m,rt << 1
#define rson m+1,r,rt<<1|1
#define SD(a) scanf("%d",&a)
#define PD(a) printf("%d",a)
#define SET(a,b) memset(a,b,sizeof(a))
#define FF(i,a) for(int i(0);i<(a);i++)
#define FD(i,a) for(int i(a);i>=(1);i--)
#define FOR(i,a,b) for(int i(a);i<=(b);i++)
#define FOD(i,a,b) for(int i(a);i>=(b);i--)
#define readf freopen("input.txt","r",stdin)
#define writef freopen("output.txt","w",stdout)
const int maxn = 1001;
const int INF = 1111;
const int dx[]={0,1,0,-1};
const int dy[]={1,0,-1,0};
const double pi = acos(-1.0);
const double eps= 1e-5;
using namespace std;
struct point{
    double x,y;
}A,B,C,D;
double P,Q,R;
double dist(point a,point b){
     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
//此时假设CD段上的点f为已知的,此时来枚举e点
double cal_ab(point a,point b,point f){
    point l=a,r=b,m,mid;
    double ans1,ans2;
    do{
        m.x=(l.x+r.x)/2.0;
        m.y=(l.y+r.y)/2.0;
        mid.x=(l.x+m.x)/2.0;
        mid.y=(l.y+m.y)/2.0;
        ans1=dist(a,m)/P+dist(m,f)/R;
        ans2=dist(a,mid)/P+dist(mid,f)/R;
        if(ans1<ans2) l=mid;
        else r=m;
    }while(fabs(ans1-ans2)>eps);
    return ans1;
}
double cal_abcd(point a,point b,point c,point d){
    point l=c,r=d,m,mid;
    double ans1,ans2;
    do{
        m.x=(l.x+r.x)/2.0;
        m.y=(l.y+r.y)/2.0;
        mid.x=(l.x+m.x)/2.0;
        mid.y=(l.y+m.y)/2.0;
        ans1=dist(m,d)/Q+cal_ab(a,b,m);
        ans2=dist(mid,d)/Q+cal_ab(a,b,mid);
        if(ans1<ans2) l=mid;
        else r=m;
    }while(fabs(ans1-ans2)>eps);
    return ans1;
}
int main()
{
    int cas;
    SD(cas);
    while(cas--){
        scanf("%lf%lf%lf%lf",&A.x,&A.y,&B.x,&B.y);
        scanf("%lf%lf%lf%lf",&C.x,&C.y,&D.x,&D.y);
        scanf("%lf%lf%lf",&P,&Q,&R);
        printf("%.2lf\n",cal_abcd(A,B,C,D));
    }
    return 0;
}


【上篇】
【下篇】

抱歉!评论已关闭.