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 A294169 #38 Apr 10 2019 10:21:25
%S A294169 65,133,301,793,2041,2413,2501,2701,3781,4699,5617,5963,7081,7991,
%T A294169 9073,9881,9937,10261,10349,12209,13213,13333,14111,14981,18721,20737,
%U A294169 24727,27133,31201,31621,35431,40321,47197,49141,49591,49601,54913,60701,64079,65869,67721,70801
%N A294169 Semiprimes k = pq such that p^k == p (mod k) and q^k == q (mod k).
%C A294169 The number k = pq is a weak pseudoprime to prime bases p and q.
%C A294169 Problem: are there infinitely many such numbers?
%C A294169 All the terms are odd squarefree semiprimes.
%C A294169 Semiprimes pq such that p^(p-1) == 1 (mod q) and q^(q-1) == 1 (mod p).
%C A294169 Odd semiprimes pq such that (q-p)^(q-p) == 1 (mod pq).
%C A294169 Semiprimes pq > 6 such that (q-p)^(q-p) == 1 (mod pq).
%C A294169 Odd semiprimes pq pseudoprime to base q-p.
%H A294169 Robert G. Wilson v, <a href="/A294169/b294169.txt">Table of n, a(n) for n = 1..2812</a>
%e A294169 65 = 5*13 is a term since 5^65 == 5 (mod 65) and 13^65 == 13 (mod 65).
%e A294169 Equivalently: 5^(5-1) == 1 (mod 13) and 13^(13-1) == 1 (mod 5).
%e A294169 Also (13-5)^(5*13-1) == 1 (mod 5*13) or (13-5)^(13-5) == 1 (mod 5*13).
%t A294169 k = 4; lst = {}; NextSemiPrime[n_, k_: 1] := Block[{c = 0, sgn = Sign[k], sp}, sp = n + sgn; While[c < Abs[k], While[ PrimeOmega@ sp != 2, If[ sgn < 0, sp--, sp++]]; If[sgn < 0, sp--, sp++]; c++]; sp + If[sgn < 0, 1, -1]]; fQ[k_] := Block[{fi = First@# & /@ FactorInteger@ k}, PowerMod[#, k, k] & /@ fi == fi]; While[k < 100000, If[ fQ@ k, AppendTo[lst, k]]; k = NextSemiPrime@ k] (* _Robert G. Wilson v_, Feb 10 2018 *)
%o A294169 (PARI) lista(nn) = {for (n=1, nn, if (bigomega(n) == 2, if (omega(n) == 2, p = factor(n)[1,1]; q = factor(n)[2,1];, p = factor(n)[1,1]; q = factor(n)[1,1];); mp = Mod(p, n); mq = Mod(q, n); if ((mp^n == mp) && (mq^n == mq), print1(n, ", "));););} \\ _Michel Marcus_, Feb 10 2018
%Y A294169 Cf. A001358.
%Y A294169 Subsequence of A046388.
%K A294169 nonn
%O A294169 1,1
%A A294169 _Thomas Ordowski_, Feb 10 2018
%E A294169 More terms from _Michel Marcus_, Feb 10 2018
%E A294169 Edited by _Thomas Ordowski_, Mar 12 2019