1. 程式人生 > >[POI 2014] Couriers

[POI 2014] Couriers

[題目連結]

          https://www.lydsy.com/JudgeOnline/problem.php?id=3524

[演算法]

         首先離線 , 將詢問按右端點排序

         如果我們知道[l , r]這個區間中[L , mid]中的數有多少個和[mid + 1 , R]中的數有多少個 , 則可以通過二分的方式求出答案

         可持久化線段樹可以完成這個任務

         時間複雜度 : O(NlogN)

[程式碼] 

          

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 5e5 + 10;

struct query
{
        int l , r;
        int id;
} q[MAXN];

int n , m , idx;
int a[MAXN] , lson[MAXN << 5
] , rson[MAXN << 5] , sum[MAXN << 5] , root[MAXN] , tmp[MAXN] , ans[MAXN]; template <typename T> inline void chkmax(T &x,T y) { x = max(x,y); } template <typename T> inline void chkmin(T &x,T y) { x = min(x,y); } template <typename T> inline void read(T &x) { T f
= 1; x = 0; char c = getchar(); for (; !isdigit(c); c = getchar()) if (c == '-') f = -f; for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + c - '0'; x *= f; } inline bool cmp(query a , query b) { return a.r < b.r; } inline void build(int &k , int l , int r) { k = ++idx; if (l == r) return; int mid = (l + r) >> 1; build(lson[k] , l , mid); build(rson[k] , mid + 1 , r); } inline void modify(int &k , int old , int l , int r , int pos , int value) { k = ++idx; lson[k] = lson[old] , rson[k] = rson[old]; sum[k] = sum[old] + value; if (l == r) return; int mid = (l + r) >> 1; if (mid >= pos) modify(lson[k] , lson[k] , l , mid , pos , value); else modify(rson[k] , rson[k] , mid + 1 , r , pos , value); } inline int query(int rt1 , int rt2 , int l , int r , int x) { if (l == r) return l; int mid = (l + r) >> 1; if (sum[lson[rt1]] - sum[lson[rt2]] > x) return query(lson[rt1] , lson[rt2] , l , mid , x); else if (sum[rson[rt1]] - sum[rson[rt2]] > x) return query(rson[rt1] , rson[rt2] , mid + 1 , r , x); else return 0; } int main() { read(n); read(m); for (int i = 1; i <= n; i++) { read(a[i]); tmp[i] = a[i]; } sort(tmp + 1 , tmp + n + 1); int len = unique(tmp + 1 , tmp + n + 1) - tmp - 1; for (int i = 1; i <= n; i++) a[i] = lower_bound(tmp + 1 , tmp + len + 1 , a[i]) - tmp; for (int i = 1; i <= m; i++) { read(q[i].l); read(q[i].r); q[i].id = i; } sort(q + 1 , q + m + 1 , cmp); build(root[0] , 1 , len); int now = 1; for (int i = 1; i <= n; i++) { modify(root[i] , root[i - 1] , 1 , n , a[i] , 1); while (now <= m && q[now].r == i) ans[q[now].id] = tmp[query(root[i] , root[q[now].l - 1] , 1 , n , (q[now].r - q[now].l + 1) / 2)] , ++now; } for (int i = 1; i <= m; i++) printf("%d\n" , ans[i]); return 0; }