HDU-1151 Air Raid

给出一个有向无环图，问最小路径覆盖是多少

最小路径覆盖=点数-最大匹配数

这个最大匹配就是，将每个点拆为入点，出点，然后形成二分图求最大匹配

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
const int N=240+10;
const int M=1e4+10;
const int INF=0x7f7f7f7f;
struct Edge
{
    int to,nxt;
}edge[M];
int tot,first[N];
void addedge(int u,int v)
{
    edge[tot].to=v;
    edge[tot].nxt=first[u];
    first[u]=tot++;
}
void init()
{
    tot=0;
    memset(first,-1,sizeof(first));
}

int uN,dis;
int Mx[N],My[N];
int dx[N],dy[N];
bool used[N];
bool SearchP()
{
    queue<int> q;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=uN;i++)
        if(Mx[i]==-1)
        {
            q.push(i);
            dx[i]=0;
        }
    dis=INF;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(dx[u]>dis) break;
        for(int i=first[u];i!=-1;i=edge[i].nxt)
        {
            int v=edge[i].to;
            if(dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(My[v]==-1) dis=dy[v];
                else
                {
                    dx[My[v]]=dy[v]+1;
                    q.push(My[v]);
                }
            }
        }
    }
    return dis!=INF;
}
bool dfs(int u)
{
    for(int i=first[u];i!=-1;i=edge[i].nxt)
    {
        int v=edge[i].to;
        if(!used[v]&&dy[v]==dx[u]+1)
        {
            used[v]=true;
            if(My[v]!=-1&&dy[v]==dis) continue;
            if(My[v]==-1||dfs(My[v]))
            {
                My[v]=u;
                Mx[u]=v;
                return true;
            }
        }
    }
    return false;
}
int MaxMatch()
{
    int res=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(My));
    while(SearchP())
    {
        memset(used,false,sizeof(used));
        for(int i=1;i<=uN;i++)
            if(Mx[i]==-1&&dfs(i)) res++;
    }
    return res;
}
int main()
{
    int T;
    int n,m,a,b;
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d%d",&n,&m);
        uN=n;
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a,b+n);
        }
        printf("%d\n",n-MaxMatch());
    }
    return 0;
}