KM算法模板

int g[maxn][maxn];int S[maxn],T[maxn];int lv[maxn],rv[maxn],leftp[maxn];int slack;int tot1,tot2;void update(){for(int i=1;i<=tot1;i++)    {        if(S[i]) lv[i]-=slack;        if(T[i]) rv[i]+=slack;     }}int match(int u){    S[u]=1;    int i;for(i=1;i<=tot1;i++)    {        if(lv[u]+rv[i]==g[u][i])        {            if(!T[i])            {                T[i]=1;                if(leftp[i]==-1||match(leftp[i]))                {                    leftp[i]=u;                    return 1;                }            }        }        else slack=min(lv[u]+rv[i]-g[u][i],slack);    }    return 0;}int  solve(){    int i,j;memset(leftp,-1,sizeof(leftp));for(i=1;i<=tot1;i++)    {        lv[i]=rv[i]=0;for(j=1;j<=tot1;j++)            lv[i]=max(lv[i],g[i][j]);    }for(i=1;i<=tot1;i++)    {        for(;;)        {for(j=1;j<=tot1;j++)                 S[j]=T[j]=0;            slack=inf;            if(match(i)) break;            else update();        }    }    int sum=0;    for(i=1;i<=tot1;i++) sum=sum+lv[i]+rv[i];    return sum;}