2022/05/02 記事中の記号の使い方やその他の軽微な修正をしました。
概要
のアフィン部分空間が の形で与えられるとき、 つのアフィン部分空間の共通部分を時間計算量 で計算する。
アフィン部分空間
ベクトル空間 の部分集合 は、部分空間 と によって と表現できるとき、アフィン部分空間であると言う。 アフィン部分空間は部分空間の一般化になっており、後述するが競プロでも出現する概念である。 に対して は一意だが は一般には一意でない。
アフィン部分空間の 通りの表現
を体として、 を のアフィン部分空間とする。 定義から明らかに、 によって と表現できる。 ここで は としたときの の次元である。
(ドット積についての) 直交補空間を考えると、 となる が存在する。 すると、 であり、諸々の記号を置きなおすと結局 という表現を得る*1。 逆に、 と表現される集合は、空であるかアフィン部分空間であることもすぐに確かめられる。
アフィン部分空間の共通部分の計算
をアフィン部分空間として、 と表現できるとする。 すると であり、容易に共通部分の表現が得られる。 この式から分かるように、これらのアフィン部分空間の共通部分は空であるかアフィン部分空間となる。 と の表現の相互の変換は、直交補空間の計算と同様に行うことができる (参考: Parity-check matrix - Wikipedia)。
実用例
注意: cp-unspoiler この問題のネタバレを含みます。
を、 を の元として見てこの順に連結したベクトルと定義する。 問題は以下のように整理される。
長さ の非負整数列 と、非負整数 が与えられる。 から つ以上選んで xor を計算し、前半 bit が 以下になるようにしつつ、後半 bit を最大化せよ。
以下の非負整数全体は、「下位 bit は自由で、残りは決まっている」という形の集合の高々 個の直和に分解される。 この分解のそれぞれについて独立に計算し、最後に最大値を取ればよい。
つ以上選ぶという条件を一旦 個以上に緩和すると、「 から 個以上選んで xor を取ることで得られる数全体」と「下位 bit は自由で、残りは決まっている数全体」はいずれも ( の) アフィン部分空間になっている。 したがって、これらの共通部分 を計算した後、 の元で後半 bit が最大となるものを計算すればよい。 実際には、 の定義を の順に連結することにして、単に を求めるのが簡単だろう。 かつ が線形独立である場合に限り最大値を に修正することで、 つ以上選ぶ条件をカバーする。
は、 の列ベクトルを掃き出して簡約化し、 に対して基底を貪欲に足していけばよい。
実装例
簡略化のため、 として以下を未満に言い換えている。
Submission #31329771 - Monoxer Programming Contest 2022(AtCoder Beginner Contest 249)
*1: の形式とどちらが綺麗なのかは分からない。