Codeforces Round #612 (Div. 1) D. LCC
異常実装になったけどなんとか自力で解けた…
問題概要
数直線上に 個の粒子がある。粒子 は初期位置 、速さ 、右に動く確率 で特徴づけられる。
この粒子を用いて以下の実験を行う。
- まず各粒子の移動方向を確率に従って決める。粒子 が右に動く確率は 、左に動く確率は である。
- 全ての粒子が時刻 にそれぞれの初期位置 からスタートし、先ほど決めた向きにそれぞれの速さ で動く。
- いずれか2つの粒子が衝突した時点で実験を終了し、その時刻を実験の経過時間とする。ただし永遠に粒子が衝突しない場合は経過時間 とする。
実験の経過時間の期待値を求めよ。この値は有理数になるので の形式で出力せよ。
制約
- は百分率で与えられ、
解法
記法
連続するいくつかの粒子の移動方向を、L
と R
を並べたものとして記します。例えば2つの粒子が向かい合って移動していることを RL
のように表記します。
以降の説明では は百分率の値ではなく確率そのものとして扱います。
衝突として考えるべきもの
いずれかの粒子が衝突した時点で終わりなので、考えるべき衝突は隣り合う2つの粒子のものだけです。粒子 と粒子 の衝突として以下のような場合があり得ます。
- 2粒子の移動方向が向かい合っている(
RL
)場合:時刻 で衝突する。 - 2粒子の移動方向がともに左(
LL
)であり、 である場合:時刻 で衝突する。 - 2粒子の移動方向がともに右(
RR
)であり、 である場合:時刻 で衝突する。
ここで、向かい合っている場合に衝突する時刻のほうが、同じ方向に動く場合に衝突する時刻より必ず早いことに注意しておきます。
全ての2粒子間においてこれらの衝突候補時刻を計算し、時刻順にソートしておきます。時刻は有理数として持っておいて、ソートのための比較関数では通分すると誤差がなくて良いです。
期待値を計算するアウトライン
まず期待値計算の概略を説明します。
時刻 で実験が継続している確率を とします。 であり、この値は先ほど計算した衝突候補時刻でのみ減少します。
を以下の図の赤線として図示すると、求めたい期待値は青く塗りつぶした領域の面積として計算できます。
衝突候補時刻 において が減少した量を と書くことにすると、この面積は を全ての衝突候補時刻について合計したものとして計算することができます。図の青領域を横に短冊状に切るイメージです。
これでこの問題は、各衝突候補時刻を順番に見ていきながら「各粒子の移動方向が、それまでに起こりうる衝突が全て発生しないような組み合わせになっている確率」を計算し続ける問題に帰着できました。
向かい合う衝突だけを考慮した確率計算
まず簡単のため、向かい合う衝突だけを考慮します。「向かい合う向きであってはいけない」という2粒子間をUnion-Find等で連結することにすると、非連結な成分同士は独立に計算できるので、結局「連続するいくつかの粒子において、どの隣り合う2粒子間も RL
の関係でない」という確率を求めれば良いことになります。
例えばインデックス の3粒子がこのような関係だとすると、衝突しないような向きの組み合わせは LLL
, LLR
, LRR
, RRR
の4通りです。その確率は
と計算できます。上記のような計算を区間 について行った結果を と表記することにします。
これを任意のインデックス区間について求めたいです。これは少しトリッキーなセグメント木を使って計算することができます。
それぞれの区間では、先ほどの とともに
として計算できるため、これをマージ演算とするセグメント木で任意区間の計算結果を求めることができます。計算のイメージは以下の図のような感じ。
Union-Findで連結成分を管理し、時刻 の時点での連結成分ごとの の値を全て掛け合わせたものが になります。
同じ向きでの衝突を考慮
同じ向きでの衝突を考慮するとさらに面倒になりますが、先ほどの延長で考えることができます。
同じ向きでの衝突が起こる時刻は向かい合う衝突が起こる時刻より遅いので、既に先ほどの連結成分となった区間に対して追加の条件が課されることになります。
以下の図を具体例として見てみます。
LL
や RR
になってはいけない箇所があることで、先ほどの の計算式のうち前の項・後ろの項がそれぞれ潰されます。つまり先ほどの連結成分に、前後どこまでの項が潰されているかの情報を追加で持たせればよいです。この生き残っている項の合計は、この図を例にすると と計算できるので、先ほどの延長で計算できます。生き残っている項がなくなった場合は衝突しない確率 、つまり実験が必ず終了することを意味します。
これで全てのケースが処理できました。あとは実装を頑張りましょう。