結果

問題 No.3073 Fraction Median
ユーザー maeshun
提出日時 2025-03-22 16:48:51
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,440 bytes
コンパイル時間 7,921 ms
コンパイル使用メモリ 468,848 KB
実行使用メモリ 15,020 KB
最終ジャッジ日時 2025-03-22 16:49:17
合計ジャッジ時間 25,603 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 14 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#include <atcoder/all>
#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>
namespace mp = boost::multiprecision;
using Bint = mp::cpp_int;
using namespace std;
using namespace atcoder;
#define rep(i, n) for(int i=0;i<(n);++i)
#define rep1(i, n) for(int i=1;i<=(n);i++)
#define ll long long
using mint = modint998244353;
using P = pair<ll,ll>;
using lb = long double;
using T = tuple<ll, ll, ll>;
#ifdef LOCAL
#  include <debug_print.hpp>
#  define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__)
#else
#  define dbg(...) (static_cast<void>(0))
#endif

int main()
{
    ll n;
    cin >> n;
    vector<ll> a(n);
    rep(i,n) cin >> a[i];
    lb high = 1e9;
    lb low = 0;
    int X = -1;
    int Y = -1;
    rep(j,60) {
        lb mid = (high+low)/2;
        int id = 0;
        ll sum = 0;
        lb d = 1e9;
        int x = -1;
        int y = -1;
        rep(i,n) {
            lb val = a[i]*mid;
            while(id<n && a[id]<=val) id++;
            if(id-1>=0 && d>abs((lb)a[id-1]/(lb)a[i]-mid)) {
                x = a[id-1];
                y = a[i];
                dbg(id,i);
                d = abs((lb)a[id-1]/(lb)a[i]-mid);
            }
            sum += id;
        }
        if(sum>=n*(n-1)/2) high = mid;
        else low = mid;
        X = x;
        Y = y;
    }
    ll g = gcd(X,Y);
    cout << X/g << " " << Y/g << endl; 
    return 0;
}
0