cp's OEIS Frontend

A272570 Inverse permutation to A175500.

%I A272570 #15 May 23 2016 17:10:39
%S A272570 1,2,3,4,5,6,7,9,10,11,14,8,16,12,18,13,20,19,22,24,28,33,35,15,25,37,
%T A272570 39,26,41,27,43,30,45,47,49,17,51,53,62,29,55,57,64,31,66,84,82,21,58,
%U A272570 70,88,74,86,59,90,60,94,96,92,23,98,102,76,32,112,68
%N A272570 Inverse permutation to A175500.
%H A272570 Max Barrentine, <a href="/A272570/b272570.txt">Table of n, a(n) for n = 1..693</a>
%t A272570 p = Last /@ Import["https://oeis.org/A175500/b175500.txt", "Table"]; n = First[ Complement[ Range[10^4], p]] - 1; ip = 0 Range@ n; Do[If[(v = p[[i]]) <= n, ip[[v]] = i], {i, Length@p}]; ip (* _Giovanni Resta_, May 21 2016 *)
%o A272570 (PARI) ok(j, va, vs, n) = {if (vecsearch(vs, j), return (0)); for (k=1, n-1, if ((numdiv(j) == numdiv(va[k])) && (numdiv(va[k-1]) == numdiv(va[n-1])), return (0));); 1;}
%o A272570 findnew(va, vs, n) = {my(j = 1); my(vs = vecsort(va)); until (ok(j, va, vs, n), j++); j;}
%o A272570 listb(nn) = {my(va = [1]); for (n=2, nn, vs = vecsort(va); newa = findnew(va, vs, n); va = concat(va, newa);); va;}
%o A272570 invp(v) = {for (i=1, vecmax(v), found = 0; for (k=1, #v, if (v[k] == i, found = k; break);); if (! found, break); print1(found, ", "););}
%o A272570 lista(nn) = invp(listb(nn)); \\ _Michel Marcus_, May 04 2016
%Y A272570 Cf. A175500.
%K A272570 nonn
%O A272570 1,2
%A A272570 _Max Barrentine_, May 03 2016