#include #include #define rep(i,b) for(int i=0;i=0;i--) #define rep1(i,b) for(int i=1;i using mpq = priority_queue, greater>; template bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (b ll sumv(const vector&a){ll res(0);for(auto&&x:a)res+=x;return res;} bool yn(bool a) { if(a) {cout << "Yes" << endl; return 1;} else {cout << "No" << endl; return 0;}} #define dame { cout << "No" << endl; return 0;} #define dame1 { cout << -1 << endl; return 0;} #define test(x) cout << "test" << x << endl; #define deb(x,y) cout << x << " " << y << endl; #define out cout << ans << endl; #define outv fore(yans , ans) cout << yans << endl; #define show(x) cerr<<#x<<" = "<; using pil = pair; using pli = pair; using pii = pair; using tp = tuple; using vi = vector; using vl = vector; using vs = vector; using vb = vector; using vpii = vector; using vpli = vector; using vpll = vector; using vpil = vector; using vvi = vector>; using vvl = vector>; using vvs = vector>; using vvb = vector>; using vvpii = vector>; using vvpli = vector>; using vvpll = vector; using vvpil = vector; using mint = modint998244353; using vm = vector; using vvm = vector>; vector dx={1,0,-1,0,1,1,-1,-1},dy={0,1,0,-1,1,-1,1,-1}; ll gcd(ll a, ll b) { return b?gcd(b,a%b):a;} ll lcm(ll a, ll b) { return a/gcd(a,b)*b;} const double eps = 1e-10; const ll LINF = 1001002003004005006ll; const int INF = 1001001001; struct mat{ vector> v; mat(long double a=1,long double b=0,long double c=0,long double d=0,long double e=1,long double f=0,long double g=0 ,long double h=0,long double i=1) { v.resize(3,vector(3)); v[0][0] = a; v[0][1] = b; v[0][2] = c; v[1][0] = d; v[1][1] = e; v[1][2] = f; v[2][0] = g; v[2][1] = h; v[2][2] = i; } mat operator*(const mat& ot) { mat ret; rep(i,3) rep(j,3){ ret.v[i][j] = 0; rep(k,3) ret.v[i][j] += v[i][k]*ot.v[k][j]; } return ret; } }; int main(){ int n; cin>>n; vector a(n),ar(n); vector sum(n),sumr(n); rep(i,n){ long double p,q,r; cin>>p>>q>>r; long double theta = (long double)r/180*M_PI; long double co=cosl(theta),si=sinl(theta); show(i); // cout << fixed << setprecision(8) << theta << endl; a[i].v[0][0] = co; a[i].v[0][1] = -si; a[i].v[0][2] = (-co+1)*p + si*q; a[i].v[1][0] = si; a[i].v[1][1] = co; a[i].v[1][2] = -si*p + (-co+1)*q; a[i].v[2][0] = 0; a[i].v[2][1] = 0; a[i].v[2][2] = 1; theta *= -1; co=cosl(theta),si=sinl(theta); ar[i].v[0][0] = co; ar[i].v[0][1] = -si; ar[i].v[0][2] = (-co+1)*p + si*q; ar[i].v[1][0] = si; ar[i].v[1][1] = co; ar[i].v[1][2] = -si*p + (-co+1)*q; ar[i].v[2][0] = 0; ar[i].v[2][1] = 0; ar[i].v[2][2] = 1; if (i == 0){ sum[i] = a[i]*sum[i]; sumr[i] = sumr[i]*ar[i]; }else{ sum[i] = a[i]*sum[i-1]; sumr[i] = sumr[i-1]*ar[i]; } // show(i); // cout << fixed << setprecision(8) << a[i].v[0][0] << endl; } int q; cin>>q; vector> ans(q); rep(i,q){ int s,t; cin>>s>>t; s--; t--; long double x,y; cin>>x>>y; mat v = sum[t]; if (s != 0) v = v*sumr[s-1]; ans[i].fi = v.v[0][0]*x + v.v[0][1]*y + v.v[0][2]; ans[i].se = v.v[1][0]*x + v.v[1][1]*y + v.v[1][2]; } rep(i,q){ cout << fixed << setprecision(20) << ans[i].fi << " " << ans[i].se << endl; } return 0; }