[BZOJ4614][极角排序]WF2016:Oil

BZOJ4614

最优解必定过两个线段的端点
枚举其中一个,然后把其他所有端点相对枚举点极角排序,前缀和扫一遍就完了

Code:

#include<bits/stdc++.h>
#define ll long long
#define db double
#define eps 1e-8
using namespace std;
inline int read(){
	int res=0,f=1;char ch=getchar();
	while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
	while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
	return res*f;
}
const int N=2005;
int n,ans,tot;
struct line{int x1,x2,y,len;}a[N];
struct slp{db k;int kd;}b[N<<1];
bool operator<(slp a,slp b){return a.k<b.k;}
inline void solve(db x,db y){
	int cnt=0;
	for(int i=1;i<=n;i++){
		if(a[i].y==y) continue;
		b[++cnt].k=(a[i].x1-x)/(a[i].y-y);
		b[++cnt].k=(a[i].x2-x)/(a[i].y-y);
		if(b[cnt].k>b[cnt-1].k){
			b[cnt-1].kd=a[i].len;
			b[cnt].kd=-a[i].len;
			b[cnt].k+=eps;
		}
		else{
			b[cnt].kd=a[i].len;
			b[cnt-1].kd=-a[i].len;
			b[cnt-1].k+=eps;
		}
	}
	sort(b+1,b+cnt+1);
	for(int i=1;i<=cnt;i++){tot+=b[i].kd;ans=max(ans,tot);}
}
int main(){
	n=read();ans=0;
	for(int i=1;i<=n;i++){
		a[i].x1=read();a[i].x2=read();a[i].y=read();
		if(a[i].x1>a[i].x2) swap(a[i].x1,a[i].x2);
		a[i].len=a[i].x2-a[i].x1;
	}
	for(int i=1;i<=n;i++){
		tot=a[i].len;solve(a[i].x1,a[i].y);
		tot=a[i].len;solve(a[i].x2,a[i].y);
		ans=max(ans,tot);
	}
	cout<<ans;
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值