結果

問題 No.620 ぐるぐるぐるりん
ユーザー 37zigen37zigen
提出日時 2020-05-11 18:20:51
言語 Java21
(openjdk 21)
結果
WA  
実行時間 -
コード長 1,904 bytes
コンパイル時間 1,979 ms
コンパイル使用メモリ 80,504 KB
実行使用メモリ 65,448 KB
最終ジャッジ日時 2024-07-18 16:45:50
合計ジャッジ時間 10,521 ms
ジャッジサーバーID
(参考情報)
judge3 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 157 ms
42,776 KB
testcase_01 AC 138 ms
42,516 KB
testcase_02 AC 147 ms
43,136 KB
testcase_03 AC 129 ms
43,036 KB
testcase_04 AC 144 ms
42,976 KB
testcase_05 AC 128 ms
42,548 KB
testcase_06 AC 154 ms
42,596 KB
testcase_07 AC 147 ms
42,904 KB
testcase_08 AC 134 ms
42,916 KB
testcase_09 AC 144 ms
42,816 KB
testcase_10 AC 132 ms
42,692 KB
testcase_11 AC 152 ms
42,928 KB
testcase_12 AC 152 ms
42,876 KB
testcase_13 AC 138 ms
43,080 KB
testcase_14 AC 154 ms
42,404 KB
testcase_15 AC 142 ms
42,904 KB
testcase_16 AC 151 ms
42,888 KB
testcase_17 AC 125 ms
42,624 KB
testcase_18 AC 155 ms
43,076 KB
testcase_19 AC 132 ms
43,096 KB
testcase_20 AC 139 ms
42,904 KB
testcase_21 AC 133 ms
42,880 KB
testcase_22 AC 133 ms
42,952 KB
testcase_23 AC 147 ms
42,936 KB
testcase_24 AC 133 ms
42,604 KB
testcase_25 AC 131 ms
42,816 KB
testcase_26 AC 576 ms
52,336 KB
testcase_27 AC 130 ms
42,752 KB
testcase_28 AC 597 ms
53,420 KB
testcase_29 WA -
testcase_30 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

import java.util.Arrays;
import java.util.Scanner;

class Main {

	double pow(double a, int n) {
		double ret = 1;
		for (; n > 0; n >>= 1, a = a * a) {
			if (n % 2 == 1) {
				ret = ret * a;
			}
		}
		return ret;
	}

	double[][] pow(double[][] a, int n) {
		double[][] ret = new double[a.length][a.length];
		for (int i = 0; i < ret.length; ++i)
			ret[i][i] = 1;
		for (; n > 0; n >>= 1, a = mul(a, a)) {
			if (n % 2 == 1) {
				ret = mul(ret, a);
			}
		}
		return ret;
	}

	double[][] mul(double[][] a, double[][] b) {
		double[][] ret = new double[a.length][b[0].length];
		for (int i = 0; i < a.length; ++i) {
			for (int j = 0; j < b[i].length; ++j) {
				for (int k = 0; k < a[i].length; ++k) {
					ret[i][j] += a[i][k] * b[k][j];
				}
			}
		}
		return ret;
	}

	void solve(int T, double w, double v, double gx, double gy) {
		double s = Math.sqrt((1 + v) * (1 + v) + w * w);
		double coe = 0;
		for (int i = 0; i <= T - 1; ++i) {
			coe += pow(s, 2 * i);
		}
		coe = 1 / coe;
		double[][] M = { { 1 + v, -w }, { w, 1 + v } };
		double[][] revM = { { 1 + v, w }, { -w, 1 + v } };
		double[][] g = { { gx }, { gy } };
		double[][] src = { { 1 }, { 0 } };
		double[][] h = mul(pow(M, T), src);
		h[0][0] = coe * (g[0][0] - h[0][0]);
		h[1][0] = coe * (g[1][0] - h[1][0]);
		for (int i = 0; i <= T - 1; ++i) {
			double[][] u = mul(pow(revM, T - 1 - i), h);
			System.out.println(u[0][0] + " " + u[1][0]);
		}
	}

	void run() {
		Scanner sc = new Scanner(System.in);
		int n = sc.nextInt();
		for (int i = 0; i < n; ++i) {
			int T = sc.nextInt();
			double p = sc.nextDouble();
			double w = sc.nextDouble();
			double v = sc.nextDouble();
			double gx = sc.nextDouble();
			double gy = sc.nextDouble();
			solve(T, w, v, gx, gy);
		}
	}

	void tr(Object... objects) {
		System.out.println(Arrays.deepToString(objects));
	}

	public static void main(String[] args) {
		new Main().run();
	}
}
0