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

A048103 Numbers not divisible by p^p for any prime p.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

If a(n) = Product p_i^e_i then p_i > e_i for all i.
Complement of A100716; A129251(a(n)) = 0. - Reinhard Zumkeller, Apr 07 2007
Density is 0.72199023441955... = Product_{p>=2} (1 - p^-p) where p runs over the primes. - Charles R Greathouse IV, Jan 25 2012
A027748(a(n),k) <= A124010(a(n),k), 1<=k<=A001221(a(n)). - Reinhard Zumkeller, Apr 28 2012
Range of A276086. Also numbers not divisible by m^m for any natural number m > 1. - Antti Karttunen, Nov 18 2024

Examples

			6 = 2^1 * 3^1 is OK but 12 = 2^2 * 3^1 is not.
625 = 5^4 is present because it is not divisible by 5^5.

Links

Crossrefs

Complement: A100716.
Positions of 0's in A129251, A342023, A376418, positions of 1's in A327936, A342007, A359550 (characteristic function).
Cf. A048102, A048104, A051674 (p^p), A054743, A054744, A377982 (a left inverse, partial sums of char. fun, see also A328402).
Cf. A276086 (permutation of this sequence, see also A376411, A376413).
Subsequences: A002110, A005117, A006862, A024451 (after its initial 0), A057588, A099308 (after its initial 0), A276092, A328387, A328832, A359547, A370114, A371083, A373848, A377871, A377992.
Disjoint union of {1}, A327934 and A358215.
Also A276078 is a subsequence, from which this differs for the first time at n=451 where a(451)=625, while that value is missing from A276078.

Programs

  • Haskell
    a048103 n = a048103_list !! (n-1)
a048103_list = filter (\x -> and $
   zipWith (>) (a027748_row x) (map toInteger $ a124010_row x)) [1..]
-- Reinhard Zumkeller, Apr 28 2012
  • Mathematica
    {1}~Join~Select[Range@ 120, Times @@ Boole@ Map[First@ # > Last@ # &, FactorInteger@ #] > 0 &] (* Michael De Vlieger, Aug 19 2016 *)
  • PARI
    isok(n) = my(f=factor(n)); for (i=1, #f~, if (f[i,1] <= f[i,2], return(0))); return(1); \\ Michel Marcus, Nov 13 2020
  • PARI
    A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); }; \\ (A359550 is the characteristic function for A048103) - Antti Karttunen, Nov 18 2024
  • Python
    from itertools import count, islice
from sympy import factorint
def A048103_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(map(lambda d:d[1]A048103_list = list(islice(A048103_gen(),30)) # Chai Wah Wu, Jan 05 2023
  • Scheme
    ;; With Antti Karttunen's IntSeq-library.
(define A048103 (ZERO-POS 1 1 A129251))
;; Antti Karttunen, Aug 18 2016

Formula

a(n) ~ kn with k = 1/Product_{p>=2}(1 - p^-p) = Product_{p>=2}(1 + 1/(p^p - 1)) = 1.3850602852..., where the product is over all primes p. - Charles R Greathouse IV, Jan 25 2012
For n >= 1, A377982(a(n)) = n. - Antti Karttunen, Nov 18 2024

Extensions

More terms from James Sellers, Apr 22 2000

A377982 a(n) is the number of terms of A048103 (numbers with no p^p divisors) less than or equal to n.

Original entry on oeis.org

0, 1, 2, 3, 3, 4, 5, 6, 6, 7, 8, 9, 9, 10, 11, 12, 12, 13, 14, 15, 15, 16, 17, 18, 18, 19, 20, 20, 20, 21, 22, 23, 23, 24, 25, 26, 26, 27, 28, 29, 29, 30, 31, 32, 32, 33, 34, 35, 35, 36, 37, 38, 38, 39, 39, 40, 40, 41, 42, 43, 43, 44, 45, 46, 46, 47, 48, 49, 49, 50, 51, 52, 52, 53, 54, 55, 55, 56, 57, 58, 58, 58, 59
Offset: 0

Views

Author

Antti Karttunen, Nov 13 2024

Keywords

Comments

The number of terms in the range of A276086 that are <= n.

Links

Crossrefs

Partial sums of A359550.
Cf. A048103, A276086, A376411.
Cf. also A377983.

Programs

  • PARI
    up_to = 10000;
A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); };
A377982list(up_to) = { my(v=vector(1+up_to), s=0); v[1] = s = A359550(0); for(n=1,up_to,s += A359550(n); v[1+n] = s); (v); };
v377982 = A377982list(up_to);
A377982(n) = v377982[1+n];

Formula

a(n) = Sum_{i=0..n} A359550(i).

A376413 a(n) = A276085(A048103(n)), where A276085 is the primorial base log-function, and A048103 is the range of the primorial base exp-function (A276086).

Original entry on oeis.org

0, 1, 2, 6, 3, 30, 4, 7, 210, 2310, 31, 8, 30030, 5, 510510, 32, 211, 9699690, 12, 2311, 223092870, 9, 6469693230, 212, 30031, 36, 200560490130, 510511, 2312, 7420738134810, 33, 304250263527210, 10, 9699691, 13082761331670030, 60, 13, 30032, 614889782588491410, 216, 510512, 223092871, 32589158477190044730
Offset: 1

Views

Author

Antti Karttunen, Nov 13 2024

Keywords

Comments

Permutation of nonnegative integers.

Links

Crossrefs

Cf. A048103, A276085, A276086.
Cf. A376411 (inverse permutation, but note the different offsets and ranges).
Cf. also A048672 (analogous permutation for base-2).

Programs

  • PARI
    A276085(n) = { my(f = factor(n), pr=1, i=1, s=0); for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); (s); };
A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); }; \\ Antti Karttunen, Nov 13 2024
for(n=1,105,if(A359550(n),print1(A276085(n),", ")));

Formula

a(n) = A276085(A048103(n)).
For all n >= 1, A376411(a(n)) = n-1.
Showing 1-3 of 3 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.