結果

問題 No.20 砂漠のオアシス
ユーザー Darjeeling0Darjeeling0
提出日時 2020-03-04 00:42:53
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 71 ms / 5,000 ms
コード長 3,184 bytes
コンパイル時間 2,038 ms
コンパイル使用メモリ 186,264 KB
実行使用メモリ 5,660 KB
最終ジャッジ日時 2023-08-04 02:20:21
合計ジャッジ時間 3,510 ms
ジャッジサーバーID
(参考情報)
judge15 / judge12
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,384 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 3 ms
4,380 KB
testcase_04 AC 3 ms
4,384 KB
testcase_05 AC 60 ms
5,640 KB
testcase_06 AC 71 ms
5,576 KB
testcase_07 AC 70 ms
5,660 KB
testcase_08 AC 70 ms
5,652 KB
testcase_09 AC 70 ms
5,588 KB
testcase_10 AC 1 ms
4,380 KB
testcase_11 AC 1 ms
4,376 KB
testcase_12 AC 4 ms
4,376 KB
testcase_13 AC 5 ms
4,376 KB
testcase_14 AC 7 ms
4,380 KB
testcase_15 AC 6 ms
4,380 KB
testcase_16 AC 15 ms
4,380 KB
testcase_17 AC 11 ms
4,376 KB
testcase_18 AC 13 ms
4,380 KB
testcase_19 AC 15 ms
4,376 KB
testcase_20 AC 3 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#define REP(i, n) for (int i = 0; i < (int)(n); i++)
#define REPR(i, n) for(int i = (int)(n); i >= 0; i--)
#define FOR(i, s, n) for (int i = (s); i < (int)(n); i++)
#define ALL(v) v.begin(), v.iner()
using namespace std;
typedef long long ll;

ll N, K, W, H, M, C, V;

const int MAX = 510000;
const int MOD = 1000000007;

long long fac[MAX], finv[MAX], inv[MAX];

vector<double> cx, cy, c;
 ll MM = 1e12;

void COMinit() {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (int i = 2; i < MAX; i++){
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}

long long COM(int n, int k){
    if (n < k) return 0;
    if (n < 0 || k < 0) return 0;
    return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}

std::vector<int> Eratosthenes(int n ){
    std::vector<bool> is_prime( n + 1 );
    for( int i = 0; i <= n; i++ ) is_prime[ i ] = true;
    std::vector<int> P;
    for( int i = 2; i <= n; i++ ){
        if( is_prime[ i ] ){
            for( int j = 2 * i; j <= n; j += i ){
                is_prime[ j ] = false;
            }
            P.emplace_back( i );
        }
    }
    return P;
}

vector<pair<int, int>> edge(pair<int, int> u){
    vector<pair<int, int>> ret;
    int i = u.first, j = u.second;
    if(i > 0) {ret.push_back(make_pair(i-1, j));}
    if(j > 0) {ret.push_back(make_pair(i, j-1));}
    if(i < N-1) {ret.push_back(make_pair(i + 1, j));}
    if(j < N-1) {ret.push_back(make_pair(i , j + 1));}
    return ret;
}

ll Dijkstra(pair<int, int> init, pair<int, int> iner, vector<vector<ll>> dist){
    vector<vector<ll>> cost(N, vector<ll> (N, MM));
    cost[init.first][init.second] = 0;
    priority_queue<pair<ll, pair<int, int>>> que;
    que.push(make_pair( 0, make_pair(init.first, init.second) ));
    REP(i, N){
        REP(j ,N){
            if( i == init.first and j == init.second) continue;
            que.push(make_pair(-MM, make_pair(i,j)));
        }
    }

    while(!que.empty()){
        pair<int, int> u = que.top().second;
        que.pop();
        for(auto v: edge(u)){

            ll cnd =  cost[u.first][u.second] + dist[v.first][v.second];

            if(cost[v.first][v.second]  > cnd){
                cost[v.first][v.second] = cnd;
                que.push(make_pair(-cnd , v));
                //cout << cnd << " ";
            }
        }
    }
   // cout << cost[iner.first][iner.second] << " ";
    return cost[iner.first][iner.second];
}

int main(){
    cin >> N >> V;
    vector<vector<ll>> dist(N, vector<ll> (N, MM));
    int ox, oy;
    cin >> oy >> ox;
    ox--; oy--;
    REP(i, N){
        REP(j, N){
            cin >> dist[i][j];
        }
    }
    if(Dijkstra(make_pair(0,0) , make_pair(N-1, N-1), dist)   < V) {cout << "YES"; return 0;}
    if(ox >= 0 and oy >= 0){
        ll to_oasis = Dijkstra(make_pair(0,0) , make_pair(ox, oy), dist) ;
        ll from_oasis = Dijkstra(make_pair(ox, oy), make_pair(N-1, N-1), dist) ;
      //  cout << 2*(V- to_oasis) << " " ;
        if(from_oasis < 2*(V- to_oasis) ) {cout << "YES";return 0;}
    }
    cout << "NO";
    return 0;
}
0