#define _USE_MATH_DEFINES #include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 1000000007; // constexpr int MOD = 998244353; constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1}; constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1}; constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1}; template inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; } template inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; template struct Matrix { explicit Matrix(const int m, const int n, const T def = 0) : data(m, std::vector(n, def)) {} int nrow() const { return data.size(); } int ncol() const { return data.front().size(); } Matrix pow(long long exponent) const { const int n = nrow(); Matrix res(n, n, 0), tmp = *this; for (int i = 0; i < n; ++i) { res[i][i] = 1; } for (; exponent > 0; exponent >>= 1) { if (exponent & 1) res *= tmp; tmp *= tmp; } return res; } inline const std::vector& operator[](const int i) const { return data[i]; } inline std::vector& operator[](const int i) { return data[i]; } Matrix& operator=(const Matrix& x) = default; Matrix& operator+=(const Matrix& x) { const int m = nrow(), n = ncol(); for (int i = 0; i < m; ++i) { for (int j = 0; j < n; ++j) { data[i][j] += x[i][j]; } } return *this; } Matrix& operator-=(const Matrix& x) { const int m = nrow(), n = ncol(); for (int i = 0; i < m; ++i) { for (int j = 0; j < n; ++j) { data[i][j] -= x[i][j]; } } return *this; } Matrix& operator*=(const Matrix& x) { const int m = nrow(), l = ncol(), n = x.ncol(); std::vector> res(m, std::vector(n, 0)); for (int i = 0; i < m; ++i) { for (int k = 0; k < l; ++k) { for (int j = 0; j < n; ++j) { res[i][j] += data[i][k] * x[k][j]; } } } data.swap(res); return *this; } Matrix operator+(const Matrix& x) const { return Matrix(*this) += x; } Matrix operator-(const Matrix& x) const { return Matrix(*this) -= x; } Matrix operator*(const Matrix& x) const { return Matrix(*this) *= x; } private: std::vector> data; }; template struct SegmentTree { using Monoid = typename T::Monoid; explicit SegmentTree(int n) : SegmentTree(std::vector(n, T::id())) {} explicit SegmentTree(const std::vector& a) : n(a.size()), p2(1) { while (p2 < n) p2 <<= 1; dat.assign(p2 << 1, T::id()); std::copy(a.begin(), a.end(), dat.begin() + p2); for (int i = p2 - 1; i > 0; --i) { dat[i] = T::merge(dat[i << 1], dat[(i << 1) + 1]); } } void set(int idx, const Monoid val) { idx += p2; dat[idx] = val; while (idx >>= 1) dat[idx] = T::merge(dat[idx << 1], dat[(idx << 1) + 1]); } Monoid get(int left, int right) const { Monoid res_l = T::id(), res_r = T::id(); for (left += p2, right += p2; left < right; left >>= 1, right >>= 1) { if (left & 1) res_l = T::merge(res_l, dat[left++]); if (right & 1) res_r = T::merge(dat[--right], res_r); } return T::merge(res_l, res_r); } Monoid operator[](const int idx) const { return dat[idx + p2]; } template int find_right(int left, const G g) { if (left >= n) return n; Monoid val = T::id(); left += p2; do { while (!(left & 1)) left >>= 1; Monoid nxt = T::merge(val, dat[left]); if (!g(nxt)) { while (left < p2) { left <<= 1; nxt = T::merge(val, dat[left]); if (g(nxt)) { val = nxt; ++left; } } return left - p2; } val = nxt; ++left; } while (__builtin_popcount(left) > 1); return n; } template int find_left(int right, const G g) { if (right <= 0) return -1; Monoid val = T::id(); right += p2; do { --right; while (right > 1 && (right & 1)) right >>= 1; Monoid nxt = T::merge(dat[right], val); if (!g(nxt)) { while (right < p2) { right = (right << 1) + 1; nxt = T::merge(dat[right], val); if (g(nxt)) { val = nxt; --right; } } return right - p2; } val = nxt; } while (__builtin_popcount(right) > 1); return -1; } private: const int n; int p2; std::vector dat; }; namespace monoid { template struct RangeMinimumQuery { using Monoid = T; static constexpr Monoid id() { return std::numeric_limits::max(); } static Monoid merge(const Monoid& a, const Monoid& b) { return std::min(a, b); } }; template struct RangeMaximumQuery { using Monoid = T; static constexpr Monoid id() { return std::numeric_limits::lowest(); } static Monoid merge(const Monoid& a, const Monoid& b) { return std::max(a, b); } }; template struct RangeSumQuery { using Monoid = T; static constexpr Monoid id() { return 0; } static Monoid merge(const Monoid& a, const Monoid& b) { return a + b; } }; } // namespace monoid int main() { struct S { using Monoid = Matrix; static Monoid id() { Matrix m(3, 3); m[0][0] = m[1][1] = m[2][2] = 1; return m; } static Monoid merge(const Monoid& a, const Monoid& b) { return b * a; } }; int n; cin >> n; vector> a(n, S::id()); REP(i, n) { double p, q, r; cin >> p >> q >> r; Matrix x = S::id(), y = S::id(); x[0][2] = -p; x[1][2] = -q; a[i] = S::merge(a[i], x); const double theta = M_PI * r / 180; y[0][0] = y[1][1] = cos(theta); y[0][1] = -sin(theta); y[1][0] = sin(theta); a[i] = S::merge(a[i], y); x[0][2] = p; x[1][2] = q; a[i] = S::merge(a[i], x); } SegmentTree seg(a); int q; cin >> q; while (q--) { int s, t; cin >> s >> t; --s; --t; Matrix ans(3, 1, 1); REP(i, 2) cin >> ans[i][0]; ans = seg.get(s, t + 1) * ans; cout << ans[0][0] << ' ' << ans[1][0] << '\n'; } return 0; }