while(1)!rem((bin2dec(fgetl(0))*8+1)^.5,1)endwhile