1,4

Curiously, sequence seems related to pentagonal (or 5-gonal) or heptagonal (or 7-gonal) numbers. Some primes follow rules in a(n) factorization. If b(n)= exponent of 13 in a(n) factorization: b(n)=0, 0, 0, 2, 4, 6, 10, 14, 18, 24, 30, 36, 44, 52, 60, 70, 80, 90...so b(3n+1)=A049450(n); b(3n+2)=A049450(n)+2*n; b(3n+3)=A049450(n)+4n. If c(n)= exponent of 11 in a(n) factorization: c(n)=4*(0, 0, 0, 0, 0, 1, 2, 3, 4, 5, 7, 9, 11, 13, 15, 18, 21, 24, 27, ..) so c(5n+1)=4*A000566(n); c(5n+2)=4*(A000566(n)+2n); c(5n+3)=4*(A000566(n)+3n); c(5n+4)=4*(A000566(n)+4n); c(5n+5)=4*(A000566(n)+5n)

(PARI) for(n=1, 15, print1((numerator(matdet(matrix(n, n, i, j, 1/(3^j+3^i))))), ", "))

Cf. A069742.

Sequence in context: A051477 A006051 A069742 this_sequence A135824 A013753 A013867

Adjacent sequences: A069740 A069741 A069742 this_sequence A069744 A069745 A069746

easy,nonn

Benoit Cloitre (benoit7848c(AT)orange.fr), Apr 21 2002