test/yuki2600.test.cpp

View this file on GitHub
Last update: 2023-04-05 23:10:22+09:00
Problem: https://yukicoder.me/problems/2600

Depends on

Code

#define PROBLEM "https://yukicoder.me/problems/2600"

#include "../template.hpp"

#include "../DataStructure/AbsoluteSum.hpp"


signed main(){
    int n;cin>>n;
    vector<ll> v(n);
    cin>>v;

    sort(ALL(v));

    if(v[0]==v.back()){
        cout<<1<<endl;
        return 0;
    }

    AbsoluteSum<ll> f;

    vector<ll> dp(n,0);
    rep(i,n){
        f.insert(v[i]);
        dp[i]=f.min();
    }

    rep(i,n){
        f.erase(v[i]);
        dp[i]+=f.min();
    }
    cout<<*min_element(ALL(dp))<<endl;
    return 0;
}

#line 1 "test/yuki2600.test.cpp"
#define PROBLEM "https://yukicoder.me/problems/2600"

#line 1 "template.hpp"
#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}

#line 4 "test/yuki2600.test.cpp"

#line 1 "DataStructure/AbsoluteSum.hpp"
template<typename T>
struct AbsoluteSum{
    multiset<T> L,R;
    T min_f;
    
    AbsoluteSum():min_f(0){}

    // add |x-a|
    void insert(T a){
        if(L.empty() or R.empty()){
            L.insert(a);
            R.insert(a);
            return ;
        }
        T pl=*prev(end(L)),pr=*begin(R);
        L.insert(a);
        R.insert(a);
        if(pl<=a and a<=pr) return ;
        // 端をswapするだけでいい
        if(*prev(end(L))>*begin(R)){
            T cl=*prev(end(L)),cr=*begin(R);
            L.erase(L.find(cl));
            R.erase(R.find(cr));
            R.insert(cl);
            L.insert(cr);
        }
        min_f+=std::min(abs(pl-a),abs(pr-a));
        return ;
    }

    // erase |x-a|
    void erase(T a){
        if(L.find(a)!=end(L) and R.find(a)!=end(R)){
            L.erase(L.find(a));
            R.erase(R.find(a));
            return ;
        }
        if(L.find(a)!=end(L)){
            L.erase(L.find(a));
            L.erase(L.find(a));
            L.insert(*begin(R));
            R.erase(begin(R));
        }else{
            R.erase(R.find(a));
            R.erase(R.find(a));
            R.insert(*prev(end(L)));
            L.erase(prev(end(L)));
        }
        T l=*prev(end(L)),r=*begin(R);
        if(l<=a and a<=r) return ;
        min_f-=std::min(abs(l-a),abs(r-a));
        return ;
    }

    T min(){return min_f;}

    T operator()(T x){
        T ret=min_f;
        for(const auto &a:L) ret+=max(a-x,T(0));
        for(const auto &b:R) ret+=max(x-b,T(0));
        return ret;
    }
};
#line 6 "test/yuki2600.test.cpp"


signed main(){
    int n;cin>>n;
    vector<ll> v(n);
    cin>>v;

    sort(ALL(v));

    if(v[0]==v.back()){
        cout<<1<<endl;
        return 0;
    }

    AbsoluteSum<ll> f;

    vector<ll> dp(n,0);
    rep(i,n){
        f.insert(v[i]);
        dp[i]=f.min();
    }

    rep(i,n){
        f.erase(v[i]);
        dp[i]+=f.min();
    }
    cout<<*min_element(ALL(dp))<<endl;
    return 0;
}