#pragma GCC target ("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx") #include #include #include #include using namespace std; #define rep(i,n) for(int i=0;itrans(matrix m, double x, double y) { vectorres(2); res[0] = m.dat[0][0] * x + m.dat[0][1] * y + m.dat[0][2]; res[1] = m.dat[1][0] * x + m.dat[1][1] * y + m.dat[1][2]; return res; } struct segtree { matrix ide; vectornode; int n; void init(int N) { n = 1; while (n < N)n *= 2; node.resize(2 * n); } void upd(double theta, double p, double q, int k) { k += n - 1; node[k].init(theta, p, q); while (k > 0) { k = (k - 1) / 2; node[k] = mul(node[k * 2 + 1], node[k * 2 + 2]); } } matrix calc(int a, int b, int k, int l, int r) { if (r <= a || b <= l)return ide; if (a <= l&&r <= b)return node[k]; matrix vl = calc(a, b, k * 2 + 1, l, (l + r) / 2); matrix vr = calc(a, b, k * 2 + 2, (l + r) / 2, r); return mul(vl, vr); } matrix get(int a, int b) { return calc(a, b, 0, 0, n); } }; segtree seg; signed main() { cin.tie(0); ios::sync_with_stdio(false); int n; cin >> n; seg.init(n); rep(i, n) { double p, q, r; cin >> p >> q >> r; seg.upd(r, p, q, i); } int q; cin >> q; rep(j, q) { int s, t; double x, y; cin >> s >> t >> x >> y; if(s>t)return 1; matrix M = seg.get(s - 1, t); vectorres = trans(M, x, y); printf("%.13lf %.13lf\n", res[0], res[1]); } }