#include #include #include #include #include #include #include #include #include #include #include #include #define debug_value(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << #x << "=" << x << endl; #define debug(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << x << endl; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } using namespace std; typedef long long ll; using P = pair; int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << setprecision(10) << fixed; int n; cin >> n; vector a(n), r(n); vector

p(n); for(int i = 0; i < n; i++) { cin >> a[i]; p[i] = P(a[i], -i); } for(int i = 0; i < n; i++) cin >> r[i]; double ans = 0.0; for(int i = 0; i < n; i++){ vector> dp(n+1, vector(n+1)); dp[0][0] = 1.0; for(int j = 0; j <= n; j++){ // cout << accumulate(dp[j].begin(), dp[j].end(), 0.0) << endl; if(j > i && a[j] > a[i]){ ans += dp[j][1]*r[j]*a[i]*0.5; } if(j == i){ // cout << i << ' ' << j << ' ' << dp[j][0] << ' ' << accumulate(dp[j].begin(), dp[j].end(), 0.0) << endl; ans += dp[j][0]*r[i]*a[i]*0.5; } if(j == n) break; // 買わない if(j < i){ if(a[j] < a[i]) dp[j+1][1] += dp[j][0]*0.5; else dp[j+1][0] += dp[j][0]*0.5; }else if(i == j){ dp[j+1][1] += dp[j][0]*0.5; }else{ dp[j+1][0] += dp[j][0]*0.5; } // 買う dp[j+1][0] += dp[j][0]*0.5; for(int k = 1; k <= j; k++){ // 買う if(p[j] <= p[i]){ dp[j+1][k] += dp[j][k]*0.5; }else{ dp[j+1][k-1] += dp[j][k]*0.5; } // 買わない if(p[j] <= p[i]){ dp[j+1][k+1] += dp[j][k]*0.5; }else{ dp[j+1][k] += dp[j][k]*0.5; } } } } cout << ans << endl; }