# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a075764 Showing 1-1 of 1 %I A075764 #23 Jun 28 2022 02:57:23 %S A075764 105,261,301,693,1605,1755,2151,2905,2907,3393,3875,4641,4833,5005, %T A075764 5655,6279,6913,7161,8883,9405,10899,11025,11289,15687,17199,19203, %U A075764 22275,27387,36855,37791,50007,50463,53493,54891,55209,55755,63327,64337 %N A075764 Schroeder pseudoprimes: Composites k that divide the k-th Schroeder number A001003(k-1). %H A075764 Amiram Eldar, Table of n, a(n) for n = 1..422 %e A075764 105 is a term because A001003(105) = 15646506064359350392347086201481965698808674470977505246623827696393838448345 which is divisible by 105. %e A075764 105 is a term because A001003(104) = 15646506064359350392347086201481965698808674470977505246623827696393838448345 which is divisible by 105. %t A075764 s = {}; k1 = k2 = 1; Do[k3 = ((6*n - 9)*k2 - (n - 3)*k1)/n; If[CompositeQ[n] && Divisible[k3, n], AppendTo[s, n]]; k1 = k2; k2 = k3, {n, 3, 10^5}]; s (* _Amiram Eldar_, Jun 28 2022 *) %o A075764 (PARI) x1 = 1; x2 = 1; for (n = 3, 100000, x = (3*(2*n - 3)*x1 - (n - 3)*x2)/n; if (!isprime(n), if (!(x%n), print(n))); x2 = x1; x1 = x); \\ _David Wasserman_, Feb 23 2005 %Y A075764 Intersection of A002808 and A075763. %Y A075764 Cf. A001003, A013998, A075762. %K A075764 nonn %O A075764 1,1 %A A075764 _Benoit Cloitre_, Oct 09 2002 %E A075764 More terms from _David Wasserman_, Feb 23 2005 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE