Submission #46484

# Submission time Handle Problem Language Result Execution time Memory
46484 2018-04-21T03:42:00 Z RockyB Building Bridges (CEOI17_building) C++17
0 / 100
3000 ms 7348 KB
/// In The Name Of God

#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

#include <bits/stdc++.h>

#define f first
#define s second

#define pb push_back
#define pp pop_back
#define mp make_pair

#define sz(x) (int)x.size()
#define sqr(x) ((x) * 1ll * (x))
#define all(x) x.begin(), x.end()

#define Kazakhstan ios_base :: sync_with_stdio(0), cin.tie(0), cout.tie(0);

#define nl '\n'
#define ioi exit(0);

typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;

const int N = (int)1e5 + 7;
const int K = 317;
const int inf = (int)1e9 + 7;
const int mod = (int)1e9 + 7;
const ll linf = (ll)1e18 + 7;

const int dx[] = {-1, 0, 1, 0, 1, -1, -1, 1};
const int dy[] = {0, 1, 0, -1, 1, -1, 1, -1};

using namespace std;

int n;
int h[N], w[N];
ll dp[N], s[N];

struct line {
	ll k, b;
	ll get(ll x) {
		return k * x + b;
	}
	bool operator < (const line &x) const {
		if (k != x.k) return k < x.k;
		return b < x.b;
	}
};

double inter(line a, line b) {
	return (b.b - a.b) / (a.k - b.k);
}
struct cht {
	vector <line> a, b;
	void build() {
		for (auto it : b) a.pb(it);
		sort(all(a));
		vector <line> nw;
		for (auto it : a) {
			if (sz(nw) && nw.back().k == it.k) continue;
			while (sz(nw) > 1 && inter(it, nw.back()) <= inter(it, nw[sz(nw) - 2])) nw.pp();
			nw.pb(it);
		}
	}
	ll get(ll x) {
		int l = 0, r = sz(a) - 2, res = 0;
		while (l <= r) {
			int mid = l + r >> 1;
			if (a[mid].get(x) >= a[mid + 1].get(x)) res = mid + 1, l = mid + 1;
			else r = mid - 1;
		}
		ll ans = sz(a) ? a[res].get(x) : linf;
		//ll ans = linf;

		for (auto it : b) {
			ans = min(ans, it.get(x));
		}
		return ans;
	}
	void add(line x) {
		b.pb(x);
		if (sz(b) >= K) build();
	}
} t;

int main() {
	#ifdef IOI2018
		freopen ("in.txt", "r", stdin);
	#endif
	Kazakhstan
	cin >> n;
	for (int i = 1; i <= n; i++) {
		cin >> h[i];
	}
	for (int i = 1; i <= n; i++) {
		cin >> w[i];
		s[i] = s[i - 1] + w[i];
	}
	memset(dp, 0x3f, sizeof(dp));
	dp[1] = 0;
	t.add({h[1], sqr(h[1]) + dp[1] - s[1]});
	for (int i = 2; i <= n; i++) {
		dp[i] = t.get(-2 * h[i]) + sqr(h[i]) + s[i - 1];
		t.add({h[i], sqr(h[i]) + dp[i] - s[i]});
		/*ll add = 0;
		for (int j = i - 1; j >= 1; j--) {
			dp[i] = min(dp[i], sqr(h[i] - h[j]) + (s[i - 1] - s[j]) + dp[j]);
			add += w[j];
		}*/
	}
	cout << dp[n];
	ioi
}

Compilation message

building.cpp: In member function 'll cht::get(ll)':
building.cpp:72:16: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
    int mid = l + r >> 1;
              ~~^~~
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1144 KB Output is correct
2 Correct 2 ms 1144 KB Output is correct
3 Correct 2 ms 1328 KB Output is correct
4 Execution timed out 3061 ms 5592 KB Time limit exceeded
5 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 3056 ms 7348 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1144 KB Output is correct
2 Correct 2 ms 1144 KB Output is correct
3 Correct 2 ms 1328 KB Output is correct
4 Execution timed out 3061 ms 5592 KB Time limit exceeded
5 Halted 0 ms 0 KB -