HDU 2485

http://acm.hdu.edu.cn/showproblem.php?pid=2485

n个车站，m条边，两边之间费用为1，问最少摧毁多少车站，使得1-n无法在k时间内到达

将2-(n-1)每个点拆成两个，并建立容量为1,费用为0的一条边，源点为1，汇点为2*n-2，这时求最小费用最大流，其中保证dis[t]<=k，求得的最大流即为摧毁的车站数量（因为1流量代表能从其中一个车站到达n点）

#include <iostream>

#include <cstdio>

#include <cstring>

#include <queue>

using namespace std ;

const int INF=0xfffffff ;

struct node{

    int s,t,cap,cost,nxt ;

}e[200005] ;

int sumflow ;

int n,m,k,cnt,head[1005],vis[1005],dis[1005],pre[1005] ;

void add(int s,int t,int cap,int cost)

{

    e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].cost=cost ;e[cnt].nxt=head[s] ;head[s]=cnt++ ;

    e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap=0 ;e[cnt].cost=-cost ;e[cnt].nxt=head[t] ;head[t]=cnt++ ;

}

int spfa(int s,int t,int N)

{

    for(int i=0 ;i<=N ;i++)

        dis[i]=INF ;

    dis[s]=0 ;

    memset(vis,0,sizeof(vis)) ;

    memset(pre,-1,sizeof(pre)) ;

    vis[s]=1 ;

    queue <int> q ;

    q.push(s) ;

    while(!q.empty())

    {

        int u=q.front() ;

        q.pop() ;

        vis[u]=0 ;

        for(int i=head[u] ;i!=-1 ;i=e[i].nxt)

        {

            int tt=e[i].t ;

            if(e[i].cap && dis[tt]>dis[u]+e[i].cost)

            {

                dis[tt]=dis[u]+e[i].cost ;

                pre[tt]=i ;

                if(!vis[tt])

                {

                    vis[tt]=1 ;

                    q.push(tt) ;

                }

            }

        }

    }

    if(dis[t]==INF || dis[t]>k)return 0 ;

    return 1 ;

}

int MCMF(int s,int t,int N)

{

    int flow,minflow,mincost ;

    mincost=flow=0 ;

    while(spfa(s,t,N))

    {

        minflow=INF ;

        for(int i=pre[t] ;i!=-1 ;i=pre[e[i].s])

            minflow=min(minflow,e[i].cap) ;

        flow+=minflow ;

        for(int i=pre[t] ;i!=-1 ;i=pre[e[i].s])

        {

            e[i].cap-=minflow ;

            e[i^1].cap+=minflow ;

        }

        mincost+=dis[t]*minflow ;

    }

    sumflow=flow ;//最大流

    return mincost ;

}

int main()

{

    while(scanf("%d%d%d",&n,&m,&k))

    {

        if(!n && !m && !k)break ;

        cnt=0 ;

        memset(head,-1,sizeof(head)) ;

        int s=1 ;

        int t=2*n-2 ;

        for(int i=1 ;i<n-1 ;i++)

            add(2*i,2*i+1,1,0) ;

        while(m--)

        {

            int a,b ;

            scanf("%d%d",&a,&b) ;

            a-- ;b-- ;

            add(a*2+1,b*2,1,1) ;

        }

        MCMF(s,t,t) ;

        printf("%d\n",sumflow) ;

    }

    return 0 ;

}

HDU 2485,码迷,mamicode.com