Description

小 T 是一名質量監督員，最近負責檢驗一批礦產的質量。這批礦產共有n個礦石，從 1 到n逐一編號，每個礦石都有自己的重量wi以及價值vi。檢驗礦產的流程是：
1. 給定 m個區間[Li,Ri]；
2. 選出一個參數W；
3. 對於一個區間[Li,Ri]，計算礦石在這個區間上的檢驗值Yi：

這批礦產的檢驗結果Y為各個區間的檢驗值之和。即：

若這批礦產的檢驗結果與所給標準值 S 相差太多，就需要再去檢驗另一批礦產。小 T 不想費時間去檢驗另一批礦產，所以他想通過調整參數 W 的值，讓檢驗結果盡可能的靠近標準值 S，即使得S?Y的絕對值最小。請你幫忙求出這個最小值。

Input

第一行包含三個整數n，m，S，分別表示礦石的個數、區間的個數和標準值。
接下來的n 行，每行2 個整數，中間用空格隔開，第i+1 行表示i 號礦石的重量wi 和價值vi 。
接下來的m 行，表示區間，每行2 個整數，中間用空格隔開，第i+n+1 行表示區間[Li,Ri]的兩個端點Li 和Ri。註意：不同區間可能重合或相互重疊。

Output

輸出只有一行，包含一個整數，表示所求的最小值。

Sample Input

5 3 15
1 5
2 5

Sample Output

10

Hint

【輸入輸出樣例說明】
當W 選4 的時候，三個區間上檢驗值分別為20、5、0，這批礦產的檢驗結果為25，此時與標準值S 相差最小為10。
【數據範圍】
對於10%的數據，有1≤n，m≤10；
對於30%的數據，有1≤n，m≤500；
對於50%的數據，有1≤n，m≤5,000；
對於70%的數據，有1≤n，m≤10,000；
對於100%的數據，有1≤n，m≤200,000，0 < wi, vi≤10^6，0 < S≤10^12，1≤Li≤Ri≤n。

思路

閑話

這道題標準解法是用前綴和+二分做，但我並沒有想到前綴和的做法，可能還是我太蒻了

所以，我就傻逼驚奇地想到了可持久化線段樹的做法...

正題

首先看題目，顯然Y是滿足單調性的，而只有在決定了W後才好統計，所以第一眼肯定知道這道題要二分

然後想怎麽統計答案，由於第x次和x的上一次la之間，對答案有貢獻有變化的點的只有滿足一下兩個條件：

　　a.Wla<=wi<Wx

　　b.Wx<=wi<Wla

其中a是對上一次有貢獻，對這一次沒貢獻;而b是對上一次沒貢獻，對這一次有貢獻。

所以我們就可以用可持久化線段樹，每一次對歷史版本進行修改，查詢

時間復雜度O(2nlogn log(maxw))

代碼

  1 #include<set>
  2 #include<map>
  3 #include<stack>
  4 #include<queue>
  5 #include<cstdio>
  6 #include<cstring>
  7 #include<iostream>
  8 #include<algorithm>
  9 #define RG register int
 10 #define rep(i,a,b)    for(RG i=a;i<=b;i++)
 11 #define per(i,a,b)    for(RG i=a;i>=b;i--)
 12 #define ll long long
 13 #define inf (1<<30)
 14 #define maxn 200005
 15 using namespace std;
 16 int n,m,tim,la=inf,id;
 17 ll S,ans;
 18 int w[maxn],v[maxn],L[maxn],R[maxn];
 19 int rt[50];
 20 struct Tree{
 21     int ls,rs,num;ll sum;
 22 }t[4000005];
 23 struct Q{
 24     ll sum;int num;
 25     Q operator + (const Q tmp)const{
 26         return (Q){sum+tmp.sum,num+tmp.num};
 27     }
 28 };
 29 inline ll read()
 30 {
 31     ll x=0,f=1;char c=getchar();
 32     while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1;c=getchar();}
 33     while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘;c=getchar();}
 34     return x*f;
 35 }
 36 
 37 void build(int &pos,int l,int r,int val)
 38 {
 39     pos=++id;
 40     if(l==r){if(w[l]>=val)    t[pos].sum=v[l],t[pos].num=1;return;}
 41     int mid=(l+r)>>1;
 42     build(t[pos].ls,l,mid,val);
 43     build(t[pos].rs,mid+1,r,val);
 44     t[pos].sum=t[t[pos].ls].sum+t[t[pos].rs].sum;
 45     t[pos].num=t[t[pos].ls].num+t[t[pos].rs].num;
 46 }
 47 
 48 void update(int form,int &pos,int l,int r,int aim,int val)
 49 {
 50     pos=++id;
 51     t[pos].ls=t[form].ls,t[pos].rs=t[form].rs;
 52     if(l==r)
 53     {    
 54         t[pos].sum=val;
 55         if(val)t[pos].num=1;else t[pos].num=0;
 56         return;
 57     }
 58     int mid=(l+r)>>1;
 59     if(aim<=mid)    update(t[form].ls,t[pos].ls,l,mid,aim,val);
 60     else            update(t[form].rs,t[pos].rs,mid+1,r,aim,val);
 61     t[pos].sum=t[t[pos].ls].sum+t[t[pos].rs].sum;
 62     t[pos].num=t[t[pos].ls].num+t[t[pos].rs].num;
 63 }
 64 
 65 Q query(int pos,int l,int r,int ql,int qr)
 66 {
 67     if(ql<=l&&r<=qr){return (Q){t[pos].sum,t[pos].num};}
 68     int mid=(l+r)>>1;
 69     Q ret={0,0};
 70     if(ql<=mid)    ret=ret+query(t[pos].ls,l,mid,ql,qr);
 71     if(qr>mid)    ret=ret+query(t[pos].rs,mid+1,r,ql,qr);
 72     return ret;
 73 }
 74 
 75 ll check(int x)
 76 {
 77     ++tim;
 78     bool flg=0;
 79     if(tim==1)    build(rt[1],1,n,x),flg=1;
 80     else
 81     {
 82         rep(i,1,n)
 83         {
 84             if(la<=w[i]&&w[i]<x)//delete
 85             {
 86                 if(!flg)    update(rt[tim-1],rt[tim],1,n,i,0),flg=1;
 87                 else        update(rt[tim],rt[tim],1,n,i,0);
 88             }
 89             else if(x<=w[i]&&w[i]<la)//add    
 90             {
 91                 if(!flg)    update(rt[tim-1],rt[tim],1,n,i,v[i]),flg=1;
 92                 else         update(rt[tim],rt[tim],1,n,i,v[i]);
 93             }
 94         }
 95     }
 96     
 97     if(!flg)    tim--;
 98     ll rans=0;Q ret;
 99     rep(i,1,m)
100         ret=query(rt[tim],1,n,L[i],R[i]),rans+=ret.sum*(ll)ret.num;
101     la=x;
102     return rans;
103 }
104 
105 int main()
106 {
107     int l=1,r=0,mid;
108     ll ans=(1ll<<40);
109     n=read(),m=read(),S=read();
110     rep(i,1,n)    w[i]=read(),v[i]=read(),r=max(r,w[i]);
111     rep(i,1,m)    L[i]=read(),R[i]=read();
112     r++;
113     while(l<=r)
114     {
115         mid=(l+r)>>1;
116         ll tmp=check(mid);
117         ans=min(ans,abs(S-tmp));
118         if(tmp<S)    r=mid-1;
119         else l=mid+1;
120     }
121     cout<<ans;
122     return 0;
123 }

聰明的質檢員 [NOIP 2011] [可持久化線段樹]