Problem:
Problem E: こたつがめを燃やさないで

Problem

こたつがめは甲羅がこたつの亀です。

こたつがめはぽきた(起床する事)途端、格子状の区画に区切られているステージに閉じ込められました。このステージには、ただ一つの出口があります。
区画は道、塀または爆弾のいずれかであり、こたつがめは上下左右に移動することで、道または爆弾の区画に侵入することができます。
斜めには移動できません。また、塀の区画は通ることができません。ステージの外にも道の区画が無限に広がっていますが亀の天敵であるワニがいる為、通る事はできません。

こたつがめはさっさと帰宅してねね(睡眠する事)を行いたいので、このステージの出口に向かおうと思いました。

こたつの赤外線を自在に操る事ができるこたつがめは、自身の周囲8区画を同時に道に変える能力を持っています。
しかし、その8区画と自身がいるいずれかの区画が爆弾の区画の場合、誘爆してこたつがめは炎上し、ねねどころか有名SNS「カメったー」すら出来なくなってしまうので能力を発動する事ができません。

例を図示します。イラストが配置されていない区画とこたつがめがいる区画は道の区画です。
こたつがめの周囲8区画とはオレンジ色の区画を指します。左の図ではこたつがめの周囲8区画、および自身がいる区画に爆弾の区画は含まれていない為、能力を発動できます。
しかし、右の図では周囲8区画、および自身がいる区画のいずれかに爆弾の区画が含まれている為、能力を発動する事ができません。

移動するコストは1区画の移動ごとに$A$、能力を発動するコストは1回あたり$B$です。道に変えた塀あたりのコストではない事に注意してください。開始地点とステージの出口は道です。
開始地点からステージの出口に着くまでの移動するコストと能力を発動するコストの合計の最小値を出力してください。ただし、ステージの出口に到達できない場合は代わりに"INF"と出力してください。

Input

入力は以下の形式で与えられる。

$H$ $W$ $A$ $B$
$c_{1,1} \cdots c_{1,W}$
$\vdots$
$c_{H,1} \cdots c_{H,W}$

入力は$H+1$行からなる。

$1$行目には上下の長さを表す整数$H$、左右の長さを表す整数$W$、移動するコストを表す整数$A$、能力を発動するコストを表す整数$B$がそれぞれ空白区切りで与えられる。

$2$行目からの$H$行には、こたつがめが閉じ込められたステージの各区画における状態$c_{i,j}$が与えられる。

$c_{i,j}$はそれぞれ's','g','.','#','*'のいずれかからなり、区画$(i,j)$が下記のような状態であることを表す。

's' : その区画が開始地点であることを表す。

'g' : その区画がステージの出口であることを表す。

'.' : その区画が道であることを表す。

'#' : その区画が塀であることを表す。

'*' : その区画が爆弾であることを表す。

Constraints

入力は以下の条件を満たす。

$1 \le H,W,A,B \le 1000$

$3 \le H+W$

$c_{i,j} \in${'s','g','.','#','*'}

's','g'はそれぞれ1回ずつ出現する

Output

開始地点からステージの出口に着くまでの移動するコストと能力を発動するコストの合計の最小値を出力してください。ただし、ステージの出口に到達できない場合は代わりに"INF"と出力してください。

Sample input 1

4 4 1 1
g#..
#...
.*..
...s

Sample output 1

7

Sample input 2

4 4 1 1
g#..
#*..
....
...s

Sample output 2

INF

Sample input 3

2 4 1 1
###g
s###

Sample output 3

6

Sample input 4

3 3 1 10
g..
##.
s..

Sample output 4

6

Sample input 5

3 3 10 1
g..
##.
s..

Sample output 5

21

Solution in Rust (judge verdict: Accepted)

#[allow(unused_imports)]
use std::cmp::*;
#[allow(unused_imports)]
use std::collections::*;
use std::io::{Write, BufWriter};
// https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8
macro_rules! input {
    ($($r:tt)*) => {
        let stdin = std::io::stdin();
        let mut bytes = std::io::Read::bytes(std::io::BufReader::new(stdin.lock()));
        let mut next = move || -> String{
            bytes
                .by_ref()
                .map(|r|r.unwrap() as char)
                .skip_while(|c|c.is_whitespace())
                .take_while(|c|!c.is_whitespace())
                .collect()
        };
        input_inner!{next, $($r)*}
    };
}

macro_rules! input_inner {
    ($next:expr) => {};
    ($next:expr, ) => {};

    ($next:expr, $var:ident : $t:tt $($r:tt)*) => {
        let $var = read_value!($next, $t);
        input_inner!{$next $($r)*}
    };
}

macro_rules! read_value {
    ($next:expr, ( $($t:tt),* )) => {
        ( $(read_value!($next, $t)),* )
    };

    ($next:expr, [ $t:tt ; $len:expr ]) => {
        (0..$len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()
    };

    ($next:expr, chars) => {
        read_value!($next, String).chars().collect::<Vec<char>>()
    };

    ($next:expr, usize1) => {
        read_value!($next, usize) - 1
    };

    ($next:expr, [ $t:tt ]) => {{
        let len = read_value!($next, usize);
        (0..len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()
    }};

    ($next:expr, $t:ty) => {
        $next().parse::<$t>().expect("Parse error")
    };
}

/*
 * Dijkstra's algorithm.
 * Verified by: AtCoder ABC035 (http://abc035.contest.atcoder.jp/submissions/676539)
 */

struct Dijkstra {
    edges: Vec<Vec<(usize, i64)>>, // adjacent list representation
}

/*
 * Code from https://doc.rust-lang.org/std/collections/binary_heap/
 */
#[derive(Copy, Clone, Eq, PartialEq)]
struct State {
    cost: i64,
    position: usize,
}

// The priority queue depends on `Ord`.
// Explicitly implement the trait so the queue becomes a min-heap
// instead of a max-heap.
impl Ord for State {
    fn cmp(&self, other: &State) -> Ordering {
        // Notice that the we flip the ordering here
        match other.cost.cmp(&self.cost) {
            std::cmp::Ordering::Equal => other.position.cmp(&self.position),
            x => x,
        }
    }
}

// `PartialOrd` needs to be implemented as well.
impl PartialOrd for State {
    fn partial_cmp(&self, other: &State) -> Option<Ordering> {
        Some(self.cmp(other))
    }
}


impl Dijkstra {
    fn new(n: usize) -> Self {
        Dijkstra { edges: vec![Vec::new(); n] }
    }
    fn add_edge(&mut self, from: usize, to: usize, cost: i64) {
        self.edges[from].push((to, cost));
    }
    /*
     * This function returns a Vec consisting of the distances from vertex source.
     */
    fn solve(&self, source: usize, inf: i64) -> Vec<i64> {
        let n = self.edges.len();
        let mut d = vec![inf; n];
        let mut que = std::collections::BinaryHeap::new();
        que.push(State {cost: 0, position: source});
        while let Some(State {cost, position: pos}) = que.pop() {
            if d[pos] <= cost {
                continue;
            }
            d[pos] = cost;
            for adj in &self.edges[pos] {
                que.push(State {cost: cost + adj.1, position: adj.0});
            }
        }
        return d;
    }
}

fn solve() {
    let out = std::io::stdout();
    let mut out = BufWriter::new(out.lock());
    macro_rules! puts {
        ($($format:tt)*) => (write!(out,$($format)*).unwrap());
    }
    input! {
        h: usize,
        w: usize,
        a: i64,
        b: i64,
        c: [chars; h],
    }
    let mut dijk = Dijkstra::new(2 * h * w);
    let dx = [0, 1, 0, -1];
    let dy = [-1, 0, 1, 0];
    let tdx = [1, 1, -1, -1];
    let tdy = [-1, 1, -1, 1];
    for i in 0..h {
        for j in 0..w {
            for d in 0..4 { 
                let nx = i as i32 + dx[d];
                let ny = j as i32 + dy[d];
                if nx < 0 || nx >= h as i32 || ny < 0 || ny >= w as i32 {
                    continue;
                }
                let nx = nx as usize;
                let ny = ny as usize;
                if c[nx][ny] == '#' {
                } else {
                    dijk.add_edge(i * w + j, nx * w + ny, a);
                }
                dijk.add_edge(h * w + i * w + j, nx * w + ny, a);
            }
            for d in 0..4 { 
                let nx = i as i32 + tdx[d];
                let ny = j as i32 + tdy[d];
                if nx < 0 || nx >= h as i32 || ny < 0 || ny >= w as i32 {
                    continue;
                }
                let nx = nx as usize;
                let ny = ny as usize;
                dijk.add_edge(h * w + i * w + j, nx * w + ny, 2 * a);
            }
            let mut nearbomb = false;
            for dx in -1..2 {
                for dy in -1..2 {
                    let nx = i as i32 + dx;
                    let ny = j as i32 + dy;
                    if nx < 0 || nx >= h as i32 || ny < 0 || ny >= w as i32 {
                        continue;
                    }
                    let nx = nx as usize;
                    let ny = ny as usize;
                    if c[nx][ny] == '*' {
                        nearbomb = true;
                        break;
                    }
                }
            }
            if !nearbomb {
                dijk.add_edge(i * w + j, h * w + i * w + j, b);
            }
        }
    }
    let mut s = 0;
    let mut g = 0;
    for i in 0..h {
        for j in 0..w {
            if c[i][j] == 's' {
                s = i * w + j;
            }
            if c[i][j] == 'g' {
                g = i * w + j;
            }
        }
    }
    const INF: i64 = 1 << 58;
    let d = dijk.solve(s, INF)[g];
    //let d = min(d, dijk.solve(s, INF)[h * w + g]);
    if d == INF {
        puts!("INF\n");
    } else {
        puts!("{}\n", d);
    }
}

fn main() {
    solve();
}