login
A372292
Numbers that occur more than once in the odd bisection of A371094.
4
21, 117, 213, 309, 341, 405, 501, 597, 693, 789, 885, 981, 1077, 1173, 1269, 1365, 1461, 1557, 1653, 1749, 1845, 1877, 1941, 2037, 2133, 2229, 2325, 2421, 2517, 2613, 2709, 2805, 2901, 2997, 3093, 3189, 3285, 3381, 3413, 3477, 3573, 3669, 3765, 3861, 3957, 4053, 4149, 4245, 4341, 4437, 4533, 4629, 4725, 4821, 4917
OFFSET
1,1
COMMENTS
Numbers that occur more than once in array A371100.
LINKS
EXAMPLE
21 is present because A371094(1) = A371094(3) = 21.
87381 is present because A371094(85) = A371094(213) = A371094(7281) = A371094(14563) = 87381.
185685 is present because A371094(469) = A371094(15473) = A371094(30947) = 185685.
PROG
(PARI)
A371094(n) = { my(m=1+3*n, e=valuation(m, 2)); ((m*(2^e)) + (((4^e)-1)/3)); };
isA372292(n) = if(!(n%2), 0, my(c=0); forstep(k=1, n, 2, if(A371094(k)==n, c++)); (c>1));
(PARI)
search_up_to = 1398101;
A371094(n) = { my(m=1+3*n, e=valuation(m, 2)); ((m*(2^e)) + (((4^e)-1)/3)); };
A372292list(up_to_n) = { my(v=vector((1+up_to_n)/2), x, lista=List([])); forstep(k=1, up_to_n, 2, x=A371094(k); if(x <= up_to_n, v[(x+1)/2]++)); for(i=1, (1+up_to_n)/2, if(v[i]>1, listput(lista, i+i-1))); Vec(lista); };
v372292 = A372292list(search_up_to);
A372292(n) = v372292[n];
CROSSREFS
Setwise difference A372290 \ A372291.
Cf. A144864 (subsequence after its initial 1), A371094, A371100.
Sequence in context: A121628 A306532 A219154 * A054257 A341570 A297316
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 26 2024
STATUS
approved