#include<bits/stdc++.h>
using namespace std;
using ll=long long;
using ld=long double;
using pll=pair<ll,ll>;
using plll=pair<ll,pll>;
using graph=vector<vector<int>>;

const ll MOD=1000000007;
const int dx[4]={1,0,-1,0};
const int dy[4]={0,1,0,-1};

ll abs1(ll a){
  if(a>=0)return a;
  else return -a;
}

ll pow1(ll a,ll b,ll p=LLONG_MAX){
  ll ret=1;
  while(b!=0){
    if((b&1)==1){
      ret*=a;
    }
    a*=a;
    if(a>=p)a%=p;
    b=b>>1;
  }
  return ret;
}

int digit(ll a,ll b){
  int i=0;
  while(a!=0){
    a/=b;
    i++;
  }
  return i;
}

ll gcd1(ll a,ll b){
  if(a==0 && b==0)return 0;
  ll n;
  while(a%b!=0){
    n=a%b;
    a=b;
    b=n;
  }
  return b;
}

ll solveLinear(ll n,ll m,ll c){
  if(n==0 && m==0)return 0;
  ll a[100],b[100];
  a[0]=n;b[0]=m;
  int i=1;
  while(b[i-1]!=0){
    a[i]=b[i-1];
    b[i]=a[i-1]%b[i-1];
    i++;
  }
  n=i;
  if(c%a[n-1]!=0){
    return 0;
  }
  ll x=c/a[n-1],y=0;
  for(i=1;i<n;i++){
    ll cnt=x;
    x=y;
    y=cnt-(a[n-i-1]/b[n-i-1])*y;
  }
  return x;
}

ll modRev(ll a,ll p=MOD){
  a%=p;
  ll n=solveLinear(a,p,1);
  n%=p;
  if(n<0)return n+p;
  return n;
}

ll fact(ll n,ll p=MOD){
  if(n<0)return 0;
  int i;
  ll ret=1;
  for(i=0;i<n;i++){
    ret=ret*(i+1);
    if(ret>=p)ret%=p;
  }
  return ret;
}

ll perm(ll n,ll k,ll p=MOD){
  if(n<0 || k<0 || n<k)return 0;
  int i;
  ll ret=1;
  for(i=0;i<k;i++){
    ret=ret*(n-i);
    if(ret>=p)ret%=p;
  }
  return ret;
}

ll binom(ll n,ll k,ll p=MOD){
  if(n<0 || k<0 || n<k)return 0;
  int i;
  ll ret=1,fact=1;
  for(i=0;i<k;i++){
    ret=ret*(n-i);
    fact*=i+1;
    if(ret>=p)ret%=p;
    if(fact>=p)fact%=p;
  }
  return (ret*modRev(fact))%p;
}

void pushHeap(vector<int>& h,int n){
  int i=h.size();
  h.push_back(n);
  while(i!=0 && h[(i-1)/2]<h[i]){
    swap(h[i],h[(i-1)/2]);
    i=(i-1)/2;
  }
  h[i]=n;
}

int popHeap(vector<int>& h){
  int n=h.size()-1;
  if(n==-1)return 0;
  int i,ret=h[0];
  swap(h[0],h[n]);
  h.pop_back();
  for(i=0;i<n;){
    int j=2*i+1;
    if(j>=n)break;
    if(j!=n-1 && h[j]<h[j+1])j++;
    if(h[i]<h[j]){
      swap(h[i],h[j]);
    }
    i=j;
  }
  return ret;
}

vector<int> ver;
void dfs(const graph &g,int v){
  ver[v]=1;
  for(auto i:g[v]){
    if(ver[i]!=1)continue;
    dfs(g,i);
  }
}

int binary_search(vector<int> v,int key){
  int left=0,right=v.size()-1;
  if(v.empty())return -1;
  while(abs1(right-left)>1){
    int mid=left+(right-left)/2;
    if(v[mid]<key)left=mid;
    else right=mid;
  }
  return left;
}

int main(){
  int i,j,k,cnt=0;
  int n,m;
  cin>>n>>m;
  int d[m+1][m+1];
  for(i=1;i<=m;i++){
    for(j=1;j<=m;j++){
      cin>>d[i][j];
    }
  }
  vector<int> v0,v1,v2;
  for(i=1;i<=n;i++){
    for(j=1;j<=n;j++){
      int c;
      cin>>c;
      int p=(i+j)%3;
      if(p==0)v0.push_back(c);
      if(p==1)v1.push_back(c);
      if(p==2)v2.push_back(c);
    }
  }
  ll sum0[m+1],sum1[m+1],sum2[m+1];
  for(i=0;i<=m;i++){
    sum0[i]=0;
    sum1[i]=0;
    sum2[i]=0;
  }
  for(i=1;i<=m;i++){
    for(j=0;j<v0.size();j++){
      sum0[i]+=d[v0[j]][i];
    }
  }
  for(i=1;i<=m;i++){
    for(j=0;j<v1.size();j++){
      sum1[i]+=d[v1[j]][i];
    }
  }
  for(i=1;i<=m;i++){
    for(j=0;j<v2.size();j++){
      sum2[i]+=d[v2[j]][i];
    }
  }
  ll mi=LLONG_MAX;
  for(i=1;i<=m;i++){
    for(j=1;j<=m;j++){
      for(k=1;k<=m;k++){
        if(i!=j && j!=k && k!=i){
          mi=min(mi,sum0[i]+sum1[j]+sum2[k]);
        }
      }
    }
  }
  cout<<mi<<endl;
}

Measured memory_bytes: 1520