[bzoj4289]PA2012 Tax——优化建图+Dijkstra

题目大意:

给出一个N个点M条边的无向图,经过一个点的代价是进入和离开这个点的两条边的边权的较大值,求从起点1到点N的最小代价。起点的代价是离开起点的边的边权,终点的代价是进入终点的边的边权
N<=100000
M<=200000

思路:

不难想到用点代表边,边代表点,然后跑最短路。
但是如果有一个菊花图的话,那个菊花就是一个完全图,边数就会有 m2 m 2 级别。
考虑这个菊花图的复杂度究竟消耗在哪里,不发现对于原图的每一条边都要便历这个菊花中的所有的边。
考虑这一题的特殊性质,即两条边权取max。我们把这个菊花图的每一条边都按照权值排序并且分为入边和出边两种边,这个菊花图中心的那个点的权值是取决于入边和出边的,如果出边的权值比入边小,那么就是取入边的权值,如果出边的权值比入边大,每当这个入边在排完序的边表中依次切换一条出边,权值增大了相邻两条出边的差。
于是可以想到每一条入边只和它对应的那条出边连边,然后对于一个排完序的出边的集合,小的向大的的连它们两个之间的差,大的向小的连0。

#include<bits/stdc++.h>

#define REP(i,a,b) for(int i=a,i##_end_=b;i<=i##_end_;++i)
#define pii pair<ll,int>
#define fi first
#define se second
#define mk make_pair
typedef long long ll;

using namespace std;

void File(){
    freopen("bzoj4289.in","r",stdin);
    freopen("bzoj4289.out","w",stdout);
}

template<typename T>void read(T &_){
    T __=0,mul=1; char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')mul=-1;
        ch=getchar();
    }
    while(isdigit(ch))__=(__<<1)+(__<<3)+(ch^'0'),ch=getchar();
    _=__*mul;
}

const int maxn=1e5+10;
const int maxm=2e5+10;
int n,m,cnt=1,ss,tt;
vector<pii>G[maxn];
int beg[maxm*2],to[maxm*6],las[maxm*6],cnte=1;
ll w[maxm*6],dis[maxm*2];

void add(int u,int v,ll val){las[++cnte]=beg[u];beg[u]=cnte;to[cnte]=v;w[cnte]=val;}

struct node{
    int u;
    ll d;
    bool operator < (const node & tm) const {
        return d>tm.d;
    }
};
priority_queue<node>qu;

void dijkstra(){
    memset(dis,63,sizeof(dis));
    dis[ss]=0;
    qu.push((node){ss,0});
    while(!qu.empty()){
        int u=qu.top().u;
        ll d=qu.top().d;
        qu.pop();
        if(d!=dis[u])continue;
        for(int i=beg[u];i;i=las[i]){
            int v=to[i];
            if(d+w[i]<dis[v]){
                dis[v]=d+w[i];
                qu.push((node){v,dis[v]});
            }
        }
    }
    printf("%lld\n",dis[tt]);
}

void init(){
    read(n); read(m);
    int u,v;
    ll val;
    ss=1; tt=1+m*2+1;
    REP(i,1,m){
        read(u); read(v); read(val);
        G[u].push_back(mk(val,++cnt));
        G[v].push_back(mk(val,++cnt));
        if(u==1)add(ss,cnt-1,val);
        if(v==1)add(ss,cnt,val);
        if(v==n)add(cnt-1,tt,val);
        if(u==n)add(cnt,tt,val);
        add(cnt-1,cnt,val);
        add(cnt,cnt-1,val);
    }
    REP(i,1,n){
        sort(G[i].begin(),G[i].end());
        REP(j,0,G[i].size()-2){
            u=G[i][j].se;
            v=G[i][j+1].se;
            val=G[i][j+1].fi-G[i][j].fi;
            add(u,v,val);
            add(v,u,0);
        }
    }
}

int main(){
    File();
    init();
    dijkstra();
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值