#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #ifdef LOCAL # include "debug_print.hpp" # define debug(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define debug(...) (static_cast(0)) #endif using namespace std; #define rep(i,n) for(int i=0; i<(n); i++) #define rrep(i,n) for(int i=(n)-1; i>=0; i--) #define FOR(i,a,b) for(int i=(a); i<(b); i++) #define RFOR(i,a,b) for(int i=(b-1); i>=(a); i--) #define ALL(v) v.begin(), v.end() #define RALL(v) v.rbegin(), v.rend() #define UNIQUE(v) v.erase( unique(v.begin(), v.end()), v.end() ); #define pb push_back using ll = long long; using D = double; using LD = long double; using P = pair; template using PQ = priority_queue>; template using minPQ = priority_queue, greater>; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b ostream &operator<<(ostream &os, const pair &p) { os << p.first << " " << p.second; return os; } template istream &operator>>(istream &is, pair &p) { is >> p.first >> p.second; return is; } template ostream &operator<<(ostream &os, const vector &v) { int s = (int)v.size(); for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i]; return os; } template istream &operator>>(istream &is, vector &v) { for (auto &x : v) is >> x; return is; } void in() {} template void in(T &t, U &...u) { cin >> t; in(u...); } void out() { cout << "\n"; } template void out(const T &t, const U &...u) { cout << t; if (sizeof...(u)) cout << sep; out(u...); } void outr() {} template void outr(const T &t, const U &...u) { cout << t; outr(u...); } template struct Fp { long long val; constexpr Fp(long long v = 0) noexcept : val(v % mod) { if (val < 0) val += mod; } constexpr int getmod() { return mod; } constexpr Fp operator - () const noexcept { return val ? mod - val : 0; } constexpr Fp operator + (const Fp& r) const noexcept { return Fp(*this) += r; } constexpr Fp operator - (const Fp& r) const noexcept { return Fp(*this) -= r; } constexpr Fp operator * (const Fp& r) const noexcept { return Fp(*this) *= r; } constexpr Fp operator / (const Fp& r) const noexcept { return Fp(*this) /= r; } constexpr Fp& operator += (const Fp& r) noexcept { val += r.val; if (val >= mod) val -= mod; return *this; } constexpr Fp& operator -= (const Fp& r) noexcept { val -= r.val; if (val < 0) val += mod; return *this; } constexpr Fp& operator *= (const Fp& r) noexcept { val = val * r.val % mod; return *this; } constexpr Fp& operator /= (const Fp& r) noexcept { long long a = r.val, b = mod, u = 1, v = 0; while (b) { long long t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v); } val = val * u % mod; if (val < 0) val += mod; return *this; } constexpr bool operator == (const Fp& r) const noexcept { return this->val == r.val; } constexpr bool operator != (const Fp& r) const noexcept { return this->val != r.val; } friend constexpr istream& operator >> (istream &is, const Fp& x) noexcept { int val; is >> val; x = Fp(val); return is; } friend constexpr ostream& operator << (ostream &os, const Fp& x) noexcept { return os << x.val; } friend constexpr Fp modpow(const Fp &a, long long n) noexcept { if (n == 0) return 1; auto t = modpow(a, n / 2); t = t * t; if (n & 1) t = t * a; return t; } }; //const int mod = 1e9+7; const int mod = 998244353; using mint = Fp; int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); ll n,m; in(n,m); vector z1(n),z2(n), w1(n),w2(n); vector zb,wb; rep(i,n){ ll x,y,h; in(x,y,h); ll zmin = x+y-m+h, zmax = x+y+m-h; ll wmin = y-x-m+h, wmax = y-x+m-h; z1[i] = zmin; z2[i] = zmax; w1[i] = wmin; w2[i] = wmax; zb.pb(zmin), zb.pb(zmax); wb.pb(wmin), wb.pb(wmax); } sort(ALL(zb)); sort(ALL(wb)); UNIQUE(zb); UNIQUE(wb); int sz = zb.size(), sw = wb.size(); vector> sum(sz, vector (sw)); rep(i,n){ int zi1 = lower_bound(ALL(zb), z1[i]) - zb.begin(); int zi2 = lower_bound(ALL(zb), z2[i]) - zb.begin(); int wi1 = lower_bound(ALL(wb), w1[i]) - wb.begin(); int wi2 = lower_bound(ALL(wb), w2[i]) - wb.begin(); sum[zi1][wi1]++; sum[zi2][wi2]++; sum[zi2][wi1]--; sum[zi1][wi2]--; } rep(i,sz){ rep(j,sw-1) sum[i][j+1] += sum[i][j]; } rep(j,sw){ rep(i,sz-1) sum[i+1][j] += sum[i][j]; } vector ans(n+1); rep(i,sz-1)rep(j,sw-1){ ans[sum[i][j]] += (mint)(zb[i+1]-zb[i])*(wb[j+1]-wb[j]); } FOR(i,1,n+1) out(ans[i]/2); }