#define _USE_MATH_DEFINES #include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define MT make_tuple #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)<; using vi = vector; using vll = vector; const double EPS = 1e-9; using Vec = vector; using Matrix = vector; Vec GaussJordan(Matrix &A) { int n = (int)A.size(); for (int i = 0; i < n; ++i) { int pivot = i; for (int j = i + 1; j < n; ++j) if (abs(A[j][i])>abs(A[pivot][i])) pivot = j; swap(A[i], A[pivot]); if (abs(A[i][i]) < EPS && abs(A[i][n]) > EPS) return Vec(); for (int j = n; j >= i; --j)A[i][j] /= A[i][i]; for (int j = 0; j < n; ++j)if (j != i) for (int k = n; k >= i; --k)A[j][k] -= A[i][k] * A[j][i]; } Vec res(n); for (int i = 0; i < n; ++i)res[i] = A[i][n]; return res; } void solve() { double xa, ya, ta, x[2][2], y[2][2], co, si; cin >> xa >> ya >> ta; rep(i, 2)rep(j, 2)cin >> x[i][j] >> y[i][j]; double xx = x[0][0] - x[0][1], yy = y[0][0] - y[0][1]; double XX = x[1][0] - x[1][1], YY = y[1][0] - y[1][1]; Matrix A(2, Vec(3)); A[0][0] = A[1][1] = xx; A[0][1] = -yy; A[1][0] = yy; A[0][2] = XX; A[1][2] = YY; auto gd = GaussJordan(A); co = gd[0], si = gd[1]; double p = x[1][0] - x[0][0] * co + y[0][0] * si; double q = y[1][0] - x[0][0] * si - y[0][0] * co; double t = ta + acos(co) / M_PI * 180; cout << xa - p << ' ' << ya - q << ' ' << t << endl; } int main() { ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(15); solve(); return 0; }