https://ningenme.hatenablog.com/entry/2020/12/12/170313

問題条件を下記のように[1]-[4]とします。
[1]. どの要素も $6$ で割ったあまりは $1$ でない
[2]. どの相異なる $2$ つの要素も、 $1$ で割ったあまりは等しい
[3]. どの要素も約数の個数は $4$ 個である
[4]. どの相異なる $2$ つの要素も、最大公約数は $1$ である

まず[1]の条件で禁止される数は全て排除しましょう。

次に[2]の条件より、 $6$ で割ったあまり $0,2,3,4,5$ に関してあまりを決め打ちして全探索を行いましょう。

以下その前提ですすめます。

まず[3]の条件より、集合の要素は全て次のいずれかです。

・ある素数 $p$ の $3$ 乗
・相異なる素数 $p$ と $q$ に対して $pq$

ここで[4]の条件より、集合内のどの $2$ 要素も同じ素因数を持たないことがわかります。

よってある素数 $p$ の $3$ 乗に関しては貪欲に全て集合に加えて良いです。またここで含んだ素数 $p$ を持つ $pq$ の形で書ける数に関してはこの時点で除外して問題ないです。

ここで $5$ 以上の素数に関しては全て $p=6n+1$ あるいは $p=6m-1$ と書くことが出来ることを考慮して、場合分けを考えましょう。

・ $6$ で割った余りが $0$ のとき

加えられる数は、 $6$ のみです。
集合の最大サイズは高々 $1$ です。

・ $6$ で割った余りが $2$ のとき

加えられる数は、 $8$ あるいは $2p$ (ここで $p$ は $p=6n+1$ の形で書ける素数)のみです
ここで $8$ と $2q$ は互いに素ではないので、集合の最大サイズは高々 $1$ です。

・ $6$ で割った余りが $3$ のとき

加えられる数は、 $27$ あるいは $3p$ (ここで $p$ は $p=6n+1$ の形で書ける素数)あるいは $3q$ (ここで $q$ は $q=6m-1$ の形で書ける素数)のみです
ここで $27,3p,3q$ はそれぞれ互いに素ではないので、集合の最大サイズは高々 $1$ です。

・ $6$ で割った余りが $4$ のとき

加えられる数は、 $2q$ (ここで $q$ は $q=6m-1$ の形で書ける素数)のみです
集合の最大サイズは高々 $1$ です。

・ $6$ で割った余りが5のとき

加えられる数は、 $p^3$ あるいは $pq$ (ここで $p$ は $p=6n+1$ の形で書ける素数、 $q$ は $q=6m-1$ の形で書ける素数)のみです
$pq$ の形で書ける数 $p_1q_1,p_2q_2$ に関しては、 $p_1 \neq p_2$ かつ $q_1 \neq q_2$ ならば、両方選択することが出来ます。
すなわち $pq$ の形で書ける素数に関して、素数を頂点としたグラフを考え頂点 $p$ と頂点 $q$ に辺を張ると、この問題は辺の張られたグラフから最大マッチングを考える問題になります。
また上記のグラフは $2$ 部グラフになるため、最大マッチングが簡単に求まります。

よってこの問題が解けました。
計算量は素因数分解の際に $O(N \sqrt{A_{max}})$ 、二部マッチングの際に $O(N^2)$ 必要で、全体で $O(N \sqrt{A_{max}}+N^2)$ となります。
素因数分解の際にポラード・ロー法などを用いることで、計算量をもう少し良くした解法もあります。