Crane POJ - 2991（线段树 + 计算几何）

The sunshine。

于 2021-04-22 13:06:08 发布

阅读量111

点赞数

分类专栏：数据结构

本文链接：https://blog.csdn.net/m0_46082460/article/details/116009923

版权

数据结构专栏收录该内容

13 篇文章 0 订阅

订阅专栏

博客详细介绍了如何利用线段树和计算几何的方法解决Crane POJ 2991题目。内容包括题意解析，解题思路的阐述，以及具体代码实现过程。通过对不同操作指令处理，计算起重机前端线段端点坐标。

摘要由CSDN通过智能技术生成

题意：

有一台起重机，我们把起重机看成由N条线段构依次首尾相接而成，第i条线段的长度为Li,最开始，所有线段都笔直朝上。

有C条操作指令，指令i给出两个整数Si和Ai,效果是使线段Si和Si+1之间的角度变成Ai度，其中角度指的是从线段Si开始沿逆时针方向旋转到Si+1所经过的角度，最开始时所有角度都是180度

按顺序执行C条指令，在每条指令执行之后，输出起重机的前端即第N根线段的端点坐标，假设起重机的支点的坐标为（0,0）。

思路：

线段树 + 计算几何

代码实现：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn = 2e5 + 5;
const double pi = acos(-1.0);

//1.数组开的太大会RE
//2.多组+memset导致超时

struct node{
	double x,y;
}stree[maxn << 2];
double cont[maxn<<2];

void pushup(int rt){//向量相加
	stree[rt].x = stree[rt << 1].x + stree[rt << 1|1].x;
	stree[rt].y = stree[rt << 1].y + stree[rt << 1|1].y;
}

void build(int l,int r,int rt){//建树
	if(l == r){
		double y;
		scanf("%lf",&y);
		stree[rt].x = 0;
		stree[rt].y = y;
		return ;
	}
	int mid = (l + r) >> 1;
	build(l,mid,rt<<1);
	build(mid + 1,r,rt<<1|1);
	pushup(rt);
}

void f(int rt,double st){
	double x = stree[rt].x;
	double y = stree[rt].y;
	stree[rt].x = x * cos(st) - y * sin(st);
	stree[rt].y = y * cos(st) + x * sin(st);
}

void pushdown(int rt){
	if(cont[rt] != 0){
		cont[rt << 1] += cont[rt];
		cont[rt << 1|1] += cont[rt];
		f(rt << 1,cont[rt]);
		f(rt << 1|1,cont[rt]);
		cont[rt] = 0;
	}
}

void update(int L,int R,double st,int l,int r,int rt){
	if(l >= L&&r <= R){
		cont[rt] += st;
		f(rt,st);
		return ;
	}
	pushdown(rt);
	int mid = (l + r) >> 1;
	if(L <= mid) update(L,R,st,l,mid,rt<<1);
	if(R > mid) update(L,R,st,mid + 1,r,rt<<1|1);
	pushup(rt);
}

int main(){
	int n,c;
	int kace = 0;
	while(~scanf("%d%d",&n,&c)){
		if(kace++) puts("");
		for(int i = 1;i <= n << 2;i++) cont[i] = 0;
		build(1,n,1);
		int ang[maxn << 2];
		for(int i = 1;i <= n;i++) ang[i] = 180;
		while(c--){
			int a,b;
			scanf("%d%d",&a,&b);
			int st1 = (b - ang[a] + 360) % 360;
			update(a + 1,n,st1*1.0/180.0*pi,1,n,1);
			ang[a] = b;
			printf("%.2lf %.2lf\n",stree[1].x,stree[1].y);
		}
	}
	return 0;
}