%I #10 Sep 23 2022 11:52:40
%S 0,1,7,19,48,112,212,383,641,988,1523,2250,3103,4364,5831,7696,9985,
%T 12945,16164,20246,24946,30145,36385,43839,51752,61610,72475,84273,
%U 97231,112733
%N a(n) is the number of distinct radii of circles passing through at least three points in a square grid of n X n points.
%C Proposed by _Ed Pegg Jr_ in Mathematics StackExchange question 273348, where terms up to 212 were given, see link.
%H Mathematics StackExchange, <a href="https://math.stackexchange.com/questions/273348/number-of-circles-in-configuration">Number of circles in configuration</a>, comment by Ed Pegg, Jan 09, 2013.
%e a(2) = 1: only one possible circle with squared radius 1/2;
%e a(3) = 7: squared radii of possible circles are {1, 2, 1/2, 5/2, 5/4, 25/16, 25/18}.
%o (PARI) \\ Function r2 determined using Mathematica
%o \\ (ArcLength [CircleThrough[{{x1, y1}, {x2, y2}, {x3, y3}}]]/(2*Pi))^2
%o radii(n) = {my (R=List(), r2(x1,y1,x2,y2,x3,y3) = ((x3 + (x2^2*y1 - x3^2*y1 - x1^2*y2 + x3^2*y2 - y1^2*y2 + y1*y2^2 + x1^2*y3 - x2^2*y3 + y1^2*y3 - y2^2*y3 - y1*y3^2 + y2*y3^2)/(2*(-x2*y1 + x3*y1 + x1*y2 - x3*y2 - x1*y3 + x2*y3)))^2 + (y3 - (-2*(x2 - x3)*(x1^2 - x3^2 + y1^2 - y3^2) + 2*(x1 - x3)*(x2^2 - x3^2 + y2^2 - y3^2))/(-4*x2*y1 + 4*x3*y1 + 4*x1*y2 - 4*x3*y2 - 4*x1*y3 + 4*x2*y3))^2)); for(x1=1, n, for(y1=1, n, for(x2=1, x1, for(y2=1, n, for(x3=1, x2, for(y3=1, n,
%o my (ax2=2*(x2-x1), ay2=2*(y2-y1), ax3=2*(x3-x1), ay3=2*(y3-y1), den=ax2*ay3 -ax3*ay2); if (den==0, next); listput (R, r2(x1,y1,x2,y2,x3,y3)))))))); Set(R)};
%o for (k=1,15, print1(#radii(k),", "))
%Y Cf. A192493, A192494, A355884.
%K nonn,more
%O 1,3
%A _Hugo Pfoertner_, Sep 23 2022, following a suggestion by _Ed Pegg Jr_, Jan 09 2013
|