The descent set $D(w)$ of a permutation $w=a_1 a_2\cdots a_n\in\frak{S}_n$ is defined by $D(w)=\{ 1\leq i\leq n-1\,:\, a_i>a_{i+1}\}$.
Given a set $S$, let $\beta_n(S)$ denote the number of permutations $w\in\frak{S}_n$ with $D(w)=S$, and let $f(n)$ be the number of subsets $S$ of $\{1,2,\dots,n-1\}$ for which $\beta_n(S)$ is odd.
QUESTION. Is there a nice formula, recurrence, etc., for $f(n)$? The values of $f(n)$ for $1\leq n\leq 15$ are (up to some glitch in my Maple code) $$ 1,2,2,8,8,16,32,128,128,256,512,1024,2048,4096,7424. $$
Note. A closely related problem is to look at permutations $w\in \frak{S}_n$ whose descent set is contained in $S$. Let $\alpha_n(S)$ denote the number of such permutations (which is equal to a simple multinomial coefficient). Let $g(n)$ be the number of subsets $S$ of $\{1,2,\dots,n-1\}$ for which $\alpha_n(S)$ is odd. It is not hard to show that $g(n)$ is equal to the number of ordered set partitions of $\{1,2,\dots,b(n)\}$, where $b(n)$ is the number of 1's in the binary expansion of $n$. See https://oeis.org/A000670 for further information on ordered set partitions.