from collections import *
import sys
import heapq
import bisect
import itertools

dxdy1 = ((0, 1), (0, -1), (1, 0), (-1, 0))
dxdy2 = ((0, 1), (0, -1), (1, 0), (-1, 0), (1, 1), (-1, -1), (1, -1), (-1, 1))
dxdy3 = ((0, 1), (1, 0))
dxdy4 = ((1, 1), (1, -1), (-1, 1), (-1, -1))
INF = float("inf")
MOD = 998244353
mod = 998244353
MOD2 = 10**9 + 7
mod2 = 10**9 + 7
# memo : len([a,b,...,z])==26

input = lambda: sys.stdin.readline().rstrip()
mi = lambda: map(int, input().split())
li = lambda: list(mi())

N = int(input())
cnt = 1
while N != 1:
    if N % 2 == 0:
        N //= 2
    else:
        N = 3 * N + 1
    cnt += 2
    if cnt > 50:
        break
cnt += 1
if cnt <= 50:
    print("Yes")
    print(cnt)
else:
    print("No")