hdu 2485(最小费用最大流)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2485

思路：题目的意思是删除最少的点使1，n的最短路大于k。将点转化为边，容量为1，费用为0，然后就是对于那些有道路的城市之间连边，若（u,v)有边，则连边(u+n)->v,容量为inf,费用为花费的时间1，然后就是跑最小费了，若dist[vt]>k，则返回false,最后输出的flow就代表要删除的点的个数。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXN 222
#define MAXM 4444444
#define inf 1<<30

struct Edge{
    int v,cap,cost,next;
}edge[MAXM];

int n,m,k,vs,vt,NE;
int head[MAXN];

void Insert(int u,int v,int cap,int cost)
{
    edge[NE].v=v;
    edge[NE].cap=cap;
    edge[NE].cost=cost;
    edge[NE].next=head[u];
    head[u]=NE++;

    edge[NE].v=u;
    edge[NE].cap=0;
    edge[NE].cost=-cost;
    edge[NE].next=head[v];
    head[v]=NE++;
}

int dist[MAXN],pre[MAXN],cur[MAXN];
bool mark[MAXN];
bool spfa(int vs,int vt)
{
    memset(mark,false,sizeof(mark));
    fill(dist,dist+2*n+1,inf);
    dist[vs]=0;
    queue<int>que;
    que.push(vs);
    while(!que.empty()){
        int u=que.front();
        que.pop();
        mark[u]=false;
        for(int i=head[u];i!=-1;i=edge[i].next){
            int v=edge[i].v,cost=edge[i].cost;
            if(edge[i].cap>0&&dist[u]+cost<dist[v]){
                dist[v]=dist[u]+cost;
                pre[v]=u;
                cur[v]=i;
                if(!mark[v]){
                    mark[v]=true;
                    que.push(v);
                }
            }
        }
    }
    if(dist[vt]>k)return false;
    return dist[vt]!=inf;
}

int MinCostFlow(int vs,int vt)
{
    int flow=0,cost=0;
    while(spfa(vs,vt)){
        int aug=inf;
        for(int u=vt;u!=vs;u=pre[u]){
            aug=min(aug,edge[cur[u]].cap);
        }
        flow+=aug,cost+=dist[vt]*aug;
        for(int u=vt;u!=vs;u=pre[u]){
            edge[cur[u]].cap-=aug;
            edge[cur[u]^1].cap+=aug;
        }
    }
    return flow;
}

int main()
{
    int u,v;
    while(~scanf("%d%d%d",&n,&m,&k)){
        if(n==0&&m==0&&k==0)break;
        NE=0;
        vs=1,vt=n;
        memset(head,-1,sizeof(head));
        for(int i=2;i<=n-1;i++)Insert(i,i+n,1,0);
        while(m--){
            scanf("%d%d",&u,&v);
            if(u==vs)Insert(vs,v,inf,1);
            else Insert(u+n,v,inf,1);
        }
        printf("%d
",MinCostFlow(vs,vt));
    }
    return 0;
}