洛谷P3195 玩具装箱
第一道斜率优化题。
首先一个基本的状态转移方程是
要使f[i]最小，即b最小。
对于每个j，可以表示为一个点。
然后我们取固定斜率时截距最小的即可，高中线性规划。
单调队列维护下凸包。
然后每次二分出j，转移。
记得给(0，L * L)赋初值。
记得开long long
++，--最好别随便用，编译器的不同会让你爆0...
1 #include <cstdio> 2 3 typedef long long LL; 4 const int N = 50010; 5 6 LL sum[N], g[N], p[N], top; 7 LL f[N], y[N]; 8 9 inline double slope(int i, int j) { 10 return ((double)(y[j] - y[i])) / (g[j] - g[i]); 11 } 12 13 inline int get(int i) { 14 if(i == 1) { 15 return 0; 16 } 17 double k = 2.0 * g[i]; 18 int l = 0, r = top, mid; 19 while(l < r) { 20 mid = (l + r) / 2; 21 //printf("%lf %lf \n", slope(p[mid], p[mid + 1]), k); 22 if(slope(p[mid], p[mid + 1]) < k) { 23 l = mid + 1; 24 } 25 else { 26 r = mid; 27 } 28 } 29 //printf("i = %d r = %d j = %d \n", i, r, p[r]); 30 return p[r]; 31 } 32 33 int main() { 34 //freopen("in.in", "r", stdin); 35 LL n, L; 36 scanf("%lld%lld", &n, &L); 37 L++; 38 for(int i = 1; i <= n; i++) { 39 LL x; 40 scanf("%lld", &x); 41 sum[i] = sum[i - 1] + x; 42 g[i] = i + sum[i]; 43 } 44 y[0] = L * L; 45 for(int i = 1; i <= n; i++) { 46 // f[i] = f[j] + (g[i] - g[j] - L) ^ 2 47 int j = get(i); 48 49 f[i] = f[j] + (g[i] - g[j] - L) * (g[i] - g[j] - L); 50 y[i] = f[i] + (g[i] + L) * (g[i] + L); 51 //printf("y[%d] = %d \n", i, y[i]); 52 53 p[++top] = i; 54 while(top > 1 && slope(p[top - 2], p[top - 1]) >= slope(p[top - 1], p[top])) { 55 p[top - 1] = p[top]; 56 top--; 57 } 58 } 59 60 /*for(int i = 1; i <= n; i++) { 61 printf("%lld ", f[i]); 62 } 63 puts("");*/ 64 printf("%lld", f[n]); 65 return 0; 66 }
[update20181208]今天又考了一次玩具装箱，发现了一个问题.......怎么能把点的坐标直接带入到斜截式里面啊！！！！
只知道y - y0 = k(x - x0)，从来没听过y0 = kx0 + b啊啊啊！！！
关于上面那个的解释:(感谢某蒋姓巨佬为我讲解)
上面那个式子化简为2gi * gj + C = F(j)
考虑有某条直线过点（gj, F(j)），且方程为kx + b = y，其中k = 2gi
那么将点带入，可得：k * gj + b = F(j)
故上面那个等式即为直线的方程。
y - F(j) = 2gi(x - gj)
y - F(j) = 2gi * x - 2gi * gj
然后反正瞎搞一搞就行了啦我也不管了啊啊啊啊阿斜率优化好难啊啊我到底在写什么东西啊
