跟POJ2318差不多,不过这次要对线进行排序,对线的排序有涉及到对点的排序
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxl 100010
#define eps 1e-8
using namespace std;
int n,m;;
double L,U,R,D;
int ans[maxl],num[maxl];
struct point
{
double x,y;
point (double a=0,double b=0)
{
x=a;y=b;
}
point operator - (const point &b) const
{
return point(x-b.x,y-b.y);
}
double operator * (const point &b) const
{
return x*b.x+y*b.y;
}
double operator ^ (const point &b) const
{
return x*b.y-y*b.x;
}
bool operator < (const point &b) const
{
return x<b.x-eps;
}
bool operator == (const point &b) const
{
return x>b.x-eps && x<b.x+eps && y>b.y-eps && y<b.y+eps;
}
};
struct line
{
point s,e;
double k;
line(point a=point(),point b=point())
{
s=a;e=b;
k=atan2(e.y-s.y,e.x-s.x);
}
bool operator < (const line &b) const
{
return s < b.s;
}
bool operator == (const line &b) const
{
return s==b.s;
}
}a[maxl];
inline void prework()
{
scanf("%d%lf%lf%lf%lf",&m,&L,&U,&R,&D);
double x1,x2;
point p1,p2;
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&x1,&x2);
p1=point(x1,U);
p2=point(x2,D);
a[i]=line(p1,p2);
}
sort(a+1,a+1+n);
for(int i=0;i<=n;i++)
ans[i]=0;
for(int i=0;i<=m;i++)
num[i]=0;
}
inline bool isleft(point p,int id)
{
point p1=a[id].s-p;
point p2=a[id].e-p;
if((p1^p2) < 0+eps)
return true;
else
return false;
}
inline void mainwork()
{
point p;int l,r,mid;
for(int i=1;i<=m;i++)
{
scanf("%lf%lf",&p.x,&p.y);
l=1,r=n;
while(l+1<r)
{
mid=(l+r)>>1;
if(isleft(p,mid))
r=mid;
else
l=mid;
}
if(isleft(p,l))
ans[l-1]++;
else if(isleft(p,r))
ans[r-1]++;
else
ans[r]++;
}
for(int i=0;i<=n;i++)
num[ans[i]]++;
}
inline void print()
{
puts("Box");
for(int i=1;i<=m;i++)
if(num[i]>0)
printf("%d: %d\n",i,num[i]);
}
int main()
{
while(~scanf("%d",&n) && n)
{
prework();
mainwork();
print();
}
return 0;
}
2019.2.11模板更新:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxl 100010
#define eps 1e-8
using namespace std;
int n,m;;
double L,U,R,D;
int ans[maxl],num[maxl];
inline int sgn(double x)
{
if(fabs(x)<eps) return 0;
if(x<0)
return -1;
else
return 1;
}
struct point
{
double x,y;
point(double a=0,double b=0)
{
x=a;y=b;
}
point operator - (const point &b)const
{
return point(x-b.x,y-b.y);
}
//按照x坐标排序
bool operator < (const point &b)const
{
return x<b.x-eps;
}
bool operator == (const point &b)const
{
return sgn(x-b.x)==0 && sgn(y-b.y)==0;
}
friend point operator * (const point &a,const double &b) //数乘
{
return point(b*a.x,b*a.y);
}
friend point operator * (const double a,const point &b)
{
return point(a*b.x,a*b.y);
}
void transxy(double sinb,double cosb)
{ //逆时针旋转,给出sin,cos
double tx=x,ty=y;
x=tx*cosb-ty*sinb;
y=tx*sinb+ty*cosb;
}
void transxy(double b)//逆时针旋转b弧度
{
double tx=x,ty=y;
x=tx*cos(b)-ty*sin(b);
y=tx*sin(b)+ty*cos(b);
}
double norm()
{
return sqrt(x*x+y*y);
}
};
inline double det(const point &a,const point &b)
{ //叉积
return a.x*b.y-a.y*b.x;
}
inline double dot(const point &a,const point &b)
{ //点积
return a.x*b.x+a.y*b.y;
}
inline double dist(const point &a,const point &b)
{
return (a-b).norm();
}
struct line
{
point s,e;
double k;
line(point a=point(),point b=point())
{
s=a;e=b;
k=atan2(e.y-s.y,e.x-s.x);
}
bool operator < (const line &b) const
{
return s < b.s;
}
bool operator == (const line &b) const
{
return s==b.s;
}
}a[maxl];
inline void prework()
{
scanf("%d%lf%lf%lf%lf",&m,&L,&U,&R,&D);
double x1,x2;
point p1,p2;
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&x1,&x2);
p1=point(x1,U);
p2=point(x2,D);
a[i]=line(p1,p2);
}
sort(a+1,a+1+n);
for(int i=0;i<=n;i++)
ans[i]=0;
for(int i=0;i<=m;i++)
num[i]=0;
}
inline bool isleft(point p,int id)
{
point p1=a[id].s-p;
point p2=a[id].e-p;
if(det(p1,p2) < 0+eps)
return true;
else
return false;
}
inline void mainwork()
{
point p;int l,r,mid;
for(int i=1;i<=m;i++)
{
scanf("%lf%lf",&p.x,&p.y);
l=1,r=n;
while(l+1<r)
{
mid=(l+r)>>1;
if(isleft(p,mid))
r=mid;
else
l=mid;
}
if(isleft(p,l))
ans[l-1]++;
else if(isleft(p,r))
ans[r-1]++;
else
ans[r]++;
}
for(int i=0;i<=n;i++)
num[ans[i]]++;
}
inline void print()
{
puts("Box");
for(int i=1;i<=m;i++)
if(num[i]>0)
printf("%d: %d\n",i,num[i]);
}
int main()
{
while(~scanf("%d",&n) && n)
{
prework();
mainwork();
print();
}
return 0;
}