I,O=input,print N=int(I()) P=[0]+list(map(int,I().split())) U=[0]+[0]*N A=[] for i in range(1,N+1): j=P[i] if U[i]or j==i:continue while j^i: U[j],k=1,i while k=i: A+=[k] k-=1 j=P[j] O("Yes") O(len(A)) for a in A:O(a)