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.

A372292 Numbers that occur more than once in the odd bisection of A371094.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Apr 26 2024

Keywords

Comments

Numbers that occur more than once in array A371100.

Examples

			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.
		

Crossrefs

Setwise difference A372290 \ A372291.
Cf. A144864 (subsequence after its initial 1), A371094, A371100.

Programs

  • 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];