定义状态 写转移方程 作差得不等式 码代码
不要用double类型的函数搞斜率 我直接T了 而且还可能出现精度问题
这道题被坑的原因是 :
忘记预处理
一直坑的原因是 网上有一份神犇题解 我确实不懂他的队尾维护 其实一般的维护方法就可以了= = 如果你们想研究就拿去看看吧
http://blog.sina.com.cn/s/blog_6e63f59e0101aeaf.html
怎么想都不明白他用 i和队尾的斜率 以及 i和队尾前一个斜率做比较的原因= = 而且两种都AC了 神奇的九寨
#include<cstdio>
#include<algorithm>
#include<cstring>
#define P(a) ((a)*(a))
using namespace std;
typedef int LL;
const int MAXN = 10000;
const int MAXM = 5000;
const int INF = 99999999;
LL d[MAXM+10][MAXN+10], Q[MAXN+10], F, R, A[MAXN+10];
int n, m;
int Slope_U(int i, int j, int p)
{
return d[p][i] + A[i+1] * A[i+1] - (d[p][j] + A[j+1] * A[j+1]);
}
int Slope_D(int i, int j, int p)
{
return 2 * (A[i+1] - A[j+1]);
}
int main()
{
int kase = 0, Case;
scanf("%d", &Case);
while(Case--)
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++) scanf("%d", &A[i]);
sort(A+1, A+1+n);
for(int i = 1; i <= n; i++) d[0][i] = INF;
for(int j = 1; j <= m; j++)
{
F = R = 0;
Q[R++] = 0;
for(int i = 1; i <= n; i++)
{
while(F < R - 1 && Slope_U(Q[F+1], Q[F], j-1) <= A[i] * Slope_D(Q[F+1], Q[F], j-1))
F++;
d[j][i] = d[j-1][Q[F]] + (A[i] - A[Q[F]+1]) * (A[i] - A[Q[F]+1]);
while(F < R - 1 && Slope_U(i, Q[R-1], j-1) * Slope_D(Q[R-1], Q[R-2], j-1) <= Slope_U(Q[R-1], Q[R-2], j-1) * Slope_D(i, Q[R-1], j-1))
R--;
Q[R++] = i;
}
}
printf("Case %d: %d\n", ++kase, d[m][n]);
}
}