题目
给定一个二维平面,平面上有 n 个点,求最多有多少个点在同一条直线上。
解题思路
来自这篇博客
思路很简单,看代码
给出2D平面中的n个坐标点,计算最多有多少个点在一条直线上
一条直线可以用斜率表示,即如果已知(x1,y1),(x2,y2)(x1,y1),(x2,y2),那么斜率为(y2−y1)/(x2−x1)(y2−y1)/(x2−x1) ,所以固定一个点,然后遍历其他点,将不同的斜率到个数的映射记录在map中,然后求最大的
思路是这样,但是由于浮点数计算有误差,导致斜率的计算不是准确的,所以不能直接记录斜率。因为斜率是由Δy/ΔxΔy/Δx计算得到的,那么可以采用<Δx,Δy><Δx,Δy>到个数的映射,不过需要将Δx,ΔyΔx,Δy化为最简的形式,这样,就可以将<Δx,Δy><Δx,Δy>这个数值对相等看成是斜率相等
化为最简的方法是找到两个数的最大公约数,使用辗转相除法,计算步骤是
用第一个数除以第二个数得到余数
如果余数为0,则第二个数就是最大公约数,返回
否则,用第二个数除以余数,回到步骤1(即将第二个数作为新的第一个数,将余数作为新的第二个数)
另外需要注意的一点是坐标点可能存在重合的情况
这里写成博客的原因是:
- 元组这种不可变的类型是可以当作字典的key
- 温习一下辗转相除法
# Definition for a point.
# class Point:
# def __init__(self, a=0, b=0):
# self.x = a
# self.y = b
def computeCom(m,n):
#辗转相除求最大公约数
if n==0:
return m
tmp=m
while True:
x=m%n
if x ==0:
return n
m=n
n=x
class Solution:
def maxPoints(self, points: List[Point]) -> int:
len_point=len(points)
if len_point<=1:
return len_point
max_v=-1
for i in range(len_point-1):
my_cal={}
samePo=1
x1=points[i].x
y1=points[i].y
for j in range(i+1,len_point):
xDiff=x1-points[j].x
yDiff=y1-points[j].y
if xDiff==0 and yDiff==0:
samePo+=1
continue
com=computeCom(xDiff,yDiff)
if (xDiff//com,yDiff//com) in my_cal:
my_cal[(xDiff//com,yDiff//com)]+=1
else:
my_cal[(xDiff//com,yDiff//com)]=1
max_s=0
for k,v in my_cal.items():
if v>max_s:
max_s=v
max_s+=samePo
if max_s>max_v:
max_v=max_s
return max_v