double f(ll N:Memoize){ll D[1d3],L=Divisor(N,D);return N>1?(N+sum[d,1,L](f(N/D[d])*EulerPhi(D[d])))/(N-1):0;}{ll@N;wt(f(N));}