本文主要是介绍POJ 1269 Intersecting Lines 直线交,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
计算几何入门
无非是平行, 重合和相交
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <map>
#include <sstream>
#include <queue>
#include <vector>
#define MAXN 100005
#define MAXM 211111
#define eps 1e-8
#define INF 500000001
using namespace std;
inline int dblcmp(double d)
{if(fabs(d) < eps) return 0;return d > eps ? 1 : -1;
}
struct point
{double x, y;point(){}point(double _x, double _y): x(_x), y(_y) {}void input(){scanf("%lf%lf", &x, &y);}bool operator ==(point a)const{return dblcmp(a.x - x) == 0 && dblcmp(a.y - y) == 0;}point sub(point p){return point(x - p.x, y - p.y);}double det(point p){return x * p.y - y * p.x;}}p[3333];
struct line
{point a, b;line(){}line(point _a, point _b){ a = _a; b = _b;}void input(){a.input();b.input();}int relation(point p){int c = dblcmp(p.sub(a).det(b.sub(a)));if (c < 0) return 1;if (c > 0) return 2;return 3;}bool parallel(line v){return dblcmp(b.sub(a).det(v.b.sub(v.a))) == 0;}int linecrossline(line v){if ((*this).parallel(v)){return v.relation(a) == 3;}return 2;}point crosspoint(line v){double a1 = v.b.sub(v.a).det(a.sub(v.a));double a2 = v.b.sub(v.a).det(b.sub(v.a));return point((a.x * a2 - b.x * a1) / (a2 - a1), (a.y * a2 - b.y * a1) / (a2 - a1));}}X, Y;
int n;
int main()
{puts("INTERSECTING LINES OUTPUT");scanf("%d", &n);for(int i = 0; i < n; i++){X.input();Y.input();int c = X.linecrossline(Y);if(c == 0) puts("NONE");else if(c == 1) puts("LINE");else{point a = X.crosspoint(Y);printf("POINT %.2f %.2f\n", a.x, a.y);}}puts("END OF OUTPUT");return 0;
}
这篇关于POJ 1269 Intersecting Lines 直线交的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!