最优贸易
题目
解析
这道题有很多种解法:暴力DP,tarjan缩点+拓扑DP,两遍spfa等等
这里我们讲分层图+spfa
分层图:
在残量网络中,满足d[y]=d[x]+1的边(x,y)构成的子图被称为分层图
将每个状态作为一层图,并在可转换的层间建图,即为分层图
像这样:(图片 by fy1234567ok)
于是直接跑一遍spfa即可
这种鬼畜图是卡不了spfa的,所以复杂度为
O
(
k
(
n
+
m
)
)
O(k(n+m))
O(k(n+m))
code:
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
int tot=0,k,n,m,head[300010],to[2000010],w[2000010],nxt[2000010],c[100010],ans[300010],x,y,z;
bool u[300010];
void add(int x,int y,int z)
{
to[++tot]=y,w[tot]=z,nxt[tot]=head[x],head[x]=tot;
}
queue <int> d;
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
ans[i]=ans[i+n]=ans[i+2*n]=-1e8;
scanf("%d",&c[i]);
add(i,i+n,-c[i]);
add(i+n,i+2*n,c[i]);
}
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
if(z==2)add(y,x,0),add(y+n,x+n,0),add(y+2*n,x+2*n,0);
add(x,y,0),add(x+n,y+n,0),add(x+2*n,y+2*n,0);
}
u[1]=1,ans[1]=0,d.push(1);
while(!d.empty())
{
k=d.front(),d.pop();
u[k]=0;
for(int i=head[k];i;i=nxt[i])
{
if(ans[to[i]]<ans[k]+w[i])
{
ans[to[i]]=ans[k]+w[i];
if(!u[to[i]])
{
u[to[i]]=1;
d.push(to[i]);
}
}
}
}
printf("%d",ans[3*n]);
return 0;
}