#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return true; } return false; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return true; } return false; }
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define mp make_pair
#define si(x) int(x.size())
const int mod=998244353,MAX=100005,INF=1<<30;

//Union Find


struct UF{
    int n;
    vector<int> par,size,edge;
    
    void init(int n_){
        n=n_;
        par.assign(n,-1);
        size.assign(n,1);
        edge.assign(n,0);
        
        for(int i=0;i<n;i++){
            par[i]=i;
        }
    }
    
    int root(int a){
        if(par[a]==a) return a;
        else return par[a]=root(par[a]);
    }
    
    void unite(int a,int b){
        edge[root(a)]++;
        if(root(a)!=root(b)){
            size[root(a)]+=size[root(b)];
            edge[root(a)]+=edge[root(b)];
            par[root(b)]=root(a);
        }
    }
    
    bool check(int a,int b){
        return root(a)==root(b);
    }
};

vector<int> G[MAX];

ll X[MAX],Y[MAX];

void solve(int u){
    for(int to:G[u]){
        solve(to);
        X[u]-=(Y[u]-X[u]);
    }
    chmin(X[u],Y[u]);
}

int main(){
    
    std::ifstream in("text.txt");
    std::cin.rdbuf(in.rdbuf());
    cin.tie(0);
    ios::sync_with_stdio(false);
    
    int N;cin>>N;
    vector<ll> A(N),B(N),C(N);
    vector<int> deg(N);
    UF uf;uf.init(N);
    for(int i=0;i<N;i++){
        cin>>A[i];
    }
    for(int i=0;i<N;i++){
        cin>>B[i];
    }
    for(int i=0;i<N;i++){
        cin>>C[i];C[i]--;
        deg[C[i]]++;
    }
    
    queue<int> Q;
    for(int i=0;i<N;i++){
        if(deg[i]==0) Q.push(i);
    }
    
    while(!Q.empty()){
        int u=Q.front();Q.pop();
        G[C[u]].push_back(u);
        deg[C[u]]--;
        if(deg[C[u]]==0) Q.push(C[u]);
    }
    
    for(int i=0;i<N;i++){
        if(deg[i]){
            uf.unite(C[i],i);
        }
    }
    
    for(int i=0;i<N;i++){
        if(deg[i]){
            if(uf.root(i)!=i){
                int ro=uf.root(i);
                A[ro]+=A[i];
                B[ro]+=B[i];
                for(int x:G[i]) G[ro].push_back(x);
                G[i].clear();
                A[i]=0;
                B[i]=0;
            }
        }
    }
    
    ll left=0,right=1LL<<60;
    while(right-left>1){
        ll mid=(left+right)/2;
        for(int i=0;i<N;i++){
            X[i]=A[i];
            if(B[i]<=(1LL<<60)/mid) Y[i]=B[i]*mid;
            else Y[i]=1LL<<60;
        }
        bool ok=true;
        for(int i=0;i<N;i++){
            if(deg[i]&&uf.root(i)==i){
                solve(i);
            }
        }
        for(int i=0;i<N;i++) ok&=(X[i]>=Y[i]);
        
        if(ok) left=mid;
        else right=mid;
    }
    cout<<left<<endl;
}