OFFSET
1,16
MAPLE
with(numtheory);
f:=proc(n) local m, i, j, c; global pi;
if n < 4 then return(0); fi;
c:=0; m:=NumberTheory[pi](n); if isprime(n) then m:=m-1; fi;
for i from 1 to m-1 do p:=ithprime(i);
for j from i+1 to m do
if ithprime(j)*p mod n =1 then c:=c+1; fi; od: od: c; end;
[seq(f(n), n=1..120)];
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 02 2018
STATUS
approved
