結果

問題 No.3483 A Forbidden Fruit
コンテスト
ユーザー tombo_
提出日時 2026-03-27 22:05:39
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
TLE  
実行時間 -
コード長 2,392 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 916 ms
コンパイル使用メモリ 110,668 KB
実行使用メモリ 9,984 KB
最終ジャッジ日時 2026-03-27 22:06:12
合計ジャッジ時間 9,878 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other TLE * 1 -- * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<iostream>
#include<iomanip>
#include<vector>
#include<queue>
#include<math.h>
using namespace std;

#define OVERLOAD_REP(_1, _2, _3, name, ...) name
#define REP1(i, n) for (auto i = std::decay_t<decltype(n)>{}; (i) != (n); ++(i))
#define REP2(i, l, r) for (auto i = (l); (i) != (r); ++(i))
#define rep(...) OVERLOAD_REP(__VA_ARGS__, REP2, REP1)(__VA_ARGS__)
#define REP(i, l, r) rep(i, l, r + 1)
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()

using ll = long long;
using ld = long double;

struct Edge {
    int to;
    ll w;
};

using Graph = vector<vector<int>>;
// using Graph = vector<vector<Edge> >;
const ll INF = 2e18;
// const int INF = 2e9;

template <class T> using vc = vector<T>;
template <class T> using vv = vector<vector<T>>;
template <class T> using vvv = vector<vector<vector<T>>>;
template <class T> using pq = priority_queue<T>;
template <class T> using pq_g = priority_queue<T, vc<T>, greater<T>>;

template <class T> istream &operator>>(istream &i, vc<T> &v) {
    rep(j, 0, v.size()) i >> v[j];
    return i;
}
template <class T> ostream &operator<<(ostream &o, vc<T> &v) {
    rep(j, 0, v.size()) o << v[j] << " ";
    return o;
}

template <class T> bool chmin(T &a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T> bool chmax(T &a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}

constexpr double EPS = 1e-9;
inline bool eq(double a, double b) { return fabs(a - b) < EPS; }
inline bool lt(double a, double b) { return a < b - EPS; }
inline bool gt(double a, double b) { return a > b + EPS; }
inline int sgn(double x) {
    if (x > EPS) return 1;
    if (x < -EPS) return -1;
    return 0;
}
int T,N,M;ll K;
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(12);

    cin>>T;
    while(T--){
        cin>>N>>M>>K;
        ll a=M-1,b=K-(N-1)*M-1;
        ld prob=0,cur=1;
        if(K>(ll)M*(N-1)+1)rep(i,0,b){
            prob+=cur/(a-i);
            cur-=cur/(a-i);
        }
        ld ans=0;
        rep(i,0,N){
            if((ll)i*M>=K)break;
            if(K<=(ll)M*(N-1)+1)ans+=1.0/N/M;
            else{
                ld a=1.0/N;
                ans+=a/M;
                ans+=a*(M-1)/M*prob;
            }
        }
        cout<<1.0-ans<<endl;
    }

    return 0;
}
0