zerocover.m

真正可以用的匈牙利算法 我校赛的时候用过 很不错 不过写的不够简洁

function [y,Row,Col]=ZeroCover(C,Row,Col,Cur)
[m,n]=size(C);
temp=C(Cur,:);      
flag=0;
Cur=[];
for i=1:m
    if temp(i)==-inf & Col(i)~=0      
        Col(i)=0;
        Cur=[Cur,i];
        flag=1;
    end
end
if flag==0                            
    y=0;
else                                 
    for i=Cur
        temp=C(:,i);
        for i=1:m
            if temp(i)==inf & Row(i)==0   
                Row(i)=1;
                y=i;
                flag=0;
                [y,Row,Col]=ZeroCover(C,Row,Col,y);
                break;
            end
        end
        if flag==1       
            y=0;
        end
    end
end