/*
* John Wick is a man of
* focus, commitment & sheer will
*/
#include<bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp> 
#include <ext/pb_ds/tree_policy.hpp>  
#include <functional>

using namespace __gnu_pbds;
using namespace std;

// -----------------<fast IO>------------------
#define thisIsSpeed		ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#define FLSH 			fflush(stdout)
#define fileIO(name) \
    freopen(name".in", "r", stdin); \
    freopen(name".out", "w", stdout);
#define PRECISION(x) 	cout << fixed << setprecision(x);
#define endl 			"\n"
// -----------------<Container manipulation / traversal macros>------------------
#define pb	 			push_back
#define sz(x) 			((int)((x).size()))
#define all(x) 			(x).begin(),(x).end()
#define rep(i,n) 		for(int (i)=0;(i)<(n);(i)++)
#define forn(i,a,b) 	for(int (i)=(a);(i)<=(b);(i)++)
#define ford(i,a,b) 	for(int (i)=(b);(i)>=(a);(i)--)
#define fill(l,a)		memset((l) , (a) , sizeof(l))
#define ss				second
#define ff 				first
// -----------------<Some common useful functions>------------------
#define checkbit(n,b)   ( (n >> b) & 1)
#define in(b,a)         ( (b).find(a) != (b).end())
// -----------------<TypeDef Start>------------------
typedef long double ld;
typedef long long ll;
typedef unsigned long long ull;
typedef vector<ll> vi;
typedef map<ll,ll> mii;
// --------------<Declaring ordered_set>---------------------- 
/* tree <data_type , null_type(set) or mapped_type(map) , 
 * 		comparison basis for two fncs , tree_type , metadata>
 * find_by_order(k):iterator to kth element
 * order_of_key(k):# of items less than k
*/
typedef tree<long long int, null_type, less<long long int>, rb_tree_tag,  
            tree_order_statistics_node_update>  
    ordered_set;  
// -----------------<Constants>------------------
const int MAXN = 1000005;
const int SQRTN = 1003;
const int LOGN = 22;
const double PI=acos(-1);
const int INF=1e9;
const int MOD = 1000000009;
const int FMOD = 998244353;
const double eps = 1e-9;
// -----------------<Modular Arithmetic>------------------
template<typename T> T gcd(T a, T b)
	{return(b?__gcd(a,b):a);}
template<typename T> T lcm(T a, T b)
	{return(a*(b/gcd(a,b)));}
template<typename T> T add(T a, T b, T c = MOD)
	{T res=a+b;return(res>=c?res-c:res);}
template<typename T> T mod_neg(T a, T b, T c = MOD)
	{T res;if(abs(a-b)<c)res=a-b;else res=(a-b)%c;return(res<0?res+c:res);}
template<typename T> T mul(T a, T b, T c = MOD)
	{ll res=(ll)a*b;return(res>=c?res%c:res);}
inline ll mulmod(ll a,ll b, ll m = MOD)
	{ll q = (ll)(((ld)a*(ld)b)/(ld)m);ll r=a*b-q*m;if(r>m)r%=m;if(r<0)r+=m;return r;}
template<typename T>T expo(T e, T n)
	{T x=1,p=e;while(n){if(n&1)x=x*p;p=p*p;n>>=1;}return x;}
template<typename T>T power(T e, T n, T m = MOD)
	{T x=1,p=e;while(n){if(n&1)x=mul(x,p,m);p=mul(p,p,m);n>>=1;}return x;}
template<typename T>T extended_euclid(T a, T b, T &x, T &y)
	{T xx=0,yy=1;y=0;x=1;while(b){T q=a/b,t=b;b=a%b;a=t;t=xx;xx=x-q*xx;x=t;t=yy;yy=y-q*yy;y=t;}return a;}
template<typename T>T mod_inverse(T a, T n = MOD)
	{T x,y,z=0;T d=extended_euclid(a,n,x,y);return(d>1?-1:mod_neg(x,z,n));}
// ------------------nCr---------------------------
/*
ll modInverse(ll n, ll p){return power(n, p-2);}
ll nCr(ll n, ll r, ll p=MOD){if (r==0) return 1;ll fac[n+1];fac[0]=1;for(lli=1;i<=n;i++)fac[i]=fac[i-1]*i%p;
	return (fac[n]*modInverse(fac[r],p)%p*modInverse(fac[n-r], p)%p)%p;}
*/
//-----------------sieve--------------------------
/*
//bool C[MAXN];
ll C[MAXN]={0};
vi V1;
void sieve(){
	ll ans=1;
    ll i,j;
    for(i=2;i<MAXN;i++){if(!C[i] && i<20){
		ans*=i;
        // V1.pb(i);
        for(j=2*i;j<MAXN;j=j+i){C[j]=true;}}}
        cout<<ans<<endl;
}
*/  
//-----------------custom structuree----------------
/*
struct cus{
	ll x,y;
	//string str;
};
*/
//-----------------custom compare---------------------
/*
bool compare(cus i1, cus i2){ 
    if(i1.x != i2.x)
		return i1.x < i2.x;
	else return i1.y < i2.y;
}
*/
// --------------------------------------------------</TEMPLATE>--------------------------------------------------

ll par[5001],c[5001];
bool vis[5001];
vector<vi> cc(5001);

void dfsu(ll v, ll p){
	if(vis[v])return;
	vis[v]=1;
	cc[p].pb(c[par[v]]);
//	cout<<par[v]<<endl;
	dfsu(par[v],p);
}


void mainn(){
	ll z=0, k=0, t=0, n=0, x=0, m=0, curr=0, num=0, n1=0, n2=0, n3=0, n4=0, i=0 , j=0;
	ll root=0, sum=0, diff=0, q=0, choice=0, d=0, len=1, beg=0, end=0, pos=0, cnt=0;
	ll lo=0, hi=0, mid=0, ans=0;
	bool flag = false;
	string s1, s2, s3, str;
	char ch='a', ch1, ch2, ch3, *ptr;
	double dub=0;
	////////////////////////////////////////////////////////////////////////////////////////////////////////////
	cin>>n>>k;
	fill(vis,0);
	rep(i,n){
		cin>>par[i+1];
		//cout<<i<<endl;
	}
	rep(i,n){
		cin>>c[i+1];
	//	cout<<i<<endl;
	}
//	cout<<"here";
	forn(i,1,n){
		if(!vis[i]){
			dfsu(i,cnt);
			cnt++;
//			cout<<"here\n";
		}
	}
	i=0,pos=cnt,ans=LLONG_MIN;
	cnt=0;
	while(i<pos){
		n4=hi=LLONG_MIN;
		curr=sz(cc[i]);
		cnt=0;
		sum=0;
		rep(j,curr){
			cc[i].pb(cc[i][j]);
			sum+=cc[i][j];
			hi=max(hi,cc[i][j]);
		}
		
		n1=k%curr;
		n2=k/curr;
		
		if(n1==0)n2--,n1=curr;
		
		if(sum>0 && n2>0)cnt+=sum*n2;
		rep(j,curr){
			n3=0;
			rep(kk,n1){
				n3+=cc[i][j+kk];
				n4=max(n4,n3);
			}
		}
		if(n4>0 && n2>0 && sum>0)cnt+=n4;
		else if(n2==0)cnt+=n4;
		n4=LLONG_MIN;
		curr=min(curr,k);
		rep(j,curr){
			n3=0;
			rep(kk,curr){
				n3+=cc[i][j+kk];
				ans=max(ans,n3);
			}
		}
		ans=max(ans,cnt);
		//cout<<i<<" "<<ans<<endl;
		i++;
	}
	
	cout<<ans<<endl;
}

int main(){
	thisIsSpeed;
	int t=1;
//	cin>>t;
	while(t--){
		mainn();
	}
	return 0;
}


Measured memory_bytes: 3956