Poj-3304 (叉积判断直线是否过线段)

题目

http://poj.org/problem?id=3304
题意大概就是给一些线段，问能不能找到一条直线使得这些线段在直线上的投影的交不为空。
线段数不大于100。

分析

• 要是投影的交不为空，那么就是至少有一条垂直于那个直线的直线能经过所有线段。
• 而且可以说明，只要有解，那么肯定可以有个解经过某两个端点。（因为直线就是被这些端点限制住的，在边界处肯定是端点）
• 然后就枚举端点弄出一条直线，再判断一下是否经过所有选段即可。

程序

#include <cstdio>#include <cmath>using namespace std;double k1,k2;int n,T,ok;struct node{double x,y;} U[1000],V[1000];bool equ(node A,node B){return abs(A.x-B.x)<=0.00000001 && abs(A.y-B.y)<=0.00000001;}bool Equ(double x,double y){return abs(x-y)<=0.00000001;}double cheng(node v1,node v2){return v1.x*v2.y-v2.x*v1.y;}bool check(node u,node v){    if (equ(u,v)) return 0;    for (int i=1; i<=n; i++){        if (equ(u,U[i]) || equ(u,V[i]) || equ(v,U[i]) || equ(v,V[i])) continue;        k1=cheng((node){U[i].x-u.x,U[i].y-u.y},(node){U[i].x-v.x,U[i].y-v.y});        k2=cheng((node){V[i].x-u.x,V[i].y-u.y},(node){V[i].x-v.x,V[i].y-v.y});        if (k1*k2>0) return 0;    }    return 1;}int main(){    for(scanf("%d",&T); T--; ok=0){        scanf("%d",&n);        for (int i=1; i<=n; i++)            scanf("%lf%lf%lf%lf",&U[i].x,&U[i].y,&V[i].x,&V[i].y);        for (int i=1; i<=n && !ok; i++)            for (int j=i; j<=n && !ok; j++){                ok=check(U[i],U[j])||check(U[i],V[j])||check(V[i],U[j])||check(V[i],V[j]);            }        puts(ok?"Yes!":"No!");    }}

提示

• 我觉得也可以枚举一个端点，假设直线过这个端点，然后利用其它线段来限制住这条直线的斜率，看看最后斜率的可取区间是否为空，可是不知道为什么过不了，代码放下面，也请帮忙看看吧。

#include <cstdio>#include <cmath>#include <algorithm>using namespace std;struct node{double x,y;} U[20000],V[20000];double Max,Min,mn,mx,k1,k2,k3,k4;int n,ok,T;double K(node A,node B){    if (A.x==B.x) return 99999999999999999999999999999.0;    return (A.y-B.y)/(A.x-B.x);}bool equ(node A,node B){return abs(A.x-B.x)<=0.00000001 && abs(A.y-B.y)<=0.00000001;}bool dy(double x,double y){return x-y>0.0001;}bool xy(double x,double y){return y-x>0.0001;}bool check(node O){    Min=-99999999999999999999999999999.0;    Max=99999999999999999999999999999.0;    for (int i=1; i<=n; i++) if(!equ(O,U[i]) && !equ(O,V[i])){        k1=K(O,U[i]);        k2=K(O,V[i]);        mn=min(k1,k2),mx=max(k1,k2);        if (dy(mn,Max) || xy(mx,Min)) return 0;        if (xy(mx,Max)) Max=mx;        if (dy(mn,Min)) Min=mn;        if (dy(Min,Max)) return 0;    }    return 1;}int main(){    for (scanf("%d",&T); T--;){        ok=0;        scanf("%d",&n);        for (int i=1; i<=n; i++) scanf("%lf%lf%lf%lf",&U[i].x,&U[i].y,&V[i].x,&V[i].y);        for (int i=1; i<=n && !ok; i++)            if (check(U[i]) || check(V[i])) ok=1;        puts(ok?"Yes!":"No!");    }}