A250469
a(1) = 1; and for n > 1, a(n) = A078898(n)-th number k for which A055396(k) = A055396(n)+1, where A055396(n) is the index of smallest prime dividing n.
Original entry on oeis.org
1, 3, 5, 9, 7, 15, 11, 21, 25, 27, 13, 33, 17, 39, 35, 45, 19, 51, 23, 57, 55, 63, 29, 69, 49, 75, 65, 81, 31, 87, 37, 93, 85, 99, 77, 105, 41, 111, 95, 117, 43, 123, 47, 129, 115, 135, 53, 141, 121, 147, 125, 153, 59, 159, 91, 165, 145, 171, 61, 177, 67, 183, 155, 189, 119, 195, 71, 201, 175, 207, 73, 213, 79, 219, 185, 225, 143, 231, 83, 237, 205, 243, 89, 249, 133, 255
Offset: 1
Cf.
A000040,
A003961,
A016945,
A046523,
A055396,
A078898,
A083140,
A083221,
A084967,
A249744,
A249810,
A249820,
A249817,
A249818,
A250471,
A266645,
A280692,
A280693,
A283465.
-
a[1] = 1; a[n_] := If[PrimeQ[n], NextPrime[n], m1 = p1 = FactorInteger[n][[ 1, 1]]; For[k1 = 1, m1 <= n, m1 += p1; If[m1 == n, Break[]]; If[ FactorInteger[m1][[1, 1]] == p1, k1++]]; m2 = p2 = NextPrime[p1]; For[k2 = 1, True, m2 += p2, If[FactorInteger[m2][[1, 1]] == p2, k2++]; If[k1+2 == k2, Return[m2]]]]; Array[a, 100] (* Jean-François Alcover, Mar 08 2016 *)
g[n_] := If[n == 1, 0, PrimePi@ FactorInteger[n][[1, 1]]]; Function[s, MapIndexed[Lookup[s, g[First@ #2] + 1][[#1]] - Boole[First@ #2 == 1] &, #] &@ Map[Position[Lookup[s, g@#], #][[1, 1]] &, Range@ 120]]@ PositionIndex@ Array[g, 10^4] (* Michael De Vlieger, Mar 08 2017, Version 10 *)
A250246
Permutation of natural numbers: a(1) = 1, a(n) = A246278(A055396(n), a(A078898(n))).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 27, 22, 23, 24, 25, 26, 21, 28, 29, 30, 31, 32, 45, 34, 35, 36, 37, 38, 33, 40, 41, 54, 43, 44, 81, 46, 47, 48, 49, 50, 75, 52, 53, 42, 125, 56, 63, 58, 59, 60, 61, 62, 39, 64, 55, 90, 67, 68, 135, 70, 71, 72, 73, 74, 51, 76, 77, 66, 79, 80, 99, 82, 83
Offset: 1
Cf.
A003961,
A005843,
A020639,
A055396,
A078898,
A246278,
A250470,
A268674,
A278524,
A302042,
A302046.
Differs from the "vanilla version"
A249818 for the first time at n=42, where a(42) = 54, while
A249818(42) = 42.
Differs from
A250250 for the first time at n=73, where a(73) = 73, while
A250250(73) = 103.
-
up_to = 16384;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
A020639(n) = if(n>1, if(n>n=factor(n, 0)[1, 1], n, factor(n)[1, 1]), 1); \\ From A020639
A055396(n) = if(1==n,0,primepi(A020639(n)));
v078898 = ordinal_transform(vector(up_to,n,A020639(n)));
A078898(n) = v078898[n];
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A250246(n) = if(1==n,n,my(k = 2*A250246(A078898(n)), r = A055396(n)); if(1==r, k, while(r>1, k = A003961(k); r--); (k))); \\ Antti Karttunen, Apr 01 2018
(Scheme, with memoizing-macro definec from Antti Karttunen's IntSeq-library, three alternative definitions)
(definec (A250246 n) (cond ((<= n 1) n) (else (A246278bi (A055396 n) (A250246 (A078898 n)))))) ;; Code for A246278bi given in A246278
(definec (A250246 n) (cond ((<= n 1) n) ((even? n) (* 2 (A250246 (/ n 2)))) (else (A003961 (A250246 (A250470 n))))))
(define (A250246 n) (A163511 (A252756 n)))
A249817
Permutation of natural numbers: a(1) = 1, a(n) = A083221(A055396(n),A246277(n)).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 27, 22, 23, 24, 25, 26, 21, 28, 29, 30, 31, 32, 39, 34, 35, 36, 37, 38, 63, 40, 41, 42, 43, 44, 33, 46, 47, 48, 49, 50, 75, 52, 53, 54, 65, 56, 99, 58, 59, 60, 61, 62, 57, 64, 95, 66, 67, 68, 111, 70, 71, 72, 73, 74, 51, 76, 77, 78, 79, 80, 45, 82, 83, 84, 155, 86, 135
Offset: 1
Other similar or related permutations:
A249815.
Cf.
A000040,
A005843,
A020639,
A055396,
A078898,
A083140,
A083221,
A246277,
A246278,
A246279,
A249821.
Differs from its inverse
A249818 for the first time at n=33, where a(33) = 39, while
A249818(33) = 45.
-
lim = 87; a083221 = Table[Take[Prime[n] Select[Range[Ceiling[lim/2]^2], GCD[# Prime@ n, Product[Prime@ i, {i, 1, n - 1}]] == 1 &], Ceiling[lim/2]], {n, Ceiling[lim/2]}]; a055396[n_] PrimePi[FactorInteger[n][[1, 1]]]; a246277[n_] := Which[n == 1, 0, EvenQ@ n, n/2, True, a246277[Times @@ Power[Which[# == 1, 1, # == 2, 1, True, NextPrime[#, -1]] & /@ First@ Transpose@ FactorInteger@ n, Last@ Transpose@ FactorInteger@ n]]]; Table[a083221[[a055396@ n, a246277@ n]], {n, 2, lim}] (* Michael De Vlieger, Jan 04 2016, after Jean-François Alcover at A055396 and Yasutoshi Kohmoto at A083140 *)
-
(define (A249817 n) (if (= 1 n) n (A083221bi (A055396 n) (A246277 n)))) ;; Code for A083221bi given in A083221
;; Alternative version:
(define (A249817 n) (if (= 1 n) n (A083221bi (A055396 n) (A249821bi (A055396 n) (A078898 n))))) ;; Code for A249821bi given in A249821.
A255407
Permutation of natural numbers: a(n) = A255127(A252460(n)).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 23, 20, 21, 22, 25, 24, 19, 26, 27, 28, 29, 30, 37, 32, 33, 34, 35, 36, 41, 38, 39, 40, 43, 42, 47, 44, 45, 46, 53, 48, 31, 50, 51, 52, 61, 54, 49, 56, 57, 58, 67, 60, 71, 62, 63, 64, 65, 66, 77, 68, 69, 70, 83, 72, 89, 74, 75, 76, 59, 78, 91, 80, 81
Offset: 1
A083221(8,1) = 19 and A255127(8,1) = 23, thus a(19) = 23.
A083221(9,1) = 23 and A255127(9,1) = 25, thus a(23) = 25.
A083221(3,2) = 25 and A255127(3,2) = 19, thus a(25) = 19.
A249822
Square array of permutations: A(row,col) = A078898(A246278(row,col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...
Original entry on oeis.org
1, 2, 1, 3, 2, 1, 4, 3, 2, 1, 5, 5, 3, 2, 1, 6, 4, 9, 3, 2, 1, 7, 8, 4, 14, 3, 2, 1, 8, 6, 12, 4, 28, 3, 2, 1, 9, 14, 5, 21, 4, 36, 3, 2, 1, 10, 13, 42, 5, 33, 4, 57, 3, 2, 1, 11, 11, 17, 92, 5, 45, 4, 67, 3, 2, 1, 12, 7, 19, 33, 305, 5, 63, 4, 93, 3, 2, 1, 13, 23, 6, 25, 39, 455, 5, 80, 4, 139, 3, 2, 1, 14, 9, 59, 6, 43, 61, 944, 5, 116, 4, 154, 3, 2, 1, 15, 17, 7, 144, 6, 52, 70, 1238, 5, 148, 4, 210, 3, 2, 1
Offset: 1
The top left corner of the array:
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, ...
1, 2, 3, 5, 4, 8, 6, 14, 13, 11, 7, 23, 9, 17, 18, ...
1, 2, 3, 9, 4, 12, 5, 42, 17, 19, 6, 59, 7, 22, 26, ...
1, 2, 3, 14, 4, 21, 5, 92, 33, 25, 6, 144, 7, 32, 39, ...
1, 2, 3, 28, 4, 33, 5, 305, 39, 43, 6, 360, 7, 48, 50, ...
1, 2, 3, 36, 4, 45, 5, 455, 61, 52, 6, 597, 7, 63, 68, ...
1, 2, 3, 57, 4, 63, 5, 944, 70, 76, 6, 1053, 7, 95, 84, ...
1, 2, 3, 67, 4, 80, 5, 1238, 96, 99, 6, 1502, 7, 106, 121, ...
...
Inverse permutations can be found from table
A249821.
Row k+1 is a right-to-left composition of the first k rows of
A251722.
Original entry on oeis.org
0, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 2, 0, -1, 0, 6, 0, 4, 0, 1, 0, -4, 0, 11, 0, -4, 4, 3, 0, 3, 0, 25, -1, -7, 0, 20, 0, -7, -1, 12, 0, 7, 0, -2, 4, -8, 0, 44, 0, 0, -2, 0, 0, 36, 0, 22, -2, -13, 0, 23, 0, -12, 8, 90, 0, 0, 0, -5, -2, 4, 0, 77, 0, -16, 4, -3, 0, 4, 0, 55, 28, -19, 0, 41, 0, -19, -4, 15, 0, 43, 0, -2, -3, -20, 0, 155, 0, 12, 5, 24, 0
Offset: 1
For n = 8 = 2*2*2, A003961(8) = 27 (3*3*3), and while 8 is on row 1 and column 4 of A083221, 27 on the next row is in column 5, thus a(8) = 5 - 4 = 1.
For n = 10 = 2*5, A003961(10) = 21 (3*7), and while 10 is on row 1 and column 5 of A083221, 21 on the next row is in column 4, thus a(10) = 4 - 5 = -1.
Cf.
A003961,
A078898,
A083221,
A083140,
A246277,
A249810,
A249817,
A249818,
A249821,
A249822,
A251721,
A251722.
A250250
Permutation of natural numbers: a(1) = 1, a(n) = A246278(a(A055396(n)),a(A078898(n))).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 27, 22, 23, 24, 25, 26, 21, 28, 29, 30, 31, 32, 45, 34, 35, 36, 37, 38, 33, 40, 41, 54, 43, 44, 81, 46, 47, 48, 49, 50, 75, 52, 53, 42, 125, 56, 63, 58, 59, 60, 61, 62, 39, 64, 55, 90, 67, 68, 135, 70, 71, 72, 103, 74, 51, 76, 77, 66, 79, 80, 99, 82, 83
Offset: 1
Differs from the "vanilla version"
A249818 for the first time at n=42, where a(42) = 54, while
A249818(42) = 42.
A266645
Permutation of natural numbers: a(n) = A064989(A250469(n)).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 10, 9, 8, 11, 14, 13, 22, 15, 12, 17, 26, 19, 34, 21, 20, 23, 38, 25, 18, 33, 16, 29, 46, 31, 58, 39, 28, 35, 30, 37, 62, 51, 44, 41, 74, 43, 82, 57, 24, 47, 86, 49, 50, 27, 52, 53, 94, 55, 42, 69, 68, 59, 106, 61, 118, 87, 40, 65, 66, 67, 122, 45, 76, 71, 134, 73, 142, 93, 36, 77, 70, 79, 146, 111, 32, 83, 158, 85, 78, 123
Offset: 1
-
f[n_] := Times @@ Power[Which[# == 1, 1, # == 2, 1, True, NextPrime[#, -1]] & /@ First@ #, Last@ #] &@ Transpose@ FactorInteger@ n; g[n_] := If[n == 1, 0, PrimePi@ FactorInteger[n][[1, 1]]]; Function[s, MapIndexed[ Function[{m, n}, f[Lookup[s, g[n] + 1][[m]] - Boole[n == 1]]][#1, First@ #2] &, #] &@ Map[Position[Lookup[s, g@ #], #][[1, 1]] &, Range@ 120]]@ PositionIndex@ Array[g, 10^4] (* Michael De Vlieger, Mar 09 2017, Version 10 *)
A255553
Permutation of natural numbers: a(n) = A255551(A252460(n)).
Original entry on oeis.org
1, 2, 3, 4, 7, 6, 9, 8, 5, 10, 13, 12, 15, 14, 11, 16, 21, 18, 25, 20, 17, 22, 31, 24, 19, 26, 23, 28, 33, 30, 37, 32, 29, 34, 39, 36, 43, 38, 35, 40, 49, 42, 51, 44, 41, 46, 63, 48, 27, 50, 47, 52, 67, 54, 61, 56, 53, 58, 69, 60, 73, 62, 59, 64, 81, 66, 75, 68, 65, 70, 79, 72, 87, 74, 71, 76, 57, 78, 93, 80, 77, 82, 99, 84, 103, 86, 83, 88, 105, 90
Offset: 1
A266646
Permutation of natural numbers: a(n) = A250470(A003961(n)).
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 10, 9, 8, 11, 16, 13, 12, 15, 28, 17, 26, 19, 22, 21, 14, 23, 46, 25, 18, 51, 34, 29, 36, 31, 82, 27, 20, 35, 76, 37, 24, 33, 64, 41, 56, 43, 40, 69, 30, 47, 136, 49, 50, 39, 52, 53, 126, 55, 100, 45, 32, 59, 106, 61, 38, 111, 244, 65, 66, 67, 58, 57, 78, 71, 226, 73, 42, 99, 70, 77, 86, 79, 190, 249, 44, 83, 166, 85
Offset: 1
Comments