Submission #1110832

#TimeUsernameProblemLanguageResultExecution timeMemory
1110832simona1230Split the sequence (APIO14_sequence)C++17
0 / 100
278 ms131072 KiB
#include <bits/stdc++.h>

using namespace std;

const int maxn=200001;

int n,k;
int a[200001];

int dp[maxn][201];
int p[maxn];

void read()
{
    cin>>n>>k;
    for(int i=1; i<=n; i++)
    {
        cin>>a[i];
        p[i]=p[i-1]+a[i];
    }
}

int idx[maxn][201];
vector<int> ans;

void rec(int i,int c)
{
    if(c==0)return;
    ans.push_back(idx[i][c]);
    rec(idx[i][c],c-1);
}

struct line
{
    int a,b,id;//a*x + b;
    line() {}
    line(int _a,int _b,int _id)
    {
        a=_a;
        b=_b;
        id=_id;
    }

    int v(int x)
    {
        return a*x+b;
    }
};

deque<line> d1,d2;

double pt(line l1,line l2)
{
    return (double)(l2.b-l1.b)/(double)(l1.a-l2.a);
}

void push_line(line l)
{
    while(d2.size()>1&&pt(l,d2[d2.size()-1])<=pt(l,d2[d2.size()-2]))d2.pop_back();
    d2.push_back(l);
}

int curr=0;

pair<int,int> query(int x,int idq)
{
    int ii=0;
    int ans=-1,id=0;
    for(int i=curr; i<d1.size(); i++)
    {
        if(d1[i].id<idq&&ans<=d1[i].v(x))
        {
            ans=d1[i].v(x);
            id=d1[i].id;
            ii=i;
        }
        else break;
    }

    curr=ii;
    return {ans,id};
}

void solve()
{
    for(int i=1; i<=n; i++)
        push_line({p[i],-p[i]*p[i],i});
    d1=d2;
    d2.clear();

    for(int j=1; j<=k; j++)
    {
        curr=0;
        //cout<<"! "<<j<<endl;
        for(int i=1; i<=n; i++)
        {
            pair<int,int> in=query(p[i],i);
            //cout<<i<<" "<<in.first<<" "<<in.second<<endl;
            dp[i][j]=in.first;
            idx[i][j]=in.second;
            push_line({p[i],dp[i][j]-p[i]*p[i],i});
            /*for(int l=0; l<i; l++)
            {
                int curr=dp1[l][j-1]+(p[i]-p[l])*p[l];
                if(dp1[i][j]<=curr)
                {
                    dp1[i][j]=curr;
                    idx[i][j]=l;
                }
            }*/
            //cout<<i<<" "<<dp1[i][j]<<endl;
        }
        d1=d2;
        d2.clear();
    }
    cout<<dp[n][k]<<'\n';
    rec(n,k);
    for(int i=k-1; i>=0; i--)
        cout<<ans[i]<<" ";
    cout<<'\n';
}

int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    cout.tie(NULL);
    read();
    solve();
    return 0;
}
/*
7 3
4 1 3 4 0 2 3

10 1
1 1 0 0 0 0 0 0 0 0
*/

Compilation message (stderr)

sequence.cpp: In function 'std::pair<int, int> query(int, int)':
sequence.cpp:69:22: warning: comparison of integer expressions of different signedness: 'int' and 'std::deque<line>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   69 |     for(int i=curr; i<d1.size(); i++)
      |                     ~^~~~~~~~~~
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...