//嘘解法 #include using namespace std; const int INF = 1000000000; vector> S; int op(int i, int j){ return S[i][j]; } //https://ei1333.github.io/library/dp/monotone-minima.cpp template< typename T, typename Compare = less< T > > vector< pair< int, T > > monotone_minima(int H, int W, const function< T(int, int) > &f, const Compare &comp = Compare()) { vector< pair< int, T > > dp(H); function< void(int, int, int, int) > dfs = [&](int top, int bottom, int left, int right) { if(top > bottom) return; int line = (top + bottom) / 2; T ma; int mi = -1; for(int i = left; i <= right; i++) { T cst = f(line, i); if(mi == -1 || comp(cst, ma)) { ma = cst; mi = i; } } dp[line] = make_pair(mi, ma); dfs(top, line - 1, left, mi); dfs(line + 1, bottom, mi, right); }; dfs(0, H - 1, 0, W - 1); return dp; } //https://ei1333.github.io/library/dp/divide-and-conquer-optimization.cpp template< typename T, typename Compare = less< T > > vector< vector< T > > divide_and_conquer_optimization(int H, int W, T INF, const function< T(int, int) > &f, const Compare &comp = Compare()) { vector< vector< T > > dp(H + 1, vector< T >(W + 1, INF)); dp[0][0] = 0; for(int i = 1; i <= H; i++) { function< T(int, int) > get_cost = [&](int y, int x) { if(x >= y) return INF; return dp[i - 1][x] + f(x, y); }; auto ret = monotone_minima(W + 1, W + 1, get_cost, comp); for(int j = 0; j <= W; j++) dp[i][j] = ret[j].second; } return dp; } int main(){ int N; cin >> N; vector P(N); for (int i = 0; i < N; i++){ cin >> P[i]; P[i]--; } S = vector>(N + 1, vector(N + 1, 0)); for (int i = 0; i < N; i++){ for (int j = i + 1; j < N; j++){ if (P[i] > P[j]){ S[i][j + 1]++; } } } for (int i = 0; i <= N; i++){ for (int j = 0; j < N; j++){ S[i][j + 1] += S[i][j]; } } for (int i = N; i >= 1; i--){ for (int j = 0; j <= N; j++){ S[i - 1][j] += S[i][j]; } } vector mx(N + 1, 0); for (int i = 0; i <= N; i++){ for (int j = i; j <= N; j++){ mx[j - i] = max(mx[j - i], S[i][j]); S[i][j] *= -1; } } vector> dp = divide_and_conquer_optimization(N + 1, N + 1, INF, op); for (int i = 1; i <= N; i++){ cout << min(-S[0][N] + dp[i][N], -S[0][N] - mx[N + 1 - i]) << endl; } }