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 C (judge verdict: Accepted)

#include <stdio.h>
#include <stdlib.h>
#define ll long long
#define rep(i,l,r)for(ll i=(l);i<(r);i++)
#define INF ((1LL<<62)-(1LL<<31))

//辺の情報を個別に持つタイプ
typedef struct edge{int s,g,c;}E;
typedef struct graph{
	int vcnt,ecnt;
	E  e[12000010];//適宜変える(ecnt)
	int id[1000010];//適宜変える(vcnt)
}G;
G g;

int esort(const void*a,const void*b){
	E*p=(E*)a,*q=(E*)b;
	if((*p).s<(*q).s)return -1;
	if((*p).s>(*q).s)return  1;
	if((*p).g<(*q).g)return -1;
	return 1;
}

#define in(ii,jj)(0<=ii&&ii<n&&0<=jj&&jj<m)
char s[1010][1010];
int  d4[]={0,1,0,-1,0};
int dd4[]={1,1,-1,-1,1};
int start,goal;
void readgraph(){
	//適宜変える
	ll n,m,a,b;
	scanf("%lld%lld%lld%lld",&n,&m,&a,&b);
	rep(i,0,n)scanf(" %s",s[i]);
	ll ecnt=0;
	rep(i,0,n)rep(j,0,m){
		if(s[i][j]=='s')start=i*m+j;
		if(s[i][j]=='g')goal=i*m+j;
		rep(k,0,4){
			int ii=i+d4[k];
			int jj=j+d4[k+1];
			if(in(ii,jj)&&s[ii][jj]!='#'){
				g.e[ecnt].s= i*m+ j;
				g.e[ecnt].g=ii*m+jj;
				g.e[ecnt].c=a;
				ecnt++;
			}
		}
		int flag=s[i][j]!='*';
		rep(k,0,4){
			int ii=i+d4[k];
			int jj=j+d4[k+1];
			flag&=!in(ii,jj)||s[ii][jj]!='*';
			ii=i+dd4[k];
			jj=j+dd4[k+1];
			flag&=!in(ii,jj)||s[ii][jj]!='*';
		}
		if(flag){
			rep(k,0,4){
				int ii=i+d4[k];
				int jj=j+d4[k+1];
				if(in(ii,jj)){
					g.e[ecnt].s= i*m+ j;
					g.e[ecnt].g=ii*m+jj;
					g.e[ecnt].c=b+a;
					ecnt++;
				}
			}
			rep(k,0,4){
				int ii=i+dd4[k];
				int jj=j+dd4[k+1];
				if(in(ii,jj)){
					g.e[ecnt].s= i*m+ j;
					g.e[ecnt].g=ii*m+jj;
					g.e[ecnt].c=b+a+a;
					ecnt++;
				}
			}
		}
	}
	g.vcnt=n*m;
	g.ecnt=ecnt;
	qsort(g.e,g.ecnt,sizeof(E),esort);
//	qsort(g.e,g.ecnt,sizeof(E),csort);
//	makeinvedge();

	int p=0;
	rep(i,0,g.vcnt){
		while(p<g.ecnt&&g.e[p].s<i)p++;
		g.id[i]=p;
	}
	g.id[g.vcnt]=g.ecnt;//番兵
}


//プラキュー（二分ヒープ）（優先度変更ありバージョン）
ll heapN,heap[1<<20],heapinv[1<<20];
int PQhikaku(int i,int j);//jの方が優先度が高いならtrueを返す
void PQchange(int n);
void heap_utod(int n){
	if(2*n>heapN)return;
	int rflag=(2*n+1<=heapN)&&(PQhikaku(2*n,2*n+1));
	if(PQhikaku(n,2*n+rflag)){
		ll temp=heap[2*n+rflag];
		heap[2*n+rflag]=heap[n];
		heap[n]=temp;
		heapinv[heap[n]]=n;
		heapinv[heap[2*n+rflag]]=2*n+rflag;	
		heap_utod(2*n+rflag);
	}
}
void heap_dtou(int n){
	if(n==1||PQhikaku(n,n/2))return;
	ll temp=heap[n];
	heap[n]=heap[n/2];
	heap[n/2]=temp;
	heapinv[heap[n]]=n;
	heapinv[heap[n/2]]=n/2;
	heap_dtou(n/2);
}
ll PQpop(){
	ll rr=heap[1];
	heapinv[heap[1]]=0;
	heap[1]=heap[heapN--];
	heapinv[heap[1]]=1;
	heap_utod(1);
	return rr;
}
void PQpush(ll n){
	heap[++heapN]=n;
	heapinv[heap[heapN]]=heapN;
	heap_dtou(heapN);
}

//早い方のダイクストラ
//グラフと始点を引いて各点への最短距離・最短経路を返す
//プラキューが必要
//O((E+V)logV)
ll daikusutorappp[1<<21];
void daikusutora2(ll s){
	rep(i,0,g.vcnt)daikusutorappp[i]=i==s?0:INF;
	rep(i,0,g.vcnt)PQpush(i);
	while(heapN){
		ll mv=PQpop();
		for(ll t=g.id[mv];t<g.ecnt&&g.e[t].s==mv;t++){
			if(daikusutorappp[g.e[t].g]>daikusutorappp[mv]+g.e[t].c){
				daikusutorappp[g.e[t].g]=daikusutorappp[mv]+g.e[t].c;
				//daikusutorappp[n+g.e[t].g]=mv;//最短経路がいる場合はこれ
				PQchange(g.e[t].g);
			}
		}
	}
}
int PQhikaku(int i,int j){return daikusutorappp[heap[i]]>daikusutorappp[heap[j]];}
void PQchange(int n){if(heapinv[n])heap_dtou(heapinv[n]);}

int main(){
	readgraph();
	daikusutora2(start);
	printf(daikusutorappp[goal]==INF?"INF\n":"%lld\n",daikusutorappp[goal]);
}