poj 2991 Crane 线段树lazy_tag

//poj 2991
//sep9
#include <iostream>
#include <cmath>
using namespace std;
const int MAXN=100024;
const double PI=acos(-1.0);
double sumx[MAXN*4],sumy[MAXN*4];
int angle[MAXN],alpha[MAXN*4];

void push_up(int k)
{
	sumx[k]=sumx[2*k]+sumx[2*k+1];
	sumy[k]=sumy[2*k]+sumy[2*k+1];	
}

void build(int l,int r,int k)
{
	alpha[k]=0;
	if(l==r){
		scanf("%lf",&sumy[k]);
		sumx[k]=0;
		return ;
	}
	int mid=(l+r)/2;
	build(l,mid,2*k);
	build(mid+1,r,2*k+1);
	push_up(k);
}

void rotate(int k,int d)
{
	double x,y,dd=d*PI/180;
	x=sumx[k]*cos(dd)-sumy[k]*sin(dd);
	y=sumx[k]*sin(dd)+sumy[k]*cos(dd);
	sumx[k]=x;
	sumy[k]=y;
}

void push_down(int k)
{
	if(alpha[k]){
		alpha[2*k]+=alpha[k];
		alpha[2*k+1]+=alpha[k];
		rotate(2*k,alpha[k]);
		rotate(2*k+1,alpha[k]);
		alpha[k]=0;
	}	
}

void update(int ql,int qr,int d,int l,int r,int k)
{
	if(ql==l&&qr==r){
		rotate(k,d);
		alpha[k]+=d;
		return ;
	}
	int mid=(r+l)/2;
	push_down(k);//lazy tag
	if(qr<=mid) update(ql,qr,d,l,mid,2*k);
	else if(ql>=mid+1) update(ql,qr,d,mid+1,r,2*k+1);
	else{
		update(ql,mid,d,l,mid,2*k);
		update(mid+1,qr,d,mid+1,r,2*k+1);
	}
	push_up(k);
}

int main()
{
	int n,q;
	while(scanf("%d%d",&n,&q)==2){
		build(1,n,1);				
		for(int i=1;i<=n;++i)
			angle[i]=180;
		while(q--){
			int s,d;
			scanf("%d%d",&s,&d);
			update(s+1,n,d-angle[s],1,n,1);
			angle[s]=d;
			printf("%.2lf %.2lf\n",sumx[1],sumy[1]);
		}
	}
	return 0;	
}