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.
%I A306451 #16 Mar 16 2019 13:17:29 %S A306451 3,6,66,561,726,7107,8205,8646,62745,100101,140097,166521,237381, %T A306451 237945,566805,656601,876129,1053426,1095186,1194285,1234806,1590513, %U A306451 1598871,1938021,2381259,2518041,3426081,4125441,5398401,5454681,5489121,5720331,5961441 %N A306451 Non-coprime pseudoprimes or primes to base 3: numbers k that are multiples of 3 and are such that k divides 3^k - 3. %C A306451 Union of {3} and (A122780 - {1} - A005935). %C A306451 Numbers of the form 3*m such that 3^(3*m-1) == 1 (mod m). %C A306451 The squarefree terms are listed in A306450. %H A306451 Jianing Song, <a href="/A306451/b306451.txt">Table of n, a(n) for n = 1..163</a> (all terms below 10^9) %F A306451 66 is a term because 66 divides 3^66 - 3 = 3*(3^65 - 1) = 3*(3^5 - 1)*(3^60 + 3^55 + ... + 3^5 + 1) and 66 is divisible by 3. %o A306451 (PARI) forstep(n=3, 1e7, 3, if(Mod(3, n)^n==3, print1(n, ", "))) %Y A306451 Cf. A005935, A006935, A122780, A306450. %Y A306451 A258801 is a subsequence. %K A306451 nonn %O A306451 1,1 %A A306451 _Jianing Song_, Feb 17 2019