Bzoj1857:[Scoi2010]传送带:三分

最新推荐文章于 2021-08-08 16:31:38 发布

原创最新推荐文章于 2021-08-08 16:31:38 发布 · 547 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#BZOJ #三分法 #OI

OI 同时被 2 个专栏收录

108 篇文章

订阅专栏

数值方法

3 篇文章

订阅专栏

本文介绍了解决Scoi2010中的传送带问题的方法，通过验证函数的下凸性，并使用三分法进行优化求解。代码中详细展示了如何计算两点之间的距离，并逐步逼近最优解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接[Scoi2010]传送带

首先由猜测法证得函数具有下凸性QwQ

然后就可以三分辣

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const double eps=1e-5;
struct point{double x,y;}A,B,C,D;
double p,q,r,ax,ay,bx,by,cx,cy,dx,dy;

double powe(double x){return x*x;}

double calc(point a,point b){
	return sqrt(powe(a.x-b.x)+powe(a.y-b.y));
}

double calc2(double x,double y,point a){
	return sqrt(powe(x-a.x)+powe(y-a.y));
}

double getval(point tmp){
	C.x=cx; C.y=cy; D.x=dx; D.y=dy;
	while (abs(C.x-D.x)>eps||abs(C.y-D.y)>eps){
		point tmpl; tmpl.x=C.x+(D.x-C.x)/3; tmpl.y=C.y+(-C.y+D.y)/3;
		point tmpr; tmpr.x=C.x+(D.x-C.x)/3*2; tmpr.y=C.y+(-C.y+D.y)/3*2;
		double ret1=calc2(ax,ay,tmp)/p+calc(tmp,tmpl)/r+calc2(dx,dy,tmpl)/q;
		double ret2=calc2(ax,ay,tmp)/p+calc(tmp,tmpr)/r+calc2(dx,dy,tmpr)/q;
		if (ret1>ret2) C=tmpl; else D=tmpr;
	}
	return calc2(ax,ay,tmp)/p+calc(tmp,C)/r+calc2(dx,dy,C)/q;
}

int main(){
	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);
	ax=A.x; ay=A.y; bx=B.x; by=B.y;
	cx=C.x; cy=C.y; dx=D.x; dy=D.y;
	scanf("%lf%lf%lf",&p,&q,&r);
	while (abs(A.x-B.x)>eps||abs(A.y-B.y)>eps){
		point tmpl; tmpl.x=A.x+(B.x-A.x)/3; tmpl.y=A.y+(B.y-A.y)/3;
		point tmpr; tmpr.x=A.x+(B.x-A.x)/3*2; tmpr.y=A.y+(B.y-A.y)/3*2;
		double ret1=getval(tmpl),ret2=getval(tmpr);
		if(ret1>ret2) A=tmpl; else B=tmpr;
	}
	printf("%.2lf\n",getval(A));
}