結果

問題 No.2520 L1 Explosion
コンテスト
ユーザー 👑 binap
提出日時 2023-10-27 23:28:50
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 2,721 bytes
コンパイル時間 4,105 ms
コンパイル使用メモリ 264,188 KB
最終ジャッジ日時 2025-02-17 16:05:12
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<atcoder/all>
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef vector<int> vi;
typedef vector<long long> vl;
typedef vector<vector<int>> vvi;
typedef vector<vector<long long>> vvl;
typedef long double ld;
typedef pair<int, int> P;

template <int m>
std::ostream& operator<<(std::ostream& os, const atcoder::static_modint<m>& a) {os << a.val(); return os;}
template<typename T>
istream& operator>>(istream& is, vector<T>& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;}
template<typename T>
ostream& operator<<(ostream& os, const vector<T>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;}
template <typename T>
ostream& operator<<(ostream& os, const vector<vector<T>>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;}

// Coodinate Compression
// https://youtu.be/fR3W5IcBGLQ?t=8550
template<typename T=int>
struct CC {
  bool initialized;
  vector<T> xs;
  CC(): initialized(false) {}
  void add(T x) { xs.push_back(x);}
  void init() {
    sort(xs.begin(), xs.end());
    xs.erase(unique(xs.begin(),xs.end()),xs.end());
    initialized = true;
  }
  int operator()(T x) {
    if (!initialized) init();
    return upper_bound(xs.begin(), xs.end(), x) - xs.begin() - 1;
  }
  T operator[](int i) {
    if (!initialized) init();
    return xs[i];
  }
  int size() {
    if (!initialized) init();
    return xs.size();
  }
};

using mint = modint998244353;

int main(){
	int n;
	ll m;
	cin >> n >> m;
	vl x(n), y(n), h(n);
	rep(i, n) cin >> x[i] >> y[i] >> h[i];
	CC<ll> cc_s, cc_t;
	rep(i, n) cc_s.add((x[i] - (m - h[i])) + y[i]);
	rep(i, n) cc_s.add((x[i] + (m - h[i])) + y[i]);
	rep(i, n) cc_t.add((x[i] - (m - h[i])) - y[i]);
	rep(i, n) cc_t.add((x[i] + (m - h[i])) - y[i]);
	int ks = cc_s.size();
	int kt = cc_t.size();
	vector<vector<tuple<int, int, int>>> event(kt);
	rep(i, n){
		event[cc_t((x[i] - (m - h[i])) - y[i])].emplace_back(+1, cc_s(x[i] - (m - h[i]) + y[i]), cc_s(x[i] + (m - h[i]) + y[i]));
		event[cc_t((x[i] + (m - h[i])) - y[i])].emplace_back(-1, cc_s(x[i] - (m - h[i]) + y[i]), cc_s(x[i] + (m - h[i]) + y[i]));
	}
	/*
	cout << cc_s.xs;
	cout << cc_t.xs;
	*/
	vi memo(ks);
	vector<mint> ans(n + 1);
	rep(t, kt){
		auto& v = event[t];
		if(t > 0){
			rep(s, ks - 1){
				ans[memo[s]] += (cc_t[t] - cc_t[t - 1]) * (cc_s[s + 1] - cc_s[s]);
			}
		}
//		cout << v.size() << "\n";
		for(auto [parity, from, to] : v){
			for(int s = from; s < to; s++){
				memo[s] += parity;
			}
		}
		/*
		cout << memo;
		cout << ans;
		cout << "\n";
		*/
	}
	for(int i = 1; i <= n; i++){
		cout << ans[i] / 2<< "\n";
	}
	return 0;
}
0