0

格子点の個数を面積で評価する

112
1

はじめに

与えられた領域内の格子点の個数を、面積を使って大まかに評価したいときってありますよね。特に高校数学では積分の応用例とかでよく登場しますね。個人的には整数論でMinkowski空間とかを考える際にこういうことをしたくなります。格子多角形に限定するとピックの定理が有名ですが、一般の領域でこの状況を包括的に扱える道具がありそうで見つからなかったので、作ってみることにしました。これを使うと、多くの図形で面積Sと格子点の個数Nの差を、周長Lを使って
|SN|=O(L)(L)
という形で評価できます。(ランダウの記号)

準備

区分的になめらかな曲線C:[a,b](x(t),y(t))に対して、そのマンハッタン長Lを、
L=ab|x(t)|dt+ab|y(t)|dt
で定める。これはL1ノルム空間(R2,1)における曲線の長さである。

馴染みない概念かもしれません。実際ふつうの弧長を使っても良かったのですが、マンハッタン長Lは通常の長さよりも簡単に計算できるので使いやすい方を選びました。たとえば、楕円x2a2+y2b2=1の長さは、
4a0π/21e2sin2θdθ,e=1b2a2
と複雑ですが、マンハッタン長は
4a+4b
となります。これは次の命題からわかります。証明は簡単なので省略します。

曲線Cを有限個の曲線に分割してC=C1+C2++Cnとし、各Ciについて条件


() Ci上でx(t),y(t)は広義単調[1]である。

を満たすとする。Ciの始点を(xi,yi)、終点を(xi+1,yi+1)とすると、Cのマンハッタン長L
L=i=1n{|xi+1xi|+|yi+1yi|}
である。

曲線Cに対し、x(t)またはy(t)が極値をとるtの個数をnとする。これを単にCの極値の個数という。

極値とか言ってますが、命題1の条件()を満たす分割がとれるようなnであればなんでもよいです。

定理

定理の主張

今回証明したいのが次の定理です。

DR2が面積Sをもつ有界領域で、その境界C:=Dが区分的になめらかな閉曲線であるとする。Cのマンハッタン長Lと極値の個数nが有限のとき、内部の格子点の数Nについて、
|SN|<L+2n
が成り立つ。

証明は後回しにして系を紹介しましょう。少し弱い評価になってしまいますが、通常の意味での弧長で評価することもできます。

DR2が面積Sをもつ有界領域で、その境界C:=Dが区分的になめらかな閉曲線であるとする。Cの弧長Lと極値の個数nが有限のとき、内部の格子点の数Nについて、
|SN|<2L+2n

L=ab|x(t)|+|y(t)|dtab2x(t)2+y(t)2dt=2L
であるから従う。

次はよくある設定ですね。

DR2が面積Sをもつ有界領域で、その境界C:=Dが区分的になめらかな閉曲線であるとする。定数r>0に対し、
rD:={(rx,ry)(x,y)D}
とする。rDに含まれる格子点の数をNrとすると、
limrNrr2=S
が成り立つ。

また、凸領域については次の評価が得られます。これについては単純に積分とかすれば証明できそうな気もします...

DR2が面積Sをもつ有界凸領域で、その境界C:=Dが区分的になめらかな曲線であるとする。
L=max(x,y)Cxmin(x,y)Cx+max(x,y)Cymin(x,y)Cy
とする。(L=横幅+縦幅)このとき、内部の格子点の数Nについて、
|SN|<2L+8

証明

証明の流れ

まずは証明の流れを説明します。

Step 1

各格子点に対して、それを中心とする一辺1の正方形のタイルを考えると、これらで平面全体を埋め尽くすことができます。そして、これらのタイルは次の3つに分けることができますね。

  1. 完全にDの内側にあるもの
  2. 部分的にDに入っているもの (すなわち曲線Cと交わるもの)
  3. 完全にDの外側にあるもの

(1),(2)の正方形の個数をN1,N2とすれば、面積はN1SN1+N2を満たします。このことから、D内にある格子点の数Nについて|SN|N2という評価ができます。

Step 2

(2)のタイルの個数N2をマンハッタン長Lを使って評価します。長さL<の曲線が踏めるタイルの数はもちろん有限なので上から評価できるはずです。ここが一番大変で、改善の余地があるところです。

Step 1. |SN|を曲線に近い格子点の個数|V(C)|で評価

各格子点pZ2に対して、pを中心とする一辺1の正方形の内部をB(p)とし、境界を含めたものをB(p)とする。また、曲線Cに対し格子点の集合V(C)を、
V(C):={pZ2B(p)C}
で定める。

C=Dとする。格子点pZ2pV(C)をみたすとき、次のいずれか一方が成り立つ。

  1. B(p)Dである。
  2. B(p)Dcである。

Dの境界Cと交わらない正方形は、Dの完全に内側にあるか外側にあるかだよ」と言っています。そりゃそうですね。


証明pV(C)すなわちB(p)D=と仮定すると、開集合によるDの分割
B(p)=(B(p)D)(B(p)(Dc))
が得られる。B(p)は連結だから、B(p)DまたはB(p)(Dc)である。後者の場合両辺の閉包を取ればB(p)Dcが得られる。

定理2の条件の下で、
|SN||V(C)|

V(C)を2つの集合X,Yに分割する。
X=V(C)D,Y=V(C)Dc
補題3より次が成り立つ。

  1. p(DZ2)Xに対して、B(p)Dである。
  2. p(DcZ2)Yに対して、B(p)Dcである。

したがって、R2=pZ2B(p)に注意すれば、
p(DZ2)XB(p)Dp(DZ2)YB(p)
であるから、次の不等式が成り立つ。
N|X|SN+|Y|
よって、
|SN||X|+|Y|=|V(C)|
である。

Step 2. |V(C)|をマンハッタン長Lで評価

Z2上の半順序を、
(x1,y1)(x2,y2):⟺x1x2y1y2
で定める。p,qZ2に対し、p<q:⟺pqpqとする。

曲線C:[a,b]c(t)=(x(t),y(t))に対し、x(t),y(t)は広義単調増加であるとする。V(C)Z2の全順序部分集合である。すなわち、
p1,p2V(C),p1p2p2p1
が成り立つ。

p1,p2V(C)をとる。c(t1)=p1,c(t2)=p2となるt1,t2[a,b]をとる。t1t2として一般性を失わない。p1=(x1,y1),p2=(x2,y2)とする。このとき、
x112<x(t1)<x1+12,x212<x(t2)<x2+12
であるから、
x112<x(t1)x(t2)<x2+12x11<x2x1x2
である。同様にして、y1y2も示される。

定理2

命題4より、
|SN||V(C)|
なので、|V(C)|Lを示せば良い。a=t0<t1<<tn=bを命題1の条件()を満たす分割とする。Ciを曲線C[ti,ti+1)上の制限とする。V(C)=i=1nV(Ci)なので、
|V(C)|i=1n|V(Ci)|
である。|V(Ci)|を上から評価しよう。Ciのパラメータ表示をx(t),y(t)とする。x(t),y(t)は広義単調増加であるとして一般性を失わない。このとき、V(Ci)Z2の全順序部分集合であるから、
V(Ci)={p1,p2,,pk},p1<p2<<pk
とおける。さて、pj=(xj,yj)に対しnj=xj+yjとおくと、n1<n2<<nkなので、
nkn1=j=1k1(nj+1nj)k1
である。また、不等式
x112x(ti)<x1+12,y112y(ti)<y1+12xk12<x(ti+1)xk+12,yk12<y(ti+1)yk+12(1)
に注意すれば、
n1=x1+y1>x(ti)+y(ti)1nk=xk+yk<x(ti+1)+y(ti+1)+1(2)
であるから、
nkn1<x(ti+1)+y(ti+1)x(ti)y(ti)+2
である。したがって、
|V(Ci)|=k<x(ti+1)+y(ti+1)x(ti)y(ti)+3
である。これより、
|V(C)|i=1n|V(Ci)|<L+3n
という評価が得られる。

もう少し強い評価を得るために、共通部分V(Ci)V(Ci+1)について考えよう。もし格子点pZ2があって、(x(ti+1),y(ti+1))B(p)であれば、pV(Ci)V(Ci+1)である。よって、
μi:={1(pZ2,x(ti+1),y(ti+1)B(p))0(otherwise)
とおけば、
|V(C)|i=1n(|V(Ci)|μi)
となる。μi=0のとき、x(ti+1),y(ti+1)のいずれかは半整数である。x(ti+1)が半整数なら、式(1)の不等式xk12<x(ti+1)xk+12x(ti+1)と修正される。したがって式(2)は
nk=xk+yk<x(ti+1)+y(ti+1)
と修正できる。したがって、
|V(C)|i=1n(|V(Ci)|μi)<L+2n
が成り立つ。

予想

最後に予想を紹介します。

DR2が面積Sの領域で、その境界C:=Dが閉曲線であるとする。Cのマンハッタン長Lが有限のとき、
D内部の格子点の数Nについて、
|SN|2L+4
が成り立つ。

nを求めるのは簡単っちゃ簡単ですが、やはりLだけで評価できたほうが便利だよねぇ...


[1]: 広義単調増加または広義単調減少


投稿日:18日前
更新日:18日前
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

dragoemon
dragoemon
143
31561
大学3年生です

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 準備
  3. 定理
  4. 定理の主張
  5. 証明
  6. 証明の流れ
  7. Step 1. |SN|を曲線に近い格子点の個数|V(C)|で評価
  8. Step 2. |V(C)|をマンハッタン長Lで評価
  9. 予想