#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using namespace atcoder; typedef long long ll; #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define repr(i, n) for (int i = (int)(n) - 1; i >= 0; i--) #define repk(i, k, n) for (int i = k; i < (int)(n); i++) #define all(v) v.begin(), v.end() #define mod1 1000000007 #define mod2 998244353 #define mod3 100000007 #define vi vector #define vs vector #define vc vector #define vl vector #define vb vector #define vvi vector> #define vvc vector> #define vvl vector> #define vvb vector> #define vvvi vector>> #define vvvl vector>> #define pii pair #define pil pair #define pli pair #define pll pair #define vpii vector> #define vpll vector> #define vvpii vector>> #define vvpll vector>> template void debug(T e) { cerr << e << endl; } template void debug(vector &v) { rep(i, v.size()) { cerr << v[i] << " "; } cerr << endl; } template void debug(vector> &v) { rep(i, v.size()) { rep(j, v[i].size()) { cerr << v[i][j] << " "; } cerr << endl; } } template void debug(vector> &v) { rep(i, v.size()) { cerr << v[i].first << " " << v[i].second << endl; } } template void debug(set &st) { for (auto itr = st.begin(); itr != st.end(); itr++) { cerr << *itr << " "; } cerr << endl; } template void debug(multiset &ms) { for (auto itr = ms.begin(); itr != ms.end(); itr++) { cerr << *itr << " "; } cerr << endl; } template void debug(map &mp) { for (auto itr = mp.begin(); itr != mp.end(); itr++) { cerr << itr->first << " " << itr->second << endl; } } void debug_out() { cerr << endl; } template void debug_out(Head H, Tail... T) { cerr << H << " "; debug_out(T...); } using mint = modint1000000007; void debug_mint1(vector &vec) { for (int i = 0; i < vec.size(); i++) { cerr << vec[i].val() << " "; } cerr << endl; } void debug_mint2(vector> &vec) { for (int i = 0; i < vec.size(); i++) { for (int j = 0; j < vec[i].size(); j++) { cerr << vec[i][j].val() << " "; } cerr << endl; } } void dfs(vector &visited, vector> &routes, vector &now_rt, ll N, ll M){ if (now_rt.size() == M){ routes.push_back(now_rt); return; } for (ll i = 0; i < N; i++){ if (!visited[i]){ visited[i] = true; now_rt.push_back(i); dfs(visited, routes, now_rt, N, M); now_rt.pop_back(); visited[i] = false; } } return; } int main() { ll N, M, K; cin >> N >> M >> K; vector A(K); vector can(N, false); rep(i, K){ cin >> A[i]; A[i]--; can[A[i]] = true; } vector> T(N, vector(N)); rep(i, N){ rep(j, N){ cin >> T[i][j]; } } vector visited(N, false); vector> routes(0, vector(0)); vector now_rt; dfs(visited, routes, now_rt, N, M); ll INF = 10000000000000023LL; ll ans = INF; for (ll i = 0; i < routes.size(); i++){ ll anssub = 0; for (ll j = 0; j < routes[i].size() - 1; j++){ anssub += T[routes[i][j]][routes[i][j + 1]]; } ll keihin_min = INF; for (ll j = 0; j < N; j++){ if (can[j]){ keihin_min = min(keihin_min, T[routes[i].back()][j]); } } ans = min(ans, anssub + keihin_min); } cout << ans << endl; }