%I #16 Sep 10 2020 09:28:24
%S 1,2,5,8,9,10,13,24,21,18,21,40,25,26,45,64,33,42,37,72,65,42,45,120,
%T 65,50,81,104,57,90,61,160,105,66,117,168,73,74,125,216,81,130,85,168,
%U 189,90,93,320,133,130,165,200,105,162,189,312,185,114,117,360,121,122,273
%N Number of solutions to x^2 == y^2 (mod n).
%H Amiram Eldar, <a href="/A062803/b062803.txt">Table of n, a(n) for n = 1..10000</a>
%H László Tóth, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL17/Toth/toth12.html">Counting Solutions of Quadratic Congruences in Several Variables Revisited</a>, J. Int. Seq. 17 (2014) # 14.11.6.
%F a(n) is multiplicative and, for an odd prime p, a(p) = 2*p - 1.
%F Multiplicative with a(2^e)=e*2^e and a(p^e)=((p-1)*e+p)*p^(e-1) for an odd prime p. - _Vladeta Jovovic_, Sep 22 2003
%t f[2, e_] := e*2^e; f[p_, e_] := ((p-1)*e+p)*p^(e-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Sep 10 2020 *)
%Y Cf. A086933.
%K nonn,mult,easy
%O 1,2
%A Ahmed Fares (ahmedfares(AT)my-deja.com), Jul 19 2001
%E More terms from _Vladeta Jovovic_, Sep 22 2003
|