您的位置:首页 > 理论基础 > 计算机网络

最小费用最大流:HDU3488

2015-08-03 20:30 459 查看
需要将每个点都拆成两个点,然后根据输入的边建立二分图,再建立超源点和超汇点,再y用最小费用最大流就搞定。

只要存在完美匹配就可以保证所有的点都在环上,本题其实就是求二分图的最小权匹配。还有一个KM算法可以更好地解决这类问题,抽时间去学学。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define maxN 205
#define maxM 30005
#define inf 0xffffff
using namespace std;

struct node
{
int c;
int f;
int c_f;
int v;
}G[maxN*2+5][maxN*2+5];
int pre[maxN*2+5];
int dist[maxN*2+5];
int inq[maxN*2+5];
int n,m;
int vertex;
int min_f;
int sum;

void init()
{
int i,j;
int u,v,w;
scanf("%d%d",&n,&m);
for(i=0;i<=2*n;i++)
{
for(j=0;j<=2*n;j++)
{
G[i][j].c=0;
G[i][j].f=0;
G[i][j].c_f=0;
G[i][j].v=inf;
}
}
for(i=0;i<m;i++)
{
scanf("%d%d%d",&u,&v,&w);
G[u][n+v].c=1;
G[u][n+v].c_f=1;
if(w<G[u][n+v].v)
G[u][n+v].v=w;
}
vertex=n*2+1;
for(i=1;i<=n;i++)
{
G[0][i].c=1;
G[0][i].c_f=1;
G[0][i].v=0;
}
for(i=1;i<=n;i++)
{
G[n+i][vertex].c=1;
G[n+i][vertex].c_f=1;
G[n+i][vertex].v=0;
}
}

void spfa(int s)
{
queue<int> q;
for(int i=0;i<=vertex;i++)
{
pre[i]=-1;
dist[i]=inf;
inq[i]=0;
}
inq[s]=1;
dist[s]=0;
q.push(s);
while(!q.empty())
{
int u=q.front();
q.pop();
inq[u]=0;
for(int i=0;i<=vertex;i++)
{
if(G[u][i].c_f==0) continue;
if(G[u][i].v==inf) G[u][i].v=-G[i][u].v;
if(dist[i]>dist[u]+G[u][i].v)
{
dist[i]=dist[u]+G[u][i].v;
pre[i]=u;
if(!inq[i])
{
inq[i]=1;
q.push(i);
}
}
}
}
}

void argument(int s,int t)
{
spfa(s);

while(pre[t]!=-1)
{
sum+=dist[t];
int u=pre[t];
int v=t;
min_f=inf;
while(u!=-1){
if(min_f>G[u][v].c_f)
{
min_f=G[u][v].c_f;
}
v=u;
u=pre[v];
}
u=pre[t];v=t;
while(u!=-1)
{
G[u][v].f+=min_f;
G[v][u].f-=min_f;
G[u][v].c_f=G[u][v].c-G[u][v].f;
G[v][u].c_f=G[v][u].c-G[v][u].f;
v=u;
u=pre[v];
}
spfa(s);
}
}

int main()
{
int T;
scanf("%d",&T);
while(T--)
{
init();
sum=0;
argument(0,vertex);
printf("%d\n",sum);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签:  网络流 算法