#451267 (C++14) No.1011 Infinite Stairs

• 問題: No.1011 Infinite Stairs
• ユーザー: ゆきのん
• 提出日時: 2020-03-29 19:42:48
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: TLE
• 実行時間: -
• コード長: 8,737 bytes
• コンパイル時間: 2,418 ms
• コンパイル使用メモリ: 186,628 KB
• 実行使用メモリ: 22,956 KB
• 最終ジャッジ日時: 2024-06-10 19:29:31
• 合計ジャッジ時間: 5,983 ms
• ジャッジサーバーID （参考情報）: judge1 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 13,756 KB
testcase_01	AC	1 ms; 6,816 KB
testcase_02	AC	200 ms; 6,944 KB
testcase_03	TLE	-
testcase_04	--	-
testcase_05	--	-
testcase_06	--	-
testcase_07	--	-
testcase_08	--	-
testcase_09	--	-
testcase_10	--	-
testcase_11	--	-
testcase_12	--	-
testcase_13	--	-
testcase_14	--	-
testcase_15	--	-
testcase_16	--	-
testcase_17	--	-
testcase_18	--	-
testcase_19	--	-
testcase_20	--	-
testcase_21	--	-
testcase_22	--	-
testcase_23	--	-
testcase_24	--	-
testcase_25	--	-
testcase_26	--	-

ソースコード

#include<bits/stdc++.h>
using namespace std;
#define fs first
#define sc second
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define ALL(A) A.begin(),A.end()
#define RALL(A) A.rbegin(),A.rend()
typedef long long LL;
typedef pair<LL,LL> P;
const LL mod=1000000007;
const LL LINF=1LL<<60;
const int INF=1<<30;
int dx[]={1,0,-1,0};
int dy[]={0,1,0,-1};


// Segment Tree Beats
// - l<=i<r について、 A_i の値を min(A_i, x) に更新
// - l<=i<r について、 A_i の値を max(A_i, x) に更新
// - l<=i<r の中の A_i の最大値を求める
// - l<=i<r の中の A_i の最小値を求める
// - l<=i<r の A_i の和を求める
// - l<=i<r について、 A_i の値に x を加える
// - l<=i<r について、 A_i の値を x に更新

class SegmentTree {
    struct Node {
        Node *left, *right;
        LL max_v, smax_v, max_c;
        LL min_v, smin_v, min_c;
        LL sum;
        LL len, ladd, lval;

        Node() : left(0), right(0), ladd(0), lval(LINF) {}

        void init(LL x) {
            max_v = min_v = sum = x;
            smax_v = -LINF; smin_v = LINF;
            max_c = min_c = 1;
        }

        void init_empty() {
            max_v = smax_v = -LINF;
            min_v = smin_v = LINF;
            max_c = min_c = 0;
        }

        void update_max(LL x) {
            sum += (x - max_v) * max_c;

            if(max_v == min_v) {
                max_v = min_v = x;
            } else if(max_v == smin_v) {
                max_v = smin_v = x;
            } else {
                max_v = x;
            }

            if(lval != LINF && x < lval) {
                lval = x;
            }
        }

        void update_min(LL x) {
            sum += (x - min_v) * min_c;

            if(max_v == min_v) {
                max_v = min_v = x;
            } else if(max_v == smin_v) {
                min_v = smax_v = x;
            } else {
                min_v = x;
            }

            if(lval != LINF && lval < x) {
                lval = x;
            }
        }

        void addall(LL x) {
            max_v += x;
            if(smax_v != -LINF) smax_v += x;
            min_v += x;
            if(smin_v != LINF) smin_v += x;

            sum += len * x;
            if(lval != LINF) {
                lval += x;
            } else {
                ladd += x;
            }
        }

        void updateall(LL x) {
            max_v = min_v = x;
            smax_v = -LINF; smin_v = LINF;
            max_c = min_c = len;

            sum = len * x;
            lval = x; ladd = 0;
        }

        void push() {

            if(lval != LINF) {
                left->updateall(lval);
                right->updateall(lval);
                lval = LINF;
                return;
            }

            if(ladd != 0) {
                left->addall(ladd);
                right->addall(ladd);
                ladd = 0;
            }

            if(max_v < left->max_v) {
                left->update_max(max_v);
            }
            if(left->min_v < min_v) {
                left->update_min(min_v);
            }

            if(max_v < right->max_v) {
                right->update_max(max_v);
            }
            if(right->min_v < min_v) {
                right->update_min(min_v);
            }
        }

        void update() {
            sum = left->sum + right->sum;

            if(left->max_v < right->max_v) {
                max_v = right->max_v;
                max_c = right->max_c;
                smax_v = max(left->max_v, right->smax_v);
            } else if(left->max_v > right->max_v) {
                max_v = left->max_v;
                max_c = left->max_c;
                smax_v = max(left->smax_v, right->max_v);
            } else {
                max_v = left->max_v;
                max_c = left->max_c + right->max_c;
                smax_v = max(left->smax_v, right->smax_v);
            }

            if(left->min_v < right->min_v) {
                min_v = left->min_v;
                min_c = left->min_c;
                smin_v = min(left->smin_v, right->min_v);
            } else if(left->min_v > right->min_v) {
                min_v = right->min_v;
                min_c = right->min_c;
                smin_v = min(left->min_v, right->smin_v);
            } else {
                min_v = left->min_v;
                min_c = left->min_c + right->min_c;
                smin_v = min(left->smin_v, right->smin_v);
            }
        }
    };

    int n, n0;
    Node *root;

    void _update_min(LL x, int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a || nd->max_v <= x) {
            return;
        }
        if(a <= l && r <= b && nd->smax_v < x) {
            nd->update_max(x);
            return;
        }

        nd->push();
        _update_min(x, a, b, nd->left, l, (l+r)/2);
        _update_min(x, a, b, nd->right, (l+r)/2, r);
        nd->update();
    }

    void _update_max(LL x, int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a || x <= nd->min_v) {
            return;
        }
        if(a <= l && r <= b && x < nd->smin_v) {
            nd->update_min(x);
            return;
        }

        nd->push();
        _update_max(x, a, b, nd->left, l, (l+r)/2);
        _update_max(x, a, b, nd->right, (l+r)/2, r);
        nd->update();
    }

    void _add_val(LL x, int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a) {
            return;
        }
        if(a <= l && r <= b) {
            nd->addall(x);
            return;
        }

        nd->push();
        _add_val(x, a, b, nd->left, l, (l+r)/2);
        _add_val(x, a, b, nd->right, (l+r)/2, r);
        nd->update();
    }

    void _update_val(LL x, int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a) {
            return;
        }
        if(a <= l && r <= b) {
            nd->updateall(x);
            return;
        }

        nd->push();
        _update_val(x, a, b, nd->left, l, (l+r)/2);
        _update_val(x, a, b, nd->right, (l+r)/2, r);
        nd->update();
    }

    LL _query_max(int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a) {
            return -LINF;
        }
        if(a <= l && r <= b) {
            return nd->max_v;
        }
        nd->push();
        LL lv = _query_max(a, b, nd->left, l, (l+r)/2);
        LL rv = _query_max(a, b, nd->right, (l+r)/2, r);
        return max(lv, rv);
    }

    LL _query_min(int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a) {
            return LINF;
        }
        if(a <= l && r <= b) {
            return nd->min_v;
        }
        nd->push();
        LL lv = _query_min(a, b, nd->left, l, (l+r)/2);
        LL rv = _query_min(a, b, nd->right, (l+r)/2, r);
        return min(lv, rv);
    }

    LL _query_sum(int a, int b, Node *nd, int l, int r) {
        if(b <= l || r <= a) {
            return 0;
        }
        if(a <= l && r <= b) {
            return nd->sum;
        }
        nd->push();
        LL lv = _query_sum(a, b, nd->left, l, (l+r)/2);
        LL rv = _query_sum(a, b, nd->right, (l+r)/2, r);
        return lv + rv;
    }

    public:
    SegmentTree(int n, vector<LL> a) : n(n) {
        n0 = 1;
        while(n0 < n) n0 <<= 1;

        Node *nds = new Node[2*n0];
        root = nds;

        nds[0].len = n0;
        for(int i=0; i<n0-1; ++i) {
            nds[i].left = &nds[2*i+1];
            nds[i].right = &nds[2*i+2];
            nds[2*i+1].len = nds[2*i+2].len = (nds[i].len >> 1);
        }

        for(int i=0; i<n; ++i) nds[n0-1+i].init(a[i]);
        for(int i=n; i<n0; ++i) nds[n0-1+i].init_empty();
        for(int i=n0-2; i>=0; i--) nds[i].update();
    }

    void update_min(int a, int b, LL x) {
        _update_min(x, a, b, root, 0, n0);
    }

    void update_max(int a, int b, LL x) {
        _update_max(x, a, b, root, 0, n0);
    }

    void add_val(int a, int b, LL x) {
        _add_val(x, a, b, root, 0, n0);
    }

    void update_val(int a, int b, LL x) {
        _update_val(x, a, b, root, 0, n0);
    }

    LL query_max(int a, int b) {
        return _query_max(a, b, root, 0, n0);
    }

    LL query_min(int a, int b) {
        return _query_min(a, b, root, 0, n0);
    }

    LL query_sum(int a, int b) {
        return _query_sum(a, b, root, 0, n0);
    }
};


int main(){
    int n,d,k;cin >> n >> d >> k;
    SegmentTree seg(k+d+1, vector<LL> (k+d+1,0));
    seg.update_val(0,1,1);
    for (int r = 0; r < n; r++) {
        for (int i = k; i >= 0; i--) {
            LL t = seg.query_sum(i,i+1)%mod;
            seg.add_val(i+1,i+d+1,t);
            seg.update_val(i,i+1,0);
        }
    }
    cout << seg.query_sum(k,k+1)%mod << endl;
    return 0;
}