kmjp's blog

競技プログラミング参加記です

FII Code 2020 Round #2 : E. Escaping Courses

ちょっと実装がめんどい問題。
https://csacademy.com/contest/fii-code-2020-round-2/task/escaping-courses/

問題

2次元グリッドが与えられる。
一部のマスは恒久的に通行不可である。
また、普段移動可能であってもある時間帯のみ移動不可になるケースがある。

上下左右の隣接マスに時間1で移動できるとき、左上マスから右下マスに至る最短時間を求めよ。

解法

移動不可の時間帯の件を除けば単なるダイクストラ法である。
ダイクストラ法に加えて、新規に通行不可になる時間と、通行不可が解除される時間を順次処理しよう。
通行不可になる時間帯になると、そのマスの最短到達時間がいったん無限大になり、以後移動不可になると考える。
また、通行不可が解除されたマスは、隣接マスがそれ以前にすでに到達可能であれば、解除の直後に再入場できると考えよう。

そうするとイベントの数はO(マス数+通行不可時間のイベント数)で収まる。

int H,W;
string C[2010];
ll D[2020][2020];
int R;
ll Y[101010],X[101010],S[101010],T[101010];
int NG[2020][2020];

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>H>>W;
	FOR(y,H) {
		cin>>C[y];
		FOR(x,W) if(C[y][x]=='#') NG[y][x]=1<<25;
	}
	
	priority_queue<pair<ll,int>> Q;
	cin>>R;
	FOR(i,R) {
		cin>>Y[i]>>X[i]>>S[i]>>T[i];
		Y[i]--,X[i]--;
		S[i]*=2;
		T[i]*=2;
		Q.push({-(S[i]-1),i+R});
		Q.push({-(T[i]-1),i});
	}
	FOR(y,H) FOR(x,W) D[y][x]=1L<<60;
	D[0][0]=0;
	Q.push({0,2*R});
	while(Q.size()) {
		ll now=-Q.top().first;
		int cur=Q.top().second;
		Q.pop();
		if(cur<R) {
			int cy=Y[cur];
			int cx=X[cur];
			NG[Y[cur]][X[cur]]--;
			int ok=0;
			FOR(i,4) {
				int dd[4]={1,0,-1,0};
				int ty=cy+dd[i];
				int tx=cx+dd[i^1];
				if(ty<0||ty>=H || tx<0||tx>=W||NG[ty][tx]>0) continue;
				if(D[ty][tx]<=now) ok=1;
			}
			//cout<<now<<" T "<<cur<<" "<<Y[cur]<<","<<X[cur]<<" "<<ok<<endl;
			if(ok) {
				D[cy][cx]=now+3;
				Q.push({-D[cy][cx],2*R+cy*2000+cx});
			}
		}
		else if(cur<2*R) {
			cur-=R;
			//cout<<now<<" S "<<cur<<" "<<Y[cur]<<","<<X[cur]<<endl;
			NG[Y[cur]][X[cur]]++;
			D[Y[cur]][X[cur]]=1LL<<60;
		}
		else {
			int cy=(cur-2*R)/2000;
			int cx=(cur-2*R)%2000;
			if(now!=D[cy][cx]) continue;
			//cout<<now<<"  "<<cy<<" "<<cx<<" "<<NG[cy][cx]<<endl;
			if(cy==H-1&&cx==W-1) {
				cout<<D[cy][cx]/2<<endl;
				return;
			}
			FOR(i,4) {
				int dd[4]={1,0,-1,0};
				int ty=cy+dd[i];
				int tx=cx+dd[i^1];
				if(ty<0||ty>=H || tx<0||tx>=W||NG[ty][tx]>0) continue;
				if(D[ty][tx]>D[cy][cx]+2) {
					D[ty][tx]=D[cy][cx]+2;
					Q.push({-D[ty][tx],2*R+ty*2000+tx});
				}
			}
			
		}
	}
	cout<<-1<<endl;
	
}

まとめ

わかってしまえばそんなにややこしくないんだけどね。