OFFSET
1,3
COMMENTS
This is to 3 as A224486 is to 2
Displayed terms complete up to 200*10^6. - Joerg Arndt, Apr 08 2013
LINKS
Joerg Arndt, Table of n, a(n) for n = 1..64 (all terms <= 10^9)
EXAMPLE
0 is a term because (3^0+1)/(3*0+1) = 2.
1 is a term because (3^1+1)/(3*1+1) = 1.
9 is a term because (3^9+1)/(3*9+1) = 703.
PROG
(PARI) for(n=0, 10^9, if((3^n+1)%(3*n+1)==0, print1(n, ", "))); /* Joerg Arndt, Apr 08 2013 */
/* the following program is significantly faster; it gives terms >=1: */
(PARI) for(n=0, 10^12, my(m=3*n+1); if( Mod(3, m)^n==Mod(-1, m), print1(n, ", ") ) ); /* Joerg Arndt, Apr 08 2013 */
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Jonathan Vos Post, Apr 07 2013
EXTENSIONS
Terms > 9 from Joerg Arndt, Apr 08 2013
STATUS
approved