#include "aliens.h"
#include <algorithm>
#include <deque>
using namespace std;
struct A
{
long long a;
long long b;
}Node[100005],ttt[100005];
bool ok[50005]={0};
long long DP[50005][4005]={0};
long long tt[50005];
deque < pair < long long , long long > > L;
bool cmp(A a,A b)
{
return a.b<b.b||a.b==b.b&&a.a>b.a;
}
long long take_photos(int n, int m, int k, std::vector<int> r, std::vector<int> c)
{
int i,j,N=0,ll,rr;
pair < long long , long long > x,y,z;
for(i=0;i<n;i++)
{
ttt[i].a=(long long) max(r[i],c[i]);
ttt[i].b=(long long) min(r[i],c[i]);
}
sort(Node,ttt+N,cmp);
ll=Node[0].b;
rr=Node[0].a;
for(i=1;i<n;i++)
{
if(ll<=Node[i].b&&Node[i].a<=rr) ok[i]=1;
else
{
ll=Node[i].b;
rr=Node[i].a;
}
}
for(i=0;i<n;i++) if(!ok[i]) Node[N++]=ttt[i];
sort(Node,Node+N,cmp);
Node[N].b=1e9;
tt[N]=1e9;
for(i=N-1;i>=0;i--) tt[i]=min(tt[i+1],Node[i].b);
for(j=0;j<N;j++) DP[j][0]=1e18;
//printf("%lld %lld\n",N,k);
for(i=1;i<=k;i++)
{
L.clear();
for(j=0;j<N;j++)
{
DP[j][i]=1e18;
if(j==0) z=make_pair(-2*Node[j].b,(Node[j].b)*Node[j].b);
else z=make_pair(-2*Node[j].b,DP[j-1][i-1]+(Node[j].b)*Node[j].b);
if(!L.empty())
{
x=L.back();
L.pop_back();
while(!L.empty()&&((double) (z.first-x.first)/(double) (x.second-z.second))>=((double) (L.back().first-x.first)/(double) (x.second-L.back().second)))
{
x=L.back();
L.pop_back();
}
L.push_back(x);
}
if(L.empty()||L.back().first!=z.first) L.push_back(z);
x=L.front();
L.pop_front();
while(!L.empty()&&x.first*(Node[j].a+1)+x.second>=L.front().first*(Node[j].a+1)+L.front().second)
{
x=L.front();
L.pop_front();
}
L.push_front(x);
DP[j][i]=x.first*(Node[j].a+1)+x.second;
DP[j][i]+=(Node[j].a+1)*(Node[j].a+1);
DP[j][i]-=max((long long) 0,Node[j].a-Node[j+1].b+1)*max((long long) 0,Node[j].a-Node[j+1].b+1);
}
}
return DP[N-1][k];
}
