题目描述:
洛谷题目链接
有向无环图,经过一条边有一定的费用,每次从1号点出发走一条路径,可以在任一点结束并返回一号点,求把每条边至少经过一次的最小费用。
题目分析:
每条边至少经过一次,即下界为1的最小费用可行流,为满足流量守恒,让每个点向一号点连一条下界为0,上界inf,费用为0的边。然后虚拟源汇补流,跑最小费用流即可。
Code:
#include<cstdio>
#include<deque>
#include<cstring>
#include<algorithm>
#define maxn 305
#define maxm 11205
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,S,T,ss,tt;
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],c[maxm],w[maxm],tot=1;
inline void line(int x,int y,int z,int v){
nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z,w[tot]=v;
nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=0,w[tot]=-v;
}
namespace Costflow{
int dis[maxn],cst;
bool vis[maxn],inq[maxn];
deque<int>q;
bool SPFA(){
memset(dis,0x3f,(T+1)<<2);
dis[T]=0,q.push_back(T);
while(!q.empty()){
int u=q.front();q.pop_front(),inq[u]=0;
for(int i=fir[u],v;i;i=nxt[i]) if(c[i^1]&&dis[v=to[i]]>dis[u]+w[i^1]){
dis[v]=dis[u]+w[i^1];
if(!inq[v]) {inq[v]=1;if(!q.empty()&&dis[v]<dis[q.front()]) q.push_front(v);else q.push_back(v);}
}
}
return dis[S]!=inf;
}
int aug(int u,int augco){
if(u==T) {cst+=dis[S]*augco;return augco;}
vis[u]=1;
int need=augco,delta;
for(int &i=cur[u];i;i=nxt[i]) if(c[i]&&!vis[to[i]]&&dis[u]==dis[to[i]]+w[i]){
delta=aug(to[i],min(need,c[i]));
c[i]-=delta,c[i^1]+=delta;
if(!(need-=delta)) break;
}
vis[u]=0;
return augco-need;
}
int solve(){
cst=0;
while(SPFA()) memcpy(cur,fir,(T+1)<<2),aug(S,inf);
return cst;
}
}
int tf[maxn],ans;
int main()
{
scanf("%d",&n);
for(int i=1,k,j,x;i<=n;i++){
scanf("%d",&k);
while(k--) scanf("%d%d",&j,&x),ans+=x,tf[i]--,tf[j]++,line(i,j,inf,x);
}
for(int i=2;i<=n;i++) line(i,1,inf,0);
ss=0,tt=n+1;
for(int i=1;i<=n;i++)
if(tf[i]>0) line(ss,i,tf[i],0);
else if(tf[i]<0) line(i,tt,-tf[i],0);
S=ss,T=tt,line(n,1,inf,0),ans+=Costflow::solve();
printf("%d\n",ans);
}