結果

問題 No.1590 Random Shopping
ユーザー tokusakuraitokusakurai
提出日時 2021-07-08 23:26:27
言語 C++17
(gcc 11.2.0 + boost 1.78.0)
結果
AC  
実行時間 334 ms / 5,000 ms
コード長 2,232 bytes
コンパイル時間 1,898 ms
使用メモリ 8,912 KB
最終ジャッジ日時 2023-02-01 19:24:03
合計ジャッジ時間 5,060 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 2 ms
3,644 KB
testcase_01 AC 2 ms
3,796 KB
testcase_02 AC 1 ms
3,620 KB
testcase_03 AC 1 ms
3,624 KB
testcase_04 AC 2 ms
3,632 KB
testcase_05 AC 2 ms
3,788 KB
testcase_06 AC 1 ms
3,832 KB
testcase_07 AC 2 ms
3,828 KB
testcase_08 AC 1 ms
3,728 KB
testcase_09 AC 2 ms
3,864 KB
testcase_10 AC 4 ms
3,896 KB
testcase_11 AC 8 ms
4,656 KB
testcase_12 AC 25 ms
5,164 KB
testcase_13 AC 90 ms
6,788 KB
testcase_14 AC 264 ms
8,872 KB
testcase_15 AC 267 ms
8,896 KB
testcase_16 AC 264 ms
8,852 KB
testcase_17 AC 259 ms
8,872 KB
testcase_18 AC 2 ms
3,632 KB
testcase_19 AC 332 ms
8,912 KB
testcase_20 AC 2 ms
3,716 KB
testcase_21 AC 334 ms
8,788 KB
testcase_22 AC 1 ms
3,672 KB
testcase_23 AC 8 ms
8,788 KB
testcase_24 AC 8 ms
8,912 KB
testcase_25 AC 2 ms
3,732 KB
testcase_26 AC 1 ms
3,804 KB
testcase_27 AC 1 ms
3,688 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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';
}
0