JavaScript (ES6), 123 118 115 111 104 96 bytes
Guardado 4 bytes gracias a @Arnauld
G=n=>n?G(n>>1,++a[n%2]):a.some(n=>(P=x=>n%--x?P(x):x)(n)-1)
F=n=>F(n-1,G(n,a=[0,0,n])||alert(n))
Una combinación de tres funciones recursivas típicas. Alerta la secuencia en orden inverso y termina en un error de "demasiada recursividad".
Fragmento de prueba
(modificado para enviar a la página)
alert = s => O.textContent += s + "\n"
G=n=>n?G(n>>1,++a[n%2]):a.some(n=>(P=x=>n%--x?P(x):x)(n)-1)
F=n=>F(n-1,G(n,a=[0,0,n])||alert(n))
F(1000)
<pre id=O></pre>
La función principal puede devolver una matriz de 104 bytes:
G=n=>n?G(n>>1,++a[n%2]):a.some(n=>(P=x=>n%--x?P(x):x)(n)-1)
F=n=>n?F(n-1).concat(G(n,a=[0,0,n])?[]:n):[]
También puede ser no recursivo a costa de otro byte:
G=n=>n?G(n>>1,++a[n%2]):a.some(n=>(P=x=>n%--x?P(x):x)(n)-1)
n=>[for(_ of Array(n))if(!G(--n,a=[0,0,n]))n]
Aquí está el que comencé con: (Guardado 6 bytes gracias a @Arnauld)
P=(n,x=n)=>n>1&--x<2||n%x&&P(n,x)
G=n=>n?G(n>>1,o+=n%2,t++):P(o)&P(t-o)
F=n=>n?F(n-1).concat(P(n)&G(n,o=t=0)?n:[]):[]
Intenté jugar más al golf y logré hacerlo en 104 bytes, luego me di cuenta de que ya había encontrado esa solución (está en la parte inferior de la respuesta). ¿No odias cuando eso sucede? :PAGS
Un intento no recursivo en la función principal (de nuevo, el mismo recuento de bytes):
n=>[for(i of Array(n))if(P(--n)&G(n,o=t=0))n]
Éste toma la ruta fácil de contar cuántos 0 y 1 hay en la representación binaria:
F=n=>n?F(n-1).concat([n,(G=x=>n.toString(2).split(x).length-1)(0),G(1)].some(n=>(P=x=>n%--x?P(x):x)(n)-1)?[]:n):[]
Lo mismo con una comprensión de matriz:
n=>[for(_ of Array(n))if(![--n,(G=x=>n.toString(2).split(x).length-1)(0),G(1)].some(n=>(P=x=>n%--x?P(x):x)(n)-1))n]
Este toma una ruta un poco más difícil para hacer lo mismo:
F=n=>n?F(n-1).concat([n,(G=(x,w=n)=>w&&G(x,w>>1)+(w%2==x))(0),G(1)].some(n=>(P=x=>n%--x?P(x):x)(n)-1)?[]:n):[]
Y esta toma otra ruta relacionada que es tan corta como la original:
F=n=>n?F(n-1).concat([n,o=(G=x=>x&&x%2+G(n>>++t))(n,t=0),t-o].some(n=>(P=x=>n%--x?P(x):x)(n)-1)?[]:n):[]
Una vez más, puedes jugar 8 bytes haciendo que alerta la secuencia en orden inverso:
F=n=>F(n-1,[n,o=(G=x=>x&&x%2+G(n>>++t))(n,t=0),t-o].some(n=>(P=x=>n%--x?P(x):x)(n)-1)||alert(n))