Codeforces 47E Cannon 询问离线,单调性尺取

题面及翻译

洛谷题面

题解

史上最水E题,用到的算法都是普及组算法.
题目隐含条件:所有炮弹速度相同,角度小于$\frac{\pi }{4}$.
可以发现,当炮弹速度相同并且角度小于$\frac{\pi }{4}$时,显然角度越大,炮弹发射的越高越远.
这样我们找到了单调性.
这样把炮弹角度从小到大排序,墙也按$x$轴从小到大排.
随着炮弹发射的越高越远,炮弹的落点的横坐标也不断变大.
按顺序去枚举炮弹,明显上一个炮弹能够过的墙,这一个炮弹也能通过.
把公式变形对炮弹打到墙上和碰到地上进行分类讨论,最后输出结果.

#include<bits/stdc++.h> //Ithea Myse Valgulious
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rec register char
#define rel register ll
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
  int x=0,f=1;char c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  return f?x:-x;
  }
template <typename mitsuha>
inline bool read(mitsuha &x){
  x=0;int f=1;char c=gc();
  for (;!isdigit(c)&&~c;c=gc()) f^=c=='-';
  if (!~c) return 0;
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  return x=f?x:-x,1;
  }
template <typename mitsuha>
inline int write(mitsuha x){
  if (!x) return 0&pc(48);
  if (x<0) x=-x,pc('-');
  int bit[20],i,p=0;
  for (;x;x/=10) bit[++p]=x%10;
  for (i=p;i;--i) pc(bit[i]+48);
  return 0;
  }
inline char fuhao(){
  char c=gc();
  for (;isspace(c);c=gc());
  return c;
  }
}using namespace chtholly;
using namespace std;
typedef double db;
const int yuzu=1e5;
const db g=9.8;

namespace gongshi{
db vxt(db v,db alp){return v*cos(alp);}
db vyt(db v,db alp,db t){return v*sin(alp)-g*t;}
db xt(db v,db alp,db t){return v*cos(alp)*t;}
db yt(db v,db alp,db t){return v*sin(alp)*t-g*t*t/2;}
db _tx(db x,db alp,db v){return x/v/cos(alp);}
/*_tx函数用于在知道横向距离的情况下求出运行的时间.*/
db _ty(db v,db alp){return 2*v*sin(alp)/g;}
/*_ty函数用于求出自然状况下炮弹落到地上运行的时间.*/
}using namespace gongshi;

struct qiang{
db x,y;
void rd(){scanf("%lf%lf",&x,&y);}
bool operator <(const qiang &b) const{
  return x!=b.x?x<b.x:y<b.y;
  }
}w[yuzu|10];

typedef pair<db,int> node;
node a[yuzu|10];

pair<db,db> llx[yuzu|10];

int main(){
int i,n=read();
db alp,v=read();
for (i=1;i<=n;++i){
  scanf("%lf",&alp);
  a[i]=node(alp,i);
  } 
sort(a+1,a+n+1);
int m=read(),r=1;
for (i=1;i<=m;++i) w[i].rd();
sort(w+1,w+m+1);
#define suan yt(v,a[i].first,_tx(w[r].x,a[i].first,v)) //炮弹在发射到某面墙的时候的高度.
for (i=1;i<=n;++i){
  for (;r<=m&&suan>w[r].y;++r);//算炮弹在第r面墙时的高度,不会被拦下来找下一面墙. 
  llx[a[i].second]=suan<=0||r>m?make_pair(xt(v,a[i].first,_ty(v,a[i].first)),0.):make_pair(w[r].x,suan); 
  /*如果落到负数高度,说明在到墙之前就落地了.*/
  }
for (i=1;i<=n;++i) printf("%.9lf %.9lf\n",llx[i].first,llx[i].second);
}

谢谢大家.