洛谷 P1314 [NOIP2011 D2T2] 聰明的質監員
題目描述
小T 是一名質量監督員,最近負責檢驗一批礦產的質量。這批礦產共有 n 個礦石,從 1到n 逐一編號,每個礦石都有自己的重量 wi 以及價值vi 。檢驗礦產的流程是:
1 、給定m 個區間[Li,Ri];
2 、選出一個引數 W;
3 、對於一個區間[Li,Ri],計算礦石在這個區間上的檢驗值Yi:
這批礦產的檢驗結果Y 為各個區間的檢驗值之和。即:Y1+Y2...+Ym
若這批礦產的檢驗結果與所給標準值S 相差太多,就需要再去檢驗另一批礦產。小T
不想費時間去檢驗另一批礦產,所以他想通過調整引數W 的值,讓檢驗結果儘可能的靠近
標準值S,即使得S-Y 的絕對值最小。請你幫忙求出這個最小值。
輸入輸出格式
輸入格式:輸入檔案qc.in 。
第一行包含三個整數n,m,S,分別表示礦石的個數、區間的個數和標準值。
接下來的n 行,每行2個整數,中間用空格隔開,第i+1 行表示 i 號礦石的重量 wi 和價值vi。
接下來的m 行,表示區間,每行2 個整數,中間用空格隔開,第i+n+1 行表示區間[Li,Ri]的兩個端點Li 和Ri。注意:不同區間可能重合或相互重疊。
輸出格式:輸出檔名為qc.out。
輸出只有一行,包含一個整數,表示所求的最小值。
輸入輸出樣例
輸入樣例#1:5 3 15 1 5 2 5 3 5 4 5 5 5 1 5 2 4 3 3
10
說明
【輸入輸出樣例說明】
當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 。
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
(這個虐狗的題號~)
被上午的考試嚇到了,做道正常題壓壓驚~
蠢蠢的我理解公式用了5分鐘......
可以用一個小小的處理降低算檢驗值的時間複雜度:Sumh[i]表示1~i號礦石中重量大於W的數量,Sumv[i]表示1~i號礦石中重量大於W的價值和,則l~r區間的檢驗值計算方法為(sumh[r]-sumh[l-1])*(sumv[r]-sumv[l-1])~[神奇]
一定要注意long long和int的使用!(我因為這個WA了3次......)
#include<cstdio>
#include<cstring>
#define ll long long
int w[200006],n,m;
int x[200006],y[200006],a[200006];
ll k,tot=12425374554373ll,mx,mi=0xfffffff,s[200006],v[200006]; //long long型
ll lmin(ll u,ll v)
{
return u<v ? u:v;
}
int min(int u,int v)
{
return u<v ? u:v;
}
int max(int u,int v)
{
return u>v ? u:v;
}
ll doo(int u)
{
ll ans=0;
memset(a,0,sizeof(a));
memset(s,0,sizeof(s));
for(int i=1;i<=n;i++) //預處理
{
if(w[i]>=u) a[i]=a[i-1]+1,s[i]=s[i-1]+v[i]; //傳說中的字首和~用了這麼久才知道這就是字首和
else a[i]=a[i-1],s[i]=s[i-1];
}
for(int i=1;i<=m;i++) ans+=(a[y[i]]-a[x[i]-1])*(s[y[i]]-s[x[i]-1]);
return ans;
}
int main()
{
scanf("%d%d%lld",&n,&m,&k);
for(int i=1;i<=n;i++)
{
scanf("%d%d",&w[i],&v[i]);
mx=max(w[i],mx);mi=min(w[i],mi);
}
for(int i=1;i<=m;i++)
{
scanf("%d%d",&x[i],&y[i]);
}
int i=mi,j=mx;
while(i<=j)
{
int mid=(i+j)/2;
ll zz=doo(mid);
if(zz>k)
{
i=mid+1;tot=lmin(tot,zz-k); //這裡i和j的順序與一般的二分不同,要反過來
}
else if(zz<k)
{
j=mid-1;tot=lmin(tot,k-zz);
}
else
{
tot=0;break;
}
}
printf("%lld",tot);
return 0;
}