传送门:Treasure Hunt
题意
有一个左下角为(0,0),右上角为(100,100)的正方形,给出n条线段,线段的起点和终点都在正方形上,给定一个宝藏的坐标,问从正方形外走到宝藏至少要经过多少条线段。每次只能走线段的中点,也就是每两个交点的中点。
题解
因为肯定是从正方形的边开始走,所以枚举每条边上的相邻两个点的中点,连接宝藏所在的点,看这条线段和几条线段有交点,并更新最小值。(为什么这么做是对的呢?因为我没找到反例(逃
存每条边上的点时记得排序 : ) ,忘记排序花式wa。
代码
#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#include<math.h>
#define eps 1e-6
using namespace std;
struct point
{
double x,y;
};
struct line
{
point s,t;
}l[1000];
double a[10010];
double cross(point a,point b)
{
return a.x*b.y-a.y*b.x;
}
point operator - (point a,point b)
{
return {a.x-b.x,a.y-b.y};
}
bool intersection(point a,point b,point c,point d)
{
//快速排斥实验
if(max(c.x,d.x)<min(a.x,b.x)||max(a.x,b.x)<min(c.x,d.x)||max(c.y,d.y)<min(a.y,b.y)||max(a.y,b.y)<min(c.y,d.y)){
return false;
}
//跨立实验
if(cross(a-d,c-d)*cross(b-d,c-d)>0||cross(d-b,a-b)*cross(c-b,a-b)>0){
return false;
}
return true;
}
int dcmp(double x)
{
return fabs(x)<eps?0:x<0?-1:1;
}
int main()
{
int n;
scanf("%d",&n);
for(int i=0;i<n;i++){
scanf("%lf%lf%lf%lf",&l[i].s.x,&l[i].s.y,&l[i].t.x,&l[i].t.y);
}
l[n++]={{0,0},{0,100}};
l[n++]={{0,0},{100,0}};
l[n++]={{100,100},{0,100}};
l[n++]={{100,100},{100,0}};
point q;
scanf("%lf%lf",&q.x,&q.y);
line p;
int ans=0x3f3f3f3f;
int k=0;
for(int i=0;i<n;i++){
if(dcmp(l[i].s.x)==0){
a[k++]=l[i].s.y;
}
if(dcmp(l[i].t.x)==0){
a[k++]=l[i].t.y;
}
}
sort(a,a+k);
for(int i=1;i<k;i++){
p.s={0,(a[i]+a[i-1])/2.0};
p.t=q;
int cnt=0;
for(int j=0;j<n;j++){
if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
}
ans=min(ans,cnt);
}
k=0;
for(int i=0;i<n;i++){
if(dcmp(l[i].s.y)==0){
a[k++]=l[i].s.x;
}
if(dcmp(l[i].t.y)==0){
a[k++]=l[i].t.x;
}
}
sort(a,a+k);
for(int i=1;i<k;i++){
p.s={(a[i]+a[i-1])/2.0,0.0};
p.t=q;
int cnt=0;
for(int j=0;j<n;j++){
if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
}
ans=min(ans,cnt);
}
k=0;
for(int i=0;i<n;i++){
if(dcmp(l[i].s.x-100)==0){
a[k++]=l[i].s.y;
}
if(dcmp(l[i].t.x-100)==0){
a[k++]=l[i].t.y;
}
}
sort(a,a+k);
for(int i=1;i<k;i++){
p.s={100,(a[i]+a[i-1])/2.0};
p.t=q;
int cnt=0;
for(int j=0;j<n;j++){
if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
}
ans=min(ans,cnt);
}
k=0;
for(int i=0;i<n;i++){
if(dcmp(l[i].s.y-100)==0){
a[k++]=l[i].s.x;
}
if(dcmp(l[i].t.y-100)==0){
a[k++]=l[i].t.x;
}
}
sort(a,a+k);
for(int i=1;i<k;i++){
p.s={(a[i]+a[i-1])/2.0,100};
p.t=q;
int cnt=0;
for(int j=0;j<n;j++){
if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
}
ans=min(ans,cnt);
}
if(q.x<0||q.x>100||q.y<0||q.y>100) ans=0;
printf("Number of doors = %d\n",ans);
return 0;
}
京公网安备 11010502036488号