Algorithm: ComputingDSS (n, q,ρ)

procedure RecPartition(n, B, N )
if n = 0
then if N=q+1 and 2*∑_i≠j a_i * a_j ≥(n-1)* ρ
                       then for Q₁ in U and | Q₁ | = a₁
do for Q₂ in U \ Q₁ and | Q₂ | = a₂
                            do for Q₃ in U \ (Q₁∪Q₂) and | Q₃ | = a₃
...
                                 do for Q_q in U \ (Q₁∪Q₂∪...∪Q_q-1 ) and | Q_q | = a_q
       doM ← { a – b ( mod n ) | a∈Q_i , b∈Q_j , i≠j}
                                                ifM is DSS
                                                then return ( [Q₁ , ..., Q_q ] )
else for i ← 1 tomin (B, n)
do a_N+1 ← i
RecPartition( n-i, i, N+1)


main
U ← {1, 2, ... , n}
    r ← Ceiling(Sqrt ( (q * n * ρ)/ (q -1) ) )
RecPartition(r, r, 0)