#682935 (C) No.1626 三角形の構築

• 問題: No.1626 三角形の構築
• ユーザー: ygussany
• 提出日時: 2021-07-23 23:14:02
• 言語: C (gcc 12.3.0)
• 結果: AC
• 実行時間: 206 ms / 2,000 ms
• コード長: 1,395 bytes
• コンパイル時間: 1,303 ms
• コンパイル使用メモリ: 32,640 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-10-12 11:48:55
• 合計ジャッジ時間: 3,358 ms
• ジャッジサーバーID （参考情報）: judge4 / judge3

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	137 ms; 5,248 KB
testcase_02	AC	1 ms; 5,248 KB
testcase_03	AC	1 ms; 5,248 KB
testcase_04	AC	1 ms; 5,248 KB
testcase_05	AC	1 ms; 5,248 KB
testcase_06	AC	1 ms; 5,248 KB
testcase_07	AC	93 ms; 5,248 KB
testcase_08	AC	140 ms; 5,248 KB
testcase_09	AC	123 ms; 5,248 KB
testcase_10	AC	161 ms; 5,248 KB
testcase_11	AC	197 ms; 5,248 KB
testcase_12	AC	70 ms; 5,248 KB
testcase_13	AC	66 ms; 5,248 KB
testcase_14	AC	74 ms; 5,248 KB
testcase_15	AC	37 ms; 5,248 KB
testcase_16	AC	78 ms; 5,248 KB
testcase_17	AC	78 ms; 5,248 KB
testcase_18	AC	68 ms; 5,248 KB
testcase_19	AC	47 ms; 5,248 KB
testcase_20	AC	102 ms; 5,248 KB
testcase_21	AC	132 ms; 5,248 KB
testcase_22	AC	206 ms; 5,248 KB
testcase_23	AC	92 ms; 5,248 KB
testcase_24	AC	65 ms; 5,248 KB
testcase_25	AC	1 ms; 5,248 KB
testcase_26	AC	17 ms; 5,248 KB

ソースコード

#include <stdio.h>
#include <math.h>

const long double eps = 1e-7;

long double area_of_triangle(long long p, long long q, long long r)
{
	if (p + q <= r) return 0.0;
	
	long double s;
	s = (p + q + r) / 2.0;
	return sqrtl(s * (s - p) * (s - q) * (s - r));
}

void naive(long long S, long long T)
{
	long long l[3], r[3], m, n;
	l[0] = 1;
	r[0] = T / 3;
	while (l[0] < r[0]) {
		m = (l[0] + r[0]) / 2;
		if (S - area_of_triangle(m, (T - m) / 2, (T - m + 1) / 2) > eps) l[0] = m + 1;
		else r[0] = m;
	}
	l[1] = 1;
	r[1] = T / 3;
	while (l[1] < r[1]) {
		m = (l[1] + r[1] + 1) / 2;
		n = m;
		if (m + n <= T - m - n) n = (T - m * 2) / 2 + 1;
		if (area_of_triangle(m, n, T - m - n) - S > eps) r[1] = m - 1;
		else l[1] = m;
	}
	
	int i, j, k = 0;
	long long ans[10000][3];
	for (i = l[0]; i <= l[1]; i++) {
		l[2] = i;
		r[2] = (T - i) / 2;
		while (l[2] < r[2]) {
			m = (l[2] + r[2]) / 2;
			if (S - area_of_triangle(i, m, T - i - m) > eps) l[2] = m + 1;
			else r[2] = m;
		}
		if (fabsl(area_of_triangle(i, l[2], T - i - l[2]) - S) < eps) {
			ans[k][0] = i;
			ans[k][1] = l[2];
			ans[k++][2] = T - i - l[2];
		}
	}
	printf("%d\n", k);
	for (i = 0; i < k; i++) printf("%lld %lld %lld\n", ans[i][0], ans[i][1], ans[i][2]);
}

int main()
{
	int i, t;
	long long S, T;
	scanf("%d", &t);
	for (i = 1; i <= t; i++) {
		scanf("%lld %lld", &S, &T);
		naive(S, T);
	}
	fflush(stdout);
	return 0;
}