hdu 2485 Destroying the bus stations（费用流）

题意：

n个点，m条有向边，问最少去掉几个点使得从1到n的路径（最短路）大于k？

思路：

这道题用费用流解，感觉好奇妙。。。

首先我们想，求费用流时，我们是每次用spfa求得最短路然后进行减流，直到找不到从汇点到源点的路径。那么对于这道题，当我们增广出一条从1到n的最短路时，我们可以将其中的某个点删去，使得这条路以及从这条路分支出来的路都使得不能从1走到n。所以我们可将点2~n进行拆点，所加的边的容量为1，代表这条路只能走一次，费用为0，因为这条边不需要花费，而其他边的容量为1，费用为1。这样，当我们跑费用流时，当dis【t】>k时，所流过的流量就是答案。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
using namespace std;
#define inf 0x3f3f3f3f
#define maxn 200
#define maxm 4050
int head[maxn];
int dis[maxn];
int pre[maxn];
int vis[maxn];
int tot;
int n,m,k;
struct Edge
{
    int to;
    int cost;
    int flow;
    int next;
} edge[maxm*2];
void addedge(int u,int v,int w,int f)
{
    edge[tot].to=v;
    edge[tot].cost=w;
    edge[tot].flow=f;
    edge[tot].next=head[u];
    head[u]=tot++;

    edge[tot].to=u;
    edge[tot].cost=0-w;
    edge[tot].flow=0;
    edge[tot].next=head[v];
    head[v]=tot++;
    return;
}
bool spfa(int s,int t)
{
    memset(dis,inf,sizeof(dis));
    memset(pre,-1,sizeof(pre));
    memset(vis,0,sizeof(vis));
    queue<int>Q;
    Q.push(s);
    dis[s]=0;
    vis[s]=1;
    int top;
    while(!Q.empty())
    {
        top=Q.front();
        Q.pop();
        vis[top]=0;
        for(int i=head[top]; i!=-1; i=edge[i].next)
        {
            if(edge[i].flow>0&&dis[edge[i].to]>dis[top]+edge[i].cost)
            {
                dis[edge[i].to]=dis[top]+edge[i].cost;
                pre[edge[i].to]=i;
                if(!vis[edge[i].to])
                {
                    vis[edge[i].to]=1;
                    Q.push(edge[i].to);
                }
            }
        }
    }
    if(dis[t]>0)return true;
    return false;
}
int solve(int s,int t)
{
    int flow=0;

    while(spfa(s,t))
    {
        //cout<<"++"<<endl;
        if(dis[t]>k)break;
        flow++;
        for(int i=pre[t]; i!=-1; i=pre[edge[i^1].to])
        {
            edge[i].flow-=1;
            edge[i^1].flow+=1;
        }


    }
    return flow;
}
int main()
{

    while(~scanf("%d%d%d",&n,&m,&k))
    {
        if(n==0&&m==0&&k==0)break;
        tot=0;
        memset(head,-1,sizeof(head));
        for(int i=2; i<n; i++)
        {
            addedge(i,i+n,0,1);
        }
        int a,b;
        for(int i=1; i<=m; i++)
        {
            scanf("%d%d",&a,&b);
            if(a!=1)a+=n;
            addedge(a,b,1,1);
        }
        printf("%d\n",solve(1,n));
    }
}