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.

A341860 a(n) is the least k such that A341857(k) = psi(k^2)/k = n, or 0 if no such k exists.

Original entry on oeis.org

1, 3, 14, 5, 22, 7, 116, 32, 38, 11, 46, 13, 212, 58, 62, 17, 618, 19, 2748, 25, 86, 23, 94, 219, 404, 106, 162, 29, 118, 31, 4476, 128, 134, 309, 142, 37, 596, 1374, 158, 41, 166, 43, 692, 115, 418, 47, 1698, 119, 788, 202, 206, 53, 214, 81, 242, 226, 916, 59
Offset: 1

Views

Author

Jianing Song, Feb 21 2021

Keywords

Comments

Conjecture: a(n) > 0 for all n.
If there exists k such that psi(k^2) | k*n and k*n+1 is prime, then A341857(k*(k*n+1)) = n, so a(n) > 0.
For odd n > 1, if a(n) > 0, then a(2n) <= a(n)/2. Proof:
Suppose k is the least integer such that A341857(k) = psi(k^2)/k = n is odd. Note that if m > 1 is odd or m is divisible by 8, then A341857(m) is even. So we must have k = 2t or k = 4t with odd t > 1.
If k = 2t, then n = psi(4t^2)/(2t) = psi(t^2)/(2t), so A341857(t) = psi(t^2)/t = 2n. This gives a(2n) <= t = k/2.
if k = 4t, then n = psi(16t^2)/(4t) = lcm(4, psi(t^2))/(4t). There are two cases: a) if psi(t^2) is divisible by 4, then psi(t^2)/(4t) = n, so A341857(2t) = psi(4t^2)/(2t) = psi(t^2)/(2t) = 2n. This gives a(2n) <= 2t = k/2. b) if psi(t^2) is not divisible by 4, then psi(t^2)/(2t) = n, so A341857(2t) = psi(4t^2)/(2t) = psi(t^2)/(2t) = n. This gives a(n) <= 2t = k/2, contradicting with minimality of k. QED.
The smallest odd n such that a(2n) < a(n)/2 is n = 71, where a(71) = 10236 and a(142) = 2276.

Examples

			A341857(116) = psi(116^2)/116 = 812/116 = 7, and there is no k < 116 such that psi(k^2)/k = 7, so a(7) = 116.
		

Crossrefs

Cf. A341857.

Programs

  • Mathematica
    A341860[n_] := Block[{k = 0}, While[CarmichaelLambda[++k^2]/k != n]; k];
    Array[A341860, 100] (* Paolo Xausa, Mar 11 2024 *)
  • PARI
    a(n) = for(k=1, oo, if(A341857(k)==n, return(k))) \\ Suppose that a(n) > 0. See A341857 for its program