[HNOI2008]水平可见直线 单调栈

题目描述:
在xoy直角坐标平面上有n条直线L1,L2,...Ln,若在y值为正无穷大处往下看,能见到Li的某个子线段,则称Li为
可见的,否则Li为被覆盖的.
例如,对于直线:
L1:y=x; L2:y=-x; L3:y=0
则L1和L2是可见的,L3是被覆盖的.
给出n条直线,表示成y=Ax+B的形式(|A|,|B|<=500000),且n条直线两两不重合.求出所有可见的直线.


题解:

一道很好的思维题。
1.简单手画一下,能被看到的直线应该是所有直线一起围成的大凸包。
2.由于是凸包,我们考虑将所有直线按斜率排序,从小到大依次加入到平面直角坐标系中。
3.我们考虑一下新加入直线的情况:

在这种情况中,我们可以看到新加入的红色直线与加入之前平面中斜率第二大的直线的交点位于先前第一大与第二大之左,显然,这就会挡住平面中斜率第二大的直线,我们就将该直线弹出,直到找到一个交点在新加入直线的交点左侧。
对于整个过程,直线的斜率单调递增,交点横坐标也单调递增,直接用单调栈维护即可。
时间复杂度为 $O(n)$

 

Code:

#include<cstdio>
#include<algorithm>
#include<string>
using namespace std;
void setIO(string a){
	freopen((a+".in").c_str(),"r",stdin);
}
const int maxn=100000+5;
struct Line{
	double slope, y;
}line[maxn];
int arr[maxn],ans[maxn],S[maxn],top;
bool cmp(int i,int j){
	if(line[i].slope==line[j].slope) return line[i].y>line[j].y;
	return line[i].slope<line[j].slope;
}
double get(int i,int j){
	return (line[i].y-line[j].y)/(line[j].slope-line[i].slope);
}
int main(){
	//setIO("input");
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;++i) {
		scanf("%lf%lf",&line[i].slope,&line[i].y);
		arr[i]=i;
	}
	sort(arr+1,arr+1+n,cmp);
    for(int i=1;i<=n;++i)
    {
    	int cur=arr[i];
  
    	if(line[cur].slope==line[arr[i-1]].slope && i!=1) continue;
    	while(top>1 && get(S[top],S[top-1])>=get(arr[i],S[top])) --top;
    	S[++top]=cur;
    	ans[top]=cur;
    }
    sort(ans+1,ans+1+top);
    for(int i=1;i<=top;++i) printf("%d ",ans[i]);
    return 0;
}

  

 

根据引用\[1\]和引用\[2\]的描述,题目中的影魔拥有n个灵魂,每个灵魂有一个战斗力ki。对于任意一对灵魂对i,j (i<j),如果不存在ks (i<s<j)大于ki或者kj,则会为影魔提供p1的攻击力。另一种情况是,如果存在一个位置k,满足ki<c<kj或者kj<c<ki,则会为影魔提供p2的攻击力。其他情况下的灵魂对不会为影魔提供攻击力。 根据引用\[3\]的描述,我们可以从左到右进行枚举。对于情况1,当扫到r\[i\]时,更新l\[i\]的贡献。对于情况2.1,当扫到l\[i\]时,更新区间\[i+1,r\[i\]-1\]的贡献。对于情况2.2,当扫到r\[i\]时,更新区间\[l\[i\]+1,i-1\]的贡献。 因此,对于给定的区间\[l,r\],我们可以根据上述方法计算出区间内所有下标二元组i,j (l<=i<j<=r)的贡献之和。 #### 引用[.reference_title] - *1* *3* [P3722 [AH2017/HNOI2017]影魔(树状数组)](https://blog.csdn.net/li_wen_zhuo/article/details/115446022)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [洛谷3722 AH2017/HNOI2017 影魔 线段树 单调栈](https://blog.csdn.net/forever_shi/article/details/119649910)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值