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.

A348018 a(n) is the index of A064549(n) = n * Product_{p prime|n} p in the sequence of powerful numbers (A001694).

This page as a plain text file.
%I A348018 #8 Sep 25 2021 04:45:10
%S A348018 1,2,4,3,6,9,10,5,7,14,16,12,20,21,24,8,28,15,31,22,35,36,39,19,17,44,
%T A348018 13,32,48,50,51,11,56,59,60,23,65,67,68,33,71,73,75,52,43,81,84,27,30,
%U A348018 37,90,64,94,29,97,46,102,104,107,74,110,111,62,18,117,119
%N A348018 a(n) is the index of A064549(n) = n * Product_{p prime|n} p in the sequence of powerful numbers (A001694).
%C A348018 A permutation of the positive integers.
%C A348018 The inverse permutation of A306458.
%H A348018 Amiram Eldar, <a href="/A348018/b348018.txt">Table of n, a(n) for n = 1..10000</a>
%F A348018 A001694(a(n)) = A064549(n).
%F A348018 A306458(a(n)) = a(A306458(n)) = n.
%F A348018 The fixed points of this permutation are 1, 2, 12, 1208, 1256, 1288 and no more below 3*10^5.
%e A348018 The sequence of powerful numbers (A001694) begins with 1, 4, 8, 9, ...
%e A348018 The position of A064549(1) = 1 in A001694 is 1, so a(1) = 1.
%e A348018 The position of A064549(3) = 9 in A001694 is 4, so a(3) = 4.
%t A348018 powQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], # > 1 &]; pows = Join[{1}, Select[Range[10^4], powQ]]; TakeWhile[Table[FirstPosition[pows, n * Times @@ (First /@ FactorInteger[n])][[1]], {n, 1, 100}], NumericQ]
%Y A348018 Cf. A001694, A064549, A306458.
%K A348018 nonn
%O A348018 1,2
%A A348018 _Amiram Eldar_, Sep 24 2021