1. 程式人生 > >bzoj1071 [SCOI2007]組隊

bzoj1071 [SCOI2007]組隊

== pan open 衡量 inline bits har 大於 long long

Description

  NBA每年都有球員選秀環節。通常用速度和身高兩項數據來衡量一個籃球運動員的基本素質。假如一支球隊裏速度最慢的球員速度為minV,身高最矮的球員高度為minH,那麽這支球隊的所有隊員都應該滿足: A * ( height – minH ) + B * ( speed – minV ) <= C 其中A和B,C為給定的經驗值。這個式子很容易理解,如果一個球隊的球員速度和身高差距太大,會造成配合的不協調。 請問作為球隊管理層的你,在N名選秀球員中,最多能有多少名符合條件的候選球員。

Input

  第一行四個數N、A、B、C 下接N行每行兩個數描述一個球員的height和speed

Output

  最多候選球員數目。

Sample Input

4 1 2 10
5 1
3 2
2 3
2 1

Sample Output

4

HINT

數據範圍: N <= 5000 ,height和speed不大於10000。A、B、C在長整型以內。
2016.3.26 數據加強 Nano_ape 程序未重測

正解:計數+單調性。

我們先把所有人按照高度從小到大排序,然後從大到小枚舉高度。

然後我們把高度$\leq$當前高度的人分別按照速度和$a*h+b*s$排序,這裏可以用歸並排序。

然後枚舉高度,我們發現一個人如果可以在這個隊,那麽他就要滿足$a*h+b*s<=c+a*minh+b*mins$,現在右邊已經是一個常數了。

然後我們可以把滿足式子且高度$\leq minh$的人加進來,同時每次把高度小於$minh$的人減去。

 1 #include <bits/stdc++.h>
 2 #define il inline
 3 #define RG register
 4 #define ll long long
 5 
 6 using namespace std;
 7 
 8 struct data{ int i; ll h,s,v; }q1[5010],q2[5010],st[5010];
 9 
10 int vis[5010],n,ans;
11 ll a,b,c;
12
13 il int gi(){ 14 RG int x=0,q=1; RG char ch=getchar(); 15 while ((ch<0 || ch>9) && ch!=-) ch=getchar(); 16 if (ch==-) q=-1,ch=getchar(); 17 while (ch>=0 && ch<=9) x=x*10+ch-48,ch=getchar(); 18 return q*x; 19 } 20 21 il ll gl(){ 22 RG ll x=0,q=1; RG char ch=getchar(); 23 while ((ch<0 || ch>9) && ch!=-) ch=getchar(); 24 if (ch==-) q=-1,ch=getchar(); 25 while (ch>=0 && ch<=9) x=x*10+ch-48,ch=getchar(); 26 return q*x; 27 } 28 29 il int cmpx(const data &a,const data &b){ return a.h<b.h; } 30 31 int main(){ 32 #ifndef ONLINE_JUDGE 33 freopen("team.in","r",stdin); 34 freopen("team.out","w",stdout); 35 #endif 36 n=gi(),a=gl(),b=gl(),c=gl(); 37 for (RG int i=1;i<=n;++i) q1[i].i=i,q1[i].h=gl(),q1[i].s=gl(),q1[i].v=a*q1[i].h+b*q1[i].s; 38 sort(q1+1,q1+n+1,cmpx); for (RG int i=1;i<=n;++i) q2[i]=q1[i]; 39 for (RG int p=n;p;--p){ 40 RG ll minx=q1[p].h; 41 if(p<n){ 42 RG int t1=p,t2=p+1; 43 for (RG int i=1;i<=n-p+1;++i) 44 if (t2>n || (t2<=n && t1<=p && q1[t1].s<q1[t2].s)) st[i]=q1[t1++]; else st[i]=q1[t2++]; 45 for (RG int i=p;i<=n;++i) q1[i]=st[i-p+1]; t1=p,t2=p+1; 46 for (RG int i=1;i<=n-p+1;++i) 47 if (t2>n || (t2<=n && t1<=p && q2[t1].v<q2[t2].v)) st[i]=q2[t1++]; else st[i]=q2[t2++]; 48 for (RG int i=p;i<=n;++i) q2[i]=st[i-p+1]; 49 } 50 for (RG int i=p,pos=p,res=0;i<=n;++i){ 51 RG ll lim=c+a*minx+b*q1[i].s; 52 while (pos<=n && q2[pos].v<=lim){ 53 if (q2[pos].s>=q1[i].s) ++res,vis[q2[pos].i]=1; ++pos; 54 } 55 ans=max(ans,res); if (vis[q1[i].i]) --res,vis[q1[i].i]=0; 56 } 57 } 58 cout<<ans; return 0; 59 }

bzoj1071 [SCOI2007]組隊