結果

問題 No.2012 Largest Triangle
コンテスト
ユーザー 👑 ygussany
提出日時 2022-07-15 23:18:42
言語 C
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 1,493 bytes
コンパイル時間 510 ms
コンパイル使用メモリ 36,924 KB
実行使用メモリ 16,640 KB
最終ジャッジ日時 2024-06-27 20:43:21
合計ジャッジ時間 7,128 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 30 TLE * 1 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("Ofast,unroll-loops")
#include <stdio.h>
#include <math.h>

typedef struct {
	double x, y;
} point;

double distance(point p, point q)
{
	return sqrt((p.x - q.x) * (p.x - q.x) + (p.y - q.y) * (p.y - q.y));
}

double area_of_triangle(point a, point b, point c)
{
	double p, q, r, s;
	p = distance(a, b);
	q = distance(b, c);
	r = distance(c, a);
	s = (p + q + r) / 2.0;
	return sqrtl(s * (s - p) * (s - q) * (s - r));
}

typedef struct {
	double key;
	int id;
} data;

void merge_sort(int n, data x[])
{
	static data y[1000001] = {};
	if (n <= 1) return;
	merge_sort(n / 2, &(x[0]));
	merge_sort((n + 1) / 2, &(x[n/2]));
	
	int i, p, q;
	for (i = 0, p = 0, q = n / 2; i < n; i++) {
		if (p >= n / 2) y[i] = x[q++];
		else if (q >= n) y[i] = x[p++];
		else y[i] = (x[p].key < x[q].key)? x[p++]: x[q++];
	}
	for (i = 0; i < n; i++) x[i] = y[i];
}

void chmax(double* a, double b)
{
	if (*a < b) *a = b;
}

int main()
{
	int i, N;
	point p[200001];
	data d[200001];
	scanf("%d", &N);
	for (i = 1, p[0].x = 0.0, p[0].y = 0.0; i <= N; i++) {
		scanf("%lf %lf", &(p[i].x), &(p[i].y));
		d[i-1].key = distance(p[0], p[i]);
		d[i-1].id = i;
	}
	merge_sort(N, d);
	
	int j;
	double ans = 0.0;
	for (i = N - 1; i >= 0; i--) {
		for (j = i - 1; j >= 0; j--) {
			if (ans >= d[i].key * d[j].key / 2.0) break;
			chmax(&ans, area_of_triangle(p[0], p[d[i].id], p[d[j].id]));
		}
	}
	printf("%lld\n", (long long)(ans * 2.0 + 0.5));
	fflush(stdout);
	return 0;
}
0