比赛 |
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;
}