最近看到了线段树,对于线段树也是有了初步的了解,还是需要时间继续研究,加深理解。
感觉线段树,个人觉得最主要的是递归过程的理解。
这一题,给定一段绳子,“ 分成 ”N段,起初,每段绳子都是垂直的。然后有C个命令,每个命令包含两个数 i , j,i 是第几段绳子,j 是 i 段绳子旋转到 i + 1 段绳子多经过的角度。也就是 i 和 i + 1 之间的角度是 j。
可以理解成向量,起初每个向量方向向上。
如果结点 i 表示的向量是vx1, vy1。角度是angle,两个儿子节点是chl,chr,那么:
vx1 = vx(chl)+ (cos(angle)* vx(chr)- sin(angle)* vy(chr));
vy1 = vy(chl)+ (sin(angle)* vx(chr)+ cos(angle)* vy(chr));
下面是AC的代码:
#include <iostream> #include <cmath> #include <cstdio> using namespace std; const double PI = acos(-1.0); //PI,直接用3.1415927过不了。 const int MAX_SIZE = (1 << 15) - 1; //线段树的最大长度 const int MAX_N = 10005; //N的最大值 const int MAX_C = 10005; //C的最大值 int N, C; int L[MAX_N], S[MAX_C], A[MAX_C]; //L是每段长度,S是每个命令中的第几段的数组,A是角度数组 double vx[MAX_SIZE], vy[MAX_SIZE]; //各个向量的数组 double angle[MAX_SIZE]; //各个向量的角度数组 double temp[MAX_SIZE]; //角度变化而保存的当前角度的数组 //k是结点编号,l,r是区间[l,r] void init(int k, int l, int r) //线段树的初始化 { angle[k] = vx[k] = 0.0; //角度都为0,x坐标为0. if(r - l == 1) //叶子结点的y坐标 { vy[k] = L[l]; } else //非叶子结点 { int chl = k * 2 + 1; int chr = k * 2 + 2; init(chl, l, (l + r) / 2); init(chr, (l + r) / 2, r); vy[k] = vy[chl] + vy[chr]; } } //将s和s+1段的角度变成a //v是结点编号 void change(int s, double a, int v, int l, int r) { if(s <= l) return; else if(s < r) //s>l & s<r { int chl = 2 * v + 1; int chr = 2 * v + 2; int m = (l + r) / 2; change(s, a, chl, l, m); change(s, a, chr, m, r); if(s <= m) angle[v] += a; double s = sin(angle[v]), c = cos(angle[v]); vx[v] = vx[chl] + (c * vx[chr] - s * vy[chr]); //改变向量 vy[v] = vy[chl] + (s * vx[chr] + c * vy[chr]); } } void solve() { init(0, 0, N); //初始化 for(int i = 0; i < N; i++) //起始每个向量的角度为180度,也就是π temp[i] = PI; for(int j = 0; j < C; j++) { int s = S[j]; double a = A[j] / 360.0 * 2 * PI; change(s, a - temp[s], 0, 0, N); temp[s] = a; //角度改变 printf("%.2lf %.2lf\n\n", vx[0], vy[0]); //输出 } } int main() { while(scanf("%d%d", &N, &C) != EOF) //输入 { for(int i = 0; i < N; i++) scanf("%d", &L[i]); for(int j = 0; j < C; j++) scanf("%d%d", &S[j], &A[j]); solve(); } return 0; }
原文地址:http://blog.csdn.net/qq_25425023/article/details/45870741