답안 #879379

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
879379 2023-11-27T08:45:43 Z marvinthang 수열 (APIO14_sequence) C++17
100 / 100
485 ms 88656 KB
/*************************************
*    author: marvinthang             *
*    created: 27.11.2023 15:25:09    *
*************************************/

#include <bits/stdc++.h>

using namespace std;

#define                  fi  first
#define                  se  second
#define                left  ___left
#define               right  ___right
#define                TIME  (1.0 * clock() / CLOCKS_PER_SEC)
#define             MASK(i)  (1LL << (i))
#define           BIT(x, i)  ((x) >> (i) & 1)
#define  __builtin_popcount  __builtin_popcountll
#define              ALL(v)  (v).begin(), (v).end()
#define           REP(i, n)  for (int i = 0, _n = (n); i < _n; ++i)
#define          REPD(i, n)  for (int i = (n); i-- > 0; )
#define        FOR(i, a, b)  for (int i = (a), _b = (b); i < _b; ++i) 
#define       FORD(i, b, a)  for (int i = (b), _a = (a); --i >= _a; ) 
#define       FORE(i, a, b)  for (int i = (a), _b = (b); i <= _b; ++i) 
#define      FORDE(i, b, a)  for (int i = (b), _a = (a); i >= _a; --i) 
#define        scan_op(...)  istream & operator >> (istream &in, __VA_ARGS__ &u)
#define       print_op(...)  ostream & operator << (ostream &out, const __VA_ARGS__ &u)
#ifdef LOCAL
    #include "debug.h"
#else
    #define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
    #define DB(...) 23
    #define db(...) 23
    #define debug(...) 23
#endif

template <class U, class V> scan_op(pair <U, V>)  { return in >> u.first >> u.second; }
template <class T> scan_op(vector <T>)  { for (size_t i = 0; i < u.size(); ++i) in >> u[i]; return in; }
template <class U, class V> print_op(pair <U, V>)  { return out << '(' << u.first << ", " << u.second << ')'; }
template <size_t i, class T> ostream & print_tuple_utils(ostream &out, const T &tup) { if constexpr(i == tuple_size<T>::value) return out << ")";  else return print_tuple_utils<i + 1, T>(out << (i ? ", " : "(") << get<i>(tup), tup); }
template <class ...U> print_op(tuple<U...>) { return print_tuple_utils<0, tuple<U...>>(out, u); }
template <class Con, class = decltype(begin(declval<Con>()))> typename enable_if <!is_same<Con, string>::value, ostream&>::type operator << (ostream &out, const Con &con) { out << '{'; for (__typeof(con.begin()) it = con.begin(); it != con.end(); ++it) out << (it == con.begin() ? "" : ", ") << *it; return out << '}'; }
template <class A, class B> bool minimize(A &a, B b)  { if (a > b) { a = b; return true; } return false; }
template <class A, class B> bool maximize(A &a, B b)  { if (a < b) { a = b; return true; } return false; }

const long long INF = 1e18;

struct Line {
    long long a, b;
    int pos;
    Line(long long a = 0, long long b = INF, int pos = -1): a(a), b(b), pos(pos) {}
    long long eval(long long x) { return a * x + b; }
};
struct ConvexHullTrick {
    vector <Line> lines; int ptr;
    ConvexHullTrick(): ptr(0) {}
    bool bad(Line a, Line b, Line c) { return (long double) (c.b - a.b) / (a.a - c.a) < (long double) (b.b - a.b) / (a.a - b.a); }
    void addLine(long long a, long long b, int c) {
        Line l(a, b, c);
        if (!lines.empty() && lines.back().a == a) {
        	if (lines.back().b <= b) return;
        	lines.pop_back();
        }
        while (lines.size() >= 2 && bad(lines.end()[-2], lines.back(), l)) lines.pop_back();
        lines.push_back(l);
    }
    Line getMin(long long x) {
        if (ptr >= (int) lines.size()) ptr = (int) lines.size() - 1;
        while (ptr < (int) lines.size() - 1 && lines[ptr].eval(x) > lines[ptr + 1].eval(x)) ++ptr;
        return lines[ptr];
    }
};
using CHT = ConvexHullTrick;

// end of template

void process(void) {

	int n, k; cin >> n >> k; ++k;
	vector <int> a(n), pref(n + 1); cin >> a;
	partial_sum(ALL(a), pref.begin() + 1);
	vector dp(2, vector(n + 1, INF));
	vector trace(n + 1, vector(k + 1, 0));
	dp[0][0] = 0;
	CHT cht;
	REP(c, k) {
		cht.lines.clear(); cht.ptr = 0;
		int cur = c & 1;
		FORE(i, c, n) {
			if (i > c) {
				Line l = cht.getMin(pref[i]);
				dp[!cur][i] = l.eval(pref[i]) + 1LL * pref[i] * pref[i];
				trace[i][c + 1] = l.pos;
			}
			if (dp[cur][i] < INF) cht.addLine(-2 * pref[i], dp[cur][i] + 1LL * pref[i] * pref[i], i);
		}
	}
	cout << (1LL * pref[n] * pref[n] - dp[k & 1][n]) / 2 << '\n';
	vector <int> res;
	int i = n;
	FORDE(c, k, 2) {
		i = trace[i][c];
		res.push_back(i);
	}
	reverse(ALL(res));
	for (int x: res) cout << x << ' ';
}

int main(void) {
	ios_base::sync_with_stdio(false); cin.tie(nullptr); // cout.tie(nullptr);
	file("sequence");
	// int t; cin >> t; while (t--)
	process();
	// cerr << "Time elapsed: " << TIME << " s.\n";
	return (0^0);
}

Compilation message

sequence.cpp: In function 'int main()':
sequence.cpp:30:61: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   30 |     #define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
      |                                                      ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
sequence.cpp:110:2: note: in expansion of macro 'file'
  110 |  file("sequence");
      |  ^~~~
sequence.cpp:30:94: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   30 |     #define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
      |                                                                                       ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~
sequence.cpp:110:2: note: in expansion of macro 'file'
  110 |  file("sequence");
      |  ^~~~
# 결과 실행 시간 메모리 Grader output
1 Correct 0 ms 348 KB contestant found the optimal answer: 108 == 108
2 Correct 0 ms 348 KB contestant found the optimal answer: 999 == 999
3 Correct 0 ms 348 KB contestant found the optimal answer: 0 == 0
4 Correct 0 ms 348 KB contestant found the optimal answer: 1542524 == 1542524
5 Correct 0 ms 348 KB contestant found the optimal answer: 4500000000 == 4500000000
6 Correct 0 ms 348 KB contestant found the optimal answer: 1 == 1
7 Correct 0 ms 348 KB contestant found the optimal answer: 1 == 1
8 Correct 0 ms 348 KB contestant found the optimal answer: 1 == 1
9 Correct 0 ms 348 KB contestant found the optimal answer: 100400096 == 100400096
10 Correct 0 ms 348 KB contestant found the optimal answer: 900320000 == 900320000
11 Correct 0 ms 348 KB contestant found the optimal answer: 3698080248 == 3698080248
12 Correct 0 ms 348 KB contestant found the optimal answer: 3200320000 == 3200320000
13 Correct 0 ms 348 KB contestant found the optimal answer: 140072 == 140072
14 Correct 1 ms 344 KB contestant found the optimal answer: 376041456 == 376041456
15 Correct 0 ms 348 KB contestant found the optimal answer: 805 == 805
16 Correct 0 ms 348 KB contestant found the optimal answer: 900189994 == 900189994
17 Correct 0 ms 348 KB contestant found the optimal answer: 999919994 == 999919994
# 결과 실행 시간 메모리 Grader output
1 Correct 0 ms 344 KB contestant found the optimal answer: 1093956 == 1093956
2 Correct 0 ms 348 KB contestant found the optimal answer: 302460000 == 302460000
3 Correct 0 ms 348 KB contestant found the optimal answer: 122453454361 == 122453454361
4 Correct 0 ms 348 KB contestant found the optimal answer: 93663683509 == 93663683509
5 Correct 0 ms 348 KB contestant found the optimal answer: 1005304678 == 1005304678
6 Correct 0 ms 348 KB contestant found the optimal answer: 933702 == 933702
7 Correct 0 ms 348 KB contestant found the optimal answer: 25082842857 == 25082842857
8 Correct 0 ms 348 KB contestant found the optimal answer: 687136 == 687136
9 Correct 0 ms 348 KB contestant found the optimal answer: 27295930079 == 27295930079
10 Correct 0 ms 348 KB contestant found the optimal answer: 29000419931 == 29000419931
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 344 KB contestant found the optimal answer: 610590000 == 610590000
2 Correct 0 ms 348 KB contestant found the optimal answer: 311760000 == 311760000
3 Correct 1 ms 604 KB contestant found the optimal answer: 1989216017013 == 1989216017013
4 Correct 0 ms 348 KB contestant found the optimal answer: 1499437552673 == 1499437552673
5 Correct 1 ms 348 KB contestant found the optimal answer: 1019625819 == 1019625819
6 Correct 1 ms 600 KB contestant found the optimal answer: 107630884 == 107630884
7 Correct 1 ms 604 KB contestant found the optimal answer: 475357671774 == 475357671774
8 Correct 0 ms 348 KB contestant found the optimal answer: 193556962 == 193556962
9 Correct 1 ms 348 KB contestant found the optimal answer: 482389919803 == 482389919803
10 Correct 0 ms 348 KB contestant found the optimal answer: 490686959791 == 490686959791
# 결과 실행 시간 메모리 Grader output
1 Correct 0 ms 348 KB contestant found the optimal answer: 21503404 == 21503404
2 Correct 0 ms 348 KB contestant found the optimal answer: 140412195 == 140412195
3 Correct 5 ms 1116 KB contestant found the optimal answer: 49729674225461 == 49729674225461
4 Correct 1 ms 348 KB contestant found the optimal answer: 37485571387523 == 37485571387523
5 Correct 4 ms 1120 KB contestant found the optimal answer: 679388326 == 679388326
6 Correct 4 ms 1316 KB contestant found the optimal answer: 4699030287 == 4699030287
7 Correct 4 ms 1348 KB contestant found the optimal answer: 12418819758185 == 12418819758185
8 Correct 3 ms 1116 KB contestant found the optimal answer: 31093317350 == 31093317350
9 Correct 1 ms 600 KB contestant found the optimal answer: 12194625429236 == 12194625429236
10 Correct 2 ms 860 KB contestant found the optimal answer: 12345131038664 == 12345131038664
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 1560 KB contestant found the optimal answer: 1818678304 == 1818678304
2 Correct 2 ms 1564 KB contestant found the optimal answer: 1326260195 == 1326260195
3 Correct 31 ms 9332 KB contestant found the optimal answer: 4973126687469639 == 4973126687469639
4 Correct 2 ms 1560 KB contestant found the optimal answer: 3748491676694116 == 3748491676694116
5 Correct 26 ms 5980 KB contestant found the optimal answer: 1085432199 == 1085432199
6 Correct 30 ms 6748 KB contestant found the optimal answer: 514790755404 == 514790755404
7 Correct 32 ms 7452 KB contestant found the optimal answer: 1256105310476641 == 1256105310476641
8 Correct 23 ms 6428 KB contestant found the optimal answer: 3099592898816 == 3099592898816
9 Correct 32 ms 6940 KB contestant found the optimal answer: 1241131419367412 == 1241131419367412
10 Correct 31 ms 8476 KB contestant found the optimal answer: 1243084101967798 == 1243084101967798
# 결과 실행 시간 메모리 Grader output
1 Correct 16 ms 12496 KB contestant found the optimal answer: 19795776960 == 19795776960
2 Correct 17 ms 11472 KB contestant found the optimal answer: 19874432173 == 19874432173
3 Correct 446 ms 88656 KB contestant found the optimal answer: 497313449256899208 == 497313449256899208
4 Correct 17 ms 13008 KB contestant found the optimal answer: 374850090734572421 == 374850090734572421
5 Correct 485 ms 88352 KB contestant found the optimal answer: 36183271951 == 36183271951
6 Correct 326 ms 66448 KB contestant found the optimal answer: 51629847150471 == 51629847150471
7 Correct 378 ms 71100 KB contestant found the optimal answer: 124074747024496432 == 124074747024496432
8 Correct 333 ms 60116 KB contestant found the optimal answer: 309959349080800 == 309959349080800
9 Correct 312 ms 65984 KB contestant found the optimal answer: 124113525649823701 == 124113525649823701
10 Correct 408 ms 80988 KB contestant found the optimal answer: 124309619349406845 == 124309619349406845