yukicoder No.1067 #いろいろな色 / Red and Blue and more various colors (Middle)
お題箱より。
No.1067 #いろいろな色 / Red and Blue and more various colors (Middle) - yukicoder
解法
箱を並び替えても答えが変わることはないので、 の昇順に並べておきます。
各箱から玉を1個ずつ取り出し、ある1つの色 の玉がちょうど 個選ばれる場合の数が、どのように計算できるか図示してみましょう。
の昇順に並べたことによって、箱は色 の玉を含まないゾーンと含むゾーンに分かれます。含まないゾーンの取り出し方は の積であり、色 の玉はもちろん 個です。
色 の玉を含むゾーンの取り出し方は以下のDPによって求めることができます。
- 箱 全てに色 の玉が1個ずつ含まれているとき、これらの箱から玉を1個ずつ取り出し、色 の玉がちょうど 個選ばれる場合の数。
この値は、箱 の全てに含まれている色であればどの色でも同じ値になります。つまりクエリで必要な色 全てについて別々に計算する必要がありません。
初期条件は で、後ろから求めていきます。遷移は
- を遷移元とし、箱 から色 の玉を選ぶ。係数 で に遷移。
- を遷移元とし、箱 から色 でない玉を選ぶ。係数 で に遷移。
となります。DPパートの計算量は です。
クエリに答えるには「ある1つの色 の玉がちょうど 個選ばれる場合の数」をたくさんの について求める必要があります。これは以下のように計算することができます。
- まず、箱が色 を含まないゾーンと含むゾーンの境界を二分探索で探す。色 を含む箱の最小のインデックスを とする(存在しない場合は )。
- 色 を含まないゾーンの選び方は なので、事前に累積「積」を計算しておいてそれを参照する。
- 色 を含むゾーンの選び方はDPテーブルを参照する。具体的には である。
これは1つの あたり で計算できます。必要な の個数は なので、ちょっと厳しいですが間に合います。
ACコード
No.1067 #いろいろな色 / Red and Blue and more various colors (Middle) - yukicoder
公式解説では多項式の係数として解説していますが、1次式を掛けた時の係数の計算と今回のDPは全く同じものなので、解法としては同じです。
University of Aizu Programming Contest 2014 (AOJ 1548) Yu-kun Likes a Directed Graph
お題箱より。
解法
重み の負辺を足してできるだけ重み合計が大きい( に近い)負閉路を作りたいので、与えられたDAGについて「1本以上の辺で構成されるパスであって、重み合計が 未満であるもののうち一番大きいものの重み」が分かれば良いです。
ここで であること、さらに元々のDAGにある辺の重みは非負であることに注目します。これらの条件から、最終的に見つけるパスも、その一部となるパスも、重みが 以下のものだけを考えれば良いことが分かります。
よって以下のようなDPテーブルを構成し、 の範囲で計算していけば良いです。
- :1本以上の辺で構成され、頂点 で終わるパスであって、重み合計が となるものが存在するか?(true/false)
トポロジカル順序に基づいて順番に計算していきます。遷移は頂点 から頂点 への重み の辺が存在する時に、
- から始まり で終わるパス: をtrueにする
- 以外から始まり で終わるパス: がtrueである について、 をtrueにする
とすれば良いです。
時間計算量は となり高速な言語ならそのままでも間に合いそうですが、bitsetを用いて高速化することもできます。
ACコード
University of Aizu Programming Contest 2012(AOJ 1508)RMQ
お題箱より。
解法
1点更新と区間最小値取得ならセグメント木でできますが、シフト操作が厄介です。シフト操作は「数列から要素を1つ削除する」「その要素を数列の別の場所に挿入する」とみなすことができますが、この操作によって間に存在する要素のインデックスがずれます。そのため、インデックスに基づいて固定された構造を持つセグメント木では対応できません。
ではどうすれば良いのかというと、インデックスに依存しない構造を持つ平衡二分木を用いれば良いです。これなら要素の削除と追加ができます。
ということで、残りはこの問題と言うより平衡二分木の解説になってしまいますが…平衡二分木については非常に有名なスライドがあるのでこれを紹介します。
プログラミングコンテストでのデータ構造 2 ~平衡二分探索木編~
平衡二分木の話は少し紛らわしくて、解説やサンプルコードなどを読む際には以下の2種類のどちらを目的としたものなのかを意識する必要があります。
- C++の
std::set
のような順序付き集合としての機能。要素を順序付きで保持する。値の追加や削除、小さい方から 番目の要素の取得、大小関係に基づく二分探索、区間クエリの処理などができる。 - 数列の拡張としての機能。要素を列としての並びに基づいて保持する。 番目の要素の取得・挿入・削除、区間クエリの処理などができる。
個人的には前者を特に「平衡二分探索木」と呼び、主に後者を指す時に「平衡二分木」という言葉を使うようにしていますが、競プロerの中でも確かな呼び分けがされていないように思います。先ほどのスライドで紹介されているTreapの実装は後者で、この問題で必要なのも後者です。
このスライドの実装に区間クエリの処理機能を追加します。やることはセグメント木と同様で、まず各ノードは自身以下の部分木の値を全てマージした計算済みの値(今回の問題では最小値)を常に持つようにしておきます。この値はスライドの実装の cnt
や sum
と同様に、自身以下の部分木が変化した時に計算し直せば良いです。
クエリ処理する時には根から再帰的に見ていって、各ノードは親に以下のような値を返すようにします。
- 対象ノード以下の部分木が受け持っているインデックス範囲が、
- クエリ対象範囲と全く重複していないならば、単位元を返す。
- クエリ対象範囲に完全に含まれているならば、計算済みの値を返す。
- そうでないならば、左の子が返した値、(もしノード自身のインデックスが範囲に含まれているならば自身の値、)右の子が返した値をマージした値を返す。
セグメント木とは異なり各ノードが受け持つインデックス範囲が固定されていませんが、部分木のサイズ cnt
は分かるので、根から再帰的に呼び出すときに各ノードのオフセット(その部分木より左にいくつ要素があったか)を渡していくことで知ることができます。
ACコード
AtCoder Beginner Contest 170 E - Smart Infants
お題箱より。実装方法に関するリクエストだったので実装中心で書きます。
解法
解法としてはシミュレーションです。園児が移動したときに、どのようなデータを管理しておけば効率的に答えを更新できるかを考えます。
今回は以下のようなデータの流れを構築することにします。それぞれのデータに対してどのような操作が必要になるか順番に見ていきます。
まず園児が移動するときに、レートいくつの園児がどこからどこに移動するのか分からないといけません。それを知るために各園児のレート A
と所属幼稚園 B
を管理します。必要な機能は数列の1点を取得または変更することなので、配列でOKです。
園児の移動は、「ある幼稚園から指定したレートの園児を1人削除」して、「ある幼稚園に指定したレートの園児を1人追加」する操作です。またその操作の前後で、幼稚園の最強園児のレートがいくつになったか分かる必要があります。つまり幼稚園ごとのレートの集合 st
に求められる機能は
- 重複を含む値の集合を管理する。
- ある値を1つ削除する。
- ある値を1つ追加する。
- 集合内の最大値を調べる。
となります。
そしてこの変更で最強園児のレートが変わった場合、それを maxes
に反映する必要があります。そして maxes
の最小値が答えとなります。つまり maxes
に求められる機能は、st
とほぼ同様の
- 重複を含む値の集合を管理する。
- ある値を1つ削除する。
- ある値を1つ追加する。
- 集合内の最小値を調べる。
となります。このようにデータの流れを整理して各段階で必要な操作を列挙することで、各データをどのようなデータ構造で持てば良いのかが分かり、操作回数も合わせて見積もることで計算量を把握することができます。
あとは、この操作を効率的に行えるデータ構造を適切に選びます。
実現方法(C++)
C++であれば言語機能として std::multiset
というものが提供されています。これを使うと、操作時点での要素数を として必要な操作は全て で実現可能です。
これは値の集合を、順序(大小関係)を保持しながら管理してくれるデータ構造です。各操作は以下のように行うことができます。
st
からある値a
を1つ削除する:st.erase(st.find(a))
st
にある値a
を1つ追加する:st.insert(a)
st
の最小値を求める:*st.begin()
st
の最大値を求める:*st.rbegin()
存在しない値を削除しようとしたり、要素数 であるときに最小値や最大値を求めようとすると異常な挙動になる可能性があるので注意です。
begin()
や rbegin()
や find()
などの関数によって得られるイテレータというものについては以下を参考にしてください。
C++のイテレータを視覚的に理解する(競プロ向け) - ARMERIA
実現方法(それ以外)
C++以外の言語の場合は少し面倒です。今回の場合、必要になる操作のうち
- 重複を含む値の集合を管理する。
- ある値を1つ追加する。
- 集合内の最小値/最大値を調べる。
という操作であれば優先度付きキューで実現できるので、これに値を削除する仕組みを入れることを考えます。(Rubyなど優先度付きキューもない言語の場合は自作しましょう…)
具体的には、メインの優先度付きキューとは別に、「削除予約」をするための優先度付きキューを用意します。そして各操作は以下のように行います。
- 値を追加する:メインキューに値を追加する。
- 値を削除する:削除予約キューに値を追加する。
- 最小値/最大値を調べる:まず、メインキューと削除予約キューの最小値/最大値が一致する限り、それらをともに除外することを繰り返す。その後、メインキューの最小値/最大値を調べる。
ACコード
C++です。
CPSCO2019 Session3 F - Flexible Permutation O(N^2)解法
お題箱より。 解法が知りたいというリクエストだったので考えてみました。
解法
後に扱うDPの都合上、個人的にはインデックスよりも値に注目したほうが理解しやすいので、条件を値中心の言葉に言い換えます。 であるような値 を「左に動いた」、 であるような値 を「右に動いた」、 であるような値 を「動かない」と呼ぶことにします。
まず、動かない値 個を固定します。この選び方は 通りあり、これらの値を除くと残りは 個の必ず動く値だけを考えれば良いことになります。
公式解説の 解法と同様に、値を順番に足していきながら以前の値とswapしていくDPを考えます。ただしそのまま適用した場合、DPの途中では相変わらず「左に動いている」「右に動いている」「(一時的に)動いていない」値がいくつあるかという情報を管理する必要があり、状態数が落ちません。
そこで以下のようなものを数える包除原理を考えます。
- 予め指定したある 個の要素の集合について、値が動かない(違反状態である)。
- 個の要素について、値が左に動く。
- 残りの 個の要素について、値が動かないか、または右に動く。
つまり「ちょうど 個の値が左に動く」という条件は絶対に満たすとして、 個の要素について包除原理を適用する、ということです。こうすると先に 個の値を除いておけば、値の扱いで区別するべきものは「左に動いている」「それ以外」の2種類になり、状態数が に落ちます。
- 既に 個の値を置いていて、そのうち左に動いている値が 個あり、それ以外の値が 個あるような並べ方の数
遷移は次に置く値を「そのまま右端に置いた」「左に動いている値とswapした」「それ以外の値とswapした」ときにそれぞれ がどのように変化するかを考えれば良いです。具体的には「それ以外の値とswapした」ときだけ が つ増えることになります。
包除原理に基づいて結果を集計します。違反状態として指定した要素の個数が である場合、その 個の選び方が 通り、残りの場所の置き方が 通りあるので、これらの積に を掛けて全て合計すれば良いです。最後に を掛けると答えになります。
ACコード
Educational Codeforces Round 92 F. Bicolored Segments
お題箱より。
問題概要
個の閉区間 が与えられ、それぞれは赤または青のいずれか1色に塗られている。
「異なる色の2つの区間が共有点を持ってはいけない」という条件のもとで、これらのうちできるだけ多くの区間を選びたい。選ぶ個数の最大値を求めよ。
制約
解法1:DP
あらかじめ区間を右端でソートしておき、右端の値が小さいものから順に見ていくDPを考えます。インデックスもソート後の順番で振り直します。
赤の区間を採用できる条件は、これまで最後に採用した青の区間の右端が、採用しようとしている赤の区間の左端よりも小さいことです。逆も然りです。なので素直に考えると、
というDPを構成することができます。区間を1-indexedで扱い、「まだその色の区間を採用していない」をインデックス で表すと都合が良いです。これで計算量 では解けますが、間に合いません。
ではどうするのかというと、さらに状態をまとめます。以下のように2色それぞれで状態をまとめ、区間を順番に見ていきながら同時に更新していきます。
これは先ほどの という2次元テーブルで管理していた状態について、それぞれの軸方向ごとに最大値を取ったものになります。このDPテーブルをどのように更新するのか考えます。
まず2次元テーブルを用いたDPがどのような遷移になるのか確認しましょう。見ようとしている区間が赤だとします。更新対象となるのは、最後に採用した青の区間の右端が見ようとしている赤の区間の左端よりも小さいものです。つまり がある値以下の範囲で更新が行われることになり、図で表すと以下のようになります。
これを両軸方向に串刺しにして最大値を取った はどうなるでしょうか。図のインデックスを用いて説明します。まず について考えると、更新対象となっている をそれぞれ 増やせば良いことが分かります。そして は増やした後のこれらの値の最大値になります。
つまり2次元テーブルで状態を管理しなくても、これらの だけで遷移を計算することができます。全ての区間を見終わったら、 のどちらか(どちらでもOK)について最大値を取ったものが答えになります。
必要な操作は区間への加算と区間の最大値取得なので、遅延伝播セグメント木などで実現することができます。
ACコード1
Submission #88638709 - Codeforces
解法2:二部グラフ
公式解説で解説されている方法です。
各区間を頂点とみなし、同時に採用してはいけない区間の間に辺を結ぶグラフを考えると、これは二部グラフになります。求めたいのはこのグラフの最大独立集合の頂点数です。
二部グラフでは最大独立集合の頂点数は全頂点数から最大マッチングの数を引いたものになります。つまり最大マッチングを求めれば答えを計算できます。
参考:二部グラフの最小点被覆、最大安定集合 (最大独立集合)、最小辺被覆を総整理! - Qiita
ということで結局、赤と青の区間の組で共有点を持つものをなるべく多くマッチングさせることを考えれば良いです。これは例えば以下のように求めることができます。
このマッチング待ち区間のプールに必要な操作は「ある値を追加する」「プールに存在する値のうち、ある値以上で最も小さい値を取り出し、それを削除する」なので、C++であれば multiset
、このようなものが言語機能にない場合は座標圧縮+BIT+二分探索などで実現することができます。
ACコード2
エイシング プログラミング コンテスト 2020 F - Two Snuke
公式解説とは違う解法です。
解法
のときは答えは なので、 とします。
まずは 解法を
まず答えを で求める方法を考えます。
の値を として、これを全通り試しましょう。そうすると、この を つに分ける分け方が 通りあります。 個並んだボールを 本の仕切りで つの領域に分ける場合の数です。
として、 も同様に定義します。そうすると「全部足して 以下」という条件は と表記することができます。これを満たす の組それぞれについて (人名ではなく積です)を計算し、全て合計したものを求めたいです。これは、以下のような場合の数と等しくなります。
- 個のボールに 本の仕切りを入れて の つの領域に分け、さらに に相当する つの領域から つずつのボールを選ぶような場合の数。
これは言い換えると、ボールと仕切りを並べる合計 箇所のスペースから、仕切りと選んだボールに相当する 箇所を選び、その 箇所には左から順に「選んだボール、仕切り、…、選んだボール、仕切り」と割り当てるような場合の数に等しいです。つまり、この値は と計算できます。
の範囲としては が全て正になれるもの、つまり である範囲を考えます。そうすると答えは
と計算することができます。階乗を前計算して についてループを回すことで で解くことができますが、間に合いません。
偶奇に分けて多項式と見る
先ほどのシグマの中身を見てみましょう。 という値は、具体的には という についての 次多項式です。同様に は についての 次多項式です。つまりシグマの中身は と についての 次多項式です。より正確には、全ての項で と の次数合計が 以下であるような多項式となっています。
ここで重要な性質として、 についての 次多項式の値を まで合計したものは についての 次多項式になります。
今回の計算式におけるループ終端は でした。このままでは扱いにくいので の偶奇で場合分けします。例えば が奇数である場合、 と置きます。そうするとループの終端は となり、 を で置き換えるとシグマの中身は と についての 次多項式になります。先ほどの性質から、このシグマを処理した結果は についての 次多項式になることが分かります。
具体的に冪乗和の公式を当てはめて大量の計算をすればこの多項式の係数が求められますが、絶対にやりたくないです。「ラグランジュ補間」という方法を使いましょう。これは「 についての 次多項式 の自由度は であり、 通りの相異なる点 から多項式を一意に決定することができる」という性質を用いるものです。特に相異なる既知の点として に対する を用いた場合、1つの未知の に対する の値は で計算することができます。
つまり、 (すなわち )に対する答えをそれぞれ で計算し、これらを用いてラグランジュ補間を行えば、入力 が奇数である場合の答えを高速に求めることができます。
が偶数である場合も と置くことで同様に計算することができます。
ACコード
Submission #15191688 - AIsing Programming Contest 2020
ラグランジュ補間のライブラリを持っていれば実装は楽です。ラグランジュ補間の実装は、長くはないですが計算式がややこしいです。
実は かどうかを場合分けせずに または の結果からラグランジュ補間しても良いのですが、一応説明に従い場合分けをしています。