被某大佬指出这是多项式板子！？
我们假设我们原始数列是(a_i, c_i)， 旋转后的数列是(a_i, b_i)，我们的增加量为x
[sum_{i = 1}^n(a_i - b_i + x)^2
]
拆开平方得：
[sum_{i = 1}^na_i^2+b_i^2+x^2+2*x*a_i-2*x*b_i-2*a_i*b_1
]
把这些东西分下类：
[x^2*n+(sum_{i=1}^na_i^2+b_i^2)+2*x*(sum_{i = 1}^n a_i+b_i)+2*(sum_{i = 1}^na_i*b_i)
]
发现(x)只有([-100, 100])，我们考虑枚举x，然后就只有最后一堆是未知的
我们考虑怎么求最后一堆：两个值乘在一起的和，是不是和多项式有关系呢？
但这并不是一个卷积的形式，但我们考虑把b反向，原式就变成了：(sum_{i=1}^na_i*b_{n-i+1})
于是我们就可与愉快的用多项式来做这道题了。不会多项式？戳戳看？
把a数组倍长，把b数组反向，于是这道题的式子就成了：(sum_{i=1}^na_{i+k}c_{n-i+1})，然后多项式的第(n+1-2*n)就分别代表(k)取(0-n-1)的值了
(Code:)
#include<bits/stdc++.h>
using namespace std;
#define il inline
#define re register
#define int long long
const double pi = acos(-1);
#define inf 12345678900000000
il int read() {
re int x = 0, f = 1; re char c = getchar();
while(c < '0' || c > '9') { if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - 48, c = getchar();
return x * f;
}
#define rep(i, s, t) for(re int i = s; i <= t; ++ i)
#define maxn 300005
struct node {
double x, y;
}a[maxn], b[maxn];
il node operator + (node a, node b) { return (node){a.x + b.x, a.y + b.y}; }
il node operator - (node a, node b) { return (node){a.x - b.x, a.y - b.y}; }
il node operator * (node a, node b) { return (node){a.x * b.x - a.y * b.y, a.x * b.y + a.y * b.x}; }
int n, m, lim, r[maxn], ans = inf, sum1, sum2, Ans = -inf;
il void FFT(node *a, int f, int len) {
rep(i, 0, len - 1) if(r[i] > i) swap(a[r[i]], a[i]);
for(re int mid = 1; mid < len; mid <<= 1) {
node base = (node){cos(pi / mid), f * sin(pi / mid)};
for(re int p = mid * 2, j = 0; j < len; j += p) {
node w = (node){1, 0};
for(re int k = 0; k < mid; ++ k, w = w * base) {
node x = a[j + k], y = a[j + k + mid] * w;
a[j + k] = x + y, a[j + k + mid] = x - y;
}
}
}
}
signed main() {
n = read(), m = read();
rep(i, 1, n) a[i].x = a[i + n].x = read();
rep(i, 1, n) b[n - i + 1].x = read();
rep(i, 1, n) sum1 += a[i].x * a[i].x + b[i].x * b[i].x, sum2 += a[i].x - b[i].x;
while((1 << lim) <= 3 * n) ++ lim;
rep(i, 0, (1 << lim)) r[i] = (r[i >> 1] >> 1) | ((i & 1) << (lim - 1));
FFT(a, 1, (1 << lim)), FFT(b, 1, (1 << lim));
rep(i, 0, (1 << lim)) a[i] = a[i] * b[i];
FFT(a, -1, (1 << lim));
rep(x, -m, m) ans = min(ans, x * x * n + sum1 + 2 * x * sum2);
rep(i, n + 1, 2 * n) Ans = max(Ans, (int)(a[i].x / (1 << lim) + 0.5));
printf("%lld", ans - 2 * Ans);
return 0;
}