comb := proc (n::integer, k::integer)
    return factorial(n)/(factorial(k)*factorial(n-k));
end proc;
perm := proc (n::integer, k::integer)
    return factorial(n)/factorial(n-k);
end proc;