比赛 Asm.Def战记之圣地亚哥“杯2015 评测结果 AAAAAAAAAA
题目名称 Asm.Def的命令 最终得分 100
用户昵称 梦那边的美好ET 运行时间 4.545 s
代码语言 C++ 内存使用 69.85 MiB
提交时间 2019-10-23 19:12:55
显示代码纯文本
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#define maxn 100010
using namespace std;
const double pi=acos(-1.0);
int n,m,op,a1,a2;
double a3,a4;
struct hs{double s[4][4];}a[maxn*4],fz,ans;
void cf(hs &A,hs s1,hs s2){
	for(int i=1;i<=3;i++)
	    for(int j=1;j<=3;j++)
		    fz.s[i][j]=(double)0.0;	
	for(int i=1;i<=3;i++)
	    for(int j=1;j<=3;j++)
	        for(int k=1;k<=3;k++)
	            fz.s[i][j]+=s1.s[i][k]*s2.s[k][j];
	for(int i=1;i<=3;i++)
	    for(int j=1;j<=3;j++)
		    A.s[i][j]=fz.s[i][j];	            
	return;
}
void maketree(int p,int l,int r){
	for(int i=1;i<=3;i++)a[p].s[i][i]=(double)1.0;
	if(l==r)return;int mid=(l+r)>>1;
	maketree(p*2,l,mid);maketree(p*2+1,mid+1,r);
	return;
}
void xc(int p){
	cf(a[p*2],a[p],a[p*2]);cf(a[p*2+1],a[p],a[p*2+1]);
	for(int i=1;i<=3;i++)for(int j=1;j<=3;j++)a[p].s[i][j]=(double)0.0;
	for(int i=1;i<=3;i++)a[p].s[i][i]=(double)1.0;
	return;
}
void insert(int p,int l,int r,int l1,int r1){
	if(l==l1&&r==r1){cf(a[p],ans,a[p]);return;}xc(p);
	int mid=(l+r)>>1;
	if(r1<=mid)insert(p*2,l,mid,l1,r1);
	if(mid<l1)insert(p*2+1,mid+1,r,l1,r1);
	if(l1<=mid&&r1>mid)insert(p*2,l,mid,l1,mid),insert(p*2+1,mid+1,r,mid+1,r1);
	return;
}
hs getans(int p,int l,int r,int x){
	if(l==r)return a[p];xc(p);
	int mid=(l+r)>>1;
	return (x<=mid)?getans(p*2,l,mid,x):getans(p*2+1,mid+1,r,x);
}
int main(){
	freopen("asm_command.in","r",stdin);
	freopen("asm_command.out","w",stdout);
    scanf("%d%d",&n,&m);
    maketree(1,1,n);
    for(int t=1;t<=m;t++){
    	scanf("%d",&op);
    	for(int i=1;i<=3;i++)for(int j=1;j<=3;j++)ans.s[i][j]=(double)0.0;
    	if(op==0){
    		scanf("%d",&a1);
    		ans=getans(1,1,n,a1);
    		printf("%.1lf ",(double)a1*ans.s[1][2]+ans.s[1][3]);
    		printf("%.1lf\n",(double)a1*ans.s[2][2]+ans.s[2][3]);
    	}
    	if(op==1){
    		scanf("%d%d%lf%lf",&a1,&a2,&a3,&a4);
    		ans.s[1][1]=ans.s[3][3]=ans.s[2][2]=(double)1.0;ans.s[1][3]=a3;ans.s[2][3]=a4;
    		insert(1,1,n,a1,a2);
    	}
    	if(op==2){
    		scanf("%d%d%lf",&a1,&a2,&a3);a3=a3/(double)180.0*pi;
    		ans.s[1][1]=ans.s[2][2]=cos(a3);ans.s[1][2]=-sin(a3);ans.s[2][1]=sin(a3);ans.s[3][3]=(double)1.0;
			insert(1,1,n,a1,a2);
    	}
    }
	return 0;
}