Access Points （单调栈）

题目大意：

平面上$n$个点对序列$a$，再选取$n$个点对序列$b$，使得在满足：$i<j,x_i\le x_j,y_i\le y_j$情况下，${\sum }_{i=1}^n(a[i].x-b[i].x{)}^2+(a[i].y-b[i].y{)}^2$最小

解题思路：

• $x$和$y$可以分开来处理
• 有一点是要清楚的是:${\sum }_{k=l}^r(x-a_k{)}^2$最小值，当$x=\frac{1}{i-j+1}{\sum }_{k=j}^i{a}_k$取到
• 依据此维护一个单调栈，当一个块均值大于栈顶块的均值时，合并这两个块，直到栈里面的均值是单调的

AC代码：

#include <bits/stdc++.h>
#define ft first
#define sd second
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define seteps(N) fixed << setprecision(N)
#define endl "\n"
const int maxn = 1e5 + 10;
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<ll, int> pli;
const ll mod = 1e9 + 7;
int n;
int x[maxn], y[maxn];
db ans;
db sqr(db x) {return x * x;}
stack <pli> sta;
void solve(int *a) {
    for (int i = 1; i <= n; i++) {
        pli res = {a[i], 1};
        while (!sta.empty() && sta.top().ft * res.sd >= sta.top().sd * res.ft) {
            res.ft += sta.top().ft;
            res.sd += sta.top().sd;
            sta.pop();
        }
        sta.push(res);
    }
    int r = n, l;
    while (!sta.empty()) {
        pli res = sta.top(); sta.pop();
        l = r - res.sd + 1;
        for (int i = l; i <= r; i++) 
            ans += sqr(a[i] - 1.0 * res.ft / res.sd);
        r = l - 1;
    }
}
int main() {
    cin >> n;
    for (int i = 1; i <= n; i++) cin >> x[i] >> y[i];
    solve(x), solve(y);
    cout << seteps(12) << ans << endl;
}