結果

問題 No.620 ぐるぐるぐるりん
ユーザー PachicobuePachicobue
提出日時 2017-12-20 22:08:54
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 56 ms / 1,000 ms
コード長 819 bytes
コンパイル時間 2,379 ms
コンパイル使用メモリ 197,608 KB
最終ジャッジ日時 2025-01-05 06:08:43
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
using ld = double;
using C = complex<ld>;

int main()
{
    int N;
    cin >> N;
    for (int q = 0; q < N; q++) {
        int T;
        cin >> T;
        ld p, omega, v, gx, gy;
        cin >> p >> omega >> v >> gx >> gy;
        const C A = C{1 + v, -omega};
        const C a = pow(A, T);
        C g{gx - a.real(), gy + a.imag()};
        const ld l = hypot(1 + v, omega);
        const ld F = (l == 1 ? T : (pow(l, 2 * T) - 1) / (l * l - 1));
        g /= F;
        vector<C> ans(T);
        ans[T - 1] = g;
        for (int i = T - 2; i >= 0; i--) {
            ans[i] = ans[i + 1] * A;
        }
        for (int i = 0; i < T; i++) {
            cout << fixed << setprecision(50) << ans[i].real() << " " << ans[i].imag() << endl;
        }
    }
    return 0;
}
0