POJ2991 线段树 区间更新 计算几何

题意：吊车由n条不同长度的线段组成，每条线段首尾相接。初始状态：每条线段垂直与x轴。每次操作改变第s条和（s+1）条的逆时针角度为a，询问每次操作后第n段末尾的坐标。
思路：
1、将每段线段都视为向量，则每次询问的结果实质上是向量和。
2、每次改变第s段和第（s+1）段的相对角度，实际上是改变了从第（s+1）段至第n段的绝对角度，则可以通过线段树进行区间更新。然而，题目给定的是第s段和第（s+1）段之间的逆时针角度，并非顺时针旋转了的角度，如何解决？通过新的degree数组解决，degree[i]代表第i段和第（i-1）段的逆时针角度，维护这个数组。
3、向量旋转公式：
http://blog.csdn.net/hjq376247328/article/details/45113563
（若知道这个公式能省不少代码量）
反思：
1、抽象能力不足。刚开始没有将线段化为向量，则必须通过 上一段末尾坐标 和 本段的绝对角度 递归计算出答案。
2、数学能力不足。刚开始不知道向量旋转公式，也没有想要自己推导，直接在向量结构体里不仅定义了坐标，还定义了长度和角度，导致代码量倍增，最后也没能成功debug。
3、思维能力不足。还没写之前也考虑到了思路2中的问题，但竟然没能想出这种简单有效的方法，化未知为已知。
4、若用G++提交，记得double对应%f，否则会TLE。

虽然耗时很长，但收获了不少新知和技巧，以及最后A出来的满足感。继续加油！

#include <iostream>#include <cmath>#include <iomanip>#include <cstdio>#include <cstring>using namespace std;const double PI = acos(-1.0);const int MAXN = 10000 + 10;int n, c;double degree[MAXN];struct Vec{    double x, y;}vec[MAXN], sum[MAXN << 2];double lazy[MAXN << 2];void PushUp(int rt){    sum[rt].x = sum[rt << 1].x + sum[rt << 1 | 1].x;    sum[rt].y = sum[rt << 1].y + sum[rt << 1 | 1].y;    return;}void PushDown(int rt){    if(lazy[rt])    {        lazy[rt << 1]  += lazy[rt];        lazy[rt << 1 | 1]  += lazy[rt];        double a = lazy[rt], xl = sum[rt << 1].x, yl = sum[rt << 1].y, xr = sum[rt << 1 | 1].x, yr = sum[rt << 1 | 1].y;        sum[rt << 1].x = xl * cos(a) + yl * sin(a);        sum[rt << 1].y = yl * cos(a) - xl * sin(a);        sum[rt << 1 | 1].x = xr * cos(a) + yr * sin(a);        sum[rt << 1 | 1].y = yr * cos(a) - xr * sin(a);        lazy[rt] = 0;    }    return;}void Build(int l, int r, int rt){    if(l == r)    {        sum[rt] = vec[l];        return;    }    int m = (l + r) >> 1;    Build(l, m, rt << 1);    Build(m + 1, r, rt << 1 | 1);    PushUp(rt);    return;}void Update(int L, int R, double c, int l, int r, int rt){    if(L <= l && r <= R)    {        double x = sum[rt].x * cos(c) + sum[rt].y * sin(c);        double y = sum[rt].y * cos(c) - sum[rt].x * sin(c);        sum[rt].x = x; sum[rt].y = y;        lazy[rt] += c;        return;    }    PushDown(rt);    int m = (l + r) >> 1;    if(L <= m) Update(L, R, c, l, m, rt << 1);    if(R > m) Update(L, R, c, m + 1, r, rt << 1 | 1);    PushUp(rt);    return;}void init(){    memset(lazy, 0, sizeof(lazy));    for(int i = 2; i <= n; i++)    {        degree[i] = PI;    }    Build(1, n, 1);    return;}int main(){    int kase = 0;    while(~scanf("%d%d", &n, &c))    {        if(kase++) printf("\n");        for(int i = 1; i <= n; i++)        {            scanf("%lf", &vec[i].y);            vec[i].x = 0;        }        init();        while(c--)        {            int s, a; scanf("%d%d", &s, &a);            double rad = (a * 1.0 / 180.0) * PI;            Update(s + 1, n, degree[s + 1] - rad, 1, n, 1);            degree[s + 1] = rad;            printf("%.2f %.2f\n", sum[1].x, sum[1].y);        }    }    return 0;}