제출 #111339

#제출 시각아이디문제언어결과실행 시간메모리
111339aleksami수열 (APIO14_sequence)C++14
100 / 100
998 ms81700 KiB
#include <bits/stdc++.h> using namespace std; typedef long long ll; #define MAXN 100005 #define MAXK 205 const ll INF = 1e18; //ll dp[MAXK][MAXN]; ll dp[2][MAXN]; ll a[MAXN]; int opt[MAXK][MAXN]; deque <int> q;//deque <line> q; //struct line {ll n;ll k;int id;}; double presekx(int x,int y,int k) { if(a[x]==a[y])return INF; return 1.00*((dp[k%2][x]-a[x]*a[x])-(dp[k%2][y]-a[y]*a[y]))/(a[y]-a[x]); //if(a.k==b.k)return INF; //return 1.00*(a.n-b.n)/(b.k-a.k); } int next_front()//line next_front() { int a = q.front();//line a = q.front(); q.pop_front(); int b = q.front();//line b = q.front(); q.push_front(a); return b; } int prev_back()//line prev_back() { int a = q.back();//line a = q.back(); q.pop_back(); int b = q.back();//line b = q.back(); q.push_back(a); return b; } ll eval(int l,int k,ll x) { return (a[l])*x+(dp[k%2][l]-a[l]*a[l]);//return l.k*x+l.n; } int main() { ios_base::sync_with_stdio(false); int n,k; cin >> n >> k; for(int i = 1; i <= n; i++)cin >> a[i],a[i]+=a[i-1]; for(int i = 1; i <= k; i++) { q.clear(); for(int j = i; j <= n; j++) { while(1) { if(q.size()<=1)break; if(eval(q.front(),i-1,a[j])<=eval(next_front(),i-1,a[j])){q.pop_front();continue;} break; } if(q.size()) { dp[i%2][j]=eval(q.front(),i-1,a[j]); opt[i][j]=q.front();//opt[i][j]=q.front().id; } int add;//line add; //add.k=a[j]; //add.n=-a[j]*a[j]+dp[i-1][j]; add=j;//add.id=j; while(1) { if(q.size()<=1)break; if(a[q.back()]==a[add]){q.pop_back();continue;}//if(q.back().k==add.k){q.pop_back();continue;} if(a[q.back()]==a[prev_back()]){q.pop_back();continue;}//if(q.back().k==prev_back().k){q.pop_back();continue;} if(presekx(prev_back(),q.back(),i-1)>=presekx(q.back(),add,i-1)){q.pop_back();continue;}//if(presekx(prev_back(),q.back())>=presekx(q.back(),add)){q.pop_back();continue;} break; } q.push_back(add); } } cout << dp[k%2][n] << "\n"; int now = n; while(k > 0) { cout << opt[k][now] << " "; now=opt[k][now]; k--; } return 0; }
#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...