noip模拟(有向无环图)题解
充斥着结构体香气~
(
因此这篇题解会其它的稍微简单些
首先呢该题的核心就是建图…这个我们就结合着图来说吧!
图中的黑色所画出的是原建图也就是题中原本所给的图,然后花里胡哨,五颜六色的是我们所建的新图。首先呢,我们需要遍历图中的所有节点(而我所给的图只是万千节点中的一个栗子),我们每到一个新节点时就需要记录一下当前的边节点和原边权(所谓边节点就是你在使用邻接链表存图时的那个变量,我们可以用存边,在本题中当作边节点),在本题中我们需要把一个边拆成两个边节点来存。之后呢?我们需要排序,为什么要排序呢?这就是这个题的精髓;我们需要在不同边上的边节点建路(这种路我们定义叫做异边路),在排序后我们可以建造使权值相邻的异边路,这样我们就可以通过权值差来构造一个类似环的通路,这样我们可以使权值相差较大的异边路通过这个类环形通路来到达又可以满足题目要求(这么干还主要是可以控制异边路的个数从而保证效率)(重申相邻权值异边路作差为该异边路的权值,如图)。在建造好正向异边路时,我们需要再建造一个反向的,而此时权值为0;(原因的话很显然不多说了)
(在处理时我们先不管1和n节点)
1节点我们在Dijkstra中预处理,而n节点我们则在输出时处理一下即可;
下面上代码…
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
#define ll long long
#define rei register int
using namespace std;
const int maxn=2e5+10;
int n,m;
struct dnode
{
int num;ll vval;
};
bool operator <(const dnode &ddx,const dnode &ddy)
{
return ddx.vval>ddy.vval;
}
priority_queue<dnode> stac;
struct hiphop
{
int st;
int nd;
}res[maxn*3];
ll d[maxn*2];
struct node
{
int ver;
int nxt;
ll val;
}edge[maxn*2];
int head[maxn],ecnt=1;
bool cmp( hiphop x, hiphop y)
{
return x.nd<y.nd;
}
void add(int x,int y,ll z)
{
edge[++ecnt].ver=y;
edge[ecnt].nxt=head[x];
edge[ecnt].val=z;
head[x]=ecnt;
}
int vis[maxn*2];
struct heapnode
{
int to,next;
ll vale;
}e[maxn*5];
int h[maxn*3];
int cnt;
void Add(int x,int y,ll z)
{
e[++cnt].next=h[x];
h[x]=cnt;
e[cnt].to=y;
e[cnt].vale=z;
}
void Dijkstra()
{
for(rei i=2;i<=ecnt;++i)
d[i]=2e18;
for(rei i=head[1];i;i=edge[i].nxt)
{
d[i]=edge[i].val;
dnode u;u.num=i,u.vval=d[i];
stac.push(u);
}
while(stac.size())
{
if(vis[stac.top().num])
{stac.pop();continue;}
int x=stac.top().num;
vis[x]=1;stac.pop();
for(rei i=h[x];i;i=e[i].next)
{
if(vis[e[i].to]) continue;
if(d[e[i].to]>d[x]+e[i].vale)
{
d[e[i].to]=d[x]+e[i].vale;
dnode v;v.num=e[i].to;v.vval=d[e[i].to];
stac.push(v);
}
}
}
}
int tot;
int main()
{
int a,b;ll c;
cin>>n>>m;
for(rei i=1;i<=m;++i)
{
scanf("%d%d%lld",&a,&b,&c);
add(a,b,c);add(b,a,c);
}
for(rei x=2;x<n;++x)
{
tot=0;
for(rei i=head[x];i;i=edge[i].nxt)
{
res[++tot].st=i,res[tot].nd=edge[i].val;
}
sort(res+1,res+tot+1,cmp);
for(rei i=1;i<=tot;++i)
{
Add(res[i].st^1,res[i].st,res[i].nd);
if(i!=1) Add(res[i].st,res[i-1].st,0);
if(i!=tot) Add(res[i].st,res[i+1].st,res[i+1].nd-res[i].nd);
}
}
Dijkstra();
ll ans=2e18;
for(rei i=head[n];i;i=edge[i].nxt)
{
ans=min(ans,d[i^1]+edge[i].val);
}
cout<<ans<<endl;
return 0;
}