1. 程式人生 > >luoguP1502 視窗的星星(掃描線)20180305 T1

luoguP1502 視窗的星星(掃描線)20180305 T1

題目

傳送門
假定天空是一個平面,每個星星都有一個座標(x,y),每顆星星都有一個亮度C,代表它的亮度。窗戶是長方形的,有固定的長和寬,邊平行於x軸和y軸。你的任務是告訴我如何擺放窗戶,才能獲得在視窗內所有星星的亮度總和最大值。注意,邊框的星星不算。視窗可以被平移,但不允許旋轉。
輸入
在輸入中有幾個測試用例。每一行的第一行包含3個整數:n,w,h,表示星形的數目,矩形視窗的水平長度和垂直高度。然後N行如下,有3個整數分別為:X,Y,C,表示星星位置(x,y)和每個恆星的亮度。沒有兩顆星星在同一點上。
天空中至少有1顆,最多有10000顆星星。1 =

題解

掃描線好題

程式碼

#include<cstdio>
#include<cstring> #include<algorithm> #include<iostream> using namespace std; const int maxn=100001; int cnt,maxnow; int n,w,h,le[maxn],ri[maxn],add[maxn<<2],v[maxn<<2]; struct Star{ int x,y,dis,id; bool f; }a[maxn],b[maxn]; bool cmp1(Star x,Star y) {return x.x<y
.x;} bool cmp2(Star x,Star y) {return x.y<y.y;} void pushdown(int rt) { if (add[rt]) { add[rt<<1]+=add[rt]; add[rt<<1|1]+=add[rt]; v[rt<<1]+=add[rt]; v[rt<<1|1]+=add[rt]; add[rt]=0; } } void change(int L,int R,int C,int l,int
r,int rt) { if (L<=l && r<=R) { add[rt]+=C; v[rt]+=C; return; } int mid=(l+r)>>1; pushdown(rt); if (L<=mid) change(L,R,C,l,mid,rt<<1); if (R>mid) change(L,R,C,mid+1,r,rt<<1|1); v[rt]=max(v[rt<<1],v[rt<<1|1]); } int ques(int L,int R,int l,int r,int rt) { if (L<=l && r<=R) return v[rt]; int mid=(l+r)>>1; pushdown(rt); int ans=0; if (L<=mid) ans=max(ans,ques(L,R,l,mid,rt<<1)); if (R>mid) ans=max(ans,ques(L,R,mid+1,r,rt<<1|1)); return ans; } int main() { freopen("stars.in","r",stdin); freopen("stars.out","w",stdout); scanf("%d%d%d",&n,&w,&h); for (int i=1; i<=n; i++) { scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].dis); a[i].id=i; a[i].f=0; b[i].y=a[i].y; b[i].dis=a[i].dis; b[i].id=i; } int tot=n; for (int i=1; i<=n; i++) { a[++tot].x=a[i].x+w-1; a[tot].id=i; a[tot].f=1; } sort(a+1,a+1+tot,cmp1); for (int i=1; i<=tot; i++) { int tt=0; if (i!=1 && a[i].x==a[i-1].x) tt=cnt; else tt=++cnt; if (!a[i].f) le[a[i].id]=tt; else ri[a[i].id]=tt; } // for (int i=1; i<=n; i++) printf("%d %d\n",le[a[i].id],ri[a[i].id]); tot=n; // for (int i=1; i<=n; i++) printf("%d %d \n",b[i].y,b[i].id); for (int i=1; i<=n; i++) { b[++tot].y=b[i].y+h-1; b[tot].id=i; b[tot].dis=-b[i].dis;// } sort(b+1,b+1+tot,cmp2); // for (int i=1; i<=tot; i++) printf("%d %d %d\n",b[i].y,b[i].id,b[i].dis); // printf("%d ",cnt); return 0; for (int i=1; i<=tot; i++) { if (b[i].dis>0) { int num=ques(le[b[i].id],ri[b[i].id],1,cnt,1); if (num+b[i].dis>maxnow) maxnow=num+b[i].dis; } change(le[b[i].id],ri[b[i].id],b[i].dis,1,cnt,1); } printf("%d\n",maxnow); return 0; }

寫過的題不能寫不對