Submission #963108

# Submission time Handle Problem Language Result Execution time Memory
963108 2024-04-14T14:21:03 Z penguin133 Aliens (IOI16_aliens) C++17
4 / 100
2 ms 4600 KB
#include <bits/stdc++.h>
using namespace std;
#include "aliens.h"
//#define int long long
typedef long long ll;
typedef long double ld;
#define pi pair<long double, ll>
#define pii pair<int, pi>
#define fi first
#define se second
#ifdef _WIN32
#define getchar_unlocked _getchar_nolock
#endif
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());

deque<pii>dq;

ld f(pi x, ll y){
	return x.fi * y + x.se;
}

long double intersect(int a, ld b, int c, ld d){
	return (long double)(d-b)/(a-c);
}

void ins(ll m, ld c, ll in){
	while(dq.size() > 1){
		int s = dq.size();
		if(dq[s-1].fi == m){
			if(dq[s-1].se.fi < c)dq.pop_back();
			else break;
			continue;
		}
		if(intersect(dq[s-1].fi, dq[s-1].se.fi, m, c) <= intersect(dq[s-2].fi, dq[s-2].se.fi, m, c))dq.pop_back();
		else break;
	}
	dq.push_back({m, {c, in}});
	
}

pi query(ll x){
	while(dq.size() > 1){
		if(f({dq[0].fi, dq[0].se.fi}, x) < f({dq[1].fi, dq[1].se.fi}, x))dq.pop_front();
		else break;
	}
	return {f({dq[0].fi, dq[0].se.fi}, x), dq[0].se.se};
}

ll mn[1000005], val[100005], N;
pi dp[100005];
pair <ll, ll> stuf[1000005];

pi cst(ld x){
	dq.clear();
	for(int i = 0; i < N; i++){
		ins(mn[i] * 2, -((i ? dp[i - 1].fi : 0) - 2 * mn[i] + mn[i] * mn[i] - (i ? max(0ll, stuf[i - 1].fi - mn[i] + 1) * (stuf[i - 1].fi - mn[i] + 1) : 0)), (i ? dp[i - 1].se : 0));
		pi res = query(stuf[i].fi);
		res.fi = -res.fi;
		res.fi += stuf[i].fi * stuf[i].fi + 2 * stuf[i].fi + 1;
		res.fi += x;
		res.se++;
		dp[i] = res;
	}
	return dp[N - 1];
}

long long take_photos(int n, int m, int k, std::vector<int> r, std::vector<int> c) {
	N = n;
	for(int i = 0; i <= m; i++)mn[i] = 1e18;
	for(int i = 0; i < n; i++){
		if(r[i] > c[i])swap(r[i], c[i]);
		mn[c[i]] = min(mn[c[i]], (ll)r[i]);
		stuf[i] = {c[i], r[i]};
	}
	sort(stuf, stuf + n);
	for(int i = 0; i < n; i++)mn[i] = stuf[i].se, val[i] = 1e18;
	val[n] = 1e18;
	for(int i = n - 2; i >= 0; i--)mn[i] = min(mn[i], mn[i + 1]);
	//return 0;
	/*
	ll dp[n + 1][k + 1];
	for(int i = 0; i <= n; i++)dp[i][0] = 1e18;

	for(int j = 1; j <= k; j++){
		for(int i = 0; i < n; i++){
			ins(mn[i] * 2, -((i ? dp[i - 1][j - 1] : 0) - 2 * mn[i] + mn[i] * mn[i] - (i ? max(0ll, stuf[i - 1].fi - mn[i] + 1) * (stuf[i - 1].fi - mn[i] + 1) : 0)));
			dp[i][j] = -query(stuf[i].fi) + stuf[i].fi * stuf[i].fi + 2 * stuf[i].fi + 1;
			
			dp[i][j] = 1e18;
			for(int a = i; a >= 0; a--){
				dp[i][j] = min(dp[i][j], (a ? dp[a - 1][j - 1] - max(0ll, stuf[a - 1].fi - mn[a] + 1) * (stuf[a - 1].fi - mn[a] + 1) : 0) + (stuf[i].fi - mn[a] + 1) * (stuf[i].fi - mn[a] + 1));
			}
			
			//cout << dp[i][j] << ' ';
		}
		//cout << '\n';
		dq.clear();
	}
	*/
	ld lo = 0, hi = 1e12;
	ld ans = 1e18;
	for(int i = 0; i < 100; i++){
		ld mid = (lo + hi) / 2;
		pi res = cst(mid);
		if(res.se <= k)ans = min(ans, res.fi - (ld)mid * res.se);
		if(res.se > k)lo = mid;
		else hi = mid;
	}
	pi brr = cst(hi);
	return (ll)round(brr.fi - hi * brr.se);
}
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 4
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4440 KB Correct answer: answer = 4
4 Correct 1 ms 4440 KB Correct answer: answer = 12
5 Correct 1 ms 4444 KB Correct answer: answer = 52
6 Correct 1 ms 4444 KB Correct answer: answer = 210
7 Correct 1 ms 4444 KB Correct answer: answer = 88
8 Correct 1 ms 4444 KB Correct answer: answer = 7696
9 Correct 1 ms 4444 KB Correct answer: answer = 1
10 Correct 1 ms 4444 KB Correct answer: answer = 2374
11 Correct 2 ms 4440 KB Correct answer: answer = 9502
12 Correct 2 ms 4444 KB Correct answer: answer = 49
13 Correct 1 ms 4600 KB Correct answer: answer = 151
14 Correct 1 ms 4440 KB Correct answer: answer = 7550
15 Correct 2 ms 4444 KB Correct answer: answer = 7220
16 Correct 1 ms 4444 KB Correct answer: answer = 7550
17 Correct 1 ms 4444 KB Correct answer: answer = 10000
18 Correct 1 ms 4444 KB Correct answer: answer = 10000
19 Correct 1 ms 4444 KB Correct answer: answer = 624
20 Correct 1 ms 4440 KB Correct answer: answer = 10000
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 1
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4444 KB Correct answer: answer = 1
4 Correct 1 ms 4440 KB Correct answer: answer = 5
5 Incorrect 1 ms 4444 KB Wrong answer: output = 40, expected = 41
6 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 4
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4440 KB Correct answer: answer = 4
4 Correct 1 ms 4440 KB Correct answer: answer = 12
5 Correct 1 ms 4444 KB Correct answer: answer = 52
6 Correct 1 ms 4444 KB Correct answer: answer = 210
7 Correct 1 ms 4444 KB Correct answer: answer = 88
8 Correct 1 ms 4444 KB Correct answer: answer = 7696
9 Correct 1 ms 4444 KB Correct answer: answer = 1
10 Correct 1 ms 4444 KB Correct answer: answer = 2374
11 Correct 2 ms 4440 KB Correct answer: answer = 9502
12 Correct 2 ms 4444 KB Correct answer: answer = 49
13 Correct 1 ms 4600 KB Correct answer: answer = 151
14 Correct 1 ms 4440 KB Correct answer: answer = 7550
15 Correct 2 ms 4444 KB Correct answer: answer = 7220
16 Correct 1 ms 4444 KB Correct answer: answer = 7550
17 Correct 1 ms 4444 KB Correct answer: answer = 10000
18 Correct 1 ms 4444 KB Correct answer: answer = 10000
19 Correct 1 ms 4444 KB Correct answer: answer = 624
20 Correct 1 ms 4440 KB Correct answer: answer = 10000
21 Correct 1 ms 4444 KB Correct answer: answer = 1
22 Correct 1 ms 4444 KB Correct answer: answer = 4
23 Correct 1 ms 4444 KB Correct answer: answer = 1
24 Correct 1 ms 4440 KB Correct answer: answer = 5
25 Incorrect 1 ms 4444 KB Wrong answer: output = 40, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 4
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4440 KB Correct answer: answer = 4
4 Correct 1 ms 4440 KB Correct answer: answer = 12
5 Correct 1 ms 4444 KB Correct answer: answer = 52
6 Correct 1 ms 4444 KB Correct answer: answer = 210
7 Correct 1 ms 4444 KB Correct answer: answer = 88
8 Correct 1 ms 4444 KB Correct answer: answer = 7696
9 Correct 1 ms 4444 KB Correct answer: answer = 1
10 Correct 1 ms 4444 KB Correct answer: answer = 2374
11 Correct 2 ms 4440 KB Correct answer: answer = 9502
12 Correct 2 ms 4444 KB Correct answer: answer = 49
13 Correct 1 ms 4600 KB Correct answer: answer = 151
14 Correct 1 ms 4440 KB Correct answer: answer = 7550
15 Correct 2 ms 4444 KB Correct answer: answer = 7220
16 Correct 1 ms 4444 KB Correct answer: answer = 7550
17 Correct 1 ms 4444 KB Correct answer: answer = 10000
18 Correct 1 ms 4444 KB Correct answer: answer = 10000
19 Correct 1 ms 4444 KB Correct answer: answer = 624
20 Correct 1 ms 4440 KB Correct answer: answer = 10000
21 Correct 1 ms 4444 KB Correct answer: answer = 1
22 Correct 1 ms 4444 KB Correct answer: answer = 4
23 Correct 1 ms 4444 KB Correct answer: answer = 1
24 Correct 1 ms 4440 KB Correct answer: answer = 5
25 Incorrect 1 ms 4444 KB Wrong answer: output = 40, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 4
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4440 KB Correct answer: answer = 4
4 Correct 1 ms 4440 KB Correct answer: answer = 12
5 Correct 1 ms 4444 KB Correct answer: answer = 52
6 Correct 1 ms 4444 KB Correct answer: answer = 210
7 Correct 1 ms 4444 KB Correct answer: answer = 88
8 Correct 1 ms 4444 KB Correct answer: answer = 7696
9 Correct 1 ms 4444 KB Correct answer: answer = 1
10 Correct 1 ms 4444 KB Correct answer: answer = 2374
11 Correct 2 ms 4440 KB Correct answer: answer = 9502
12 Correct 2 ms 4444 KB Correct answer: answer = 49
13 Correct 1 ms 4600 KB Correct answer: answer = 151
14 Correct 1 ms 4440 KB Correct answer: answer = 7550
15 Correct 2 ms 4444 KB Correct answer: answer = 7220
16 Correct 1 ms 4444 KB Correct answer: answer = 7550
17 Correct 1 ms 4444 KB Correct answer: answer = 10000
18 Correct 1 ms 4444 KB Correct answer: answer = 10000
19 Correct 1 ms 4444 KB Correct answer: answer = 624
20 Correct 1 ms 4440 KB Correct answer: answer = 10000
21 Correct 1 ms 4444 KB Correct answer: answer = 1
22 Correct 1 ms 4444 KB Correct answer: answer = 4
23 Correct 1 ms 4444 KB Correct answer: answer = 1
24 Correct 1 ms 4440 KB Correct answer: answer = 5
25 Incorrect 1 ms 4444 KB Wrong answer: output = 40, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Correct answer: answer = 4
2 Correct 1 ms 4444 KB Correct answer: answer = 4
3 Correct 1 ms 4440 KB Correct answer: answer = 4
4 Correct 1 ms 4440 KB Correct answer: answer = 12
5 Correct 1 ms 4444 KB Correct answer: answer = 52
6 Correct 1 ms 4444 KB Correct answer: answer = 210
7 Correct 1 ms 4444 KB Correct answer: answer = 88
8 Correct 1 ms 4444 KB Correct answer: answer = 7696
9 Correct 1 ms 4444 KB Correct answer: answer = 1
10 Correct 1 ms 4444 KB Correct answer: answer = 2374
11 Correct 2 ms 4440 KB Correct answer: answer = 9502
12 Correct 2 ms 4444 KB Correct answer: answer = 49
13 Correct 1 ms 4600 KB Correct answer: answer = 151
14 Correct 1 ms 4440 KB Correct answer: answer = 7550
15 Correct 2 ms 4444 KB Correct answer: answer = 7220
16 Correct 1 ms 4444 KB Correct answer: answer = 7550
17 Correct 1 ms 4444 KB Correct answer: answer = 10000
18 Correct 1 ms 4444 KB Correct answer: answer = 10000
19 Correct 1 ms 4444 KB Correct answer: answer = 624
20 Correct 1 ms 4440 KB Correct answer: answer = 10000
21 Correct 1 ms 4444 KB Correct answer: answer = 1
22 Correct 1 ms 4444 KB Correct answer: answer = 4
23 Correct 1 ms 4444 KB Correct answer: answer = 1
24 Correct 1 ms 4440 KB Correct answer: answer = 5
25 Incorrect 1 ms 4444 KB Wrong answer: output = 40, expected = 41
26 Halted 0 ms 0 KB -