POJ 2991 Crane

线段树+计算几何，区间更新，区间求和，向量旋转。

/* ***********************************************
Author        :Zhou Zhentao
Email         :774388357@qq.com
Created Time  :2015/11/27 9:58:30
File Name     :main.cpp
************************************************ */

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;

const int maxn=+;
const double pi= acos(-1.0);
struct SegTree
{
    double SumX,SumY;//记录向量
    int add;//延迟更新的标记，记录角度转变量
}st[maxn*];
int f[maxn];
int N,Q;

double H(double x)
{
    return x*pi/180.0;
}

void Rot(int rt,double A)
{
    double NewX,NewY;
    NewX=st[rt].SumX*cos(H(A))-st[rt].SumY*sin(H(A));
    NewY=st[rt].SumX*sin(H(A))+st[rt].SumY*cos(H(A));
    st[rt].SumX=NewX;
    st[rt].SumY=NewY;
}

void pushUp(int rt)
{
    st[rt].SumX=st[rt*].SumX+st[rt*+].SumX;
    st[rt].SumY=st[rt*].SumY+st[rt*+].SumY;
}

void pushDown(int rt)
{
    if(st[rt].add)
    {
        st[*rt].add+=st[rt].add;
        st[*rt+].add+=st[rt].add;
        Rot(*rt,(double)st[rt].add);
        Rot(*rt+,(double)st[rt].add);
        st[rt].add=;
    }
}

void build(int l,int r,int rt)
{
    st[rt].add=;
    if(l==r)
    {
        scanf("%lf",&st[rt].SumY);
        st[rt].SumX=;
        return;
    }
    int m=(l+r)/;
    build(l,m,*rt);
    build(m+,r,*rt+);
    pushUp(rt);
}

void update(int L,int R,int add,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        st[rt].add+=add;
        Rot(rt,(double)add);
        return;
    }

    pushDown(rt);
    int m=(l+r)/;
    if(L<=m) update(L,R,add,l,m,*rt);
    if(R>=m+) update(L,R,add,m+,r,*rt+);
    pushUp(rt);
}

int main()
{

    int tot=;
    while(~scanf("%d%d",&N,&Q))
    {
        if(tot) printf("\n");tot=;

        build(,N,);

        for(int i=;i<=N;i++) f[i]=;

        for(int i=;i<=Q;i++)
        {
            int x;int y;
            scanf("%d%d",&x,&y);
            int A=y-f[x];
            update(x+,N,A,,N,);
            f[x]=y;

            printf("%.2lf %.2lf\n",st[].SumX,st[].SumY);
        }
    }
    return ;
}