如果可以当然要走直线
如果不行呢
那就是几段相连的折线
折线。为了叙述方便,考虑它跟
y
=
m
x
n
y=\frac{mx}{n}
y=nmx 围成的多边形
P(n,m) O(0,0)
先不考虑会不会交到格点,那么你会发现取多边形必然没有取三角形好
但是实际上三角形OAP可能会经过格点……比如说B。
然而你会发现这样的话取OBP比取OAP更好。
假如再经过格点,那就继续取,显然能够取完?
直觉上是这样子的。
到这里,我们知道最优的三角形边上不会有格点。
这还不够,继续考虑一下:三角形内部又怎样呢?
不妨再考虑三角形内部的格点
那么你会发现,假如三角形内部有格点
那当然是取里面的好啊
所以最优的三角形内部不会有格点
那么实际上
我们对
∀
x
∈
{
0
,
1
,
⋯
,
n
}
\forall x\in\{0,1,\cdots,n\}
∀x∈{0,1,⋯,n}
取
y
1
=
⌈
m
x
n
⌉
,
y
2
=
⌊
m
x
n
⌋
y_1=\lceil \frac{mx}{n}\rceil,\quad y_2=\lfloor\frac{mx}{n}\rfloor
y1=⌈nmx⌉,y2=⌊nmx⌋
枚举所有的
(
x
,
y
1
)
(x,y_1)
(x,y1)
(
x
,
y
2
)
(x,y_2)
(x,y2)
找到最优的合法答案即可
以上是
O
(
n
)
O(n)
O(n)
稍微优化一下,即可得到
O
(
min
{
n
,
m
}
)
O(\min\{n,m\})
O(min{n,m}) 的做法
还可以优化
合法三角形边界上的格点数目为 s=gcd(n,m)+2
由Pick’s Theorem
最优三角形面积 S=0+(s/2)-1=gcd(n,m)/2
那么,三角形顶点A在直线
y
=
m
x
±
(
m
,
n
)
n
y=\frac{mx\pm(m,n)}{n}
y=nmx±(m,n)上
不妨
n
=
a
(
m
,
n
)
,
m
=
b
(
m
,
n
)
n=a(m,n),\quad m=b(m,n)
n=a(m,n),m=b(m,n)
反正不需要求出所有的 A,那直接取
±
\pm
± 为
+
+
+
则直线解析式
y
=
b
x
+
1
a
y=\frac{bx+1}{a}
y=abx+1
就有
a
y
−
b
x
=
1
=
(
a
,
b
)
ay-bx=1=(a,b)
ay−bx=1=(a,b)
exgcd
然后用中垂线剖开,二分
当然直接三分也行吧
O
(
log
n
)
O(\log n)
O(logn)