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-5 of 5 results.

A258409 Greatest common divisor of all (d-1)'s, where the d's are the positive divisors of n.

Original entry on oeis.org

1, 2, 1, 4, 1, 6, 1, 2, 1, 10, 1, 12, 1, 2, 1, 16, 1, 18, 1, 2, 1, 22, 1, 4, 1, 2, 1, 28, 1, 30, 1, 2, 1, 2, 1, 36, 1, 2, 1, 40, 1, 42, 1, 2, 1, 46, 1, 6, 1, 2, 1, 52, 1, 2, 1, 2, 1, 58, 1, 60, 1, 2, 1, 4, 1, 66, 1, 2, 1, 70, 1, 72, 1, 2, 1, 2, 1, 78, 1
Offset: 2

Views

Author

Ivan Neretin, May 29 2015

Keywords

Comments

a(n) = 1 for even n; a(p) = p-1 for prime p.
a(n) is even for odd n (since all divisors of n are odd).
It appears that a(n) = A052409(A005179(n)), i.e., it is the largest integer power of the smallest number with exactly n divisors. - Michel Marcus, Nov 10 2015
Conjecture: GCD of all (p-1) for prime p|n. - Thomas Ordowski, Sep 14 2016
Conjecture is true, because the set of numbers == 1 (mod g) is closed under multiplication. - Robert Israel, Sep 14 2016
Conjecture: a(n) = A289508(A328023(n)) = GCD of the differences between consecutive divisors of n. See A328163 and A328164. - Gus Wiseman, Oct 16 2019

Examples

			65 has divisors 1, 5, 13, and 65, hence a(65) = gcd(1-1,5-1,13-1,65-1) = gcd(0,4,12,64) = 4.
		

Crossrefs

Cf. A084190 (similar but with LCM).
Looking at prime indices instead of divisors gives A328167.
Partitions whose parts minus 1 are relatively prime are A328170.

Programs

  • Haskell
    a258409 n = foldl1 gcd $ map (subtract 1) $ tail $ a027750_row' n
    -- Reinhard Zumkeller, Jun 25 2015
  • Maple
    f:= n -> igcd(op(map(`-`,numtheory:-factorset(n),-1))):
    map(f, [$2..100]); # Robert Israel, Sep 14 2016
  • Mathematica
    Table[GCD @@ (Divisors[n] - 1), {n, 2, 100}]
  • PARI
    a(n) = my(g=0); fordiv(n, d, g = gcd(g, d-1)); g; \\ Michel Marcus, May 29 2015
    
  • PARI
    a(n) = gcd(apply(x->x-1, divisors(n))); \\ Michel Marcus, Nov 10 2015
    
  • PARI
    a(n)=if(n%2==0, return(1)); if(n%3==0, return(2)); if(n%5==0 && n%4 != 1, return(2)); gcd(apply(p->p-1, factor(n)[,1])) \\ Charles R Greathouse IV, Sep 19 2016
    

A328026 Number of divisible pairs of consecutive divisors of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 03 2019

Keywords

Comments

The number m = 2^n, n >= 0, is the smallest for which a(m) = n. - Marius A. Burtea, Nov 20 2019

Examples

			The divisors of 500 are {1,2,4,5,10,20,25,50,100,125,250,500}, with consecutive divisible pairs {1,2}, {2,4}, {5,10}, {10,20}, {25,50}, {50,100}, {125,250}, {250,500}, so a(500) = 8.
		

Crossrefs

Positions of 1's are A000040.
Positions of 0's and 2's are A328028.
Positions of terms > 2 are A328189.
Successive pairs of consecutive divisors are counted by A129308.

Programs

  • Magma
    f:=func;  g:=func; [g(n):n in [1..100]]; // Marius A. Burtea, Nov 20 2019
  • Mathematica
    Table[Length[Split[Divisors[n],!Divisible[#2,#1]&]]-1,{n,100}]
  • PARI
    a(n) = {my(d=divisors(n), nb=0); for (i=2, #d, if ((d[i] % d[i-1]) == 0, nb++)); nb;} \\ Michel Marcus, Oct 05 2019
    

Formula

a(p^k) = k for any prime number p and k >= 0. - Rémy Sigrist, Oct 05 2019

Extensions

Data section extended up to a(105) by Antti Karttunen, Feb 23 2023

A060766 Least common multiple of differences between consecutive divisors of n (ordered by size).

Original entry on oeis.org

1, 2, 2, 4, 3, 6, 4, 6, 15, 10, 6, 12, 35, 10, 8, 16, 9, 18, 10, 28, 99, 22, 12, 20, 143, 18, 42, 28, 60, 30, 16, 88, 255, 28, 18, 36, 323, 130, 60, 40, 21, 42, 154, 60, 483, 46, 24, 42, 75, 238, 234, 52, 27, 132, 84, 304, 783, 58, 60, 60, 899, 84, 32, 104, 165, 66, 442
Offset: 2

Views

Author

Labos Elemer, Apr 24 2001

Keywords

Examples

			For n=98, divisors={1,2,7,14,49,98}; differences={1,5,7,35,49}; a(98) = LCM of differences = 245.
		

Crossrefs

The GCD version appears to be A258409.
The LCM of the prime indices of n is A290103(n).
The differences between consecutive divisors of n are row n of A193829.

Programs

  • Mathematica
    a[n_ ] := LCM@@(Drop[d=Divisors[n], 1]-Drop[d, -1])
    Table[LCM@@Differences[Divisors[n]],{n,2,70}] (* Harvey P. Dale, Oct 08 2012 *)

Formula

a(n) = A290103(A328023(n)). - Gus Wiseman, Oct 16 2019

Extensions

Edited by Dean Hickerson, Jan 22 2002

A328024 Heinz numbers of multisets representing the differences between some positive integer's consecutive divisors.

Original entry on oeis.org

1, 2, 3, 6, 7, 13, 20, 29, 37, 39, 42, 53, 61, 79, 107, 110, 113, 151, 173, 181, 199, 239, 261, 271, 281, 312, 317, 349, 359, 374, 397, 421, 457, 497, 503, 541, 557, 577, 593, 613, 701, 733, 769, 787, 798, 857, 863, 903, 911, 953, 983, 1021, 1061, 1069, 1151
Offset: 1

Views

Author

Gus Wiseman, Oct 02 2019

Keywords

Comments

The Heinz number of an integer partition or multiset {y_1,...,y_k} is prime(y_1)*...*prime(y_k).
There is exactly one entry with any given sum of prime indices A056239.

Examples

			The sequence of terms together with their prime indices begins:
     1: {}
     2: {1}
     3: {2}
     6: {1,2}
     7: {4}
    13: {6}
    20: {1,1,3}
    29: {10}
    37: {12}
    39: {2,6}
    42: {1,2,4}
    53: {16}
    61: {18}
    79: {22}
   107: {28}
   110: {1,3,5}
   113: {30}
   151: {36}
   173: {40}
   181: {42}
   199: {46}
   239: {52}
   261: {2,2,10}
   271: {58}
   281: {60}
   312: {1,1,1,2,6}
For example, the divisors of 8 are {1,2,4,8}, with differences {1,2,4}, with Heinz number 42, so 42 belongs to the sequence.
		

Crossrefs

A permutation of A328023.
Also the set of possible Heinz numbers of rows of A193829, A328025, or A328027.

Programs

  • Mathematica
    nn=1000;
    Select[Union[Table[Times@@Prime/@Differences[Divisors[n]],{n,nn}]],#<=nn&]

A328025 Irregular triangle read by rows where row n gives the differences between consecutive divisors of n in weakly decreasing order.

Original entry on oeis.org

1, 2, 2, 1, 4, 3, 1, 1, 6, 4, 2, 1, 6, 2, 5, 3, 1, 10, 6, 2, 1, 1, 1, 12, 7, 5, 1, 10, 2, 2, 8, 4, 2, 1, 16, 9, 3, 3, 1, 1, 18, 10, 5, 2, 1, 1, 14, 4, 2, 11, 9, 1, 22, 12, 4, 2, 2, 1, 1, 1, 20, 4, 13, 11, 1, 18, 6, 2, 14, 7, 3, 2, 1, 28, 15, 5, 4, 2, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 02 2019

Keywords

Examples

			Triangle begins:
   {}
   1
   2
   2  1
   4
   3  1  1
   6
   4  2  1
   6  2
   5  3  1
  10
   6  2  1  1  1
  12
   7  5  1
  10  2  2
   8  4  2  1
  16
   9  3  3  1  1
  18
  10  5  2  1  1
For example, the divisors of 18 are {1,2,3,6,9,18}, with differences {1,1,3,3,9}, so row 18 is {9,3,3,1,1}.
		

Crossrefs

Same as A193829 with rows sorted in weakly decreasing order.
Same as A328027 with rows reversed.
Row sums are A001477.
Row lengths are A000005.
First column is A060681.
Heinz numbers of rows are A328023.

Programs

  • Mathematica
    Table[Sort[Differences[Divisors[n]],Greater],{n,30}]
Showing 1-5 of 5 results.