#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define int long long #define ii pair #define app push_back #define all(a) a.begin(), a.end() #define bp __builtin_popcountll #define ll long long #define mp make_pair #define x first #define y second #define Time (double)clock()/CLOCKS_PER_SEC #define debug(x) std::cerr << #x << ": " << x << '\n'; #define FORI(i,a,b) for (int i = (a); i < (b); ++i) #define FOR(i,a) FORI(i,0,a) #define ROFI(i,a,b) for (int i = (b)-1; i >= (a); --i) #define ROF(i,a) ROFI(i,0,a) #define rep(a) FOR(_,a) #define each(a,x) for (auto& a: x) #define FORN(i, n) FORI(i, 1, n + 1) using vi = vector; template std::istream& operator >>(std::istream& input, std::pair & data) { input >> data.x >> data.y; return input; } template std::istream& operator >>(std::istream& input, std::vector& data) { for (T& x : data) input >> x; return input; } template std::ostream& operator <<(std::ostream& output, const pair & data) { output << "(" << data.x << "," << data.y << ")"; return output; } template std::ostream& operator <<(std::ostream& output, const std::vector& data) { for (const T& x : data) output << x << " "; return output; } std::ostream& operator <<(std::ostream& output, const __int128 &x) { __int128 n = x; if (n == 0) { output << "0"; return output; } if (n < 0) { n = -n; output << "-"; } string s; while (n) { s += '0' + (n%10); n /= 10; } reverse(all(s)); cout << s; return output; } ll div_up(ll a, ll b) { return a/b+((a^b)>0&&a%b); } // divide a by b rounded up ll div_down(ll a, ll b) { return a/b-((a^b)<0&&a%b); } // divide a by b rounded down ll math_mod(ll a, ll b) { return a - b * div_down(a, b); } #define tcT template using V = vector; tcT> bool ckmin(T& a, const T& b) { return b < a ? a = b, 1 : 0; } // set a = min(a,b) tcT> bool ckmax(T& a, const T& b) { return a < b ? a = b, 1 : 0; } ll gcd(ll a, ll b) { while (b) { tie(a, b) = mp(b, a % b); } return a; } int Bit(int mask, int bit) { return (mask >> bit) & 1; } #define ld long double const int N = 507; ld dp[N][N]; signed main() { #ifdef LOCAL #else #define endl '\n' ios_base::sync_with_stdio(0); cin.tie(0); #endif int n; cin >> n; vi a(n); cin >> a; vi r(n); cin >> r; ld ans = 0; FOR (who, n) { FOR (i, n) { FOR (j, n) { dp[i][j] = 0; } } dp[0][0] = 1; FOR (i, n) { FOR (j, i + 1) { if (dp[i][j] > 0) { if (mp(a[i], i) < mp(a[who], who)) { dp[i + 1][j + 1] += dp[i][j]/2; dp[i + 1][j] += dp[i][j]/2; } else { if (j == 0) { if (i >= who) { //debug(who); //debug(mp(i, j)); //debug(dp[i][j]); ans += dp[i][j] * a[who] * r[i]; dp[i + 1][0] += dp[i][j]/2; } else { dp[i + 1][0] += dp[i][j]; } } else { dp[i + 1][j] += dp[i][j]/2; dp[i + 1][j - 1] += dp[i][j]/2; } } } } } } cout.precision(20); cout << fixed << ans/2 << endl; }