結果

問題 No.3653 Space-Time Courier
コンテスト
ユーザー t98slider
提出日時 2026-08-28 23:13:57
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 2,727 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,291 ms
コンパイル使用メモリ 349,744 KB
実行使用メモリ 6,272 KB
最終ジャッジ日時 2026-08-28 23:14:39
合計ジャッジ時間 41,366 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 22 WA * 5 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;

template<class T> istream& operator >> (istream& is, vector<T>& vec) {
    for(T& x : vec) is >> x;
    return is;
}

template<class T> ostream& operator << (ostream& os, const vector<T>& vec) {
    if(vec.empty()) return os;
    os << vec[0];
    for(auto it = vec.begin(); ++it != vec.end(); ) os << ' ' << *it;
    return os;
}

unsigned long long xor64() {
    static unsigned long long x
        = (unsigned long long)(chrono::duration_cast<chrono::nanoseconds>(
            chrono::high_resolution_clock::now().time_since_epoch()).count())
            * 10150724397891781847ULL;
    x ^= x << 7;
    return x ^= x >> 9;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	int n, m;
	cin >> n >> m;
	vector<int> p(n);
	for(auto &&v : p) cin >> v;
	vector<tuple<int,int,int>> E(m);
	for(auto &&[u, v, w] : E){
		cin >> u >> v >> w;
		u--, v--;
	}
    if(n <= 800){
        vector<vector<ll>> dp(n, vector<ll>(n, 1ll << 60));
        for(auto [u, v, w] : E) dp[u][v] = min(dp[u][v], (ll)w);
        for(int k = 0; k < n; k++){
            for(int i = 0; i < n; i++){
                for(int j = 0; j < n; j++){
                    dp[i][j] = min(dp[i][j], dp[i][k] + dp[k][j]);
                }
            }
        }
        ll mn = 1ll << 60;
        int ans = 0;
        for(int i = 0; i < n; i++){
            for(int j = 0; j < n; j++){
                if(i == j) continue;
                ll v = dp[i][j] + p[i] + p[j];
                if(v < mn){
                    mn = v;
                    ans = 1;
                }else if(v == mn) ans++;
            }
        }
        cout << mn << " " << ans << '\n';
        return 0;
    }
    for(int i = 0; i + 1 < m; i++){
        swap(E[i], E[xor64() % m]);
    }
	ll mn = 1ll << 60;
	vector<ll> dp(n);
    constexpr int BL = 2500;
    vector<bitset<BL>> dp2(n);
	for(int i = 0; i < n; i++){
        dp[i] = p[i];
        dp2[i][i] = 1;
    }
    vector<ll> ndp(n, 1ll << 60);
    vector<bitset<BL>> ndp2(n);
    for(auto [u, v, w] : E){
        if(dp[u] + w < ndp[v]){
            ndp[v] = dp[u] + w;
            ndp2[v].reset();
        }
        if(dp[u] + w == ndp[v]) ndp2[v] |= dp2[u];
    }
    swap(dp, ndp);
    swap(dp2, ndp2);
	for(int i = 0; i < 20 * n; i++){
		for(auto [u, v, w] : E){
			if(dp[u] + w < dp[v]){
                dp[v] = dp[u] + w;
                dp2[v].reset();
            }
            if(dp[u] + w == dp[v]) dp2[v] |= dp2[u];
		}
	}
    for(int i = 0; i < n; i++) mn = min(mn, dp[i] + p[i]);
    int ans = 0;
    for(int i = 0; i < n; i++){
        if(mn != dp[i] + p[i]) continue;
        ans += dp2[i].count();
    }
	cout << mn << " " << ans << '\n';
}
0