i=int(input())
x=1

def df(n):
    return 2*n

def f(n):
    return 0.5*(n**2)+0.5*n-i
def s(n):
    return 0.5*(n**2)+0.5*n
op=1
while True:
    prev=x
    x=x-(f(x)/df(x)) # ニュートン法
    if abs(prev - x) < 0.01:
        break
x=int(x-3)
for j in range(6):
    x=x+1
    if s(x)==float(i):
        print("YES")
        print(x)
        exit()
print("NO")