제출 #737829

#제출 시각아이디문제언어결과실행 시간메모리
737829amirhoseinfar1385수열 (APIO14_sequence)C++17
60 / 100
126 ms131072 KiB
#include<bits/stdc++.h>
using namespace std;
vector<vector<long long>>dpa,dp;
vector<long long>ps,suf;
int n,k;

void solve(int l,int r,int tl,int tr,int fj){
	//cout<<l<<" "<<r<<" "<<tl<<" "<<tr<<" "<<fj<<endl;
	if(l>r||tl>tr){
		return ;
	}
	if(tl==tr){
		for(int i=l;i<=r;i++){
			dpa[i][fj]=tl;
		}
		return ;
	}
	int mid=(l+r)>>1;
	for(int i=tl;i<=min(tr,mid-1);i++){
		long long fake=dp[i][fj-1]+(ps[mid]-ps[i])*suf[mid+1];
		if(fake>dp[mid][fj]){
			dp[mid][fj]=fake;
			dpa[mid][fj]=i;
		}
	}
	solve(l,mid-1,tl,dpa[mid][fj],fj);
	solve(mid+1,r,dpa[mid][fj],tr,fj);
}

int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>k;
	vector<long long>all(n+1);
	for(int i=1;i<=n;i++){
		cin>>all[i];
	}
	dp.resize(n+3,vector<long long>(k+1));
	dpa.resize(n+3,vector<long long>(k+1));
	ps.resize(n+3);
	suf.resize(n+3);
	for(int i=1;i<=n;i++){
		ps[i]=ps[i-1]+all[i];
	}
	for(int i=n;i>=1;i--){
		suf[i]=all[i]+suf[i+1];
	}
	for(int i=1;i<=n;i++){
		dp[i][1]=ps[i]*suf[i+1];
		dpa[i][1]=-1;
	}
	for(int j=2;j<=k;j++){
		solve(1,n,1,n,j);
		for(int i=1;i<=n;i++){
			dp[i][j]=dp[dpa[i][j]][j-1]+(ps[i]-ps[dpa[i][j]])*suf[i+1];
		}
	}
	long long maxa=-1,w=0;
	for(int i=1;i<=n;i++){
		if(dp[i][k]>maxa){
			w=i;
			maxa=dp[i][k];
		}	
	}
	cout<<maxa<<"\n";
	int now=k;
	vector<int>res;
	while(now>0){
		res.push_back(w);
		w=dpa[w][now];
		now--;
	}
	reverse(res.begin(),res.end());
	res.resize(unique(res.begin(),res.end())-res.begin());
	for(auto x:res){
		cout<<x<<" ";
	}
	cout<<endl;
}
#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...