1. 程式人生 > >bzoj2216

bzoj2216

++i const std zoj sed one bzoj max hid

決策單調性+整體二分

這裏就是j<k且kj劣於j,j不會再選,所以我們整體二分

pos是因為從L->R中這個是最優點,所以對於mid+1->r選pos之前肯定不優,l->mid-1不會選>pos,因為每個位置都小於mid,並且pos->mid-1這段區間的決策點沒有pos優,因為當前f[i]的i小於mid,選的決策的位置大於pos,由於i小於mid,所以sqrt(i-j),j越大,下降越快,所以pos+1->mid-1肯定沒pos優

技術分享
#include<bits/stdc++.h>
using namespace std;
const int
N = 500010; int n; int a[N], id[N]; double f1[N], f2[N]; double calc(int i, int j) { return (double)a[j] - (double)a[i] + sqrt(abs((double)i - (double)j)); } void solve(int l, int r, int L, int R, double *f) { if(l > r) return; int mid = (l + r) >> 1, lim = min(mid, R), pos = lim;
double mx = 0; for(int i = L; i <= lim; ++i) if(calc(mid, i) > mx) { mx = calc(mid, i); pos = i; } f[id[mid]] = mx; solve(l, mid - 1, L, pos, f); solve(mid + 1, r, pos, R, f); } int main() { scanf("%d", &n); for(int i = 1; i <= n; ++i) { id[i]
= i; scanf("%d", &a[i]); } solve(1, n, 1, n, f1); reverse(a + 1, a + n + 1); reverse(id + 1, id + n + 1); solve(1, n, 1, n, f2); for(int i = 1; i <= n; ++i) printf("%d\n", (int)ceil(max(f1[i], f2[i]))); return 0; }
View Code

bzoj2216