제출 #517078

#제출 시각아이디문제언어결과실행 시간메모리
517078Gurban수열 (APIO14_sequence)C++17
100 / 100
772 ms81536 KiB
#include <bits/stdc++.h> using namespace std; using ll = long long; const int N = 205; const int maxn=1e5+5; int n,k; ll p[maxn]; ll dp[maxn],now[maxn]; int go[N][maxn]; struct line { ll m, c; int idx; long long eval(long long x) { return m * x + c; } long double intersectX(line l) { if(l.m == m){ cout<<"HAHAHAHA\n"; exit(0); } return (long double) (c - l.c) / (l.m - m); } }; int main(){ ios::sync_with_stdio(false); cin.tie(0); cin >> n >> k; for(int i = 1;i <= n;i++){ int x; cin >> x; p[i] = x + p[i - 1]; } for(int i = 1;i <= k;i++){ deque<line>dq; dq.push_front({p[i],dp[i] - p[i] * p[i],i}); for(int j = i + 1;j <= n;j++){ while(dq.size() >= 2 and dq[0].eval(p[j]) <= dq[1].eval(p[j])) dq.pop_front(); now[j] = dq[0].eval(p[j]); go[i][j] = dq[0].idx; line cur = {p[j],dp[j] - p[j] * p[j],j}; while(dq.size() >= 2 and cur.intersectX(dq[(int)dq.size()-2]) <= dq.back().intersectX(dq[(int)dq.size()-2])) dq.pop_back(); dq.push_back(cur); } for(int j = i + 1;j <= n;j++) dp[j] = now[j]; } cout<<dp[n]<<'\n'; int idx = go[k][n]; for(int i = k-1;i >= 0;i--){ cout<<idx<<' '; // assert(idx >= 0 and idx <= n); idx = go[i][idx]; } }
#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...