关闭

【bzoj4061】[Cerc2012]Farm and factory

876人阅读 评论(0) 收藏 举报
分类:

  这题不错
  反正我想了一天是没想出来hhhh
  后来看了糖老师的题解 恍然大悟
  设dis[i][j],dis[i][j]分别为原图和新图中i到j的最短路径。
  设新加的一个点为x
  若要使在新图中最短路不经过x,应要满足
  u,dis[u][x]>=max{dis[u][1]dis[x][1],dis[u][2]dis[x][2]}
  这个还是能想得到的
  得到了这个以后,因为要使dis[u][x]最小,故应该取等号。把这个式子写出来,就会是
  max{dis[u][1]dis[x][1],dis[u][2]dis[x][2]}
  如果把每个(dis[u][1],dis[u][2])看做一个点,答案就是是这些点到点(dis[x][1],dis[x][2])的切比雪夫距离和的最小值。
  于是就乱搞咯2333
  用dijkstra,时间复杂度O(m+nlogn)。亲测SPFA会TLE,不知道是不是姿势问题。

#include <bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,u)  for(int i=head[u];i;i=nxt[i])
#define cr(x) memset(x , 0 , sizeof x)
#define pii pair<ll , int>
#define mp make_pair
#define maxn 100003
#define maxm 600007

inline int rd() {
    char c = getchar();
    while (!isdigit(c)) c = getchar() ; int x = c - '0';
    while (isdigit(c = getchar())) x = x * 10 + c - '0';
    return x;
}

typedef long long ll;
typedef int arr_int[maxn];
typedef ll  arr_ll [maxn];
typedef int adj[maxm];

const ll inf = 1e12;

arr_int head , vis;
arr_ll  dis[2] , x , y;
adj to , val , nxt;
int n , m , ett;

priority_queue<pii , vector<pii> , greater<pii> > Q;

inline void ins(int u , int v , int w) {
    to[++ ett] = v , val[ett] = w , nxt[ett] = head[u] , head[u] = ett;
}

void input() {
    n = rd() , m = rd();
    cr(head) , ett = 0;
    rep (i , 1 , m) {
        int u = rd() , v = rd() , w = rd();
        ins(u , v , w) , ins(v , u , w);
    }
}

void dijk(int s) {
    int c = s - 1;
    rep (i , 1 , n) dis[c][i] = inf;
    cr(vis);
    dis[c][s] = 0;
    Q.push(mp(0 , s));
    while (!Q.empty()) {
        int u = Q.top().second ; Q.pop();
        if (vis[u]) continue;
        vis[u] = 1;
        fore (i , u) {
            int v = to[i] , w = val[i];
            if (dis[c][v] > dis[c][u] + w) {
                dis[c][v] = dis[c][u] + w;
                Q.push(mp(dis[c][v] , v));
            }
        }
    }
}

void solve() {
    dijk(1) , dijk(2);
    rep (i , 1 , n) {
        x[i] = dis[0][i] + dis[1][i];
        y[i] = dis[0][i] - dis[1][i];
    }
    sort(x + 1 , x + n + 1) , sort(y + 1 , y + n + 1);
    ll nx = x[(n + 1) / 2] , ny = y[(n + 1) / 2];
    double ans = 0;
    rep (i , 1 , n)
        ans += abs(x[i] - nx) + abs(y[i] - ny);
    printf("%.12lf\n" , ans / (2.0 * n));
}

int main() {
    #ifndef ONLINE_JUDGE
        freopen("data.txt" , "r" , stdin);
    #endif
    int T = rd();
    rep (i , 1 , T) {
        input();
        solve();
    }
    return 0;
}
0
0

猜你在找
【直播】机器学习&数据挖掘7周实训--韦玮
【套餐】系统集成项目管理工程师顺利通关--徐朋
【直播】3小时掌握Docker最佳实战-徐西宁
【套餐】机器学习系列套餐(算法+实战)--唐宇迪
【直播】计算机视觉原理及实战--屈教授
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之矩阵--黄博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之凸优化--马博士
【套餐】Javascript 设计模式实战--曾亮
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:76729次
    • 积分:2169
    • 等级:
    • 排名:第17409名
    • 原创:136篇
    • 转载:0篇
    • 译文:0篇
    • 评论:44条
    小伙伴
    大神犇
    最新评论