cache = [1]
A000961(n) = {
    if (#cache < n,
        for (v=cache[#cache]+1, oo,
           if (omega(v)==1,
              cache = concat(cache, v);
              if (#cache==n, break);
           );
        );
    );
    return (cache[n])
}

\\ c(t,n) = (n+1)-th tuple of t pairwise coprime numbers whose factorizations depend on the base t+1 representation of n 
\\ in particular:
\\     c(1,n) == [A289272(n)]
\\     c(2,n) == [A289815(n), A289816(n)]
c(t,n) = {
    my (v=vector(t, i, 1));
    my (x=1);
    for (i=2, oo,
        my (o=A000961(i));
        if (n==0,
            return (v)
        );
        if (gcd(x,o)==1 && omega(o)==1,
            my (d = n % (t+1));
            if (d>0,
               v[d] *= o;
               x    *= o;
            );
            n \= t+1;
        );
    );
}