/** * author: otera **/ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; //#define int long long typedef long long ll; typedef unsigned long long ul; typedef unsigned int ui; typedef long double ld; const int inf=1e9+7; const ll INF=1LL<<60 ; const ll mod=1e9+7 ; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef complex Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair P; typedef pair LDP; typedef pair LP; #define fr first #define sc second #define all(c) c.begin(),c.end() #define pb push_back #define debug(x) cerr << #x << " = " << (x) << endl; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } template struct LazySegmentTree{ private: int n; vector node; vector lazy; T E0; S E1; inline void updatef(S& lazy,S& value){ lazy.fr += value.fr; if(value.sc > 0) lazy.sc = value.sc; //lazy += value; //lazy = max(lazy, value); //lazy = min(lazy, value); } inline void calculatef(T& node,S& lazy,int len){ //node += lazy * len; //区間sumはこっち //node += lazy ; //区間maxとか if(lazy.sc > 0) { ld now = sqrt(node.fr * node.fr + node.sc * node.sc); node.fr *= lazy.sc / now; node.sc *= lazy.sc / now; } if(lazy.fr > eps) { ld nxtx = node.fr * cos(2.0 * pi * lazy.fr / 360.0) - node.sc * sin(2.0 * pi * lazy.fr / 360.0); ld nxty = node.fr * sin(2.0 * pi * lazy.fr / 360.0) + node.sc * cos(2.0 * pi * lazy.fr / 360.0); node.fr = nxtx, node.sc = nxty; } } inline T queryf(T& x,T& y){ //return x + y; //return x * y; //return max(x, y); return P{x.fr + y.fr, x.sc + y.sc}; } public: LazySegmentTree(int sz,T nodeE,S lazyE ):E0(nodeE), E1(lazyE){ n=1; while(n& v,T E0,S E1 ):E0(E0),E1(E1){ n=1; int sz=v.size(); while(n=0; --i){ node[i] = queryf(node[2*i+1],node[2*i+2]); } } void eval(int k,int l,int r){ if(lazy[k]==E1)return ; calculatef(node[k], lazy[k], r-l); if(r-l>1){ updatef(lazy[2*k+1], lazy[k]); updatef(lazy[2*k+2], lazy[k]); } lazy[k]=E1; } void update(int a, int b, S x,int k=0,int l=0,int r=-1){ if(r<0)r=n; eval(k,l,r); if(r<=a||b<=l)return; if(a<=l&&r<=b){ updatef(lazy[k], x); eval(k,l,r); } else { update(a,b,x,2*k+1,l,(l+r)/2); update(a,b,x,2*k+2,(l+r)/2,r); node[k]=queryf(node[2*k+1], node[2*k+2]); } } T query(int a,int b,int k=0,int l=0,int r=-1){ if(r<0)r=n; eval(k,l,r); if(r<=a||b<=l)return E0; if(a<=l&&r<=b)return node[k]; T xl=query(a,b,2*k+1,l,(l+r)/2); T xr=query(a,b,2*k+2,(l+r)/2,r); return queryf(xl, xr); } }; void solve() { int n, q; cin >> n >> q; vector> a(n, P{1.0, 0.0}); LazySegmentTree, pair> seg(a, P{1.0, 0.0}, P{0.0, -1}); rep(_, q) { int t; cin >> t; int i; cin >> i; if(t == 0) { int x; cin >> x; -- i; seg.update(i, n, P{(ld)x ,-1}); } else if(t == 1) { int x; cin >> x; -- i; seg.update(i, i + 1, P{0.0, x}); } else { auto res = seg.query(0, i + 1); cout << res.fr << " " << res.sc << endl; } } } signed main() { ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(10); //int t; cin >> t; rep(i, t)solve(); solve(); return 0; }