当前位置: 代码迷 >> 综合 >> codevs 1128 导弹拦截
  详细解决方案

codevs 1128 导弹拦截

热度:44   发布时间:2023-12-15 07:54:59.0

题目:

codevs 1128 导弹拦截

题解:

这题简略看了一下真没啥思路,而且是2010年的普及题目啊!!瞬间感觉智商回到西元前(雾)。然后就想只有两个,每一个导弹只有两种选择,给一号拦截或者二号拦截,枚举?2^n,好像当年给了20分,数据范围如下:


对于10%的数据,N = 1
对于20%的数据,1 ≤ N ≤ 2
对于40%的数据,1 ≤ N ≤ 100
对于70%的数据,1 ≤ N ≤ 1000
对于100%的数据,1 ≤ N ≤ 100000,且所有坐标分量的绝对值都不超过1000。

不明白前两个点是干什么的。然后看到1e5想一想是nlogn的,那就八成是排序,先按离一号拦截点距离从小到大排序,在数组中选定一个一号最多能拦截到的导弹,那么左边的都能拦截到,右边都不能,这时二号拦截系统的半径就是这些导弹中离它最远的那个到它的距离,感觉是n^2的,居然就这样有了70分?!。。。n^3是什么神奇算法。然后发现找最大值的操作是对连续区间,那就上线段树吧,区间最大值。好像常数有点大,加了优读都要跑200+ms。

代码:

//QWsin
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100000+10;
const int INF=(1<<30);int x[maxn],y[maxn],id[maxn];
int n,x1,x2,y1,y2;inline int read()
{int ret=0,ok=1;char ch=getchar();while(ch<'0'||ch>'9') {
   if(ch=='-')ok=-1;ch=getchar();}for(;ch>='0'&&ch<='9';ch=getchar()) ret=ret*10+ch-'0';return ret*ok;
}int f(int x){
   return x*x;}
int dis(int X1,int Y1,int X2,int Y2){
   return f(X1-X2)+f(Y1-Y2);}
int cmp(int a,int b)
{return dis(x[a],y[a],x1,y1) < dis(x[b],y[b],x1,y1);
}struct Node{int l,r,m;Node *lc,*rc;Node(int l=0,int r=0):l(l),r(r){
   m=0;lc=rc=NULL;}
}*root;#define mid ((l+r)>>1)
void UP(Node* &p){p->m=max(p->lc->m,p->rc->m);}
void build(Node* &p,int l,int r)
{p=new Node(l,r);if(l==r){ p->m=dis(x[id[l]],y[id[l]],x2,y2);return ;}build(p->lc,l,mid);build(p->rc,mid+1,r);UP(p);
}int query(Node* p,int L,int R)
{if(L>R) return 0;int l=p->l,r=p->r;  if(L<=l&&r<=R) return p->m;int ret=0;if(L<=mid) ret=max(ret,query(p->lc,L,R));if(mid <R) ret=max(ret,query(p->rc,L,R));return ret;
}int main()
{x1=read();y1=read();x2=read();y2=read();n=read();for(int i=1;i<=n;i++) x[i]=read(),y[i]=read();  for(int i=1;i<=n;i++) id[i]=i;sort(id+1,id+n+1,cmp);build(root,1,n);int ans=INF;for(int i=1;i<=n;i++)ans=min(ans,dis(x[id[i]],y[id[i]],x1,y1)+query(root,i+1,n));cout<<ans<<endl;return 0;
}