結果
問題 | No.574 正多面体サイコロ |
ユーザー |
|
提出日時 | 2017-10-07 04:12:33 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 39 ms / 2,000 ms |
コード長 | 1,496 bytes |
コンパイル時間 | 1,443 ms |
コンパイル使用メモリ | 166,996 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-17 03:39:54 |
合計ジャッジ時間 | 2,302 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include <bits/stdc++.h>#define show(x) cerr << #x << " = " << x << endlusing namespace std;using ll = long long;using pii = pair<int, int>;using vi = vector<int>;template <typename T>ostream& operator<<(ostream& os, const vector<T>& v){os << "sz=" << v.size() << "\n[";for (const auto& p : v) {os << p << ",";}os << "]\n";return os;}template <typename S, typename T>ostream& operator<<(ostream& os, const pair<S, T>& p){os << "(" << p.first << "," << p.second<< ")";return os;}constexpr ll MOD = 1e9 + 7;template <typename T>constexpr T INF = numeric_limits<T>::max() / 100;int main(){cin.tie(0);ios::sync_with_stdio(false);ll F, N, K;cin >> F >> N >> K;using ld = long double;ld sum = 0;for (ll i = 1; i <= F; i++) {for (ll k = 0; k < K; k++) {for (ll l = K - k; l <= N - k; l++) {ld prob = 1;for (ll j = 0; j < k; j++) {prob *= (long double)(N - j) / (k - j);}for (ll j = 0; j < l; j++) {prob *= (long double)(N - k - j) / (l - j);}prob *= pow(((long double)(i - 1) / F), N - k - l);prob *= pow(((long double)(F - i) / F), k);prob /= pow(F, l);sum += i * prob;}}}cout << fixed << setprecision(10) << sum << endl;return 0;}