OFFSET
1,2
COMMENTS
No more terms through 10^6. - Ryan Propper, May 05 2006
The first 8 terms are all in A032863, all known subsequent terms, i.e., at least up to a(21), are not in A032863. - M. F. Hasler, Oct 05 2018
MATHEMATICA
t = 1; Do[t = 3t/2; If[ Mod[ Floor[t], n] == 0, Print[n]], {n, 500000}] (* Robert G. Wilson v, Apr 06 2006 *)
PROG
(PARI) a=1; for(n=1, 10^6, a*=3; b=shift(a, -n); if(b%n==0, print1(n, ", "))) \\ Robert Gerbicz, Aug 23 2006
(PARI) P=1; for(n=1, oo, (P*=3)>>n%n||print1(n", ")) \\ M. F. Hasler, Oct 05 2018
(Python)
from gmpy2 import mpz, t_div_2exp, t_mod
A073633_list, m = [], mpz(1)
for n in range(1, 10**9):
m *= 3
if t_mod(t_div_2exp(m, n), n) == 0:
A073633_list.append(n) # Chai Wah Wu, Mar 30 2020
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Benoit Cloitre, Aug 29 2002
EXTENSIONS
More terms from Michel ten Voorde Jun 20 2003
2 more terms from Ryan Propper, May 05 2006
More terms from Robert Gerbicz, Aug 23 2006
a(22)-a(24) from Chai Wah Wu, Mar 30 2020
STATUS
approved