◆19999から20139の範囲に
素数は15個

20011 20021 20023 20029 20047 20051
20063 20071 20089 20101 20107 20113
20117 20123 20129

◆superPCM関数
Table[Product[(2n-1)^(C(0,3-a))C(0,C(0,((n-a)^(2a-2)mod(2a-1)))),{a,3,100}],{n,10000,10070}]

{0, 0, 0, 0, 0, 0, 20011, 0, 0, 0, 0, 20021,
20023, 0, 0, 20029, 0, 0, 0, 0, 0, 0, 0, 0,
20047, 0, 20051, 0, 0, 0, 0, 0, 20063, 0, 0,
0, 20071, 0, 0, 0, 0, 0, 0, 0, 0, 20089, 0, 0,
0, 0, 0, 20101, 0, 0, 20107, 0, 0, 20113, 0,
20117, 0, 0, 20123, 0, 0, 20129, 0, 0, 0, 0, 0}


◆的中率100％