>>436
xをx=k(a[1]+a[p]√p+a[q]√q+a[pq]√(pq))/d^2のような標準形で表す。
a[1],a[p],...は整数で、最小公倍数が1になるようにする(kとしてくくり出す)。
p,q,...は素数で、xに出現する根号の中を素因数分解して出現するもの全てを考える。
すると√x∈Sならばx=√k(b[1]+b[p]√p+b[q]√q+b[pq]√(pq))/dの形で表現できるはず。
これで未知変数b[1],b[p],...の数と同じ本数の連立2次方程式が立つから、
それを解けるような数値計算パッケージがもしあれば解決する。