正题
题目链接:https://www.ybtoj.com.cn/contest/114/problem/1
题目大意
给出nnn个点mmm条边的一张无向图,对于每个点iii求不经过i∼1i\sim 1i∼1的最短路的第一条边的情况下iii到111的最短路
数据保证这条边唯一
n∈[1,105],m∈[1,2×105],c∈[1,103]n\in[1,10^5],m\in[1,2\times 10^5],c\in[1,10^3]n∈[1,105],m∈[1,2×105],c∈[1,103]
解题思路
因为保证的那个东西,所以图的最短路树真的是一棵树了,所以先跑出最短路树考虑在最短路树上面搞。
然后题目限制了我们不能从树上的祖先那条边过来,这样就分为了两种情况。一种是从该点的子树外面连过来的边,另一种是从子树中走上来的边。第二种很麻烦,因为子树的最短路是用该节点的最短路扩展的,所以不能直接使用。
考虑一条非树边(x,y)(x,y)(x,y),这条边会扩展一条disy+wdis_y+wdisy+w到xxx的路径。(disxdis_xdisx表示1∼x1\sim x1∼x的最短路)。
并且这条边可以使用到LCA(x,y)LCA(x,y)LCA(x,y)处,此时xxx的祖先们都不包含yyy在子树内,可以直接用yyy的子树扩展。
所以可以维护一个左偏树,每次合并两个儿子的信息,如果堆顶的边需要被删除就删除。需要写一个lazylazylazy标记来修改整棵树
时间复杂度O(nlogn)O(n\log n)O(nlogn)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#define mp(x,y) make_pair(x,y)
using namespace std;
const int N=4e5+10;
struct point{int val,x,y;point(int v=0,int xx=0,int yy=0){val=v;x=xx;y=yy;return;}
};
bool operator<(point x,point y)
{return x.val<y.val;}
struct Heap{point val[N];int t[N][2],lazy[N],dis[N];void Downdata(int x){if(!lazy[x])return;int ls=t[x][0],rs=t[x][1];lazy[ls]+=lazy[x];lazy[rs]+=lazy[x];val[ls].val+=lazy[x];val[rs].val+=lazy[x];lazy[x]=0;return;}int Merge(int x,int y){Downdata(x);Downdata(y);if(!x||!y)return x+y;if(val[y]<val[x])swap(x,y);int &ls=t[x][0],&rs=t[x][1];rs=Merge(rs,y);if(dis[rs]>dis[ls])swap(ls,rs);dis[x]=dis[rs]+1;return x;}int Del(int x){int &ls=t[x][0],&rs=t[x][1];val[x]=0;return Merge(ls,rs);}
}T;
struct node{int to,next,w;
}a[N];
int n,m,tot,cnt,num,ls[N],f[N];
int rfn[N],p[N],ans[N];
bool v[N];vector<int> G[N];
priority_queue<pair<int,int> >q;
void addl(int x,int y,int w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return;
}
void dij(){memset(f,0x3f,sizeof(f));q.push(mp(0,1));f[1]=0;while(!q.empty()){int x=q.top().second;q.pop();if(v[x])continue;v[x]=1;for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(f[x]+a[i].w<f[y]){f[y]=f[x]+a[i].w;q.push(mp(-f[y],y));}}}return;
}
void dfs(int x){rfn[x]=++cnt;for(int i=0;i<G[x].size();i++){int y=G[x][i];dfs(y);T.val[p[y]].val+=f[y]-f[x];T.lazy[p[y]]+=f[y]-f[x];p[x]=T.Merge(p[x],p[y]);}for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(f[x]+a[i].w==f[y])continue;if(f[y]+a[i].w==f[x])continue;T.val[++num]=point(f[y]+a[i].w,x,y);p[x]=T.Merge(p[x],num);}while(1){if(!p[x]){ans[x]=-1;break;}point w=T.val[p[x]];if(rfn[w.y]>=rfn[x]){p[x]=T.Del(p[x]);continue;}ans[x]=w.val;break;}return;
}
int main()
{freopen("pal.in","r",stdin);freopen("pal.out","w",stdout);scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int x,y,w;scanf("%d%d%d",&x,&y,&w);addl(x,y,w);addl(y,x,w);}dij();for(int x=1;x<=n;x++)for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(f[x]+a[i].w==f[y])G[x].push_back(y);}dfs(1);for(int i=2;i<=n;i++)if(!ans[i])puts("-1");else printf("%d\n",ans[i]);return 0;
}