zoj zju 3597 Hit the Target! 线段树 扫描线
这题很锻炼模型转换的能力http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3597
题意:有一排的枪编号依次为1~n 有一排靶子编号依次为1~m
告诉你哪些枪能打中哪些靶子,然后如果每次只能选连续的P把枪,连续的Q个靶子,每次能打中的靶子的最大值为多少
答案是把每次打中的最大值相加再除以总的次数
即选择 1~P 的枪能打中的最多的靶子的数量 + 2~p+1 的枪能打中的最多的靶子的数量 +。。。n-p+1~n的枪能打中的最多的靶子的数量 /(n-p+1)
注意,每把枪最多只能打一个靶子
解法:将题目中的关系转换为坐标 以枪为纵坐标 a 能打中 b ,在坐标系中为(b,a),然后连续的P把枪和连续的Q个靶子则可以表示为
用一个P x Q的矩形去覆盖,最多能覆盖的总的点数,就是poj 2482 了,但要注意一点,如果两根同一高度的水平线的距离小于Q,则重叠的部分职能算一次,因为每把枪只能打一个靶子
我的代码
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define lson l,m,x<<1 #define rson m+1,r,x<<1|1 const int maxn = 100010; int Max[maxn<<2],col[maxn<<2]; struct point{ int x,y; bool operator <(const point &cmp) const { return y<cmp.y || (y==cmp.y&&x<cmp.x); } }pp[maxn]; inline int max(int a,int b){ return a>b?a:b; } inline int min(int a,int b){ return a<b?a:b; } void pushdown(int x){ if(col[x]){ Max[x<<1]+=col[x]; Max[x<<1|1]+=col[x]; col[x<<1]+=col[x]; col[x<<1|1]+=col[x]; col[x]=0; } } void build(int l,int r,int x){ Max[x]=col[x]=0; if(l==r) return ; int m=(l+r)>>1; build(lson); build(rson); } void update(int L,int R,int val,int l,int r,int x){ if(L<=l&&r<=R){ col[x]+=val; Max[x]+=val; return ; } pushdown(x); int m=(l+r)>>1; if(L<=m) update(L,R,val,lson); if(R>m) update(L,R,val,rson); Max[x]=max(Max[x<<1],Max[x<<1|1]); } void solve(int n,int m,int p,int q,int K) { int i,j,y0; double ret=0; build(1,m,1); i=j=0; for(y0=1;y0+p-1<=n;y0++) { for(;j<K && pp[j].y-y0 < p;j++) { if(j+1 < K && pp[j+1].y==pp[j].y && pp[j+1].x-pp[j].x < q) update(pp[j].x,pp[j+1].x-1,1,1,m,1); else update(pp[j].x,min(m,pp[j].x+q-1),1,1,m,1); } ret+=Max[1]; for(;i<j&&pp[i].y==y0;i++) { if(i+1<j && pp[i+1].y==y0 && pp[i+1].x-pp[i].x < q) update(pp[i].x,pp[i+1].x-1,-1,1,m,1); else update(pp[i].x,min(pp[i].x+q-1,m),-1,1,m,1); } } printf("%.2lf\n",ret/(n-p+1)); } int main() { int t,n,m,p,q,k,i; scanf("%d",&t); while(t--) { scanf("%d%d%d%d%d",&n,&m,&p,&q,&k); for(i=0;i<k;i++) scanf("%d%d",&pp[i].y,&pp[i].x); sort(pp,pp+k); solve(n,m,p,q,k); } return 0; }