OFFSET
1,1
COMMENTS
MAPLE
nn:=10^20:
for n from 1 to 10 do:
ii:=0:
for k from 1 to nn while(ii=0) do:
it0:=0:it1:=0:m:=k:
for i from 1 to nn while(m<>1) do:
if irem(m, 2)=0
then
m:=m/2:it0:=it0+1:
else
m:=3*m+1:it1:=it1+1:
fi:
od:
if it1<>0 and it0/it1 = ithprime(n)
then
ii:=1:printf(`%d %d \n`, n, k):
else
fi:
od:
od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Sep 10 2018
STATUS
approved