const INF: usize = 1usize << 60; struct UnionFind { n: usize, parents: Vec, depths: Vec, chains: Vec, } impl UnionFind { fn new(n: usize) -> Self { UnionFind { n: n, parents: (0..n).collect(), depths: vec![0; n], chains: vec![1; n], } } fn equiv(&mut self, a: usize, b: usize) -> bool { self.find(a) == self.find(b) } fn unite(&mut self, a: usize, b: usize) { if self.equiv(a, b) { return; } let x = self.parents[a]; let y = self.parents[b]; if self.depths[x] > self.depths[y] { self.parents[x] = self.parents[y]; self.chains[y] += self.chains[x]; } else { self.parents[y] = self.parents[x]; self.chains[x] += self.chains[y]; if self.depths[x] == self.depths[y] { self.depths[x] += 1; } } } fn find(&mut self, a: usize) -> usize { if self.parents[a] == a { return a; } let p = self.find(self.parents[a]); self.parents[a] = p; p } fn chain_cnt(&mut self, i: usize) -> usize { let idx = self.find(i); self.chains[idx] } } fn combinations(n: usize, r: usize) -> Vec> { fn dfs(cur: usize, limit: usize, size_limit: usize, stack: &mut Vec, result: &mut Vec>) { if cur == limit { result.push(stack.clone().to_owned()); return; } if limit - cur > size_limit - stack.len() { dfs(cur+1, limit, size_limit, stack, result); } if size_limit > stack.len() { stack.push(cur); dfs(cur+1, limit, size_limit, stack, result); stack.pop(); } } let mut stack = vec![]; let mut result = vec![]; dfs(0, n, r, &mut stack, &mut result); result } fn main() { let mut nmk = String::new(); std::io::stdin().read_line(&mut nmk).ok(); let nmk: Vec = nmk.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let n = nmk[0]; let m = nmk[1]; let k = nmk[2]; let mut a = String::new(); std::io::stdin().read_line(&mut a).ok(); let a: Vec = a.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let mut paths = vec![vec![INF; n]; n]; for _ in 0..m { let mut temp = String::new(); std::io::stdin().read_line(&mut temp).ok(); let temp: Vec = temp.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let x = temp[0]-1; let y = temp[1]-1; let z = temp[2]; paths[x][y] = z; paths[y][x] = z; } for k in 0..n { for i in 0..n { for j in 0..n { paths[i][j] = paths[i][j].min(paths[i][k] + paths[k][j]); } } } let mut result = INF; for combs in combinations(n, k).iter() { let mut uf = UnionFind::new(k); let mut idxmap = vec![n; n]; let mut costs = 0usize; combs.iter().enumerate().for_each(|(idx, &i)| { idxmap[i] = idx; costs += a[i]; }); let mut lines = vec![]; for i in 0..k { for j in i+1..k { if idxmap[combs[i]] == n || idxmap[combs[j]] == n { continue; } lines.push((combs[i], combs[j], paths[combs[i]][combs[j]])); } } lines.sort_by_key(|&(_, _, cost)| cost); for &(u, v, cost) in lines.iter() { let u = idxmap[u]; let v = idxmap[v]; if u == n || v == n { continue; } if uf.find(u) == uf.find(v) { continue; } uf.unite(u, v); costs += cost; } result = result.min(costs); } println!("{}", result); }