算法竞赛模板

Development

题目大意是给一个无向图，只有500个点代表城市，一些双向边代表路，边权代表通行时间。
还有一些城市有飞机场，给定飞机场之间的通行时间都一样为t。
给q次询问，三种op,op=1,加边;op=2,加机场,op=3,查询所有城市两两通行的时间和，无法通行记为0
q=1000

数据范围很容易想到floyd O^2维护新加边，关键在于机场的维护
有暴力做法，每次新增机场相当于和所有的机场都连接一条边，然后跑单边的floyd更新，这样的复杂度是O^3会T
增加一个虚点，所有飞机场到虚点的时间为t,虚点到机场的时间为0，这样就可以统一成第一种更新方法了

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define int long long
const int inf = 1e18;
const int N = 510;
int f[N][N];
int n, m, k, q, t;
vector<int> pl;
int pp[N];
struct nd {
    int u, v;
};

signed main() {
    ios::sync_with_stdio(false), cin.tie(nullptr);
    cin >> n >> m;
    for (int i = 0; i <= n + 5; i++) {
        for (int j = 0; j <= n + 5; j++) {
            f[i][j] = inf;
        }
    }
    
    for (int i = 0; i <= n; i++) {
        f[i][i] = 0;
    }

    for (int i = 1; i <= m; i++) {
        int u, v, dis;
        cin >> u >> v >> dis;
        f[u][v] = min(f[u][v], dis);
        f[v][u] = f[u][v];
    }

    cin >> k >> t;
    for (int i = 1; i <= k; i++) {
        int x;
        cin >> x;
        f[x][0] = t;
        f[0][x] = 0;
    }
    for (int k = 0; k <= n; k++) {
        for (int i = 0; i <= n; i++) {
            for (int j = 0; j <= n; j++) {
                f[i][j] = min(f[i][j], f[i][k] + f[k][j]);
            }
        }
    }

    cin >> q;
    int op, x, y, dis;
    while (q--) {
        cin >> op;
        if (op == 1) {
            cin >> x >> y >> dis;
            if (f[x][y] <= dis) {
                continue;
            }
            f[x][y] = f[y][x] = dis;
            for (int i = 0; i <= n; i++) {
                for (int j = 0; j <= n; j++) {
                    f[i][j] = min(f[i][j], f[i][x] + f[y][j] + dis);
                    f[i][j] = min(f[i][j], f[i][y] + f[x][j] + dis);
                }
            }
        }
        else if (op == 2) {
            cin >> x;
            f[x][0] = t;
            f[0][x] = 0;
            for (int i = 0; i <= n; i++) {
                for (int j = 0; j <= n; j++) {
                    f[i][j] = min(f[i][j], f[i][x] + f[0][j] + t);
                    f[i][j] = min(f[i][j], f[i][0] + f[x][j]);
                }
            }
        }
        else {
            int ans = 0;
            for (int i = 1; i <= n; i++) {
                for (int j = 1; j <= n; j++) {
                    if (f[i][j] >= inf) continue;
                    ans += f[i][j];
                }
            }
            cout << ans << '\n';
        }
    }
    return 0;
}