login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A064124
Number of divisors of 7^n - 1 that are relatively prime to 7^m - 1 for all 0 < m < n.
0
4, 1, 2, 3, 2, 2, 4, 2, 4, 4, 4, 4, 2, 4, 4, 4, 4, 2, 4, 4, 2, 4, 8, 8, 4, 4, 16, 2, 8, 2, 8, 4, 8, 2, 4, 2, 8, 4, 8, 4, 8, 4, 4, 8, 4, 2, 4, 2, 16, 2, 16, 4, 8, 2, 2, 16, 8, 2, 16, 4, 16, 8, 2, 8, 32, 4, 16, 8, 8, 16, 4, 8, 32, 4, 4, 2, 8, 8, 4, 8, 16, 16, 128, 2, 8
OFFSET
1,1
MATHEMATICA
a = {1}; Do[ d = Divisors[ 7^n - 1 ]; l = Length[ d ]; c = 0; k = 1; While[ k < l + 1, If[ Union[ GCD[ a, d[ [ k ] ] ] ] == {1}, c++ ]; k++ ]; Print[ c ]; a = Union[ Flatten[ Append[ a, Transpose[ FactorInteger[ 7^n - 1 ] ][ [ 1 ] ] ] ] ], {n, 1, 70} ]
PROG
(PARI) a(n) = sumdiv(7^n-1, d, vecsum(vector(n-1, k, gcd(d, 7^k-1) == 1)) == n-1); \\ Michel Marcus, Jun 23 2018
CROSSREFS
Cf. A063982.
Sequence in context: A378626 A321091 A375576 * A106654 A252094 A196138
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Sep 10 2001
EXTENSIONS
a(71)-a(85) from Giovanni Resta, Jun 26 2018
STATUS
approved