#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=5009;
int n,s,w[N],p[N],a[N],vis[N],cyc[N],m;
int tag[N],lim;ll val[N][N],dp[N][N];vector<int> G[N];
void dfs(int u){
    ++tag[u];for(int i=0;i<=lim;i++)dp[u][i]=val[u][i]=0;
    for(int v:G[u]){
        if(tag[v])continue;dfs(v);
        for(int i=1;i<=lim;i++)
            dp[u][i]+=val[v][i];
    }
    for(int i=1;i<=lim;i++){
        dp[u][i]+=(i&1)*w[u]-1ll*i*p[u];
        val[u][i]=max(val[u][i-1],dp[u][i]);
    }
    --tag[u];
}
int main(){
    // freopen("cheese.in","r",stdin);
    // freopen("cheese.out","w",stdout);
    scanf("%d%d",&n,&s);
    for(int i=1;i<=n;i++)scanf("%d",w+i);
    for(int i=1;i<=n;i++)scanf("%d",p+i);
    for(int i=1;i<=n;i++)scanf("%d",a+i),G[a[i]].push_back(i);
    for(int u=s;;u=a[u]){cyc[++m]=u;if(vis[u])break;vis[u]=1;}
    if(cyc[m]!=s){
        lim=5;dfs(s);
        printf("%lld\n",max(dp[s][1],dp[s][2])+p[s]);return 0;
    }
    else if(m==3){
        ll ans=w[s];tag[a[s]]=1;lim=5;dfs(s);
        ans=max(ans,max(dp[s][1],dp[s][2])+p[s]);
        tag[a[s]]=0;dfs(s);ans=max(ans,dp[s][1]+p[s]);
        printf("%lld\n",ans);return 0;
    }
    else{
        --m;reverse(cyc+2,cyc+m+1);lim=n+2;
        for(int i=1;i<=m;i++)tag[cyc[i]]=1;
        for(int i=1;i<=m;i++)dfs(cyc[i]);
        for(int i=0;i<=lim;i++){
            val[i][0]=0;
            for(int j=1;j<=m;j++)
                val[i][j]=val[i][j-1]+dp[cyc[j]][i];
        }
        ll ans=dp[s][2];
        for(int i=0;i<=n;i++){
            ll pre=0;
            for(int j=1;j<=m;j++){
                pre=max(pre,val[i+2][j-1]-val[i+1][j-1]);
                ans=max(ans,pre+val[i+1][j]-val[i][j]+val[i][m]);
            }
            if(i>=1)ans=max(ans,val[i][m]);
        }
        printf("%lld\n",ans+p[s]);
    }
    return 0;
}