| ||||||||||
Online Judge | Problem Set | Authors | Online Contests | User | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest |
套模板。。。#include <algorithm> #include <iostream> #include <vector> #include <string> #include <queue> #include <cstring> #include <cstdlib> #include <cstdio> #include <cmath> #define eps 1e-8 #define zero(x) (((x)>0?(x):-(x))<eps) #define N 40 using namespace std; struct point { double x, y; } p1, p2, p3, p4; // 判两直线平行 bool Parallel(point u1, point u2, point v1, point v2) { return zero((u1.x-u2.x)*(v1.y-v2.y)-(v1.x-v2.x)*(u1.y-u2.y)); } // 计算cross product (P1-P0)x(P2-P0) double Xmult(point p1, point p2, point p0) { return (p1.x-p0.x)*(p2.y-p0.y) - (p2.x-p0.x)*(p1.y-p0.y); } // 两点距离 double Distance(point p1,point p2) { return sqrt((p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y)); } // 点到直线距离 double Disptoline(point p, point l1, point l2) { return fabs(Xmult(p, l1, l2)) / Distance(l1, l2); } // 计算两直线交点,注意事先判断直线是否平行 point Intersection(point u1, point u2, point v1, point v2) { point ret = u1; double t = ((u1.x-v1.x)*(v1.y-v2.y) - (u1.y-v1.y)*(v1.x-v2.x)) / ((u1.x-u2.x)*(v1.y-v2.y) - (u1.y-u2.y)*(v1.x-v2.x)); ret.x += (u2.x-u1.x) * t; ret.y += (u2.y-u1.y) * t; return ret; } int main() { int n, x1, y1, x2, y2, x3, y3, x4, y4; scanf("%d", &n); printf("INTERSECTING LINES OUTPUT\n"); while(n--) { scanf("%d %d %d %d %d %d %d %d", &x1, &y1, &x2, &y2, &x3, &y3, &x4, &y4); p1.x = x1, p1.y = y1; p2.x = x2, p2.y = y2; p3.x = x3, p3.y = y3; p4.x = x4, p4.y = y4; if(Parallel(p1, p2, p3, p4)) { if(zero(Disptoline(p2, p3, p4))) printf("LINE\n"); else printf("NONE\n"); } else { point p = Intersection(p1, p2, p3, p4); printf("POINT %.2lf %.2lf\n", p.x, p.y); } } printf("END OF OUTPUT\n"); return 0; } Followed by: Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator