#include <iostream>
#include <fstream>
#include <vector>
#include <set>
#include <map>
#include <cstring>
#include <string>
#include <cmath>
#include <cassert>
#include <ctime>
#include <algorithm>
#include <sstream>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <cstdlib>
#include <cstdio>
#include <iterator>
#include <functional>
#include <bitset>
#include <random>

#ifdef LOCAL
#define eprintf(...) fprintf(stderr,__VA_ARGS__)
clock_t __my_start_clock = clock();
#else
#define eprintf(...)
#endif

#define TIMESTAMP(x) eprintf("["#x"] Time : %.3lf s.\n", (clock() - __my_start_clock)*1.0/CLOCKS_PER_SEC)
#define TIMESTAMPf(x,...) eprintf("[" x "] Time : %.3lf s.\n", __VA_ARGS__, (clock() - __my_start_clock)*1.0/CLOCKS_PER_SEC)

#define db(x) cerr << #x << " = " << x << endl
#define db2(x, y) cerr << "(" << #x << ", " << #y << ") = (" << x << ", " << y << ")\n";
#define db3(x, y, z) cerr << "(" << #x << ", " << #y << ", " << #z << ") = (" << x << ", " << y << ", " << z << ")\n"
#define dbv(a) cerr << #a << " = "; for (auto xxxx: a) cerr << xxxx << " "; cerr << endl


using namespace std;

#ifdef LOCAL
static struct __timestamper {
    string what;
    __timestamper(const char* what) : what(what){};
    __timestamper(const string& what) : what(what){};
    ~__timestamper(){
        TIMESTAMPf("%s", what.data());
    }
} __TIMESTAMPER("end");
#else
struct __timestamper {};
#endif

typedef long long ll;
typedef long double dbl;

const int MOD = (int)1e9 + 7;

void add(int &x, ll y) {
    x = (x + y) % MOD;
}

int mul(int x, int y) {
    return 1LL * x * y % MOD;
}

int solve(int h, vector<int> p) {
    int k = 1 << (h - 1);
    assert((int)p.size() == k);

    int n = 1 << h;
    vector<int> L(n), R(n);
    for (int i = 0; i < k; i++) L[k + i] = R[k + i] = i;
    for (int i = k - 1; i >= 1; i--) {
        L[i] = L[2 * i];
        R[i] = R[2 * i + 1];
    }

    vector<int> sum(n);
    auto enable = [&](int x, int prod) {
        x += k;
        while (x > 0) {
            prod = mul(prod, x);
            add(sum[x], prod);
            x >>= 1;
        }
    };

    vector<int> cprod(k);
    for (int i = 0; i < k; i++) cprod[i] = k + i;

    int ans = 0;
    for (int v = k - 1; v >= 1; v--) {
        for (int i = L[v * 2 + 1]; i <= R[v * 2 + 1]; i++) {
            enable(p[i], cprod[i]);
        }

        for (int i = L[v * 2]; i <= R[v * 2]; i++) {
            int oth = k + p[i];

            int rprod = mul(mul(v, cprod[i]), oth);
            while (oth > 1) {
                int pr = oth >> 1;
                rprod = mul(rprod, pr);
                add(ans, mul(rprod, sum[oth ^ 1]));
                oth = pr;
            }
        }

        for (int i = L[v * 2 + 1]; i <= R[v  * 2 + 1]; i++) {
            enable(p[i], mul(cprod[i], MOD - 1));
        }

        for (int i = L[v]; i <= R[v]; i++) cprod[i] = mul(cprod[i], v);
    }
    return ans;
}

int main() {
#ifdef LOCAL
    freopen("in", "r", stdin);
#endif
    ios_base::sync_with_stdio(0);
    cin.tie(0);

    int h;
    while (cin >> h) {
        int k = 1 << (h - 1);
        vector<int> p(k);
        for (int i = 0; i < k; i++) {
            cin >> p[i];
            p[i]--;
        }

        cout << solve(h, p) << endl;
    }

    return 0;
}


Measured memory_bytes: 7896