Submission #759255

# Submission time Handle Problem Language Result Execution time Memory
759255 2023-06-16T01:04:54 Z NK_ UFO (IZhO14_ufo) C++17
100 / 100
681 ms 102468 KB
#include <bits/stdc++.h>

using namespace std;

/////////////////////// MACROS ////////////////////////////////////////////
using ll = long long;
using ld = long double;
using db = double;
using str = string;

using pi = pair<int,int>;
using pl = pair<ll,ll>;

using vi = vector<int>;
using vl = vector<ll>;
using vs = vector<str>;
using vc = vector<char>;
using vpi = vector<pi>;
using vpl = vector<pl>;

#define tcT template<class T
#define tcTU tcT, class U
tcT> using V = vector<T>;
tcT, size_t SZ> using AR = array<T,SZ>;
tcTU> using PR = pair<T,U>;
tcTU> using umap = unordered_map<T, U>;
tcT> using uset = unordered_set<T>;
tcT> using mset = multiset<T>;

#define mp make_pair
#define f first
#define s second

#define sz(x) int((x).size())
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define rsz resize
#define ins insert
#define ft front()
#define bk back()
#define ppb pop_back()
#define ppf pop_front()
#define pb push_back
#define eb emplace_back
#define pf push_front

#define lb lower_bound
#define ub upper_bound

// LOOPS
#define FOR(i,a,b) for (int i = (a); i < (b); ++i)
#define F0R(i,a) FOR(i,0,a)
#define ROF(i,a,b) for (int i = (b)-1; i >= (a); --i)
#define R0F(i,a) ROF(i,0,a)
#define rep(a) F0R(_, a)
#define each(a,x) for (auto& a: x)

/////////////////////// IMPORANT VARS /////////////////////////////////////

const int MOD = 1e9+7; // 998244353;
const int MX = 2e5+5;
const ll INFL = ll(3e18)+10;
const int INF = int(1e9)+10;
const ld PI = acos((ld)-1);
const int dx[4] = {1,0,-1,0}, dy[4] = {0,1,0,-1};
tcT> using pqg = priority_queue<T,vector<T>,greater<T>>;
tcT> using pql = priority_queue<T,vector<T>,less<T>>;
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
#define nl '\n'

constexpr int pct(int x) { return __builtin_popcount(x); } // # of bits set
constexpr int bits(int x) { 
	return x == 0 ? 0 : 31-__builtin_clz(x); }
constexpr int p2(int x) { return 1<<x; }
constexpr int msk2(int x) { return p2(x)-1; }
 
ll cdiv(ll a, ll b) { return a/b+((a^b)>0&&a%b); } // divide a by b rounded up
ll fdiv(ll a, ll b) { return a/b-((a^b)<0&&a%b); } // divide a by b rounded down
 
tcT> bool ckmin(T& a, const T& b) {
	return b < a ? a = b, 1 : 0; } // set a = min(a,b)
tcT> bool ckmax(T& a, const T& b) {
	return a < b ? a = b, 1 : 0; }
 
tcTU> T fstTrue(T lo, T hi, U f) {
	hi ++; assert(lo <= hi); // assuming f is increasing
	while (lo < hi) { // find first index such that f is true 
		T mid = lo+(hi-lo)/2;
		f(mid) ? hi = mid : lo = mid+1; 
	} 
	return lo;
}
tcTU> T lstTrue(T lo, T hi, U f) {
	lo --; assert(lo <= hi); // assuming f is decreasing
	while (lo < hi) { // find first index such that f is true 
		T mid = lo+(hi-lo+1)/2;
		f(mid) ? lo = mid : hi = mid-1;
	} 
	return lo;
}
tcT> void remDup(vector<T>& v) { // sort and remove duplicates
	sort(all(v)); v.erase(unique(all(v)),end(v)); }
tcTU> void erase(T& t, const U& u) { // don't erase
	auto it = t.find(u); assert(it != end(t));
	t.erase(it); } // element that doesn't exist from (multi)set


// #include <ext/pb_ds/assoc_container.hpp>
// using namespace __gnu_pbds;

// tcT> using iset = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>; 
// #define ook order_of_key
// #define fbo find_by_order

// struct chash { 
// 	const uint64_t C = ll(2e18*PI)+71;
// 	const int RANDOM = rng();
// 	ll operator()(ll x) const {
// 		return __builtin_bswap64((x^RANDOM)*C); }
// };
// template<class K,class V> using um = unordered_map<K,V,chash>;
// template<class K,class V> using ht = gp_hash_table<K,V,chash>;
// template<class K,class V> V get(ht<K,V>& u, K x) {
// 	auto it = u.find(x); return it == end(u) ? 0 : it->s; }

/////////////////////// OUPUT /////////////////////////////////////////////
#define ts to_string
str ts(char c) { return str(1,c); }
str ts(const char* s) { return (str)s; }
str ts(str s) { return s; }
str ts(bool b) {
	#ifdef LOCAL
		return b ? "true" : "false";
	#else
		return ts((int)b);
	#endif
}
tcTU> str ts(pair<T,U> p) {
	#ifdef LOCAL
		return "("+ts(p.f)+", "+ts(p.s)+")";
	#else
		return ts(p.f)+" "+ts(p.s);
	#endif
}

tcTU> str ts(V<pair<T, U>> v) {
	#ifdef LOCAL
		bool fst = 1; str res = "{";
		for (const auto& x: v) {
			if (!fst) res += ", ";
			fst = 0; res += ts(x);
		}
		res += "}"; return res;
	#else
		bool fst = 1; str res = "";
		for (const auto& x: v) {
			if (!fst) res += " ";
			fst = 0; res += ts(x);
		}
		return res;
	#endif
}

tcT> str ts(T v) {
	#ifdef LOCAL
		bool fst = 1; str res = "{";
		for (const auto& x: v) {
			if (!fst) res += ", ";
			fst = 0; res += ts(x);
		}
		res += "}"; return res;
	#else
		bool fst = 1; str res = "";
		for (const auto& x: v) {
			if (!fst) res += " ";
			fst = 0; res += ts(x);
		}
		return res;

	#endif
}

///////////////////////// DEBUG ///////////////////////////////////////////
#define tcTUU tcT, class ...U
void DBG() { cerr << "]" << "\e[0m" << endl; }
tcTUU> void DBG(const T& t, const U&... u) {
	cerr << ts(t); if (sizeof...(u)) cerr << ", ";
	DBG(u...); }
#ifdef LOCAL
	#define dbg(...) cerr << "\e[1m" << "Line(" << __LINE__ << ") -> [" << #__VA_ARGS__ << "]: [", DBG(__VA_ARGS__);
	#define chk(...) if (!(__VA_ARGS__)) cerr << "Line(" << __LINE__ << ") -> function(" \
		 << __FUNCTION__  << ") -> CHK FAILED: (" << #__VA_ARGS__ << ")" << "\n", exit(0);
#else
	#define dbg(...) 0
	#define chk(...) 0
#endif	

///////////////////////// FILE I/O ////////////////////////////////////////
void unsyncIO() { cin.tie(0)->sync_with_stdio(0); }
void setPrec() { cout << fixed << setprecision(15); }
void setIn(str s) { freopen(s.c_str(),"r",stdin); }
void setOut(str s) { freopen(s.c_str(),"w",stdout); }
void setIO(str s = "") {
	unsyncIO(); setPrec();
	#ifndef LOCAL	
		if (sz(s)) setIn(s+".in"), setOut(s+".out"); // for USACO
	#endif
}

///////////////////////// TEMPLATE ABOVE //////////////////////////////////

// REMEMBER
// - Don't Focus On Only One Approach
// - Read And Understand Problem Fully
// - Think Of Edges Cases
// - Implement Carefully
// - Always Check For Overflows
// - Reset Global Variables
// - Look At The Bigger Picture
// - Don't Get Discouraged, You Can Pull It Back

tcT> struct Seg { // comb(ID,b) = b
	const T ID = -INF; T comb(T a, T b) { return max(a, b); } 
	int n; V<T> seg;
	void init(int _n) { // upd, query also work if n = _n
		for (n = 1; n < _n; ) n *= 2; 
		// n = _n;
		seg.assign(2*n,ID); }
	void pull(int p) { seg[p] = comb(seg[2*p],seg[2*p+1]); }
	void set(int p, T val) { // set val at position p
		seg[p += n] = val; for (p /= 2; p; p /= 2) pull(p); }
	void upd(int p, T val) { // set val at position p
		seg[p += n] += val; for (p /= 2; p; p /= 2) pull(p); }

	int get(int p) { return seg[p += n]; }

	int first_at_least(int lo, int val, int ind, int l, int r) { // if seg stores max across range
		if (r <= lo || val > seg[ind]) return -1;
		if (l == r) return l;
		int m = (l+r)/2;
		int res = first_at_least(lo,val,2*ind,l,m); if (res != -1) return res;
		return first_at_least(lo,val,2*ind+1,m+1,r);
	}

	int get_first(int lo, int v) { return first_at_least(lo, v, 1, 0, n-1); }

	int last_at_least(int hi, int val, int ind, int l, int r) { // if seg stores max across range
		if (hi <= l || val > seg[ind]) return -1;
		if (l == r) return l;
		int m = (l+r)/2;
		int res = last_at_least(hi,val,2*ind+1,m+1,r); if (res != -1) return res;
		return last_at_least(hi,val,2*ind,l,m);
	}
	
	int get_last(int hi, int v) { return last_at_least(hi, v, 1, 0, n-1); }
};

void solve() {
	int N, M, T, K, X; cin >> N >> M >> T >> K >> X;

	V<Seg<int>> R(N), C(M); 
	F0R(i, N) R[i].init(M);
	F0R(i, M) C[i].init(N);

	F0R(i, N) {
		F0R(j, M) {
			int x; cin >> x;
			R[i].set(j, x);
			C[j].set(i, x);
		}
	}

	F0R(l_, K) {
		char dir; cin >> dir;
		dbg(l_, dir);
		if (dir == 'N') {
			int c, h; cin >> c >> h; --c;
			dbg(c, h);
			int cur = -1;

			rep(T) {
				int nxt = C[c].get_first(cur, h);
				if (nxt == -1) break;

				dbg(nxt, c, C[c].get(nxt), R[nxt].get(c));
				C[c].upd(nxt, -1);
				R[nxt].upd(c, -1);
				cur = nxt;
			}	
		}

		if (dir == 'S') {
			int c, h; cin >> c >> h; --c;
			dbg(c, h);
			int cur = N + 1;

			rep(T) {
				int nxt = C[c].get_last(cur, h);
				if (nxt == -1) break;

				dbg(nxt, c, C[c].get(nxt), R[nxt].get(c));
				C[c].upd(nxt, -1);
				R[nxt].upd(c, -1);
				cur = nxt;
			}
		}

		if (dir == 'W') {
			int r, h; cin >> r >> h; --r;
			dbg(r, h);
			int cur = -1;

			rep(T) {
				int nxt = R[r].get_first(cur, h);
				if (nxt == -1) break;

				dbg(r, nxt, R[r].get(nxt), C[nxt].get(r));
				R[r].upd(nxt, -1);
				C[nxt].upd(r, -1);
				cur = nxt;
			}
		}	

		if (dir == 'E') {
			int r, h; cin >> r >> h; --r;
			dbg(r, h);
			int cur = M + 1;

			rep(T) {
				int nxt = R[r].get_last(cur, h);
				if (nxt == -1) break;

				dbg(r, nxt, R[r].get(nxt), C[nxt].get(r));
				R[r].upd(nxt, -1);
				C[nxt].upd(r, -1);
				cur = nxt;
			}
		}
	}	

	V<vl> A(N, vl(M));
	F0R(i, N) F0R(j, M) {
		A[i][j] = R[i].get(j);
		assert(R[i].get(j) == C[j].get(i));
	}

	// F0R(i, N) cout << ts(A[i]) << nl;


	V<vl> P(N+1, vl(M+1));
	FOR(i, 1, N+1) FOR(j, 1, M+1) {
		P[i][j] = P[i-1][j] + P[i][j-1] - P[i-1][j-1] + A[i-1][j-1];
	}

	auto qry = [&](int r1, int c1, int r2, int c2) {
		return P[r2][c2] - P[r1-1][c2] - P[r2][c1-1] + P[r1-1][c1-1];
	};

	ll ans = 0;
	FOR(r, 1, N+1) FOR(c, 1, M+1) if (r+X-1 <= N && c+X-1 <= M) {
		ckmax(ans, qry(r, c, r+X-1, c+X-1));
	}

	cout << ans << nl;
}

int main() {
	setIO(); 

	int TT = 1;
	// cin >> TT;

	rep(TT) solve();
	
	exit(0-0);
}

Compilation message

ufo.cpp: In function 'void solve()':
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:275:3: note: in expansion of macro 'dbg'
  275 |   dbg(l_, dir);
      |   ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:278:4: note: in expansion of macro 'dbg'
  278 |    dbg(c, h);
      |    ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:285:5: note: in expansion of macro 'dbg'
  285 |     dbg(nxt, c, C[c].get(nxt), R[nxt].get(c));
      |     ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:294:4: note: in expansion of macro 'dbg'
  294 |    dbg(c, h);
      |    ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:301:5: note: in expansion of macro 'dbg'
  301 |     dbg(nxt, c, C[c].get(nxt), R[nxt].get(c));
      |     ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:310:4: note: in expansion of macro 'dbg'
  310 |    dbg(r, h);
      |    ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:317:5: note: in expansion of macro 'dbg'
  317 |     dbg(r, nxt, R[r].get(nxt), C[nxt].get(r));
      |     ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:326:4: note: in expansion of macro 'dbg'
  326 |    dbg(r, h);
      |    ^~~
ufo.cpp:194:19: warning: statement has no effect [-Wunused-value]
  194 |  #define dbg(...) 0
      |                   ^
ufo.cpp:333:5: note: in expansion of macro 'dbg'
  333 |     dbg(r, nxt, R[r].get(nxt), C[nxt].get(r));
      |     ^~~
ufo.cpp: In function 'void setIn(str)':
ufo.cpp:201:28: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  201 | void setIn(str s) { freopen(s.c_str(),"r",stdin); }
      |                     ~~~~~~~^~~~~~~~~~~~~~~~~~~~~
ufo.cpp: In function 'void setOut(str)':
ufo.cpp:202:29: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  202 | void setOut(str s) { freopen(s.c_str(),"w",stdout); }
      |                      ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 0 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
4 Correct 8 ms 656 KB Output is correct
5 Correct 45 ms 2596 KB Output is correct
6 Correct 115 ms 19268 KB Output is correct
7 Correct 223 ms 44980 KB Output is correct
8 Correct 167 ms 44996 KB Output is correct
9 Correct 306 ms 41912 KB Output is correct
10 Correct 360 ms 44976 KB Output is correct
11 Correct 243 ms 38560 KB Output is correct
12 Correct 372 ms 44976 KB Output is correct
13 Correct 442 ms 51260 KB Output is correct
14 Correct 319 ms 38592 KB Output is correct
15 Correct 370 ms 44980 KB Output is correct
16 Correct 184 ms 38572 KB Output is correct
17 Correct 525 ms 51288 KB Output is correct
18 Correct 157 ms 43208 KB Output is correct
19 Correct 235 ms 51300 KB Output is correct
20 Correct 681 ms 102468 KB Output is correct