社区讨论

有没有巨佬能帮忙看一看我这只洛谷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 条回复,欢迎继续交流。

正在加载回复...