結果

問題 No.1932 動く点 P / Moving Point P
コンテスト
ユーザー kabipoyo
提出日時 2022-05-06 23:57:14
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,281 ms / 6,000 ms
コード長 2,852 bytes
コンパイル時間 5,001 ms
コンパイル使用メモリ 261,028 KB
最終ジャッジ日時 2025-01-29 04:15:44
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 11
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
typedef int64_t lint;
#define rep(i, n) for(int i=0; i<n; i++)
#define repx(i, l, n) for(int i=l; i<n; i++)
#define all(v) v.begin(), v.end()
#define show(x) cout << #x << ": " << x << endl;
#define list(x) cout << #x << ": " << x << "  ";
#define pb push_back
using vi = vector<double>;
using vvi = vector<vector<double>>;
template<class T> inline void vin(vector<T>& v) { rep(i, v.size()) cin >> v.at(i); }
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
template<class T> inline void drop(T x) { cout << x << endl; exit(0); }
template<class T> void vout(vector<T> v) { rep(i, v.size()) { cout << v.at(i) << ' '; } cout << endl; }
constexpr lint LINF = LLONG_MAX/2;
constexpr double PI=3.14159265358979323846;
#define deg_to_rad(deg) (((deg)/360)*2*PI)
#define rad_to_deg(rad) (((rad)/2/PI)*360)
vvi matrix_mul(vvi A, vvi B) {
	vvi C(A.size(), vi(B[0].size()));
	rep(i, A.size()) {
		rep(k, B.size()) {
			rep(j, B[0].size()) {
				C[i][j] += A[i][k] * B[k][j];
			}
		}
	}
	return C;
}

int main() {
	lint N, Q;
	cin >> N;
	vi p(N), q(N), r(N);
	rep(i, N) cin >> p[i] >> q[i] >> r[i];

	lint a=0, b=0, c=0;
	double x, y, z;
	std::vector<vvi> v, w;
	vvi e(3, vi(3));
	e[0][0] = 1;
	e[1][1] = 1;
	e[2][2] = 1;
	v.pb(e);
	w.pb(e);
	rep(i, N) {
		vvi t(3, vi(3));
		x = sin(deg_to_rad(r[i]));
		y = cos(deg_to_rad(r[i]));
		t[0][0] = y;
		t[0][1] = -x;
		t[0][2] = p[i]-p[i]*y+q[i]*x;
		t[1][0] = x;
		t[1][1] = y;
		t[1][2] = q[i]-p[i]*x-q[i]*y;
		t[2][2] = 1;
		v.pb(t);
		x = sin(deg_to_rad(-r[i]));
		y = cos(deg_to_rad(-r[i]));
		t[0][0] = y;
		t[0][1] = -x;
		t[0][2] = p[i]-p[i]*y+q[i]*x;
		t[1][0] = x;
		t[1][1] = y;
		t[1][2] = q[i]-p[i]*x-q[i]*y;
		t[2][2] = 1;
		w.pb(t);
	}
	// cout << fixed << setprecision(10);
	// rep(i, N) {
	// 	rep(j, 3) {
	// 		rep(k, 3) {
	// 			std::cout << v[i][j][k] << ' ';
	// 		}
	// 		std::cout << '\n';
	// 	}
	// 	std::cout << '\n';
	// }
	rep(i, N) {
		v[i+1] = matrix_mul(v[i+1], v[i]);
		w[i+1] = matrix_mul(w[i], w[i+1]);
	}

	cin >> Q;
	rep(_, Q) {
		cin >> a >> b >> x >> y;
		a--;
		vvi t(3, vi(1));
		t[0][0] = x;
		t[1][0] = y;
		t[2][0] = 1;
		// t = matrix_mul(w[a], t);
		// cout << fixed << setprecision(10) << t[0][0] << ' ' << t[1][0] << endl;
		// t = matrix_mul(v[b], t);
		// cout << fixed << setprecision(10) << t[0][0] << ' ' << t[1][0] << endl;
		// std::cout << '\n';
		t = matrix_mul(matrix_mul(v[b], w[a]), t);
		cout << fixed << setprecision(10) << t[0][0] << ' ' << t[1][0] << endl;
	}
	// rep(i, N) {
	// 	rep(j, 3) {
	// 		rep(k, 3) {
	// 			std::cout << v[i][j][k] << ' ';
	// 		}
	// 		std::cout << '\n';
	// 	}
	// 	std::cout << '\n';
	// }
}
0