oj.uz

Submission #788507

# Submission time Handle Problem Language Result Execution time Memory
788507 2023-07-20T09:59:51 Z hugo_pm Chorus (JOI23_chorus) C++17
61 / 100
 7000 ms 5756 KB 
#include <bits/stdc++.h>
#define int long long
using namespace std;
	
#define all(v) (v).begin(), (v).end()
#define rall(v) (v).rbegin(), (v).rend()
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define sz(v) ((int)((v).size()))
	
template<typename T>
void chmax(T &x, const T &v) { if (x < v) x = v; }
template<typename T>
void chmin(T &x, const T &v) { if (x > v) x = v; }
	
using pii = pair<int, int>;
using vi = vector<int>;
	
string to_string(string s) { return s; }
template <typename T> string to_string(T v) {
	bool first = true;
	string res = "[";
	for (const auto &x : v) {
		if (!first)
			res += ", ";
		first = false;
		res += to_string(x);
	}
	res += "]";
	return res;
}
	
void dbg_out() { cout << endl; }
template <typename Head, typename... Tail> void dbg_out(Head H, Tail... T) {
	cout << ' ' << to_string(H);
	dbg_out(T...);
}
	
#ifdef DEBUG
#define dbg(...) cout << "(" << #__VA_ARGS__ << "):", dbg_out(__VA_ARGS__)
#else
#define dbg(...)
#endif
	
const int INF = 1e12;
struct _Line {
	mutable int k, m, p;
	bool operator<(const _Line& o) const { return k < o.k; }
	bool operator<(int x) const { return p < x; }
};

struct _LineContainer : multiset<_Line, less<>> {
	// (for doubles, use inf = 1/.0, div(a,b) = a/b)
	static const int inf = LLONG_MAX;
	int div(int a, int b) { // floored division
		return a / b - ((a ^ b) < 0 && a % b); }
	bool isect(iterator x, iterator y) {
		if (y == end()) return x->p = inf, 0;
		if (x->k == y->k) x->p = x->m > y->m ? inf : -inf;
		else x->p = div(y->m - x->m, x->k - y->k);
		return x->p >= y->p;
	}
	void add(int k, int m) {
		auto z = insert({k, m, 0}), y = z++, x = y;
		while (isect(y, z)) z = erase(z);
		if (x != begin() && isect(--x, y)) isect(x, y = erase(y));
		while ((y = x) != begin() && (--x)->p >= y->p)
			isect(x, erase(y));
	}
	int queryMax(int x) {
		assert(!empty());
		auto l = *lower_bound(x);
		return l.k * x + l.m;
	}
};

struct LineContainer {
	_LineContainer fs;
	void add (int k, int m) { // kx+m
		// -f(x) = -(kx+m) = (-k)x + (-m)
		fs.add(-k, -m);
	}
	int min(int x) {
		if (fs.empty()) return INF;
		return -fs.queryMax(x);
	}
};

int N, K;
string S;
// ferme[i] = nombre d'ouvrants avant Fi
vector<int> ouvre, ferme;
vector<int> cumuFerme;
int coutExclus(int left, int right) {
	if (left >= right) {
		return 0;
	}
	int projExclus = min(right, ouvre[right-1]);
	if (projExclus <= left) return 0;
	return (projExclus - left)*(right) - (cumuFerme[projExclus] - cumuFerme[left]);
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	
	cin >> N >> K >> S;
	{
		int curOuvre = 0, curFerme = 0;
		cumuFerme.push_back(0);
		for (char c : S) {
			if (c == 'A') {
				++curOuvre;
				ouvre.push_back(curFerme);
			} else {
				++curFerme;
				ferme.push_back(curOuvre);
				cumuFerme.push_back(cumuFerme.back() + curOuvre);
			}
		}
	}
	
	vector<int> old(N+1, INF), now;
	old[N] = 0;
	rep(nbCoup, 1, K+1) {
		now.assign(N+1, INF);
		LineContainer cont;
		// [J EXCLUS]
		// gright(left) = (projExclus - left)*(right) - (cumuFerme[projExclus] - cumuFerme[left])
		// gj(x) = (pj - x)*j - (cf[pj] - cf[left])
		// gj(x) = -jx + (j*pj - cf[pj]) + cf[x]
		//         mx  + p               + terme fixe
		auto getWithFix = [&] (int x) {
			return cont.min(x) + cumuFerme[x];
		};
		auto proj = [&] (int j) {
			return min(j, ouvre[j-1]);
		};
		auto add = [&] (int j) {
			int pj = proj(j);
			cont.add(-j, old[j] + j*pj - cumuFerme[pj]);
		};
		int jToAdd = N;
		for (int left = N-1; left >= 0; --left) {
			while (jToAdd > 0 && left <= min(proj(jToAdd), jToAdd-1)) {
				add(jToAdd--);
			}
			now[left] = getWithFix(left);
			if (jToAdd > 0) {
				chmin(now[left], old[jToAdd] + coutExclus(left, jToAdd));
			}
		}
		old = now;
	}
	cout << old[0] << '\n';
}

Subtask #1
16.0 / 16.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 212 KB Output is correct
4 Correct 0 ms 212 KB Output is correct
5 Correct 0 ms 212 KB Output is correct
6 Correct 1 ms 212 KB Output is correct
7 Correct 1 ms 212 KB Output is correct
8 Correct 1 ms 212 KB Output is correct
9 Correct 1 ms 324 KB Output is correct
10 Correct 1 ms 212 KB Output is correct
11 Correct 1 ms 212 KB Output is correct
12 Correct 0 ms 212 KB Output is correct
13 Correct 1 ms 212 KB Output is correct
14 Correct 0 ms 320 KB Output is correct
15 Correct 1 ms 212 KB Output is correct
16 Correct 0 ms 316 KB Output is correct

Subtask #2
24.0 / 24.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 212 KB Output is correct
4 Correct 0 ms 212 KB Output is correct
5 Correct 0 ms 212 KB Output is correct
6 Correct 1 ms 212 KB Output is correct
7 Correct 1 ms 212 KB Output is correct
8 Correct 1 ms 212 KB Output is correct
9 Correct 1 ms 324 KB Output is correct
10 Correct 1 ms 212 KB Output is correct
11 Correct 1 ms 212 KB Output is correct
12 Correct 0 ms 212 KB Output is correct
13 Correct 1 ms 212 KB Output is correct
14 Correct 0 ms 320 KB Output is correct
15 Correct 1 ms 212 KB Output is correct
16 Correct 0 ms 316 KB Output is correct
17 Correct 1 ms 340 KB Output is correct
18 Correct 4 ms 340 KB Output is correct
19 Correct 10 ms 340 KB Output is correct
20 Correct 1 ms 340 KB Output is correct
21 Correct 1 ms 340 KB Output is correct
22 Correct 17 ms 372 KB Output is correct
23 Correct 17 ms 340 KB Output is correct
24 Correct 1 ms 324 KB Output is correct
25 Correct 17 ms 376 KB Output is correct
26 Correct 15 ms 380 KB Output is correct
27 Correct 5 ms 352 KB Output is correct
28 Correct 5 ms 356 KB Output is correct

Subtask #3
21.0 / 21.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 212 KB Output is correct
4 Correct 0 ms 212 KB Output is correct
5 Correct 0 ms 212 KB Output is correct
6 Correct 1 ms 212 KB Output is correct
7 Correct 1 ms 212 KB Output is correct
8 Correct 1 ms 212 KB Output is correct
9 Correct 1 ms 324 KB Output is correct
10 Correct 1 ms 212 KB Output is correct
11 Correct 1 ms 212 KB Output is correct
12 Correct 0 ms 212 KB Output is correct
13 Correct 1 ms 212 KB Output is correct
14 Correct 0 ms 320 KB Output is correct
15 Correct 1 ms 212 KB Output is correct
16 Correct 0 ms 316 KB Output is correct
17 Correct 1 ms 340 KB Output is correct
18 Correct 4 ms 340 KB Output is correct
19 Correct 10 ms 340 KB Output is correct
20 Correct 1 ms 340 KB Output is correct
21 Correct 1 ms 340 KB Output is correct
22 Correct 17 ms 372 KB Output is correct
23 Correct 17 ms 340 KB Output is correct
24 Correct 1 ms 324 KB Output is correct
25 Correct 17 ms 376 KB Output is correct
26 Correct 15 ms 380 KB Output is correct
27 Correct 5 ms 352 KB Output is correct
28 Correct 5 ms 356 KB Output is correct
29 Correct 933 ms 836 KB Output is correct
30 Correct 2217 ms 800 KB Output is correct
31 Correct 559 ms 792 KB Output is correct
32 Correct 1 ms 596 KB Output is correct
33 Correct 2 ms 836 KB Output is correct
34 Correct 3186 ms 864 KB Output is correct
35 Correct 3255 ms 864 KB Output is correct
36 Correct 7 ms 720 KB Output is correct
37 Correct 553 ms 724 KB Output is correct
38 Correct 2645 ms 860 KB Output is correct
39 Correct 1213 ms 724 KB Output is correct
40 Correct 627 ms 660 KB Output is correct
41 Correct 624 ms 660 KB Output is correct
42 Correct 483 ms 716 KB Output is correct

Subtask #4
0 / 26.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 212 KB Output is correct
4 Correct 0 ms 212 KB Output is correct
5 Correct 0 ms 212 KB Output is correct
6 Correct 1 ms 212 KB Output is correct
7 Correct 1 ms 212 KB Output is correct
8 Correct 1 ms 212 KB Output is correct
9 Correct 1 ms 324 KB Output is correct
10 Correct 1 ms 212 KB Output is correct
11 Correct 1 ms 212 KB Output is correct
12 Correct 0 ms 212 KB Output is correct
13 Correct 1 ms 212 KB Output is correct
14 Correct 0 ms 320 KB Output is correct
15 Correct 1 ms 212 KB Output is correct
16 Correct 0 ms 316 KB Output is correct
17 Correct 1 ms 340 KB Output is correct
18 Correct 4 ms 340 KB Output is correct
19 Correct 10 ms 340 KB Output is correct
20 Correct 1 ms 340 KB Output is correct
21 Correct 1 ms 340 KB Output is correct
22 Correct 17 ms 372 KB Output is correct
23 Correct 17 ms 340 KB Output is correct
24 Correct 1 ms 324 KB Output is correct
25 Correct 17 ms 376 KB Output is correct
26 Correct 15 ms 380 KB Output is correct
27 Correct 5 ms 352 KB Output is correct
28 Correct 5 ms 356 KB Output is correct
29 Correct 933 ms 836 KB Output is correct
30 Correct 2217 ms 800 KB Output is correct
31 Correct 559 ms 792 KB Output is correct
32 Correct 1 ms 596 KB Output is correct
33 Correct 2 ms 836 KB Output is correct
34 Correct 3186 ms 864 KB Output is correct
35 Correct 3255 ms 864 KB Output is correct
36 Correct 7 ms 720 KB Output is correct
37 Correct 553 ms 724 KB Output is correct
38 Correct 2645 ms 860 KB Output is correct
39 Correct 1213 ms 724 KB Output is correct
40 Correct 627 ms 660 KB Output is correct
41 Correct 624 ms 660 KB Output is correct
42 Correct 483 ms 716 KB Output is correct
43 Execution timed out 7062 ms 5756 KB Time limit exceeded
44 Halted 0 ms 0 KB -

Subtask #5
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 212 KB Output is correct
4 Correct 0 ms 212 KB Output is correct
5 Correct 0 ms 212 KB Output is correct
6 Correct 1 ms 212 KB Output is correct
7 Correct 1 ms 212 KB Output is correct
8 Correct 1 ms 212 KB Output is correct
9 Correct 1 ms 324 KB Output is correct
10 Correct 1 ms 212 KB Output is correct
11 Correct 1 ms 212 KB Output is correct
12 Correct 0 ms 212 KB Output is correct
13 Correct 1 ms 212 KB Output is correct
14 Correct 0 ms 320 KB Output is correct
15 Correct 1 ms 212 KB Output is correct
16 Correct 0 ms 316 KB Output is correct
17 Correct 1 ms 340 KB Output is correct
18 Correct 4 ms 340 KB Output is correct
19 Correct 10 ms 340 KB Output is correct
20 Correct 1 ms 340 KB Output is correct
21 Correct 1 ms 340 KB Output is correct
22 Correct 17 ms 372 KB Output is correct
23 Correct 17 ms 340 KB Output is correct
24 Correct 1 ms 324 KB Output is correct
25 Correct 17 ms 376 KB Output is correct
26 Correct 15 ms 380 KB Output is correct
27 Correct 5 ms 352 KB Output is correct
28 Correct 5 ms 356 KB Output is correct
29 Correct 933 ms 836 KB Output is correct
30 Correct 2217 ms 800 KB Output is correct
31 Correct 559 ms 792 KB Output is correct
32 Correct 1 ms 596 KB Output is correct
33 Correct 2 ms 836 KB Output is correct
34 Correct 3186 ms 864 KB Output is correct
35 Correct 3255 ms 864 KB Output is correct
36 Correct 7 ms 720 KB Output is correct
37 Correct 553 ms 724 KB Output is correct
38 Correct 2645 ms 860 KB Output is correct
39 Correct 1213 ms 724 KB Output is correct
40 Correct 627 ms 660 KB Output is correct
41 Correct 624 ms 660 KB Output is correct
42 Correct 483 ms 716 KB Output is correct
43 Execution timed out 7062 ms 5756 KB Time limit exceeded
44 Halted 0 ms 0 KB -