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.

A286324 a(n) is the number of bi-unitary divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 6, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 6, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 8, 4, 4, 2, 8, 4, 4, 4, 8, 2, 8
Offset: 1

Views

Author

Michel Marcus, May 07 2017

Keywords

Comments

a(n) is the number of terms of the n-th row of A222266.

Examples

			From _Michael De Vlieger_, May 07 2017: (Start)
a(1) = 1 since 1 is the empty product; all divisors of 1 (i.e., 1) have a greatest common unitary divisor that is 1. 1 is a unitary divisor of all numbers n.
a(p) = 2 since 1 and p have greatest common unitary divisor 1.
a(6) = 4 since the divisor pairs {1, 6} and {2, 3} have greatest common unitary divisor 1.
a(24) = 8 since {1, 24}, {2, 12}, {3, 8}, {4, 6} have greatest unitary divisors {1, {1, 3, 8, 24}}, {{1, 2}, {1, 3, 4, 12}}, {{1, 3}, {1, 8}}, {{1, 4}, {1, 2, 3, 6}}: 1 is the greatest common unitary divisor among all 4 pairs. (End)

Links

Crossrefs

Cf. A222266, A188999, A293185 (indices of records), A340232, A350390.
Cf. A000005, A034444 (unitary), A037445 (infinitary).

Programs

  • Mathematica
    f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; Table[DivisorSum[n, 1 &, Last@ Intersection[f@ #, f[n/#]] == 1 &], {n, 90}] (* Michael De Vlieger, May 07 2017 *)
f[p_, e_] := If[OddQ[e], e + 1, e]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 120] (* Amiram Eldar, Dec 19 2018 *)
  • PARI
    udivs(n) = {my(d = divisors(n)); select(x->(gcd(x, n/x)==1), d); }
gcud(n, m) = vecmax(setintersect(udivs(n), udivs(m)));
biudivs(n) = select(x->(gcud(x, n/x)==1), divisors(n));
a(n) = #biudivs(n);
  • PARI
    a(n)={my(f=factor(n)[,2]); prod(i=1, #f, my(e=f[i]); e + e % 2)} \\ Andrew Howroyd, Aug 05 2018
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (X^3 - X^2 + X + 1) / ((X-1)^2 * (X+1)))[n], ", ")) \\ Vaclav Kotesovec, Jan 11 2024

Formula

Multiplicative with a(p^e) = e + (e mod 2). - Andrew Howroyd, Aug 05 2018
a(A340232(n)) = 2*n. - Bernard Schott, Mar 12 2023
a(n) = A000005(A350390(n)) (the number of divisors of the largest exponentially odd number dividing n). - Amiram Eldar, Sep 01 2023
From Vaclav Kotesovec, Jan 11 2024: (Start)
Dirichlet g.f.: zeta(s)^2 * Product_{p prime} (1 - (p^s - 1)/((p^s + 1)*p^(2*s))).
Let f(s) = Product_{p prime} (1 - (p^s - 1)/((p^s + 1)*p^(2*s))).
Sum_{k=1..n} a(k) ~ f(1) * n * (log(n) + 2*gamma - 1 + f'(1)/f(1)), where
f(1) = Product_{p prime} (1 - (p-1)/((p+1)*p^2)) = A306071 = 0.80733082163620503914865427993003113402584582508155664401800520770441381...,
f'(1) = f(1) * Sum_{p prime} 2*(p^2 - p - 1) * log(p) /(p^4 + 2*p^3 + 1) = f(1) * 0.40523703144422392508596509911218523410441417240419849262346362977537989... = f(1) * A306072
and gamma is the Euler-Mascheroni constant A001620. (End)

A340233 a(n) is the least number with exactly n exponential divisors.

Original entry on oeis.org

1, 4, 16, 36, 65536, 144, 18446744073709551616, 576, 1296, 589824
Offset: 1

Views

Author

Amiram Eldar, Jan 01 2021

Keywords

Comments

a(11) = 2^(2^10) has 309 digits and is too large to be included in the data section.
See the link for more values of this sequence.

Examples

			a(2) = 4 since 4 is the least number with 2 exponential divisors, 2 and 4.

Links

Crossrefs

Subsequence of A025487.
Cf. A049419, A318278, A322791.
Similar sequences: A005179 (all divisors), A038547 (odd divisors), A085629 (coreful divisors), A309181 (nonunitary), A340232 (bi-unitary).

Programs

  • Mathematica
    f[p_, e_] := DivisorSigma[0, e]; d[1] = 1; d[n_] := Times @@ (f @@@ FactorInteger[n]);  max = 6; s = Table[0, {max}]; c = 0; n = 1;  While[c < max, i = d[n]; If[i <= max && s[[i]] == 0, c++; s[[i]] = n]; n++]; s (* ineffective for n > 6 *)

Formula

A049419(a(n)) = n and A049419(k) != n for all k < a(n).

A358252 a(n) is the least number with exactly n non-unitary square divisors.

Original entry on oeis.org

1, 8, 32, 128, 288, 864, 1152, 2592, 4608, 13824, 10368, 20736, 28800, 41472, 64800, 279936, 115200, 331776, 345600, 663552, 259200, 1679616, 518400, 1620000, 1166400, 4860000, 1036800, 17915904, 2073600, 15552000, 6998400, 26873856, 4147200, 53747712, 8294400
Offset: 0

Views

Author

Amiram Eldar, Nov 05 2022

Keywords

Comments

a(n) is the least number k such that A056626(k) = n.
Since A056626(k) depends only on the prime signature of k, all the terms of this sequence are in A025487.

Examples

			a(1) = 8 since 8 is the least number that has exactly one non-unitary square divisor, 4.

Links

Crossrefs

Cf. A056626, A025487, A358253.
Similar sequences: A005179 (all divisors), A038547 (odd divisors), A085629 (coreful divisors), A130279 (square), A187941 (even), A309181 (non-unitary), A340232 (bi-unitary), A340233 (exponential), A357450 (odd square).

Programs

  • Mathematica
    f1[p_, e_] := 1 + Floor[e/2]; f2[p_, e_] := 2^(1 - Mod[e, 2]); f[1] = 0; f[n_] := Times @@ f1 @@@ (fct = FactorInteger[n]) - Times @@ f2 @@@ fct; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[21, 10^6]
  • PARI
    s(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + floor(f[i,2]/2)) - 2^sum(i = 1, #f~, 1 - f[i,2]%2);}
lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};

A358262 a(n) is the least number with exactly n noninfinitary square divisors.

Original entry on oeis.org

1, 16, 144, 256, 3600, 1296, 2304, 65536, 129600, 16777216, 32400, 20736, 57600, 331776, 589824, 4294967296, 6350400, 1099511627776, 150994944, 810000, 1587600, 1679616, 518400, 5308416, 2822400, 84934656, 8294400, 26873856, 14745600, 21743271936, 38654705664
Offset: 0

Views

Author

Amiram Eldar, Nov 06 2022

Keywords

Comments

a(n) is the least number k such that A358261(k) = n.
Since A358261(k) depends only on the prime signature of k, all the terms of this sequence are in A025487.

Examples

			a(1) = 16 since 16 is the least number with exactly one noninfinitary divisor, 4.

Links

Crossrefs

Cf. A025487, A358261, A358263.
Similar sequences: A005179 (all divisors), A038547 (odd divisors), A085629 (coreful divisors), A130279 (square), A187941 (even), A309181 (non-unitary), A340232 (bi-unitary), A340233 (exponential), A357450 (odd square), A358252 (non-unitary square).

Programs

  • Mathematica
    f1[p_, e_] := 1 + Floor[e/2]; f2[p_, e_] := 2^DigitCount[If[OddQ[e], e - 1, e], 2, 1]; f[1] = 0; f[n_] := Times @@ f1 @@@ (fct = FactorInteger[n]) - Times @@ f2 @@@ fct; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[15, 2*10^7]
  • PARI
    s(n) = {my(f = factor(n));  prod(i=1, #f~, 1+f[i,2]\2) - prod(i=1, #f~, 2^hammingweight(if(f[i,2]%2, f[i,2]-1, f[i,2])))};
lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};

A361418 a(n) is the least number with exactly n noninfinitary divisors.

Original entry on oeis.org

1, 4, 12, 16, 60, 36, 48, 256, 360, 4096, 180, 144, 240, 576, 768, 65536, 2520, 1048576, 12288, 900, 1260, 1296, 720, 2304, 1680, 9216, 2880, 5184, 3840, 147456, 196608, 36864, 27720, 46656, 3145728, 4398046511104, 61440, 3600, 6300, 18014398509481984, 10080, 20736
Offset: 0

Views

Author

Amiram Eldar, Mar 11 2023

Keywords

Comments

a(n) is the least number k such that A348341(k) = n.
Since A348341(k) depends only on the prime signature of k, all the terms of this sequence are in A025487.

Examples

			a(1) = 4 since 4 is the least number with exactly one noninfinitary divisor, 2.

Links

Crossrefs

Cf. A025487, A348341, A348342.
Similar sequences: A005179 (all divisors), A038547 (odd divisors), A085629 (coreful divisors), A130279 (square), A187941 (even), A309181 (non-unitary), A340232 (bi-unitary), A340233 (exponential), A357450 (odd square), A358252 (non-unitary square).

Programs

  • Mathematica
    f[1] = 0; f[n_] := DivisorSigma[0, n] - Times @@ Flatten[2^DigitCount[#, 2, 1] & /@ FactorInteger[n][[;; , 2]]];
seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s];
seq[35, 10^7]
  • PARI
    s(n) = {my(f = factor(n)); numdiv(f) - prod(i = 1, #f~, 2^hammingweight(f[i,2]));}
lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.