P:=proc(q, h) local a, b, d, n, t; a:=1/h; b:=ilog10(h)+1; d:=h;
print(d); t:=-2; for n from 1 to q do if trunc(evalf(a+t/n, 100)*10^(b+ilog10(n)+1))=d*10^(ilog10(n)+1)+n then b:=b+ilog10(n)+1; d:=d*10^(ilog10(n)+1)+n;
a:=a+t/n; t:=(-1)^abs(t)*(abs(t)+1); print(n); fi; od; end: P(10^20, 3);
