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.

A261890 Second differences of A055615, first differences of A261869.

Original entry on oeis.org

2, 4, -8, 16, -24, 20, -7, 10, -31, 32, -24, 40, -26, -16, -2, 34, -36, 38, 2, -20, -46, 68, -23, 26, -52, 26, -29, 28, 0, 32, 2, -32, 0, -36, -2, 112, -74, -40, -2, 40, 0, 44, -43, 46, -139, 140, -47, 0, 51, -102, -2, 106, 2, -110, 112, -56, -118, 176, -120
Offset: 1

Views

Author

Reinhard Zumkeller, Sep 05 2015

Keywords

Comments

a(n) = (n+2)*mu(n+2) - 2*(n+1)*mu(n+1) + n*mu(n), where mu = A008683 (Moebius function);
a(A114180(n)) = 0;
abs(a(n)) <= 4*(n+1).

Crossrefs

Programs

  • Haskell
    import Data.List (genericIndex)
    a261890 n = genericIndex a261890_list (n - 1)
    a261890_list = zipWith (-) (tail a261869_list) a261869_list
  • Mathematica
    Differences[Table[n MoebiusMu[n],{n,80}],2] (* Harvey P. Dale, Mar 12 2023 *)

Formula

a(n) = A055615(n+2) - 2*A055615(n+1) + A055615(n).

A055615 a(n) = n * mu(n), where mu is the Möbius function A008683.

Original entry on oeis.org

1, -2, -3, 0, -5, 6, -7, 0, 0, 10, -11, 0, -13, 14, 15, 0, -17, 0, -19, 0, 21, 22, -23, 0, 0, 26, 0, 0, -29, -30, -31, 0, 33, 34, 35, 0, -37, 38, 39, 0, -41, -42, -43, 0, 0, 46, -47, 0, 0, 0, 51, 0, -53, 0, 55, 0, 57, 58, -59, 0, -61, 62, 0, 0, 65, -66, -67, 0, 69, -70, -71, 0
Offset: 1

Views

Author

Michael Somos, Jun 04 2000

Keywords

Comments

Dirichlet inverse of n (A000027).
Absolute values give n if n is squarefree, otherwise 0.
a(n) is multiplicative because both mu(n) and n are. - Mitch Harris, Jun 09 2005
a(n) is multiplicative with a(p^1) = -p, a(p^e) = 0 if e > 1. - David W. Wilson, Jun 12 2005
Negative of the Moebius number of the dihedral group of order 2n. - Eric M. Schmidt, Jul 28 2013

Examples

			G.f. = x - 2*x^2 - 3*x^3 - 5*x^5 + 6*x^6 - 7*x^7 + 10*x^10 - 11*x^11 - 13*x^13 + ...
		

Crossrefs

Moebius transform of A023900.
Cf. A000027 (Dirichlet inverse), A061669 (sum with it).
Cf. A062004.
Cf. A013929 (positions of 0's), A068340 (partial sums), A261869 (first differences), A261890 (second differences).

Programs

  • Haskell
    a055615 n = a008683 n * n  -- Reinhard Zumkeller, Sep 04 2015
    
  • Magma
    [n*MoebiusMu(n): n in [1..80]]; // Vincenzo Librandi, Nov 19 2014
    
  • Maple
    with(numtheory): A055615:=n->n*mobius(n): seq(A055615(n), n=1..100); # Wesley Ivan Hurt, Nov 18 2014
  • Mathematica
    Table[n MoebiusMu[n], {n,80}] (* Harvey P. Dale, May 26 2011 *)
  • PARI
    {a(n) = if( n<1, 0, n * moebius(n))};
    
  • PARI
    {a(n) = if( n<1, 0, direuler(p=2, n, 1 - p*X)[n])};
    
  • Python
    from sympy import mobius
    def A055615(n): return n*mobius(n) # Chai Wah Wu, Apr 01 2023
  • SageMath
    [n*moebius(n) for n in (1..100)] # G. C. Greubel, May 24 2022
    

Formula

a(n) = n * A008683(n).
Dirichlet g.f.: 1/zeta(s-1).
Multiplicative with a(p^e) = -p*0^(e-1), e>0 and p prime. - Reinhard Zumkeller, Jul 17 2003
Conjectures: lim b->1+ Sum n=1..inf a(n)*b^(-n) = -12 and lim b->1- Sum n=1..inf a(n)*b^n = -12 (+ indicates that b decreases to 1, - indicates it increases to 1), both considering that zeta(-1) = -1/12 and calculations (more generally mu(n)*n^s is Abel summable to zeta(-s)). - Gerald McGarvey, Sep 26 2004
Dirichlet generating function for the absolute value: zeta(s-1)/zeta(2s-2). - Franklin T. Adams-Watters, Sep 11 2005
G.f. A(x) satisfies: A(x) = x - Sum_{k>=2} k*A(x^k). - Ilya Gutkovskiy, May 11 2019
Sum_{k=1..n} abs(a(k)) ~ 3*n^2/Pi^2. - Amiram Eldar, Feb 02 2024

A068781 Lesser of two consecutive numbers each divisible by a square.

Original entry on oeis.org

8, 24, 27, 44, 48, 49, 63, 75, 80, 98, 99, 116, 120, 124, 125, 135, 147, 152, 168, 171, 175, 188, 207, 224, 242, 243, 244, 260, 275, 279, 288, 296, 315, 324, 332, 342, 343, 350, 351, 360, 363, 368, 375, 387, 404, 423, 424, 440, 459, 475, 476, 495, 507, 512
Offset: 1

Views

Author

Robert G. Wilson v, Mar 04 2002

Keywords

Comments

Also numbers m such that mu(m)=mu(m+1)=0, where mu is the Moebius-function (A008683); A081221(a(n))>1. - Reinhard Zumkeller, Mar 10 2003
The sequence contains an infinite family of arithmetic progressions like {36a+8}={8,44,80,116,152,188,...} ={4(9a+2)}. {36a+9} provides 2nd nonsquarefree terms. Such AP's can be constructed to any term by solution of a system of linear Diophantine equation. - Labos Elemer, Nov 25 2002
1. 4k^2 + 4k is a member for all k; i.e., 8 times a triangular number is a member. 2. (4k+1) times an odd square - 1 is a member. 3. (4k+3) times odd square is a member. - Amarnath Murthy, Apr 24 2003
The asymptotic density of this sequence is 1 - 2/zeta(2) + Product_{p prime} (1 - 2/p^2) = 1 - 2 * A059956 + A065474 = 0.1067798952... (Matomäki et al., 2016). - Amiram Eldar, Feb 14 2021
Maximum of the n-th maximal anti-run of nonsquarefree numbers (A013929) differing by more than one. For runs instead of anti-runs we have A376164. For squarefree instead of nonsquarefree we have A007674. - Gus Wiseman, Sep 14 2024

Examples

			44 is in the sequence because 44 = 2^2 * 11 and 45 = 3^2 * 5.
From _Gus Wiseman_, Sep 14 2024: (Start)
Splitting nonsquarefree numbers into maximal anti-runs gives:
  (4,8)
  (9,12,16,18,20,24)
  (25,27)
  (28,32,36,40,44)
  (45,48)
  (49)
  (50,52,54,56,60,63)
  (64,68,72,75)
  (76,80)
  (81,84,88,90,92,96,98)
  (99)
The maxima are a(n). The corresponding pairs are (8,9), (24,25), (27,28), (44,45), etc.
(End)
		

Crossrefs

Subsequence of A261869.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
A053797 gives lengths of runs of nonsquarefree numbers, firsts A373199.

Programs

  • Haskell
    a068781 n = a068781_list !! (n-1)
    a068781_list = filter ((== 0) . a261869) [1..]
    -- Reinhard Zumkeller, Sep 04 2015
    
  • Mathematica
    Select[ Range[2, 600], Max[ Transpose[ FactorInteger[ # ]] [[2]]] > 1 && Max[ Transpose[ FactorInteger[ # + 1]] [[2]]] > 1 &]
    f@n_:= Flatten@Position[Partition[SquareFreeQ/@Range@2000,n,1], Table[False,{n}]]; f@2 (* Hans Rudolf Widmer, Aug 30 2022 *)
    Max/@Split[Select[Range[100], !SquareFreeQ[#]&],#1+1!=#2&]//Most (* Gus Wiseman, Sep 14 2024 *)
  • PARI
    isok(m) = !moebius(m) && !moebius(m+1); \\ Michel Marcus, Feb 14 2021

Formula

A261869(a(n)) = 0. - Reinhard Zumkeller, Sep 04 2015
Showing 1-3 of 3 results.