Codeforces 321E Ciel and Gondolas

传送门：http://codeforces.com/problemset/problem/321/E

【题解】

首先有一个$O(n^2k)$的dp。

# include <stdio.h>
# include <string.h>
# include <iostream>
# include <algorithm>
// # include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
const int M = 5e5 + , N =  + ;
const int mod = 1e9+;
const ll inf = 1e17;

inline int getint() {
    int x = ; char ch = getchar();
    while(!isdigit(ch)) ch = getchar();
    while(isdigit(ch)) x = (x<<) + (x<<) + ch - '', ch = getchar();
    return x;
}

int n, K, a[N][N], d[N][N];
ll f[][N];

int main() {
    n = getint(), K = getint();
    for (int i=; i<=n; ++i)
        for (int j=; j<=n; ++j)
            a[i][j] = a[i-][j] + a[i][j-] - a[i-][j-] + getint();
    for (int i=; i<=n; ++i) {
        d[i][i] = ;
        for (int j=i+; j<=n; ++j)
            d[i][j] = (a[j][j] - a[j][i-] - a[i-][j] + a[i-][i-])/;
    }
    for (int i=; i<=n; ++i) f[][i] = inf;
    for (int i=; i<=K; ++i) {
        for (int j=; j<=n; ++j) {
            f[i][j] = inf;
            for (int k=; k<j; ++k)
                f[i][j] = min(f[i][j], f[i-][k] + d[k+][j]);
        }
    }
    cout << f[K][n];

    return ;
}

然后发现每层当n向右移动的时候，决策点一定向右移动

（可能可以观察+证明）

然后用整体二分trick就可以了。复杂度$O(Knlogn)$。

# include <stdio.h>
# include <string.h>
# include <iostream>
# include <algorithm>
// # include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
const int M = 5e5 + , N =  + ;
const int mod = 1e9+;
const ll inf = 1e17;

inline int getint() {
    int x = ; char ch = getchar();
    while(!isdigit(ch)) ch = getchar();
    while(isdigit(ch)) x = (x<<) + (x<<) + ch - '', ch = getchar();
    return x;
}

int n, K, a[N][N], d[N][N];
ll f[][N], *F, *G;

inline void solve(int l, int r, int al, int ar) {
    if(l > r) return ;
    int mid = l+r>>; ll mx = inf, t; int pos = ;
    for (int j=al; j<=ar && j<mid; ++j)
        if((t = G[j] + d[j+][mid]) < mx) mx = t, pos = j;
    F[mid] = mx;
    solve(l, mid-, al, pos);
    solve(mid+, r, pos, ar);
}

int main() {
    n = getint(), K = getint();
    for (int i=; i<=n; ++i)
        for (int j=; j<=n; ++j)
            a[i][j] = a[i-][j] + a[i][j-] - a[i-][j-] + getint();
    for (int i=; i<=n; ++i) {
        d[i][i] = ;
        for (int j=i+; j<=n; ++j)
            d[i][j] = (a[j][j] - a[j][i-] - a[i-][j] + a[i-][i-])/;
    }
    for (int i=; i<=n; ++i) f[][i] = inf;
    for (int i=; i<=K; ++i) {
        F = f[i], G = f[i-];
        solve(, n, , n);
    }
    cout << f[K][n];

    return ;
}