結果
問題 | No.1590 Random Shopping |
ユーザー |
|
提出日時 | 2021-07-08 23:26:27 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 272 ms / 5,000 ms |
コード長 | 2,232 bytes |
コンパイル時間 | 2,335 ms |
コンパイル使用メモリ | 200,788 KB |
最終ジャッジ日時 | 2025-01-22 19:29:11 |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i, n) for(int i = 0; i < n; i++) #define rep2(i, x, n) for(int i = x; i <= n; i++) #define rep3(i, x, n) for(int i = x; i >= n; i--) #define each(e, v) for(auto &e: v) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define sz(x) (int)x.size() using ll = long long; using pii = pair<int, int>; using pil = pair<int, ll>; using pli = pair<ll, int>; using pll = pair<ll, ll>; const int MOD = 1000000007; //const int MOD = 998244353; const int inf = (1<<30)-1; const ll INF = (1LL<<60)-1; template<typename T> bool chmax(T &x, const T &y) {return (x < y)? (x = y, true) : false;}; template<typename T> bool chmin(T &x, const T &y) {return (x > y)? (x = y, true) : false;}; struct io_setup{ io_setup(){ ios_base::sync_with_stdio(false); cin.tie(NULL); cout << fixed << setprecision(15); } } io_setup; int main(){ int N; cin >> N; vector<int> A(N), R(N); rep(i, N) cin >> A[i]; rep(i, N) cin >> R[i]; vector<bool> used(1001, false); rep(i, N) used[A[i]] = true; vector<vector<double>> P(N, vector<double>(1001, 0.0)); vector<vector<double>> dp(N+1, vector<double>(N+1, 0.0)); rep(i, 1001){ if(!used[i]) continue; rep(j, N+1){ rep(k, N+1){ dp[j][k] = 0.0; } } dp[0][0] = 1.0; rep(j, N){ rep(k, N){ if(A[j] < i){ dp[j+1][k] += dp[j][k]*0.5; dp[j+1][k+1] += dp[j][k]*0.5; } else{ if(k == 0){ P[j][i] += dp[j][k]*0.5; dp[j+1][k] += dp[j][k]; } else{ dp[j+1][k] += dp[j][k]*0.5; dp[j+1][k-1] += dp[j][k]*0.5; } } } } } double ans = 0; int pre = 0; rep(i, 1001){ if(used[i]){ rep(j, N){ ans += P[j][i]*R[j]*(i-pre); } pre = i; } } cout << ans << '\n'; }