Submission #90037

#TimeUsernameProblemLanguageResultExecution timeMemory
90037KLPPSplit the sequence (APIO14_sequence)C++14
71 / 100
2041 ms83940 KiB
#include<stdio.h> #include<iostream> #include<deque> #include<queue> using namespace std; typedef long long int lld; #define INF 1000000000000000000 typedef pair<lld,lld> line; int pre[100001][202]; lld sq(lld x){ return x*x; } class CH{ deque<pair<line,int> > d; public: void print(){ queue<pair<line,int> >q; while(!d.empty()){ q.push(d.front()); cout<<d.front().first.first<<" "<<d.front().first.second<<" "<<d.front().second<<endl; d.pop_front(); } while(!q.empty()){ d.push_back(q.front()); q.pop(); } } double intersect(pair<line,int> a, pair<line,int> b){ return (double)(a.first.first-b.first.first)/(b.first.second-a.first.second); } void add(line l, int index){ pair<line,int> P; P.first=l; P.second=index; while(d.size()>0 && P.first.second==d.back().first.second){ if(P.first.first<d.back().first.first){ d.pop_back(); }else{ return; } } if(d.size()<2){ d.push_back(P); return; } pair<line,int> A; pair<line,int> B; do{ A=d.back(); d.pop_back(); B=d.back(); //cout<<intersect(P,A)<<" "<<intersect(P,B)<<endl; }while(d.size()>1 && intersect(P,A)<=intersect(P,B)); if(intersect(P,A)>intersect(P,B))d.push_back(A); d.push_back(P); } int size(){ return d.size(); } int query(lld c){ if(d.size()==1){ return d.front().second; } pair<line,int> A; pair<line,int> B; do{ A=d.front(); d.pop_front(); B=d.front(); }while(d.size()>1 && intersect(A,B)<=c); if(intersect(A,B)>c)d.push_front(A); return d.front().second; } void clear(){ while(!d.empty())d.pop_back(); } }; int main(){ int n,k; scanf("%d %d",&n,&k); lld arr[n]; lld acc[n+1]; acc[0]=0; for(int i=0;i<n;i++){ scanf("%lld",&arr[i]); acc[i+1]=acc[i]+arr[i]; } lld DP[n+1]; for(int i=0;i<=n;i++){ DP[i]=-1; } DP[0]=0; CH *C=new CH(); for(int parts=0;parts<=k;parts++){ for(int i=parts;i<=n;i++){ if(i>parts){ pre[i][parts+1]=C->query(acc[i]); //N[i]=DP[pre[i][parts+1]]+sq(acc[i]-acc[pre[i][parts+1]]); } if(DP[i]!=-1){ C->add(line(sq(acc[i])+DP[i],-2*acc[i]),i); } } C->clear(); for(int i=n;i>-1;i--)DP[i]=DP[pre[i][parts+1]]+sq(acc[i]-acc[pre[i][parts+1]]); DP[parts]=-1; /*for(int i=0;i<=n;i++)cout<<DP[i]<<" "; cout<<endl;*/ } lld ans=sq(acc[n])-DP[n]; ans/=2; printf("%lld\n",ans); int index=n; int partition=k+1; while(partition>0){ index=pre[index][partition]; partition--; if(index>0){ printf("%d ",index); } } printf("\n"); return 0; }

Compilation message (stderr)

sequence.cpp: In function 'int main()':
sequence.cpp:81:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d %d",&n,&k);
  ~~~~~^~~~~~~~~~~~~~~
sequence.cpp:86:12: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
       scanf("%lld",&arr[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...