結果

問題 No.3437 [Cherry 8th Tune C] Silhouette
コンテスト
ユーザー abc864197532
提出日時 2026-01-23 21:41:15
言語 C++23
(gcc 15.2.0 + boost 1.89.0)
結果
AC  
実行時間 215 ms / 2,000 ms
コード長 2,076 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,116 ms
コンパイル使用メモリ 339,124 KB
実行使用メモリ 7,848 KB
最終ジャッジ日時 2026-01-23 21:43:18
合計ジャッジ時間 6,899 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 11
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define pb push_back
#define all(a) a.begin(), a.end()
#define sz(a) ((int)a.size())
#ifdef Doludu
template <typename T>
ostream& operator << (ostream &o, vector <T> vec) {
    o << "{"; int f = 0;
    for (T i : vec) o << (f++ ? " " : "") << i;
    return o << "}";
}
void bug__(int c, auto ...a) {
    cerr << "\e[1;" << c << "m";
    (..., (cerr << a << " "));
    cerr << "\e[0m" << endl;
}
#define bug_(c, x...) bug__(c, __LINE__, "[" + string(#x) + "]", x)
#define bug(x...) bug_(32, x)
#define bugv(x...) bug_(36, vector(x))
#define safe bug_(33, "safe")
#else
#define bug(x...) void(0)
#define bugv(x...) void(0)
#define safe void(0)
#endif
const int mod = 998244353, N = 1 << 20;

int add(int a, int b) { a += b; if (a >= mod) a -= mod; return a; }
int sub(int a, int b) { a -= b; if (a < 0) a += mod; return a; }
int mul(int a, int b) { return 1ll * a * b % mod; }
int Pow(int a, int b) { int ans = 1; for (; b; b >>= 1, a = mul(a, a)) if (b & 1) ans = mul(ans, a); return ans; }

int main() {
    ios::sync_with_stdio(false), cin.tie(0);
    int t;
    cin >> t;
    while (t--) {
        vector <array <ll, 3>> a(4);
        for (int i = 0; i <= 3; ++i) for (int j = 0; j < 3; ++j) cin >> a[i][j];
        vector <array <ll, 3>> vec;
        for (int i = 0; i < 3; ++i) {
            auto [x, y, z] = a[i];
            auto [sx, sy, sz] = a[3];
            x -= sx, y -= sy, z -= sz;
            z = abs(z);
            vec.pb({-x * sz + sx * z, -y * sz + sy * z, z});
        }
        __int128 ansa = 0, ansb = vec[0][2] * vec[1][2] * vec[2][2];
        for (int i = 0; i < 3; ++i) {
            auto [x1, y1, m1] = vec[i];
            auto [x2, y2, m2] = vec[(i + 1) % 3];
            ll m3 = vec[(i + 2) % 3][2];
            ansa += ((__int128)x1 * y2 - (__int128)y1 * x2) * m3;
        }
        if (ansa < 0) ansa = -ansa;
        ansa %= mod, ansb %= mod;
        ansa = mul(ansa, Pow(2, mod - 2));
        cout << mul(ansa, Pow(ansb, mod - 2)) << "\n";
    }
}
0