#include "aliens.h"
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define MASK(i) (1LL << (i))
#define GETBIT(mask, i) (((mask) >> (i)) & 1)
#define ALL(v) (v).begin(), (v).end()
#define block_of_code if(true)
ll max(ll a, ll b){return (a > b) ? a : b;}
ll min(ll a, ll b){return (a < b) ? a : b;}
ll gcd(ll a, ll b){return __gcd(a, b);}
ll lcm(ll a, ll b){return a / gcd(a, b) * b;}
ll LASTBIT(ll mask){return (mask) & (-mask);}
int pop_cnt(ll mask){return __builtin_popcountll(mask);}
int ctz(ull mask){return __builtin_ctzll(mask);}
int logOf(ull mask){return 63 - __builtin_clzll(mask);}
mt19937_64 rng(chrono::high_resolution_clock::now().time_since_epoch().count());
ll rngesus(ll l, ll r){return l + (ull) rng() % (r - l + 1);}
double rngesus_d(double l, double r){
double cur = rngesus(0, MASK(60) - 1);
cur /= MASK(60) - 1;
return l + cur * (r - l);
}
template <class T1, class T2>
bool maximize(T1 &a, T2 b){
if (a < b) {a = b; return true;}
return false;
}
template <class T1, class T2>
bool minimize(T1 &a, T2 b){
if (a > b) {a = b; return true;}
return false;
}
template <class T>
void printArr(T container, string separator = " ", string finish = "\n", ostream &out = cout){
for(auto item: container) out << item << separator;
out << finish;
}
template <class T>
void remove_dup(vector<T> &a){
sort(ALL(a));
a.resize(unique(ALL(a)) - a.begin());
}
const ll INF = 1e18 + 69;
ll sqr(ll x){return x * x;}
vector<pair<int, int>> compress(vector<pair<int, int>> a){
vector<pair<int, int>> tmp = a;
a.clear();
for(pair<int, int> i: tmp){
if (a.empty() || a.back().second != i.second) a.push_back(i);
else minimize(a.back().first, i.first);
}
return a;
}
struct CHT{
#define Node pair<ll, ll>
vector<Node> a;
vector<double> bubble;
bool flag;
// if flag = 0: if inserted in decreasing order, return min, otherwise return max
// flag = 1: the opposite
CHT(bool _flag){
flag = _flag;
}
double getInter(Node a, Node b){
double x = (double) (b.second - a.second) / (a.first - b.first);
return x;
}
void add(Node x, bool isMin = 1){
if (a.empty()) {a.push_back(x); return;}
if (a.back().first == x.first){
if (isMin) minimize(a[a.size() - 1].second, x.second);
else maximize(a[a.size() - 1].second, x.second);
return;
}
while(a.size() >= 2){
double x1 = getInter(a.back(), x), x2 = bubble.back();
if (flag){
if (x1 >= x2) break;
}
else{
if (x1 <= x2) break;
}
a.pop_back(); bubble.pop_back();
}
bubble.push_back(getInter(a.back(), x));
a.push_back(x);
}
ll get(ll x){
if (a.empty()) return INF;
long idx;
if (flag) idx = lower_bound(ALL(bubble),x) - bubble.begin();
else idx = lower_bound(ALL(bubble), x, greater<double>()) - bubble.begin();
return a[idx].first * x + a[idx].second;
}
};
ll take_photos(int n, int m, int k, vector<int> r, vector<int> c) {
vector<pair<int, int>> point;
for(int i = 0; i<n; ++i) {
r[i]++; c[i]++;
if (r[i] > c[i]) swap(r[i], c[i]);
point.push_back({r[i], c[i]});
}
point.push_back({0, 0});
sort(ALL(point), [](pair<int, int> x, pair<int, int> y){return x.second < y.second;});
point = compress(point);
n = point.size() - 1;
vector<int> suff_min(n+2, m+1);
for(int i = n; i>=0; --i)
suff_min[i] = min(point[i].first, suff_min[i+1]);
vector<int> nxt(n+2);
for(int i = 1; i<=n+1; ++i){
if (suff_min[i] != suff_min[i-1]) nxt[i] = i-1;
else nxt[i] = nxt[i-1];
}
vector<ll> dp(n+1, INF);
dp[0] = 0;
ll ans = INF;
while(k--){
vector<ll> cur(n+1, INF);
CHT hull_min(1);
int j = 1;
for(int i = 1; i<=n; ++i) {
while(j <= nxt[i+1]){
if (dp[j-1] != INF)
hull_min.add({-2 * suff_min[j], sqr(suff_min[j]) + dp[j-1]});
j++;
}
int L = suff_min[i+1];
ll lmao = 0;
if (L <= point[i].second) lmao = sqr(point[i].second - L + 1);
cur[i] = hull_min.get(point[i].second + 1);
if (cur[i] != INF){
cur[i] += sqr(point[i].second + 1) - lmao;
}
}
// printArr(cur);
dp = cur;
minimize(ans, dp.back());
}
return ans;
}
// int main(void){
// ios::sync_with_stdio(0);cin.tie(0); cout.tie(0);
// clock_t start = clock();
// int n, m, k; cin >> n >> m >> k;
// vector<int> r(m), c(m);
// for(int i = 0; i<n; ++i) cin >> r[i] >> c[i];
// cout << take_photos(n, m, k, r, c) << "\n";
// cerr << "Time elapsed: " << clock() - start << " ms\n";
// return 0;
// }
