社区讨论
有没有巨佬能帮忙看一看我这只洛谷AC,bzojWA的代码
P3648[APIO2014] 序列分割参与者 2已保存回复 1
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 1 条
- 当前快照
- 1 份
- 快照标识符
- @mi6nhq2e
- 此快照首次捕获于
- 2025/11/20 07:45 4 个月前
- 此快照最后确认于
- 2025/11/20 07:45 4 个月前
CPP
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=1e5+1;
ll n,k,q[N],sum[N],dp[2][N],x,l,r;//dp[i][j]表示前j个数切i次
int to[205][N];//这玩意儿要开int,否则MLE
bool t;
il ll gi()
{
re ll x=0,t=1;
re char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if(ch=='-') t=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*t;
}
il ll pf(re ll x){return x*x;}
il double xzy(re int j,re int k)
{
//分母有可能为0
if(sum[j]==sum[k]) return -1e18;
return 1.0*(pf(sum[j])-pf(sum[k])-dp[t^1][j]+dp[t^1][k])/(sum[j]-sum[k]);
}
int main()
{
x=n=gi();k=gi();
fp(i,1,n) sum[i]=sum[i-1]+gi();
fp(o,1,k)
{
l=r=0; t^=1;//优先队列搞切点(WC,为什么把t打在i循环内。。。)
fp(i,1,n)
{
while(l<r&&xzy(q[l],q[l+1])<=sum[i]) ++l;//舍
dp[t][i]=dp[t^1][q[l]]+sum[q[l]]*(sum[i]-sum[q[l]]);
to[o][i]=q[l];
while(l<r&&xzy(q[r-1],q[r])>=xzy(q[r],i)) --r;
q[++r]=i;//注意这两只while很套路啊
}
}
printf("%lld\n",dp[t][n]);
fq(i,k,1)
{
x=to[i][x];
printf("%lld ",x);
}
return 0;
}
回复
共 1 条回复,欢迎继续交流。
正在加载回复...