当前位置:编程学习 > C/C++ >>

hdu_4170 Supply Mission (计算几何)

题意:
飞机在位置(x0,y0),飞行速度为v km/h,有N(0<N<8)艘船分别为(xi,yi)速度向量为(vxi,vyi) km/h,坐标单位为km
飞机必须在每艘船上要一小时卸载货物,最后飞回原来的位置(x0,y0),求最少时间花费,用时分秒输出。
思路:
想了会儿,没什么思路,看了看数据范围,N最大才8,时限10秒,完全可以暴力嘛。
直接把到达船的顺序做个全排列,取最小的一次时间就行了,算相遇时要一定要细心。
还有个小trick,秒进位后等于60,这里搞了我好久好久。。。
[cpp] 
/*
program:hdu_4170 & hunnu_11215
author:BlackAndWhite
*/ 
#include<stdio.h> 
#include<math.h> 
#include<algorithm> 
#define eps 1e-6 
using namespace std; 
struct point 

    int x,y; 
}a[10],v[10],s; 
double t,ans; 
int N,sv,i,p[10],T[3],js=1; 
double crosstime(double x0,double y0,double v0,double x1,double y1,double vx,double vy) 
{//表示飞机在(x0,y0),速度为v0,要与船(x1,y1)速度向量为(vx,vy)相遇所花费的时间。  
    double a,b,c,t; 
    a=vx*vx+vy*vy-v0*v0; 
    b=2*(vx*(x1-x0)+vy*(y1-y0)); 
    c=(x1-x0)*(x1-x0)+(y1-y0)*(y1-y0); 
    if(a<eps&&a>-eps) return -c/b; 
    if(a<eps) {a=-a;b=-b;c=-c;} 
    return (-b+sqrt(b*b-4*a*c))/(2*a); 

double back(double x0,double y0,double x1,double y1,double sv)//最后回到起点的时间花费  

    return sqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0))/sv; 

void trs(double ans,int T[3])//小数形式转换为时分秒  

    T[0]=(int)ans; 
    ans-=T[0]; 
    T[1]=(int)(ans*=60.0);ans-=T[1]; 
    T[2]=(int)(ans*60.0+0.9999999); 
    if(T[2]==60) {T[2]=0;T[1]++;}//trick,进位啊。。。  
    if(T[1]==60) {T[1]=0;T[0]++;} 

int main() 

    while(scanf("%d",&N),N) 
    { 
        for(i=0;i<N;i++) scanf("%d%d%d%d",&a[i].x,&a[i].y,&v[i].x,&v[i].y); 
        scanf("%d%d%d",&s.x,&s.y,&sv); 
        for(i=0;i<=N;i++) p[i]=i; 
        ans=1e9; 
        do 
        { 
            t=crosstime(s.x*1.0,s.y*1.0,sv*1.0,a[p[0]].x,a[p[0]].y,v[p[0]].x,v[p[0]].y)+1.0; 
            for(i=1;i<N;i++)//当前在p[i-1]船上准备出发瞬间,p[i-1]船已经行驶t小时 ,p[i]船同样行驶t小时  
                t+=1.0+crosstime(a[p[i-1]].x*1.0+t*v[p[i-1]].x,a[p[i-1]].y*1.0+t*v[p[i-1]].y,sv*1.0,a[p[i]].x*1.0+t*v[p[i]].x,a[p[i]].y*1.0+t*v[p[i]].y,v[p[i]].x*1.0,v[p[i]].y*1.0); 
            t+=back(a[p[N-1]].x*1.0+t*v[p[N-1]].x,a[p[N-1]].y*1.0+t*v[p[N-1]].y,s.x,s.y,sv); 
            if(t<ans) ans=t; 
        } 
        while(next_permutation(p,p+N)); 
        trs(ans,T); 
        printf("Case %d: %d hour(s) %d minute(s) %d second(s)\n",js++,T[0],T[1],T[2]); 
    } 
    return 0; 

补充:软件开发 , C++ ,
CopyRight © 2012 站长网 编程知识问答 www.zzzyk.com All Rights Reserved
部份技术文章来自网络,