結果

問題 No.1626 三角形の構築
コンテスト
ユーザー KoD
提出日時 2021-07-23 22:58:39
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 2,416 bytes
コンパイル時間 1,821 ms
コンパイル使用メモリ 202,788 KB
最終ジャッジ日時 2025-01-23 08:37:41
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using i32 = std::int32_t;
using u32 = std::uint32_t;
using i64 = std::int64_t;
using u64 = std::uint64_t;
using i128 = __int128_t;
using u128 = __uint128_t;
using isize = std::ptrdiff_t;
using usize = std::size_t;

class rep {
    struct Iter {
        usize itr;
        constexpr Iter(const usize pos) noexcept : itr(pos) {}
        constexpr void operator++() noexcept { ++itr; }
        constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; }
        constexpr usize operator*() const noexcept { return itr; }
    };
    const Iter first, last;

  public:
    explicit constexpr rep(const usize first, const usize last) noexcept : first(first), last(std::max(first, last)) {}
    constexpr Iter begin() const noexcept { return first; }
    constexpr Iter end() const noexcept { return last; }
};

template <class T> using Vec = std::vector<T>;

void main_() {
    i64 S, T;
    std::cin >> S >> T;
    if (16 * S * S % T != 0) {
        std::cout << 0 << '\n';
        return;
    }
    const auto U = 16 * S * S / T;
    std::set<std::array<i64, 3>> set;
    for (i64 x = 1; x * x * x <= U and 3 * x <= T; ++x) {
        if (U % x == 0 and (T - x) % 2 == 0) {
            const auto sum = T - x;
            const auto prd = U / x;
            // t * t - sum * t  + prd = 0 の 2 解
            const auto D = sum * sum - 4 * prd;
            if (D < 0) {
                continue;
            }
            const i64 sqrt = std::sqrt(D);
            if (sqrt * sqrt != D) {
                continue;
            }
            if (sum % 2 != sqrt % 2) {
                continue;
            }
            const auto y = (sum - sqrt) / 2;
            const auto z = (sum + sqrt) / 2;
            if ((T - y) % 2 != 0 or (T - z) % 2 != 0 or y <= 0 or z <= 0) {
                continue;
            }
            std::array<i64, 3> arr;
            arr[0] = (T - x) / 2;
            arr[1] = (T - y) / 2;
            arr[2] = (T - z) / 2;
            std::sort(arr.begin(), arr.end());
            set.insert(std::move(arr));
        }
    }
    std::cout << set.size() << '\n';
    for (const auto& [a, b, c] : set) {
        std::cout << a << ' ' << b << ' ' << c << '\n';
    }
}

int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);
    usize cases;
    std::cin >> cases;
    while (cases--) main_();
    return 0;
}
0