import sys
input = sys.stdin.readline

def main():
    N = int(input().strip())
    P_list = list(map(int, input().split()))
    P = [0] + P_list
    C = [0] * (N + 1)
    for i in range(1, N + 1):
        C[i] = P[i] + 1 if P[i] < N else 1
    visited = [False] * (N + 1)
    cbl = {}
    for i in range(1, N + 1):
        if not visited[i]:
            cur = i
            cyc = []
            while not visited[cur]:
                visited[cur] = True
                cyc.append(cur)
                cur = C[cur]
            L = len(cyc)
            cbl.setdefault(L, []).append(cyc)

    B = [0] * (N + 1)
    for L, lists in list(cbl.items()):
        if L % 2 == 1:
            s = (L + 1) // 2
            for cyc in lists:
                for j in range(L):
                    B[cyc[j]] = cyc[(j + s) % L]

    for L, lists in list(cbl.items()):
        if L % 2 == 0:
            if len(lists) % 2 == 1:
                print("No")
                return
            for k in range(0, len(lists), 2):
                x = lists[k]
                y = lists[k + 1]
                for j in range(L):
                    B[x[j]] = y[j]
                    B[y[j]] = x[(j + 1) % L]

    Q = [0] * (N + 1)
    for i in range(1, N + 1):
        bi = B[i]
        if bi == 1:
            Q[i] = N
        else:
            Q[i] = bi - 1

    print("Yes")
    print(" ".join(map(str, Q[1:])))
if __name__ == "__main__":
    main()