Problem:
Score : 400 points

Problem Statement

There are N towns in the State of Atcoder, connected by M bidirectional roads.

The i-th road connects Town A_i and B_i and has a length of C_i.

Joisino is visiting R towns in the state, r_1,r_2,..,r_R (not necessarily in this order).

She will fly to the first town she visits, and fly back from the last town she visits, but for the rest of the trip she will have to travel by road.

If she visits the towns in the order that minimizes the distance traveled by road, what will that distance be?

Constraints

2≤N≤200

1≤M≤N×(N-1)/2

2≤R≤min(8,N) (min(8,N) is the smaller of 8 and N.)

r_i≠r_j (i≠j)

1≤A_i,B_i≤N, A_i≠B_i

(A_i,B_i)≠(A_j,B_j),(A_i,B_i)≠(B_j,A_j) (i≠j)

1≤C_i≤100000

Every town can be reached from every town by road.

All input values are integers.

Input

Input is given from Standard Input in the following format:

N M R
r_1 ... r_R
A_1 B_1 C_1
:
A_M B_M C_M

Output

Print the distance traveled by road if Joisino visits the towns in the order that minimizes it.

Sample Input 1

3 3 3
1 2 3
1 2 1
2 3 1
3 1 4

Sample Output 1

2

For example, if she visits the towns in the order of 1, 2, 3, the distance traveled will be 2, which is the minimum possible.

Sample Input 2

3 3 2
1 3
2 3 2
1 3 6
1 2 2

Sample Output 2

4

The shortest distance between Towns 1 and 3 is 4. Thus, whether she visits Town 1 or 3 first, the distance traveled will be 4.

Sample Input 3

4 6 3
2 3 4
1 2 4
2 3 3
4 3 1
1 4 1
4 2 2
3 1 6

Sample Output 3

3

Solution in Rust (1.15.1) (judge verdict: Time Limit Exceeded)

#[allow(unused_imports)]
use std::cmp::{Ordering};

#[allow(unused_imports)]
use std::collections::{BTreeMap, BTreeSet, BinaryHeap, HashMap, HashSet, VecDeque};

#[allow(unused_imports)]
use std::io::{stdin, stdout, BufWriter, Write};

#[allow(unused_imports)]
use std::iter::FromIterator;

#[allow(unused_macros)]
macro_rules! min {
    ($x: expr) => ($x);
    ($x: expr, $($z: expr),+) => {{
        std::cmp::min($x, min!($($z),*))
    }}
}

#[allow(unused_macros)]
macro_rules! max {
    ($x: expr) => ($x);
    ($x: expr, $($z: expr),+) => {{
        std::cmp::max($x, max!($($z),*))
    }}
}

#[allow(unused_macros)]
macro_rules! get { 
    ([$t: ty]) => { 
        { 
            let mut line = String::new(); 
            stdin().read_line(&mut line).unwrap(); 
            line.split_whitespace().map(|t|t.parse::<$t>().unwrap()).collect::<Vec<_>>()
        }
    };
    ([$t: ty];$n: expr) => {
        (0..$n).map(|_|get!([$t])).collect::<Vec<_>>()
    };
    ($t: ty) => {
        {
            let mut line = String::new();
            stdin().read_line(&mut line).unwrap();
            line.trim().parse::<$t>().unwrap()
        }
    };
    ($($t: ty),*) => {
        { 
            let mut line = String::new();
            stdin().read_line(&mut line).unwrap();
            let mut iter = line.split_whitespace();
            ($(iter.next().unwrap().parse::<$t>().unwrap(),)*)
        }
    };
    ($t: ty; $n: expr) => {
        (0..$n).map(|_|get!($t)).collect::<Vec<_>>()
    };
    ($($t: ty),*; $n: expr) => {
        (0..$n).map(|_|get!($($t),*)).collect::<Vec<_>>()
    };
}

#[allow(unused_macros)]
macro_rules! debug {
    ($($a:expr),*) => {
        #[cfg(debug_assertions)]
        writeln!(&mut std::io::stderr(), concat!("[DEBUG] ", $(stringify!($a), "={:?} "),*), $($a),*).unwrap();
    }
}

const BIG_STACK_SIZE: bool = true;

#[allow(dead_code)]
fn main() {
    use std::thread;
    if BIG_STACK_SIZE {
        thread::Builder::new()
            .stack_size(32 * 1024 * 1024)
            .name("solve".into())
            .spawn(solve)
            .unwrap()
            .join()
            .unwrap();
    } else {
        solve();
    }
}

#[derive(Eq, PartialEq, Clone, Debug)]
pub struct Rev<T>(pub T);
impl<T: PartialOrd> PartialOrd for Rev<T> {
    fn partial_cmp(&self, other: &Rev<T>) -> Option<Ordering> {
        other.0.partial_cmp(&self.0)
    }
}
impl<T: Ord> Ord for Rev<T> {
    fn cmp(&self, other: &Rev<T>) -> Ordering {
        other.0.cmp(&self.0)
    }
}

#[derive(Clone)]
pub struct ElementSwaps {
    sdir: Vec<SizeDirection>,
    emit_reset: bool,
    swaps_made: usize,
}
impl ElementSwaps {
    pub fn new(length: usize) -> ElementSwaps {
        ElementSwaps {
            emit_reset: true,
            sdir: (0..length)
                .map(|i| SizeDirection {
                    size: i,
                    dir: Direction::Neg,
                })
            .collect(),
            swaps_made: 0,
        }
    }
}
#[derive(Copy, Clone)]
enum Direction {
    Pos,
    Neg,
}
#[derive(Copy, Clone)]
struct SizeDirection {
    size: usize,
    dir: Direction,
}
impl Iterator for ElementSwaps {
    type Item = (usize, usize);
    fn next(&mut self) -> Option<(usize, usize)> {
        fn new_pos_wrapping(i: usize, s: Direction) -> usize {
            i.wrapping_add(match s {
                Direction::Pos => 1,
                Direction::Neg => !0,
            })
        }
        fn new_pos(i: usize, s: Direction) -> usize {
            match s {
                Direction::Pos => i + 1,
                Direction::Neg => i - 1,
            }
        }
        let max = self
            .sdir
            .iter()
            .cloned()
            .enumerate()
            .filter(|&(i, sd)| {
                new_pos_wrapping(i, sd.dir) < self.sdir.len()
                    && self.sdir[new_pos(i, sd.dir)].size < sd.size
            })
        .max_by(|&(_, sd1), &(_, sd2)| sd1.size.cmp(&sd2.size));
        match max {
            Some((i, sd)) => {
                let j = new_pos(i, sd.dir);
                self.sdir.swap(i, j);
                for x in &mut self.sdir {
                    if x.size > sd.size {
                        x.dir = match x.dir {
                            Direction::Pos => Direction::Neg,
                            Direction::Neg => Direction::Pos,
                        };
                    }
                }
                self.swaps_made += 1;
                Some((i, j))
            }
            None => {
                if self.emit_reset {
                    self.emit_reset = false;
                    if self.sdir.len() > 1 {
                        self.swaps_made += 1;
                        Some((0, 1))
                    } else {
                        self.swaps_made += 1;
                        Some((0, 0))
                    }
                } else {
                    None
                }
            }
        }
    }
#[inline]
#[allow(dead_code)]
    fn size_hint(&self) -> (usize, Option<usize>) {
        let n: usize = (2..self.sdir.len() + 1).product();
        (n - self.swaps_made, Some(n - self.swaps_made))
    }
}
pub struct Permutations<T> {
    swaps: ElementSwaps,
    v: Vec<T>,
}
impl<T: Clone> Permutations<T> {
    pub fn new(v: Vec<T>) -> Self {
        Permutations {
            swaps: ElementSwaps::new(v.len()),
            v: v,
        }
    }
}
impl<T: Clone> Iterator for Permutations<T> {
    type Item = Vec<T>;
#[inline]
    fn next(&mut self) -> Option<Vec<T>> {
        match self.swaps.next() {
            None => None,
            Some((0, 0)) => Some(self.v.clone()),
            Some((a, b)) => {
                let elt = self.v.clone();
                self.v.swap(a, b);
                Some(elt)
            }
        }
    }
#[inline]
#[allow(dead_code)]
    fn size_hint(&self) -> (usize, Option<usize>) {
        self.swaps.size_hint()
    }
}

fn solve() {
    let (n,m,rl) = get!(usize, usize, usize);
    let r = get!([usize]);
    let abc = get!(usize, usize, usize; m);
    let mut vec = vec![Vec::new();n+1];
    for (a,b,c) in abc {
        vec[a].push((b,c));
        vec[b].push((a,c));
    }

    let perm = Permutations::new(r);
    let mut mn = std::usize::MAX;
    for r in perm {
        let mut costsum = 0;
        for i in 1..rl {
            let mut heap = BinaryHeap::new();
            let mut seen = BTreeSet::new();
            heap.push(Rev((0, r[i-1])));
            seen.insert(r[i-1]);
            let mut tcost = 0;
            while let Some(Rev((cost, node))) = heap.pop() {
                if node == r[i] {
                    tcost = cost;
                    break;
                }
                for &(next, ncost) in &vec[node] {
                    if seen.contains(&next) {
                        continue;
                    }
                    heap.push(Rev((cost+ncost, next)));
                }
            }
            costsum += tcost;
        }
        mn = min!(costsum, mn);
    }
    println!("{}", mn);
}