#37923 (C++11) No.245 貫け！

• 問題: No.245 貫け！
• ユーザー: airis
• 提出日時: 2015-07-19 20:27:33
• 言語: C++11 (gcc 11.4.0)
• 結果: WA
• 実行時間: -
• コード長: 2,359 bytes
• コンパイル時間: 687 ms
• コンパイル使用メモリ: 83,620 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-07-08 10:23:35
• 合計ジャッジ時間: 1,596 ms
• ジャッジサーバーID （参考情報）: judge4 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	1 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	1 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	1 ms; 5,376 KB
testcase_07	WA	-
testcase_08	AC	2 ms; 5,376 KB
testcase_09	WA	-
testcase_10	WA	-
testcase_11	WA	-
testcase_12	WA	-
testcase_13	WA	-
testcase_14	WA	-
testcase_15	AC	41 ms; 5,376 KB
testcase_16	WA	-
testcase_17	AC	42 ms; 5,376 KB
testcase_18	WA	-
testcase_19	WA	-

ソースコード

#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <numeric>
#include <bitset>
#include <complex>
#define rep(x, to) for (int x = 0; x < (to); x++)
#define REP(x, a, to) for (int x = (a); x < (to); x++)
#define foreach(itr, x) for (typeof((x).begin()) itr = (x).begin(); itr != (x).end(); itr++)
#define EPS (1e-14)

using namespace std;

typedef long long ll;
typedef pair<int, int> PII;
typedef pair<ll, ll> PLL;
typedef complex<double> Complex;
typedef vector< vector<int> > Mat;

struct P {
	double x, y;
	P() : x(0), y(0) {}
	P(double x, double y) : x(x), y(y) {};
	double add(double a, double b) { // 誤差を考慮して加算
		if (abs(a + b) < EPS * (abs(a) + abs(b))) return 0;
		return a + b;
	}
	P operator+(P p) {
		return P(add(x, p.x), add(y, p.y));
	}
	P operator-(P p) {
		return P(add(x, -p.x), add(y, -p.y));
	}
	P operator*(double d) {
		return P(x * d, y * d);
	}
	double dot(P p) { // 内積
		return add(x * p.x, y * p.y);
	}
	double cross(P p) { // 外積
		return add(x * p.y, -y * p.x);
	}
};

bool intersect(P s1, P e1, P s2, P e2) { // 線分交差
	//x座標的に離れて位置している
	if (max(s1.x, e1.x) < min(s2.x, e2.x) || max(s2.x, e2.x) < min(s1.x, e1.x)) return false;
	//y座標的に離れて位置している
	if (max(s1.y, e1.y) < min(s2.y, e2.y) || max(s2.y, e2.y) < min(s1.y, e1.y)) return false;
	//交差している
	double res1 = (e1 - s1).cross(s2 - s1) * (e1 - s1).cross(e2 - s1);
	double res2 = (e2 - s2).cross(e1 - s2) * (e2 - s2).cross(s1 - s2);
	return (res1 < EPS) && (res2 < EPS);
}

bool on_seg(P p1, P p2, P q) { // 線分p1-p2上にある点qがあるか判定
	return (p1 - q).cross(p2 - q) == 0 && (p1 - q).dot(p2 - q) <= 0;
}

int N;
P s[105], e[105];
int ans;

int calc(P p, P q) {
	int res = 0;
	rep(i, N) {
		res += intersect(p, q, s[i], e[i]);
	}
	return res;
}

int main() {

	cin >> N;
	rep(i, N) {
		cin >> s[i].x >> s[i].y >> e[i].x >> e[i].y;
	}

	for (int i = 0; i < N; i++) {
		for (int j = i + 1; j < N; j++) {
			ans = max(ans, calc(s[i], s[j]));
			ans = max(ans, calc(s[i], e[j]));
			ans = max(ans, calc(e[i], s[j]));
			ans = max(ans, calc(e[i], e[j]));
		}
	}

	cout << ans << endl;
	
	return 0;
}