function n = Calc(p)
    N = 1e5;
    for pessoas = 2 : 366
        a = randi(365, pessoas, N);
        for i = 1:N
            res(i) = length(unique(a(:, i))) < pessoas;
        end
        prob = sum(res)/N;
        if (prob > p)
            n = pessoas;
            break
        end
    end
end

fprintf("a) = %.4f\n", Calc(0.5))
fprintf("b) = %.4f\n", Calc(0.9))