首页 > 代码库 > 计算几何习题集
计算几何习题集
poj 1269 (直线交点):
这道题是给两条直线(输入直线上的两个点),然后问你两条直线是相交、重合还是平行。
很简单一道题。不过用g++测试的话double要用%f输出。不知道为什么。
(求交点的模板感觉很科学)
1 #include <cstdio> 2 #include <iostream> 3 #define FOR(i,l,r) for(int i=(l);i<=(r);i++) 4 #define FE(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++) 5 #define rep(i,n) for(int i=0;i<(n);i++) 6 #define zero(x) ((x>=0?x:-x)<1e-10) 7 #define debug(x) cout<<#x<<" = "<<x<<endl 8 using namespace std; 9 10 struct pot{11 double x,y;12 pot() { x=0; y=0; }13 pot(double _x, double _y):x(_x),y(_y){};14 };15 double cross(pot a,pot b) { return a.x*b.y-a.y*b.x; }16 pot to_v(pot a,pot b) { pot c=pot(b.x-a.x,b.y-a.y); return c; }17 pot intersection_line(pot u1,pot u2,pot v1,pot v2)18 {19 pot ret = u1;20 double t = ((u1.x-v1.x)*(v1.y-v2.y) - (u1.y-v1.y)*(v1.x-v2.x))21 / ((u1.x-u2.x)*(v1.y-v2.y) - (u1.y-u2.y)*(v1.x-v2.x));22 ret.x += (u2.x-u1.x) * t;23 ret.y += (u2.y-u1.y) * t;24 return ret;25 }26 27 pot a[5];28 int main()29 {30 puts("INTERSECTING LINES OUTPUT");31 int n; scanf("%d",&n);32 for (pot ans;n--;)33 {34 FOR (i,1,4) scanf("%lf%lf",&a[i].x,&a[i].y);35 if (zero(cross(to_v(a[1],a[2]),to_v(a[3],a[4]))))36 if (zero(cross(to_v(a[1],a[2]),to_v(a[1],a[3])))) puts("LINE");37 else puts("NONE");38 else ans=intersection_line(a[1],a[2],a[3],a[4]),printf("POINT %.2f %.2f\n",ans.x,ans.y);39 }40 puts("END OF OUTPUT");41 return 0;42 }43 /*pot intersection_line(pot u1,pot u2,pot v1,pot v2)44 {45 pot ret;46 double a1=u1.y-u2.y,b1=u2.x-u1.x,c1=u1.y*(u1.x-u2.x)-u1.x*(u1.y-u2.y);47 double a2=v1.y-v2.y,b2=v2.x-v1.x,c2=v1.y*(v1.x-v2.x)-v1.x*(v1.y-v2.y);48 ret.y=(a1*c2-a2*c1)/(a2*b1-a1*b2);49 ret.x=(b1*c2-b2*c1)/(a1*b2-a2*b1);50 return ret;51 }*/
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。