はまやんはまやんはまやん

hamayanhamayan's blog

注文の多い順列 [yukicoder 1001]

https://yukicoder.me/problems/no/1001

前提知識

解説

https://yukicoder.me/submissions/436553

難しい。
実装のため、数は全部0~N-1に直しておく。
考え方としては、挿入DPが近いかもしれない。
そっちを学んでから、以下解説を見るとわかりやすさが違うかも。

まず、条件がt=1だけだった場合を考えて、以下のようなDPを考えてみる。
dp[i] := 数i以前まで配置済みのときの順列の組み合わせ
このDPができる大切な要素として、単調性がある。
例えば条件として、全部t=1として{1,4,3,0,2}だったとする。
数0が置けるのは、4番目だけで、数1が置けるのは1番目と4番目。
数2が置けるのは、1,4,5番目…といったように数が大きくなると、置ける範囲が単調増加していく。
つまり、何が言いたいかというと、ある数iが置ける選択肢がk通りあったとして、1つを選んだ場合、
残りのk-1個はそのあとの数を自由に置くことができるということ。
なので、更新の選択肢は常に「置ける個数-既に置いた個数」となる。
よって、ある数iが置ける個数を以下のように前計算しておけば、簡単に計算可能。
cnt1[i] := 数iが置ける場所の個数
dp[i + 1] = dp[i] * (cnt[i + 1] - i)

だが、今回はt=0もある。
基本的な考え方は同じである。
dp[done0][done1] := 数を0から順番に配置して、t=0の部分をdone0個配置済み、t=1の部分をdone1個配置済みのときの組み合わせ
配列cnt1と同様に配列cnt0も作っておこう。
t=1に配置する操作は上と同じようにやればいい。

問題はt=0に配置する操作である(正直こっちあんまし分かってない。間違っていたらすみません)。
これはt=1と異なり、数を増やしていくと候補が単調減少していく。
こういう範囲が狭まっていくものの考え方として、それより後が構築可能な状態を保ちながら貪欲というテクがある。
これを応用するような感じで解いていく。
結論から言うと、更新の選択肢は(cnt0[nxt] - (n0 - (done0 + 1)))個となる。
cnt0[nxt]が次置く数が置ける場所の個数
n0がt=0の全体個数
done0が既に置いているt=0の個数
今置ける組み合わせはcnt0[nxt]通りであるが、全部使える訳じゃない。
置ける範囲はどんどん減っていくので、ここで置ける場所というのは後半のために一部残しておく必要がある。
残しておくべき個数は今置く予定のものを含めると、n0 - (done0 + 1)となる。
done0+1で今置く予定を含めた置いた個数。全体からそれを引くと、後々置くべき個数となる。
この分はあとで置けるようにとっておく必要がある。
なので、それを置ける場所cnt0[nxt]から引くと、置ける場所の組み合わせが求められる。

int N;
vector<int> X[2];
mint dp[3010][3010];
int cnt0[3010], cnt1[3010];
//---------------------------------------------------------------------------------------------------
void _main() {
    cin >> N;
    rep(i, 0, N) {
        int t, x; cin >> t >> x;
        X[t].push_back(x - 1);
    }

    int n0 = X[0].size(), n1 = X[1].size();

    rep(i, 0, N) fore(x, X[0]) if (i <= x) cnt0[i]++;
    rep(i, 0, N) fore(x, X[1]) if (x <= i) cnt1[i]++;

    dp[0][0] = 1;
    rep(done0, 0, n0 + 1) rep(done1, 0, n1 + 1) {
        int nxt = done0 + done1;

        if (done0 < n0) dp[done0 + 1][done1] += dp[done0][done1] * (cnt0[nxt] - (n0 - (done0 + 1)));
        if (done1 < n1) dp[done0][done1 + 1] += dp[done0][done1] * (cnt1[nxt] - done1);
    }

    cout << dp[n0][n1] << endl;
}