bzoj2765 [JLOI2010]铁人双项比赛

题目大意:
铁人双项比赛由长跑和骑自行车组成。现在给定总赛程s,以及每个选手长跑和骑车的平均速度,请你求出对于某个指定的选手最有利的k和r。所谓最有利,是指选择了这个k和r后,该选手可以获得冠军,且领先第2名尽量地多。

分析:
很早就知道这道题的解法是半平面交
实际上我们可以很容易的得到一个方程:
这里写图片描述
我们可以看作是以k为x轴,以t为y轴的一条直线
如果我们对于1..n-1这些直线求交点,那么最优解一定出现在交点中。
那么相当于用n所在的直线去卡这些交点求最优值

于是我就开始码。。。
但是总感觉有点麻烦,而且好像有点不对劲。。。

然而发现网上的代码,码量都是出人意料的少
实际上:数据范围比较小直接暴力就可以了
感觉貌似可以求出交点在再求凸包,然后直接判断凸包上的点
好像也可以用半平面交直接求交点,然后判断

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>

using namespace std;

const double INF=1e13;
const int N=100010;
double a[N],b[N],s,pos,ans;
int n;
struct node{
    double x,y;
    node (double xx=0,double yy=0)
    {
        x=xx; y=yy;
    }
};
node p[N],q[N],po[N];

node operator +(const node &a,const node &b){return node(a.x+b.x,a.y+b.y);}
node operator -(const node &a,const node &b){return node(a.x-b.x,a.y-b.y);}
node operator *(const node &a,const double &b){return node(a.x*b,a.y*b);}
node operator /(const node &a,const double &b){return node(a.x/b,a.y/b);}

double Dot(node a,node b){return a.x*b.x+a.y*b.y;}
double Len(node a){return sqrt(Dot(a,a));}
double Cross(node a,node b){return a.x*b.y-a.y*b.x;}

void solve(double x)
{
    if (x<0 || x>s) return;
    double ti=INF;
    for (int i=1;i<n;i++)
        ti=min(ti,a[i]*x+b[i]);
    ti-=a[n]*x+b[n];
    if (ti>ans) ans=ti,pos=x;
}

int main()
{
    scanf("%lf%d",&s,&n);
    for (int i=1;i<=n;i++)
    {
        double x,y;
        scanf("%lf%lf",&x,&y);
        a[i]=1.0/x-1.0/y;
        b[i]=s/y;
    }
    ans=-INF;   
    for (int i=1;i<n;i++)
        for (int j=i+1;j<n;j++)
            solve((b[j]-b[i])/(a[i]-a[j]));    //两直线交点 
    solve(0.0);
    solve(s);
    if (ans>=0) printf("%.2lf %.2lf %.0lf\n",pos,s-pos,ans*3600);
    else printf("NO\n");
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值