デンソークリエイトプログラミングコンテスト2022 Winter(AtCoder Beginner Contest 280) G問題メモ

デンソークリエイトプログラミングコンテスト2022 Winter(AtCoder Beginner Contest 280)

G - Do Use Hexagon Grid 2

問題

無限に広がるHEXグリッドがある（問題ページ参照）
辺で隣接するマス間を移動できる
マス $X→Y$ へ移動するのに必要な最小移動回数を、距離とする
$N$ 個のマス $(X_1,Y_1),...,(X_N,Y_N)$ が与えられる
この中からマスを1つ以上選ぶ方法は全部で $2^N-1$ 個あるが、そのうち「選んだマスのうちどの2マスも距離が $D$ 以下」という条件を満たす方法の個数を $\mod{998244353}$ で求めよ
$1 \le N \le 300$

解法

本問題のような座標の決め方をしたHEXグリッドにおいて、$(0,0)$ と $(x,y)$ との距離は $\max(|x|,|y|,|x-y|)$ で求めることができるらしい。なるほど。

まずはこれを知ってるか気付くかしないと難しい。

で、各マスを改めて $(x,y,x-y)$ の3次元空間上の点として置きなおした時、ある $D \times D \times D$ の立方体の中に入る点同士は、どのように選んでも距離が $D$ 以下となる、ということができる。

これなら、立方体の中の点の個数を数えて、それぞれを選ぶ/選ばないを2の何乗かで計算できるので、かなり考えやすくなった。

組み合わせの数え方

とはいえ、じゃあどのような立方体を考えればいいのか？

中途半端な位置にある立方体を考える必要は無いので、3軸方向のそれぞれ一方の面には点があると考えてよい。
このように、端に来る点を固定するのはよくある考え方である。
ここでは、各軸方向で小さい座標の方の面に点がある立方体を考える。

ただ、固定する3点を全探索し、さらにそれぞれで立方体の中に入る点を数えると、$O(N^4)$ となり、TLEする。
固定するのはまず2点だけでよく、最後に3点目を固定すると $O(N^3)$ や $O(N^3 \log{N})$ となる。

まず、x座標最小の点 $P=(X_P,Y_P)$ と、y座標最小の点 $Q=(X_Q,Y_Q)$ を固定する。
$P$ と $Q$ が同じ点である場合もあり得る。

$P,Q$ は必ず選ぶようにする前提なので、以下の場合は処理しない。

$P$ と $Q$ の距離が $D$ を超える
$P=Q$ である場合を除き、以下のいずれか
- $X_Q \le X_P$（$P$ がx座標最小という前提に違反）
- $Y_P \le Y_Q$（$Q$ がy座標最小という前提に違反）
- 不等号にイコールを含む理由については、「重複の除き方」参照

その上で、x座標が $[X_P,X_P+D]$、y座標が $[Y_Q,Y_Q+D]$ の正方形に含まれ、かつ $P$ とも $Q$ とも距離が $D$ 以下である点をリストアップする。

D ┼・─────┼
  │    ・      │
  Ｐ  ・    ・  │
  │            ・
  │      ・    │
  ┼─Ｑ────┼
                 D

リストアップされた点は、個別には $P,Q$ とともに選ぶことができるが、しかし、残る1軸（x-y軸）によっては、リストアップされた点同士は一緒に選べないかも知れない。

ただ、ここまで来れば残りは1次元なので、$x-y$ でソートし、ここで3軸目の最小座標の点を固定してやれば、どこまで選ぶことができるかは尺取りや二分探索等で簡単にわかる。

固定した3点目を $R$、その $x-y$ の値を $Z_R$ とする。

リストアップされた点のx-yのソート結果例

          R
Zi: [ -3 -1 -1  1  2  4 ]

$[Z_R,Z_R+D]$ 内に、$R$ を除いてたとえば4点が含まれるとすると、 $\{P,Q,R\}$ を全て含む選び方に対して、その4点は自由に選べる。よってこの3点を固定した場合は $2^4$ 通りとなる。

これを、$R$ を全通り（上記の例なら6点）ずらしながら、結果を合計する。

最後に、$\{P,Q\}$ の2点（$P=Q$ の場合は1点）だけしか選ばれないという1通りを入れれば、 $P,Q$ を固定した場合の数え上げが完了する。

重複の防ぎ方

x座標やy座標が同じ点が複数あるとき、特に重複がややこしくなる。

主に2つの状況で、重複を防ぐよう注意が必要となる。

x座標最小とy座標最小の2点を固定するとき
2点を固定し、正方形に含まれる点をリストアップするとき

2点を固定するとき

x座標最小 $P$ とy座標最小 $Q$ は、同じ点である場合もある。

y │・
  │    ・
  ②
  │  ・
─①─③── x
  │

「x座標最小を①、y座標最小も①」と固定した場合の処理では、 ②や③も正方形に含まれる点としてリストアップされ、{①,②} や {①,②,③} などが数え上げられる。

ここで「x座標最小を②、y座標最小を①」と固定した場合も処理してしまうと、{①,②} などが重複して数え上げられてしまう。
「x座標最小を①、y座標最小を③」とした場合も同様。

よって、「2点が選ばれたとき、その座標のいずれかが同じ2点」である場合は、既にx軸最小・y座標最小とも同じ点とした場合で数えられているため、処理をしてはいけない。

正方形に含まれる点をリストアップするとき

y │
  ③
・│          ・
  │    ・
  ②
  │                ・
  ①
  ┼─④─⑤───⑥─ x
  ・

たとえば上記の例で、$(①,②,④,⑤)$ は、

①をx座標最小、④をy座標最小と固定して、含まれる点として②,⑤が選ばれた
②をx座標最小、⑤をy座標最小と固定して、含まれる点として①,④が選ばれた

など、計4通りで重複して数えてしまう。

これを防ぐ方法としては、正方形に含まれる点をリストアップするとき、

①をx座標最小と固定したら、②も③も含む
②をx座標最小と固定したら、x座標が同じでそれよりy座標が小さい①は含めず、③のみ含む
③をx座標最小と固定したら、①も②も含まない

というように、同じのx座標の点の中では、固定した点がy座標でも最小であるように選べばよい。

とにかく、様々な局面で「最小である1点を固定」という考え方が重宝する問題。

Python3

n, d = map(int, input().split())
points = [tuple(map(int, input().split())) for _ in range(n)]
MOD = 998244353


def honeycomb_dist(x1, y1, x2, y2):
    dx = x2 - x1
    dy = y2 - y1
    return max(abs(dx), abs(dy), abs(dx - dy))


points.sort()

ans = 0
for i in range(n):
    x1, y1 = points[i]
    xy1 = x1 - y1

    for j in range(i, n):
        x2, y2 = points[j]
        xy2 = x2 - y2

        if i != j and y1 <= y2:
            continue
        if honeycomb_dist(x1, y1, x2, y2) > d:
            continue

        ans += 1

        cur_points = []
        for k in range(i + 1, n):
            if j == k:
                continue
            x3, y3 = points[k]
            xy3 = x3 - y3
            if (x1 <= x3 <= x1 + d
                    and ((j > k and y2 < y3 <= y2 + d) or (j < k and y2 <= y3 <= y2 + d))
                    and abs(xy1 - xy3) <= d
                    and abs(xy2 - xy3) <= d):
                cur_points.append(xy3)

        if len(cur_points) == 0:
            continue

        cur_points.sort()
        m = len(cur_points)
        l = 0
        r = 0
        while l < m:
            while r < m and cur_points[l] + d >= cur_points[r]:
                r += 1
            ans += pow(2, r - l - 1, MOD)
            l += 1
        ans %= MOD

print(ans % MOD)

目次