斜率dp

要点：

1. 列出dp关系式，找出哪些量是变量，哪些是常量，换成y=kx+b的形式
2. 分析凸包是上凸还是下凸
3. 初始化一定思考全面

题目1：P3195 [HNOI2008]玩具装箱

分析：

参考博客
有一点注意，一开始其实有元素0在队列中，要注意其代表的意义，并进行初始化

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//typedef __int128 lll;
#define print(i) cout << "debug: " << i << endl
#define close() ios::sync_with_stdio(0), cin.tie(0), cout.tie(0)
#define mem(a, b) memset(a, b, sizeof(a))
const ll mod = 1e9 + 7;
const int maxn = 5e4 + 10;
const int inf = 0x3f3f3f3f;

ll dp[maxn], v[maxn], sum[maxn], a[maxn], b[maxn];
ll n, t;
ll q[maxn]; 

ll dy(int i, int j)
{
    return dp[j] + b[j] * b[j] - dp[i] - b[i] * b[i];
}

ll dx(int i, int j)
{
    return b[j] - b[i];
}

int main()
{
    cin >> n >> t;
    for(int i = 1; i <= n; i++)
    {
        cin >> v[i];
        sum[i] = sum[i - 1] + v[i];
        a[i] = sum[i] + i, b[i] = sum[i] + i + 1 + t;
    }
    a[0] = 0, b[0] = 1 + t;
    int l, r;
    l = r = 0;
    for(int i = 1; i <= n; i++)
    {
        while(l < r && dy(q[l], q[l + 1]) < dx(q[l], q[l + 1]) * 2 * a[i])
            l++;
        dp[i] = dp[q[l]] + (a[i] - b[q[l]]) * (a[i] - b[q[l]]);
        while(l < r && dy(q[r - 1], q[r]) * dx(q[r], i) > dy(q[r], i) * dx(q[r - 1], q[r]))
            r--;
        q[++r] = i;
    }
    //for(int i = 1; i <= n; i++) print(dp[i]);
    cout << dp[n] << endl;
}

题目2：Lawrence HDU - 2829

分析：

这个题是也是斜率优化dp，只不过常量除了i还有j。在这里着重讲一下单调队列在这里的含义：
根据dp状态方程可以得出，dp[i][j]跟dp[k][j-1]的状态有关，因此队列中维护的值是相对于j-1来说的，这样第二层循环之前先把j入队就显得很合理了。
个人认为斜率dp的入队初始化是重中之重，方程谁都能推出来，但是推出方程之后能否细心的分析清楚如何初始化才是最关键的。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//typedef __int128 lll;
#define print(i) cout << "debug: " << i << endl
#define close() ios::sync_with_stdio(0), cin.tie(0), cout.tie(0)
#define mem(a, b) memset(a, b, sizeof(a))
const ll mod = 1e9 + 7;
const int maxn = 1010;
const int inf = 0x3f3f3f3f;
ll dp[maxn][maxn], sum[maxn], s[maxn], v[maxn];
int q[maxn];
int n, m, l, r;

ll dx(int k1, int k2)
{
    return sum[k2] - sum[k1];
}

ll dy(int k1, int k2, int j)
{
    return dp[k2][j - 1] - s[k2] + sum[k2] * sum[k2] - (dp[k1][j - 1] - s[k1] + sum[k1] * sum[k1]);
}

int main()
{
    while(cin >> n >> m && (n || m))
    {
        for(int i = 1; i <= n; i++)
        {
            cin >> v[i];
            sum[i] = sum[i - 1] + v[i];
            s[i] = s[i - 1] + v[i] * sum[i - 1];
        }
        for(int i = 1; i <= n; i++)
            dp[i][0] = s[i];
        for(int j = 1; j <= m; j++)
        {
            q[0] = j;
            l = r = 0;
            for(int i = j + 1; i <= n; i++)
            {
                while(l < r && dy(q[l], q[l + 1], j) < dx(q[l], q[l + 1]) * sum[i])
                    l++;
                dp[i][j] = dp[q[l]][j - 1] + s[i] - s[q[l]] - sum[q[l]] * (sum[i] - sum[q[l]]);
                while(l < r && dy(q[r - 1], q[r], j) * dx(q[r], i) > dx(q[r - 1], q[r]) * dy(q[r], i, j))
                    r--;
                q[++r] = i;
            }
        }
        cout << dp[n][m] << endl;
    }
}