nhom(n1,18). nhom(n2,18).
phong(p30,30). phong(p40,40).
si_so(T) :- T = #sum { N,M : nhom(M,N) }.
1 { chon(P) : phong(P,_) } 1.
:- chon(P), phong(P,C), si_so(T), T>C.
#show si_so/1.
#show chon/1.