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.

A291167 Numbers k such that psi(k) is a perfect square where psi(k) = A001615(k).

This page as a plain text file.
%I A291167 #16 Jul 23 2021 02:08:25
%S A291167 1,3,18,20,22,27,60,66,70,72,80,88,92,94,99,115,119,162,170,210,212,
%T A291167 214,217,240,243,252,264,265,276,280,282,288,308,310,315,320,322,345,
%U A291167 352,357,368,376,382,385,423,497,500,510,517,527,540,594,596,612,636,637,642,648,651,679,680,710,725,742
%N A291167 Numbers k such that psi(k) is a perfect square where psi(k) = A001615(k).
%C A291167 The product of an even number of distinct members of A066436 is in the sequence. - _Robert Israel_, Aug 22 2017
%H A291167 Robert Israel, <a href="/A291167/b291167.txt">Table of n, a(n) for n = 1..10000</a>
%e A291167 60 is a term because psi(60) = 144 is a perfect square.
%p A291167 filter:= proc(n) issqr(n*mul(1+1/p,p=numtheory:-factorset(n))) end proc:
%p A291167 select(filter, [$1..1000]); # _Robert Israel_, Aug 22 2017
%t A291167 Select[Range@ 750, IntegerQ@ Sqrt[# Sum[MoebiusMu[d]^2/d, {d, Divisors@ #}]] &] (* _Michael De Vlieger_, Aug 19 2017 *)
%o A291167 (PARI) a001615(n) = n*sumdivmult(n, d, issquarefree(d)/d);
%o A291167 is(n) = issquare(a001615(n));
%Y A291167 Cf. A001615, A006532, A066436.
%K A291167 nonn,easy
%O A291167 1,2
%A A291167 _Altug Alkan_, Aug 19 2017