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.

Previous Showing 11-13 of 13 results.

A286154 Compound filter: a(n) = T(A055396(n), A000010(n)), where T(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

0, 1, 5, 2, 18, 2, 40, 7, 23, 7, 96, 7, 142, 16, 38, 29, 238, 16, 308, 29, 80, 46, 444, 29, 234, 67, 173, 67, 676, 29, 791, 121, 212, 121, 328, 67, 1093, 154, 302, 121, 1339, 67, 1499, 191, 302, 232, 1785, 121, 994, 191, 530, 277, 2227, 154, 864, 277, 668, 379, 2718, 121, 2944, 436, 668, 497, 1228, 191, 3505, 497, 992, 277, 3936, 277, 4207, 631, 822, 631
Offset: 1

Views

Author

Antti Karttunen, May 04 2017

Keywords

Links

Crossrefs

Cf. A000010, A000027, A055396, A286142, A286143, A286144, A286152, A286160, A286161, A286162, A286163, A286164.

Programs

  • Mathematica
    Table[(2 + (#1 + #2)^2 - #1 - 3 #2)/2 & @@ {If[n == 1, 0, PrimePi[ FactorInteger[n][[1, 1]] ]], EulerPhi@ n}, {n, 76}] (* Michael De Vlieger, May 04 2017 *)
  • PARI
    A000010(n) = eulerphi(n);
A055396(n) = if(n==1, 0, primepi(factor(n)[1, 1])); \\ This function from Charles R Greathouse IV, Apr 23 2015
A286154(n) = (2 + ((A055396(n)+A000010(n))^2) - A055396(n) - 3*A000010(n))/2;
for(n=1, 10000, write("b286154.txt", n, " ", A286154(n)));
  • Python
    from sympy import primepi, isprime, primefactors, totient
def a049084(n): return primepi(n)*(1*isprime(n))
def a055396(n): return 0 if n==1 else a049084(min(primefactors(n)))
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def a(n): return T(a055396(n), totient(n)) # Indranil Ghosh, May 05 2017
  • Scheme
    (define (A286154 n) (* (/ 1 2) (+ (expt (+ (A055396 n) (A000010 n)) 2) (- (A055396 n)) (- (* 3 (A000010 n))) 2)))

Formula

a(n) = (1/2)*(2 + ((A055396(n)+A000010(n))^2) - A055396(n) - 3*A000010(n)).

A286253 Compound filter: a(n) = P(A055396(n), A001511(1+n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

0, 1, 8, 1, 9, 1, 25, 1, 5, 1, 26, 1, 27, 1, 17, 1, 35, 1, 53, 1, 5, 1, 75, 1, 9, 1, 8, 1, 65, 1, 131, 1, 5, 1, 13, 1, 90, 1, 12, 1, 104, 1, 134, 1, 5, 1, 186, 1, 14, 1, 8, 1, 152, 1, 18, 1, 5, 1, 188, 1, 189, 1, 30, 1, 9, 1, 229, 1, 5, 1, 273, 1, 252, 1, 8, 1, 14, 1, 347, 1, 5, 1, 323, 1, 9, 1, 12, 1, 324, 1, 19, 1, 5, 1, 31, 1, 350, 1, 8, 1, 377, 1, 462, 1, 5
Offset: 1

Views

Author

Antti Karttunen, May 07 2017

Keywords

Links

Crossrefs

Cf. A000027, A001511, A055396, A286164, A286251, A286252, A286254.

Programs

  • PARI
    A001511(n) = (1+valuation(n,2));
A055396(n) = if(n==1, 0, primepi(factor(n)[1, 1])); \\ This function from Charles R Greathouse IV, Apr 23 2015
A286253(n) = (2 + ((A055396(n)+A001511(1+n))^2) - A055396(n) - 3*A001511(1+n))/2;
for(n=1, 10000, write("b286253.txt", n, " ", A286253(n)));
  • Python
    from sympy import primepi, isprime, primefactors
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def a049084(n): return primepi(n)*(1*isprime(n))
def a055396(n): return 0 if n==1 else a049084(min(primefactors(n)))
def a001511(n): return 2 + bin(n - 1)[2:].count("1") - bin(n)[2:].count("1")
def a(n): return T(a055396(n), a001511(n + 1)) # Indranil Ghosh, May 07 2017
  • Scheme
    (define (A286253 n) (* (/ 1 2) (+ (expt (+ (A055396 n) (A001511 (+ 1 n))) 2) (- (A055396 n)) (- (* 3 (A001511 (+ 1 n)))) 2)))

Formula

a(n) = (1/2)*(2 + ((A055396(n)+A001511(1+n))^2) - A055396(n) - 3*A001511(1+n)).

A286149 Compound filter: a(n) = T(A046523(n), A109395(n)), where T(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

1, 5, 8, 14, 17, 34, 30, 44, 19, 51, 68, 103, 93, 72, 196, 152, 155, 103, 192, 132, 72, 126, 278, 349, 32, 159, 53, 165, 437, 976, 498, 560, 709, 237, 786, 739, 705, 282, 159, 402, 863, 660, 948, 243, 337, 384, 1130, 1273, 49, 132, 1546, 288, 1433, 349, 126, 459, 282, 567, 1772, 2761, 1893, 636, 165, 2144, 2421, 1921, 2280, 390, 2707, 2046, 2558, 2773, 2703
Offset: 1

Views

Author

Antti Karttunen, May 04 2017

Keywords

Links

Crossrefs

Cf. A000027, A046523, A109395, A285729, A286142, A286143, A286144, A286152, A286154, A286160, A286161, A286162, A286163, A286164.

Programs

  • Mathematica
    Table[(2 + (#1 + #2)^2 - #1 - 3 #2)/2 & @@ {Times @@ MapIndexed[ Prime[First@ #2]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]] - Boole[n == 1], Denominator[EulerPhi[n]/n]}, {n, 73}] (* Michael De Vlieger, May 04 2017 *)
  • PARI
    A109395(n) = n/gcd(n, eulerphi(n));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ This function from Charles R Greathouse IV, Aug 17 2011
A286149(n) = (1/2)*(2 + ((A046523(n)+A109395(n))^2) - A046523(n) - 3*A109395(n));
for(n=1, 10000, write("b286149.txt", n, " ", A286149(n)));
  • Python
    from sympy import factorint, totient, gcd
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def P(n):
    f = factorint(n)
    return sorted([f[i] for i in f])
def a046523(n):
    x=1
    while True:
        if P(n) == P(x): return x
        else: x+=1
def a(n): return T(a046523(n), n/gcd(n, totient(n))) # Indranil Ghosh, May 05 2017
  • Scheme
    (define (A286149 n) (* (/ 1 2) (+ (expt (+ (A046523 n) (A109395 n)) 2) (- (A046523 n)) (- (* 3 (A109395 n))) 2)))

Formula

a(n) = (1/2)*(2 + ((A046523(n)+A109395(n))^2) - A046523(n) - 3*A109395(n)).
Previous Showing 11-13 of 13 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.