[线性规划 对偶 凸包 三分] Codeforces 605C #335 (Div. 1) C. Freelancer's Dreams

很显然的线性规划
对偶一下

maximum px+qyaix+biy<=1

这个肯定是半平面交出一个凸包然后用直线去切
直接三分就好了
有点需要精度

//ai*x+bi*y<=1 max p*x+q*y
#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long double ld;

const int N=100005;

int n; ld p,q;
ld a[N],b[N];

inline ld F(ld x){
  ld y=1e130;
  for (int i=1;i<=n;i++)
    y=min(y,(1.0-a[i]*x)/b[i]);
  return p*x+q*y;
}

int main(){
  double _p,_q,_a,_b;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d%lf%lf",&n,&_p,&_q),p=_p,q=_q;
  for (int i=1;i<=n;i++) scanf("%lf%lf",&_a,&_b),a[i]=_a,b[i]=_b;
  ld L=0,R=1e130,m1,m2;
  for (int i=1;i<=n;i++) R=min(R,1.0/a[i]);
  while (R-L>1e-15)
    if (F(m1=L+(R-L)/3.0)<F(m2=R-(R-L)/3.0))
      L=m1;
    else
      R=m2;
  printf("%.15lf\n",(double)F((R+L)/2));
  return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值