cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A348999 a(n) = A348929(A276086(n)), where A348929(n) = gcd(n, A003959(n)), A003959 is multiplicative with a(p^e) = (p+1)^e, and A276086 gives the prime product form of primorial base expansion of n.

Original entry on oeis.org

1, 1, 1, 6, 1, 6, 1, 2, 3, 6, 3, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 1, 6, 1, 6, 1, 2, 3, 6, 3, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 1, 6, 1, 6, 1, 2, 3, 6, 3, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 3, 6, 9, 18, 1, 2, 1, 6, 1, 6, 1, 2, 3, 6, 3, 18
Offset: 0

Views

Author

Antti Karttunen, Nov 07 2021

Keywords

Comments

After each primorial number (A002110), the apparent periodicity grows more complex.

Crossrefs

Programs

  • PARI
    A348999(n) = { my(m1=1, m2=1, p=2); while(n, m1 *= (p^(n%p)); m2 *= ((1+p)^(n%p)); n = n\p; p = nextprime(1+p)); gcd(m1,m2); };

Formula

a(n) = A348929(A276086(n)).
a(n) = gcd(A276086(n), A348949(n)) = gcd(A276086(n), A348950(n)).

A348928 a(n) = gcd(n, A003958(n)), where A003958 is multiplicative with a(p^e) = (p-1)^e.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 4, 3, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 3, 4, 1, 6, 1, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 5, 2, 3, 2, 1, 4, 1, 2, 3, 1, 1, 2, 1, 4, 1, 2, 1, 4, 1, 2, 1, 2, 1, 6, 1, 4, 1, 2, 1, 12, 1, 2, 1, 2, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 4, 3
Offset: 1

Views

Author

Antti Karttunen, Nov 07 2021

Keywords

Crossrefs

Differs from similar A126864 for the first time at n=36, where a(36) = 4, while A126864(36) = 2.

Programs

  • Mathematica
    f[p_, e_] := (p - 1)^e; a[n_] := GCD[n, Times @@ f @@@ FactorInteger[n]]; Array[a, 100] (* Amiram Eldar, Nov 07 2021 *)
  • PARI
    A003958(n) = if(1==n,n,my(f=factor(n)); for(i=1,#f~,f[i,1]--); factorback(f));
    A348928(n) = gcd(n, A003958(n));

Formula

a(n) = gcd(n, A003958(n)) = gcd(n, A322582(n)) = gcd(A003958(n), A322582(n)).
Showing 1-2 of 2 results.