POJ 2991 Crane (선분 수 + 계산 기하학)

기중기 한 대가 설명 되 어 있다.우 리 는 기 중 기 를 N 개의 선분 이 순서대로 머리 와 꼬리 를 연결 하여 만 든 것 으로 간주한다.제 i 라인 의 길 이 는 Li 입 니 다.처음에는 모든 라인 이 똑바로 연결 되 어 위 를 가 리 켰 다.기 존 C 조 크레인 조종 명령명령 i 는 두 개의 정수 Si 와 Ai 를 주 고 효 과 는 선분 Si 와 Si + 1 사이 의 각 도 를 Ao 도로 바 꾸 는 것 이다.그 중에서 각 도 는 선분 인 Si 부터 시계 반대 방향 으로 Si + 1 이 지나 가 는 각 도 를 말한다.처음에는 모든 각도 가 180 도 였 다.이 C 조 명령 을 순서대로 집행 하 다.모든 명령 이 실 행 된 후에 크레인 의 앞부분 (N 라인 의 점) 의 좌 표를 출력 한다.크레인 의 지점 좌표 가 (0, 0) Input 여러 조 의 사례 라 고 가정 하면 각 조 의 사례 는 첫 번 째 행위 두 개의 정수 N 과 C 는 크레인 의 선분 수 와 명령 수 를 나타 내 고 두 번 째 행위 N 개의 정 수 는 각 선분 의 길 이 를 나타 낸다. 그 다음 에 C 행 은 각 줄 의 두 개의 정수 로 명령 을 표시 하고 파일 의 끝 에 Output 을 입력 하여 각 조 의 사례 에 대해출력 이 명령 을 실행 한 후 크레인 앞부분 의 좌 표를 빈 줄 로 두 그룹 으로 나 누 어 출력 Sample Input 2 1 10 5 1 90 3 2 5 5 5 1 270 2 90 Sample Output 5.00 10.00
- 10.00 5.00 - 5.00 10.00 Solution 본 문 제 는 선분 트 리 로 해결한다. 각 노드 는 연속 적 인 선분 의 집합 을 표시 하고 아래 의 두 값 1 을 유지한다. 대응 하 는 선분 집합 중의 첫 번 째 선분 을 수직 방향 으로 전환 한 후 첫 번 째 선분 의 기점 에서 마지막 선분 의 종점 을 가리 키 는 벡터 2 (이 노드 에 잎 노드 가 있다 면)두 아들 노드 가 대응 하 는 부분 을 연결 한 후에 오른쪽 아들 이 돌려 야 하 는 각 도 는 바로 노드 i 가 표시 하 는 벡터 가 Xi, Yi 이 고 각 도 는 angi 이 며 두 아들 노드 가 chl 과 chr 이면 Xi = Xchl + (cos (angi) * Xchr - sin (angi) * Ychr) Yi = Ychl + (sin (angi) * Xchr + cos (angi) * Ychr) 가 있 습 니 다. 그러면 매번 업데이트 가 O (logn) 시간 안에 완 성 됩 니 다.출력 값 은 루트 노드 가 벡터 에 대응 하 는 값 Code 입 니 다.

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
using namespace std;
#define maxn 10001
const int ST_SIZE=(1<<15)-1;
int N,C;
int L[maxn];
int S[maxn],A[maxn];
//          
double vx[ST_SIZE],vy[ST_SIZE];//       
double ang[ST_SIZE];//       
double prv[maxn];//                    
void init(int k,int l,int r)//      ,k     ,l,r          [l,r)   
{
    ang[k]=vx[k]=0.0;//    x       0 
    if(r-1==l)//     
        vy[k]=L[l];
    else//      
    {
        int chl=2*k+1,chr=2*k+2;
        init(chl,l,(l+r)/2);
        init(chr,(l+r)/2,r);
        vy[k]=vy[chl]+vy[chr];
    }
}
void change(int s,double a,int v,int l,int r)// s s+1     a,v     ,l,r          [l,r)   
{
    if(s<=l)//s      
        return ;
    else if(s<r)
    {
        int chl=v*2+1,chr=v*2+2;
        int m=(l+r)/2;
        change(s,a,chl,l,m);
        change(s,a,chr,m,r);
        if(s<=m)
            ang[v]+=a;
        double s=sin(ang[v]),c=cos(ang[v]);
        vx[v]=vx[chl]+(c*vx[chr]-s*vy[chr]);
        vy[v]=vy[chl]+(s*vx[chr]+c*vy[chr]);
    }
}
int main()
{
    int res=0;
    const double PI=acos((double)-1);
    while(scanf("%d%d",&N,&C)!=EOF)
    {
        res++;
        if(res!=1)//             
            printf("
");
        memset(vx,0,sizeof(vx));//    
        memset(vy,0,sizeof(vy));
        memset(prv,0,sizeof(prv));
        memset(ang,0,sizeof(ang));
        for(int i=0;i<N;i++)
            scanf("%d",&L[i]);
        for(int i=0;i<C;i++)
            scanf("%d%d",&S[i],&A[i]);
        init(0,0,N);//    
        for(int i=1;i<N;i++)
            prv[i]=PI;
        for(int i=0;i<C;i++)//     
        {
            int s=S[i];
            double a=A[i]/360.0*2*PI;//         
            change(s,a-prv[s],0,0,N);
            prv[s]=a;
            printf("%.2lf %.2lf
",vx[0],vy[0]);
        }
    }
    return 0;
}