最小费用最大流: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;
}
只要存在完美匹配就可以保证所有的点都在环上,本题其实就是求二分图的最小权匹配。还有一个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;
}
相关文章推荐
- 动易2006序列号破解算法公布
- Ruby实现的矩阵连乘算法
- C#插入法排序算法实例分析
- 超大数据量存储常用数据库分表分库算法总结
- C#数据结构与算法揭秘二
- C#冒泡法排序算法实例分析
- 算法练习之从String.indexOf的模拟实现开始
- C#算法之关于大牛生小牛的问题
- C#实现的算24点游戏算法实例分析
- c语言实现的带通配符匹配算法
- 浅析STL中的常用算法
- 算法之排列算法与组合算法详解
- C++实现一维向量旋转算法
- Ruby实现的合并排序算法
- C#折半插入排序算法实现方法
- 基于C++实现的各种内部排序算法汇总
- C++线性时间的排序算法分析
- C++实现汉诺塔算法经典实例
- PHP实现克鲁斯卡尔算法实例解析
- C#获取关键字附近文字算法实例