greedy_min_facility_location.m

function list = greedy_min_facility_location(similarity_mat, K)
N = size(similarity_mat, 1);
remain_list = 1:N;
list = [];
for idx = 1:K
    max_value = -inf;
    
    for jdx = 1:length(remain_list)
        new_list = [list, remain_list(jdx)];
        if length(new_list) == 1
            value = min(similarity_mat(new_list,:));
        else
            value = min(max(similarity_mat(new_list,:)));
        end
        if value > max_value
            max_value = value;
            max_index = remain_list(jdx);
        end
    end
    list = [list, max_index];
    remain_list = setdiff(remain_list, max_index);
end
end