結果

問題 No.124 門松列(3)
コンテスト
ユーザー koyumeishi
提出日時 2015-03-12 02:35:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 9 ms / 5,000 ms
コード長 3,058 bytes
コンパイル時間 861 ms
コンパイル使用メモリ 87,464 KB
実行使用メモリ 5,504 KB
最終ジャッジ日時 2024-06-24 17:18:53
合計ジャッジ時間 1,849 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <cstdio>
#include <sstream>
#include <map>
#include <string>
#include <algorithm>
#include <queue>
#include <cmath>
#include <deque>

using namespace std;

typedef struct{
	int to;
	int cost;
}edge;


//using some techniques writen in wikipedia, SLF and LLL
void SPFA(vector<vector<edge>> &G, vector<int> &dist, int s){
	const int INF = 1e8;
	fill(dist.begin(), dist.end(), INF);

	vector<int> cnt(dist.size(), 0);

	dist[s] = 0;
	deque<int> Q;

	long long sum = 0;

	auto my_push = [&](int node){
		//SLF
		if(Q.size() == 0 || dist[Q.front()] > dist[node]){
			Q.push_front(node);
		}else{
			Q.push_back(node);
		}
		cnt[node]++;
	};
	auto my_pop = [&]() -> int{
		int ret = Q.front();
		cnt[ret]--;
		Q.pop_front();
		return ret;
	};

	auto adjust = [&](){
		//LLL
		double ave = 1.0*sum / Q.size();
		if(Q.size() <= 1) return;
		while( dist[Q.front()] > ave ){
			Q.push_back( Q.front() );
			Q.pop_front();
		}
	};
	
	my_push(s);

	while(!Q.empty()){
		int pos = my_pop();
		sum -= dist[pos];

		for(edge& E : G[pos]){
			if(dist[ E.to ] > dist[ pos ] + E.cost ){

				if(cnt[E.to] != 0){
					sum -= dist[ E.to ];
				}

				dist[ E.to ] = dist[ pos ] + E.cost;
				
				sum += dist[E.to];

				if(cnt[ E.to ] == 0){
					my_push( E.to );

					adjust();
				}
			}

		}
		
	}

}
 
void add_edge(vector<vector<edge> > &G, int from, int to, int cost){
	G[from].push_back( (edge){to, cost} );
	//G[to].push_back( (edge){from, cost} );
}


bool is_kadomatsu_sec(int a, int b, int c){
	if(a==b || b==c || c==a) return false;
	if(b<a && b<c) return true;
	if(b>a && b>c) return true;
	return false;
}

int x[] = {0,1,0,-1};
int y[] = {-1,0,1,0};

int main(){
	int w,h;
	cin >> w >> h;
	vector<vector<int> > v(h, vector<int>(w));
	for(int i=0; i<h; i++){
		for(int j=0; j<w; j++){
			cin >> v[i][j];
		}
	}

	vector<vector<edge> > G(w*h*4);
	for(int i=0; i<h; i++){
		for(int j=0; j<w; j++){
			
			for(int k=0; k<4; k++){
				int xx = j+x[k];
				int yy = i+y[k];
				if(xx<0 || xx>=w) continue;
				if(yy<0 || yy>=h) continue;

				for(int l=0; l<4; l++){
					if(l==k) continue;
					int xxx = j+x[l];
					int yyy = i+y[l];
					if(xxx<0 || xxx>=w) continue;
					if(yyy<0 || yyy>=h) continue;

					if(is_kadomatsu_sec(v[yy][xx], v[i][j], v[yyy][xxx])){
						add_edge(G, i*w+j + w*h*k, yyy*w+xxx + w*h*((l+2)%4), 1);
					}
				}
				
			}
			
		}
	}

	const int INF = 10000000;

	int ans = INF;

	vector<int> dist(w*h*4);
	
	for(int i=0; i<4; i++){
		int x1 = 0+x[i];
		int y1 = 0+y[i];
		if(x1<0 || x1>=w) continue;
		if(y1<0 || y1>=h) continue;
		
		for(int j=0; j<4; j++){
			int x2 = x1+x[j];
			int y2 = y1+y[j];
			if(x2<0 || x2>=w) continue;
			if(y2<0 || y2>=h) continue;
			if(x2==0 && y2==0) continue;
			
			
			if(is_kadomatsu_sec(v[0][0], v[y1][x1], v[y2][x2])){
				
				SPFA(G, dist, y2*w+x2 + w*h*((j+2)%4));
				for(int k=0; k<4; k++){
					ans = min(ans, dist[ (h-1)*w+(w-1) + w*h*k ]+2);
				}
			}
			
		}
	}

	if(ans < INF) cout << ans << endl;
	else cout << -1 << endl;
	return 0;
}
0