login
A323249
Sequence lists numbers k > 1 such that k^2 == d(k) (mod sigma(k)), where d = A000005 and sigma = A000203.
3
8, 9, 14, 26, 38, 62, 74, 86, 122, 134, 146, 158, 194, 206, 218, 254, 278, 302, 314, 326, 362, 386, 398, 422, 446, 458, 482, 542, 554, 566, 614, 626, 662, 674, 698, 734, 746, 758, 794, 818, 842, 866, 878, 914, 926, 974, 998, 1046, 1082, 1094, 1142, 1154, 1202, 1214
OFFSET
1,1
COMMENTS
A112772 is a subsequence. After 8 and 9, the following numbers not in A112772 are 5438, 5462, 5504, 11096, ...
LINKS
FORMULA
Solutions of k^2 mod sigma(k) = d(k).
EXAMPLE
sigma(8) = 15 and 8^2 mod 15 = 4 = d(8).
MAPLE
with(numtheory): op(select(n->n^2 mod sigma(n)=tau(n), [$1..1214]));
MATHEMATICA
Select[Range[1225], PowerMod[#1, 2, #3] == #2 & @@ Prepend[DivisorSigma[{0, 1}, #], #] &] (* Michael De Vlieger, Jan 18 2019 *)
PROG
(PARI) for(k=1, 2000, x=sigma(k); if(Mod(k, x)^2==Mod(numdiv(k), x), print1(k, ", "))) \\ Jinyuan Wang, Feb 03 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Jan 08 2019
STATUS
approved