結果

問題 No.1222 -101
ユーザー KudeKude
提出日時 2020-09-05 01:18:41
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 370 ms / 2,000 ms
コード長 5,314 bytes
コンパイル時間 2,080 ms
コンパイル使用メモリ 199,084 KB
実行使用メモリ 20,424 KB
最終ジャッジ日時 2024-05-05 04:01:10
合計ジャッジ時間 6,822 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 7 ms
7,936 KB
testcase_11 AC 5 ms
7,936 KB
testcase_12 AC 220 ms
13,584 KB
testcase_13 AC 220 ms
13,396 KB
testcase_14 AC 221 ms
13,480 KB
testcase_15 AC 117 ms
10,856 KB
testcase_16 AC 165 ms
12,284 KB
testcase_17 AC 199 ms
13,000 KB
testcase_18 AC 115 ms
10,808 KB
testcase_19 AC 142 ms
11,408 KB
testcase_20 AC 179 ms
12,328 KB
testcase_21 AC 169 ms
12,260 KB
testcase_22 AC 2 ms
5,376 KB
testcase_23 AC 2 ms
5,376 KB
testcase_24 AC 2 ms
5,376 KB
testcase_25 AC 1 ms
5,376 KB
testcase_26 AC 1 ms
5,376 KB
testcase_27 AC 2 ms
5,376 KB
testcase_28 AC 2 ms
5,376 KB
testcase_29 AC 2 ms
5,376 KB
testcase_30 AC 1 ms
5,376 KB
testcase_31 AC 2 ms
5,376 KB
testcase_32 AC 363 ms
20,424 KB
testcase_33 AC 370 ms
20,296 KB
testcase_34 AC 148 ms
9,552 KB
testcase_35 AC 148 ms
9,672 KB
testcase_36 AC 158 ms
9,548 KB
testcase_37 AC 166 ms
9,544 KB
testcase_38 AC 155 ms
9,548 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define chmax(a, b) a = max(a, b)
#define chmin(a, b) a = min(a, b)
#define all(x) (x).begin(), (x).end()
using namespace std;
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
#define MOD 1000000007
//#define MOD 998244353
struct mint {
	int i;
	mint() : i(0) {}
	mint(int x) {
		i = int(x % MOD);
		if (i < 0) i += MOD;
	}
	template<class T> mint(T x) {
		i = int(x % MOD);
		if (i < 0) i += MOD;
	}
	mint operator+(const mint x) const {return i + x.i;}
	mint operator-(const mint x) const {return i - x.i;}
	mint operator*(const mint x) const {return (long long)i * x.i;}
	mint operator/(const mint x) const {return (long long)i * x.pow(MOD - 2).i;}
	mint inv() {return pow(MOD - 2);}
	template<class T> mint pow(T p) const {
		long long r = 1;
		long long t = i;
		for(; p > 0; p >>= 1) {
			if (p & 1) r = r * t % MOD;
			t = t * t % MOD;
		}
		return r;
	}
	template<class T1, class T2> static mint pow(T1 a, T2 b) {
		long long r = 1;
		long long t = (long long)(a % MOD);
		for(; b > 0; b >>= 1) {
			if (b & 1) r = r * t % MOD;
			t = t * t % MOD;
		}
		return r;
	}
	mint& operator+=(const mint x) {
		i = (i + x.i) % MOD;
		return *this;
	}
	mint& operator-=(const mint x) {
		i = i - x.i;
		if (i < 0) i += MOD;
		return *this;
	}
	mint& operator*=(const mint x) {
		i = (int)((long long)i * x.i % MOD);
		return *this;
	}
	mint& operator/=(const mint x) {
		i = (long long)i * x.pow(MOD - 2).i % MOD;
		return *this;
	}
};

template<typename T>
struct BIT {
  int n;
  vector<T> d;
  BIT(int n=0):n(n),d(n+1) {}
  void add(int i, T x=1) {
    for (i++; i <= n; i += i&-i) {
      d[i] += x;
    }
  }
  T sum(int i) {
    T x = 0;
    for (i++; i; i -= i&-i) {
      x += d[i];
    }
    return x;
  }
  T sum(int l, int r) {
    return sum(r-1) - sum(l-1);
  }
  T operator[](int x) {
    return sum(x) - sum(x-1);
  }
};

std::ostream& operator<<(std::ostream& os, const mint& m) {
	return os << m.i;
}

struct COND {
    int l;
    int r;
    COND(int l, int r):l(l),r(r){}
};

int main() {
    int n, m;
    cin >> n >> m;
    vector<COND> c0, c1;
    rep(i, m) {
        int l, r, p;
        cin >> l >> r >> p;
        l--; r--;
        if (p == 0) c0.emplace_back(l, r);
        else c1.emplace_back(l, r);
    }
    map<int, P> cover;
    for(auto x: c1) {
        int l = x.l, r = x.r, cnt = 1;
        auto it1 = cover.lower_bound(l);
        if (it1 != cover.begin()) {
            auto it0 = it1;
            it0--;
            if (l <= it0->second.first) {
                l = it0->first;
                cnt += it0->second.second;
                chmax(r, it0->second.first);
                cover.erase(it0);
            }
        }
        if (it1 != cover.end() && it1->first <= r) {
            chmax(r, it1->second.first);
            cnt += it1->second.second;
            cover.erase(it1);
        }
        cover[l] = {r, cnt};
    }
    vector<bool> filled(n);
    mint ans = 1;
    for(auto x: cover) {
        int l = x.first, r = x.second.first, cnt = x.second.second;
        ans *= mint::pow(2, r - l + 1 - cnt);
        for(int i = l; i <= r; i++) filled[i] = 1;
    }

    int sz = c0.size();
    {
        int idx = 0;
        vector<vector<P>> q(n);
        set<int> alive;
        rep(i, sz) {
            q[c0[i].l].emplace_back(0, i);
            q[c0[i].r].emplace_back(1, i);
        }
        rep(i, n) {
            if (filled[i]) {
                for(auto p: q[i]) {
                    if (p.first == 0) {
                        c0[p.second].l = idx;
                        alive.insert(p.second);
                    } else {
                        c0[p.second].r = idx - 1;
                        alive.erase(p.second);
                    }
                }
            } else {
                if (!alive.empty()) filled[i] = 1;
                for(auto p: q[i]) {
                    filled[i] = 1;
                    if (p.first == 0) {
                        c0[p.second].l = idx;
                        alive.insert(p.second);
                    } else {
                        c0[p.second].r = idx;
                        alive.erase(p.second);
                    }
                }
                idx++;
            }
        }
    }
    for(auto p: c0) {
        if (p.l > p.r) {
            cout << 0 << endl;
            return 0;
        }
    }

    sort(all(c0), [](COND x, COND y){return x.l < y.l;});
    rep(i, sz) {
        int l = c0[i].l, r = c0[i].r;
        vector<COND> t;
        t.emplace_back(0, r - l);
        while(i + 1 < sz && c0[i+1].l <= r) {
            i++;
            t.emplace_back(c0[i].l - l, c0[i].r - l);
            chmax(r, c0[i].r);
        }
        int num = r - l + 1;
        sort(all(t), [](COND x, COND y){return x.r < y.r;});
        BIT<mint> dp(num + 1);
        int lb = 0, nxtcid = 0;
        mint tinv = mint(1) / 2;
        dp.add(0, mint::pow(2, num));
        rep(i, num) {
            dp.add(i + 1, dp.sum(lb, i + 1) * tinv);
            if (i == t[nxtcid].r) {
                chmax(lb, t[nxtcid].l + 1);
                nxtcid++;
            }
        }
        ans *= dp.sum(lb, num + 1);
    }
    rep(i, n) if (!filled[i]) ans *= 3;
    cout << ans << endl;
}
0