結果

問題 No.2395 区間二次変換一点取得
ユーザー 👑 NachiaNachia
提出日時 2023-07-28 21:43:59
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 3,108 bytes
コンパイル時間 985 ms
コンパイル使用メモリ 84,248 KB
実行使用メモリ 10,496 KB
最終ジャッジ日時 2024-04-16 05:35:35
合計ジャッジ時間 4,486 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#line 1 "..\\Main.cpp"

#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <atcoder/modint>
#line 1 "D:\\Programming\\VSCode\\competitive-cpp\\nachia\\array\\dual-segment-tree.hpp"

#line 3 "D:\\Programming\\VSCode\\competitive-cpp\\nachia\\array\\dual-segment-tree.hpp"

namespace nachia{

template<
    class F,
    F composition(F f, F x)
>
struct DualSegtree {

    struct Node { F f; bool propagated; };
    int N;
    int logN;
    std::vector<Node> A;

    void mapf(Node& a, F f) {
        a.propagated = false;
        a.f = composition(f, a.f);
    }
    void spread(int i) {
        if(A[i].propagated || !(i < N)) return;
        mapf(A[i*2], A[i].f);
        mapf(A[i*2+1], A[i].f);
        A[i] = A[0];
    }

    DualSegtree(int n, F id) {
        N=1; logN=0;
        while(N<n){ N *= 2; logN++; }
        A.assign(N*2, { id, true });
    }
    DualSegtree(const std::vector<F>& a) : DualSegtree(a.size()) {
        for(int i=0; i<a.size(); i++){ A[i+N].f = a[i]; A[i+N].propagated = false; }
    }

    void clear(int p) {
        p += N;
        for(int d=logN; d; d--) spread(p >> d);
        A[p] = A[0];
    }
    F get(int p){
        p += N;
        for(int d=logN; d; d--) spread(p >> d);
        return A[p].f;
    }
    void apply(int l, int r, F f){
        if(!(l < r)) return;
        if(l == 0 && r == N){ mapf(A[1], f); return; }
        l += N; r += N;
        for(int d=logN; d; d--){
            if((l >> d) << d != l) spread(l >> d);
            if((r >> d) << d != r) spread(r >> d);
        }
        while(l < r){
            if(l&1){ mapf(A[l++], f); } l /= 2;
            if(r&1){ mapf(A[--r], f); } r /= 2;
        }
    }
    void apply(int p, F f){
        p += N;
        for(int d=logN; d; d--) spread(p >> d);
        mapf(A[p], f);
    }
};

} // namespace nachia;
#line 8 "..\\Main.cpp"
using namespace std;
using i32 = int;
using u32 = unsigned int;
using i64 = long long;
using u64 = unsigned long long;
#define rep(i,n) for(int i=0; i<(int)(n); i++)
const i64 INF = 1001001001001001001;

//using Modint = atcoder::static_modint<998244353>;
using Modint = atcoder::modint;

struct S{
    Modint z=1, xz=1, y=1;
};
struct F{
    Modint YY=1, XY=0, ZY=0;
    Modint XX=1, ZX=0;
    Modint ZZ=1;
};

F composition(F f, F x){
    return {
        x.YY*f.YY,
        x.XX*f.XY + x.XY*f.YY,
        x.ZZ*f.ZY + x.ZX*f.XY + x.ZY*f.YY,
        x.XX*f.XX,
        x.ZZ*f.ZX + x.ZX*f.XX,
        x.ZZ*f.ZZ
    };
}


int main(){
    int N; cin >> N;
    int B; cin >> B;
    Modint::set_mod(B);
    int Q; cin >> Q;
    auto rq = nachia::DualSegtree<F, composition>(N, F());
    rep(q,Q){
        int l,m,r; cin >> l >> m >> r; l--; m--;
        rq.apply(l, r, {3,2,2,3,3,3});
        auto g = rq.get(m);
        auto x = g.XX + g.ZX;
        auto y = g.XY + g.YY + g.ZY;
        auto z = g.ZZ;
        cout << x.val() << ' ' << y.val() << ' ' << z.val() << '\n';
    }
    return 0;
}


struct ios_do_not_sync{
    ios_do_not_sync(){
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
    }
} ios_do_not_sync_instance;
0