HDU 3435 费用流

141 篇文章 0 订阅
34 篇文章 0 订阅

点击打开链接

题意:给个无向的图,问你删除任意边后,使这个图是哈密顿图,若有多个,输出路径上的所有权值和最小,没有就输出NO

思路:今天开始从10年多校开始刷题,敌人留给我们的时间不多了,这题看完题意后,看了看样例,自己yy了一下,写了一发,交了ac,看了样例后我是这样想的,因为是哈密顿图,那么每个点肯定是走了两次,而图是无向图,和求最大匹配有些类似,然后就瞎YY的过了

#include <queue>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=2010;
typedef pair<int,int> P;
struct edge{
    int to,cap,rev,cost;
    edge();
    edge(int a,int b,int c,int d){to=a,cap=b,cost=c,rev=d;};
};
vector<edge>G[maxn];
int h[maxn],dis[maxn];
int prevv[maxn],preve[maxn];
void add_edge(int st,int en,int cap,int cost){
    G[st].push_back(edge(en,cap,cost,G[en].size()));
    G[en].push_back(edge(st,0,-cost,G[st].size()-1));
}
int min_cost_flow(int st,int en,int f){
    int ans=0;
    memset(h,0,sizeof(h));
    while(f>0){
        priority_queue<P,vector<P>,greater<P> >que;
        for(int i=0;i<maxn;i++) dis[i]=inf;
        dis[st]=0;que.push(P(0,st));
        while(!que.empty()){
            P p=que.top();que.pop();
            int v=p.second;
            if(dis[v]<p.first) continue;
            for(unsigned int i=0;i<G[v].size();i++){
                edge &e=G[v][i];
                if(e.cap>0&&dis[e.to]>dis[v]+e.cost+h[v]-h[e.to]){
                    dis[e.to]=dis[v]+e.cost+h[v]-h[e.to];
                    prevv[e.to]=v;
                    preve[e.to]=i;
                    que.push(P(dis[e.to],e.to));
                }
            }
        }
        if(dis[en]==inf) return -1;
        for(int i=0;i<maxn;i++) h[i]+=dis[i];
        int d=f;
        for(int i=en;i!=st;i=prevv[i]){
            d=min(d,G[prevv[i]][preve[i]].cap);
        }
        f-=d;
        ans+=d*h[en];
        for(int i=en;i!=st;i=prevv[i]){
            edge &e=G[prevv[i]][preve[i]];
            e.cap-=d;
            G[i][e.rev].cap+=d;
        }
    }
    return ans;
}
int main(){
    int T,n,m,u,v,cost,cas=1;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=0;i<maxn;i++) G[i].clear();
        for(int i=1;i<=n;i++) add_edge(0,i,1,0);
        for(int i=1+n;i<=2*n;i++) add_edge(i,2*n+1,1,0);
        for(int i=0;i<m;i++){
            scanf("%d%d%d",&u,&v,&cost);
            add_edge(u,v+n,1,cost);
            add_edge(v,u+n,1,cost);
        }
        int ans=min_cost_flow(0,2*n+1,n);
        if(ans==-1) printf("Case %d: NO\n",cas++);
        else printf("Case %d: %d\n",cas++,ans);
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值