| # | Time | Username | Problem | Language | Result | Execution time | Memory |
|---|---|---|---|---|---|---|---|
| 173692 | red1108 | Aliens (IOI16_aliens) | C++17 | 0 ms | 0 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define fi first
#define se second
#define fastio ios_base::sync_with_stdio(false);cin.tie(0)
#define fopen freopen("input.txt", "r", stdin)
#define pb push_back
#define prec(a) cout<<fixed;cout.precision(a);
#define all(a) (a).begin(), (a).end()
typedef long long ll;
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef tuple<int,int,int> tiii;
const ll INF = 1e18;
const int inf = 2e9;
template<class T>
void pr(T t) {cout << t << " ";}
template<class T, class ...Args>
void pr(T a, Args ...args) {cout << a << " ";pr(args...);}
template<class ...Args>
void prl(Args ...args) {pr(args...);cout << endl;}
ll dp[101010];
vector<pll> tmp, p;
int cur=0;
struct Line{
ll a, b;
int cnt;
Line(){}
Line(ll aa, ll bb, int cc){a=aa;b=bb;cnt=cc;}
};
vector<Line> line;
double cross(Line A, Line B){
return (double)(B.b-A.b)/(double)(A.a-B.a);
}
pll get(ll x){
while(cur+1<line.size()&&cross(line[cur],line[cur+1])<=(double)x) cur++;
//printf("why?? %d\n", cur);
//for(auto i:line) printf("deug (%lld,%lld)\n", i.a,i.b);
if(cur>=line.size()) return pll(0,0);
return pll(line[cur].a*x+line[cur].b,line[cur].cnt);
}
void add(Line nl){
while(line.size()>=2&&cross(line[line.size()-2],line[line.size()-1])>=cross(line[line.size()-2],nl)) line.pop_back();
line.push_back(nl);
}
int get_take(ll c){
//dp초기화
cur=0;
line.clear();
add(Line(-4*(p[1].fi-1),c+2LL*(p[1].fi-1)*(p[1].fi-1),1));
pll t;
for(int i=1;i<(int)p.size();i++){
if(i>1){
if(p[i].fi<=p[i-1].se) add(Line(-4*(p[i].fi-1),c+t.fi+2LL*(p[i].fi-1)*(p[i].fi-1) - 2LL*(p[i-1].se-p[i].fi+1)*(p[i-1].se-p[i].fi+1),t.se+1));
else add(Line(-4*(p[i].fi-1),c+t.fi+2LL*(p[i].fi-1)*(p[i].fi-1),t.se+1));
}
t = get(p[i].se);
t.fi+=2LL*p[i].se*p[i].se;
dp[i]=t.fi;
}
return (int)t.se;
}
ll take_photos(int n, int m, int k, vector<int>x, vector<int>y){
k=min(n,k);
for(int i=0;i<n;i++){
if(x[i]>y[i]) swap(x[i], y[i]);
tmp.pb({x[i],y[i]});
}
sort(tmp.begin(), tmp.end(), [](pii a, pii b){if(a.se==b.se) return a.fi>b.fi;return a.se<b.se;});
p.pb({-1,-1});
for(auto i:tmp){
while(!p.empty()&&p[p.size()-1].fi>=i.fi) p.pop_back();
p.pb(i);
}
k=min(k,(int)p.size()-1);
ll s = 0,e=2*(ll)m*m;
ll t = INF;
while(s<e){
ll mid = (s+e)/2;
if(get_take(mid*2+1)<=k) e = mid;
else s = mid+1;
}
int num = get_take(s*2);
return dp[p.size()-1]/2-s*num;
}
int main(){
int n, m, k;
int inx[100010], iny[100010];
cin>>n>>m>>k;
vector<int> x,y;
for(int i=0;i<n;i++){
int a;
cin>>a;
x.pb(a);
}
for(int i=0;i<n;i++){
int a;
cin>>a;
y.pb(a);
}
printf("%lld\n",take_photos(n,m,k,x,y));
}