結果

問題 No.199 星を描こう
コンテスト
ユーザー asugen0402
提出日時 2019-06-06 16:54:51
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 1 ms / 2,000 ms
コード長 4,034 bytes
コンパイル時間 911 ms
コンパイル使用メモリ 32,768 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-25 11:31:48
合計ジャッジ時間 2,078 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <float.h>
#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

// 内部定数
#define D_POINT_CNT		5										// 座標数

// 内部構造体 - 座標情報
typedef struct Point {
	double mdX, mdY;											// 座標
} Point;

// 内部変数
static FILE *szpFpI;											// 入力
static Point sz1Point[D_POINT_CNT];								// 座標

// 内部変数 - テスト用
#ifdef D_TEST
	static int siRes;
	static FILE *szpFpA;
	static int siTNo;
#endif

// 出力
int
fOut(
	char *pcpLine				// <I> １行
)
{
	char lc1Buf[1024];

#ifdef D_TEST
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, pcpLine)) {
		siRes = -1;
	}
#else
	printf("%s", pcpLine);
#endif

	return 0;
}

// 対象点の三角形の内外
int
fInPoint(
	double pdX1					// <I> 三角形頂点１Ｘ
	, double pdY1				// <I> 三角形頂点１Ｙ
	, double pdX2				// <I> 三角形頂点２Ｘ
	, double pdY2				// <I> 三角形頂点２Ｙ
	, double pdX3				// <I> 三角形頂点３Ｘ
	, double pdY3				// <I> 三角形頂点３Ｙ
	, double pdX				// <I> 対象点Ｘ
	, double pdY				// <I> 対象点Ｙ
)
{
	double ldArea = 0.5 * (-pdY2 * pdX3 + pdY1 * (-pdX2 + pdX3) + pdX1 * (pdY2 - pdY3) + pdX2 * pdY3);
	double ldS = 1 / (2 * ldArea) * (pdY1 * pdX3 - pdX1 * pdY3 + (pdY3 - pdY1) * pdX + (pdX1 - pdX3) * pdY);
	double ldT = 1 / (2 * ldArea) * (pdX1 * pdY2 - pdY1 * pdX2 + (pdY1 - pdY2) * pdX + (pdX2 - pdX1) * pdY);
	double ldVal = 1 - ldS - ldT;

	if ((0 < ldS && ldS < 1) && (0 < ldT  && ldT < 1) && (0 < ldVal && ldVal < 1)) {
		return 0;
	}
	else {
		return -1;
	}
}

// 実行メイン
int
fMain(
)
{
	int i, j, k, l, liRet;
	char lc1Buf[1024];

	// 座標 - 取得
	for (i = 0; i < D_POINT_CNT; i++) {
		fgets(lc1Buf, sizeof(lc1Buf), szpFpI);
		sscanf(lc1Buf, "%lf%lf", &sz1Point[i].mdX, &sz1Point[i].mdY);
	}

	// １・２点目 - 選択
	for (i = 0; i < D_POINT_CNT - 2; i++) {
		for (j = i + 1; j < D_POINT_CNT - 1; j++) {

			// 傾き
			int liA1;
			double ldA1;
			double ldX = sz1Point[i].mdX - sz1Point[j].mdX;
			if (ldX == 0.0) {
				liA1 = -1;
			}
			else {
				liA1 = 0;
				ldA1 = (sz1Point[i].mdY - sz1Point[j].mdY) / ldX;
			}

			// ３点目 - 選択
			for (k = j + 1; k < D_POINT_CNT; k++) {

				// 傾き
				int liA2;
				double ldA2;
				ldX = sz1Point[j].mdX - sz1Point[k].mdX;
				if (ldX == 0.0) {
					liA2 = -1;
				}
				else {
					liA2 = 0;
					ldA2 = (sz1Point[j].mdY - sz1Point[k].mdY) / ldX;
				}

				// １直線
				if (liA1 == liA2) {
					if (liA1 == 0) {
						if (ldA1 == ldA2) {
							return -1;
						}
					}
					else {
						return -1;
					}
				}

				// 点の三角形の内外
				for (l = 0; l < D_POINT_CNT; l++) {
					if (i != l && j != l && k != l) {
						liRet = fInPoint(sz1Point[i].mdX, sz1Point[i].mdY, sz1Point[j].mdX, sz1Point[j].mdY, sz1Point[k].mdX, sz1Point[k].mdY, sz1Point[l].mdX, sz1Point[l].mdY);
						if (liRet == 0) {
							return -1;
						}
					}
				}
			}
		}
	}

	return 0;
}

// １回実行
int
fOne(
)
{
	int liRet;
	char lc1Buf[1024];

	// 入力 - セット
#ifdef D_TEST
	sprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo);
	szpFpI = fopen(lc1Buf, "r");
	sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo);
	szpFpA = fopen(lc1Buf, "r");
	siRes = 0;
#else
	szpFpI = stdin;
#endif

	// 実行メイン
	liRet = fMain();

	// 出力
	if (liRet == 0) {
		sprintf(lc1Buf, "YES\n");
	}
	else {
		sprintf(lc1Buf, "NO\n");
	}
	fOut(lc1Buf);

	// 残データ有無
#ifdef D_TEST
	lc1Buf[0] = '\0';
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, "")) {
		siRes = -1;
	}
#endif

	// テストファイルクローズ
#ifdef D_TEST
	fclose(szpFpI);
	fclose(szpFpA);
#endif

	// テスト結果
#ifdef D_TEST
	if (siRes == 0) {
		printf("OK %d\n", siTNo);
	}
	else {
		printf("NG %d\n", siTNo);
	}
#endif

	return 0;
}

// プログラム開始
int
main()
{

#ifdef D_TEST
	int i;
	for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {
		siTNo = i;
		fOne();
	}
#else
	fOne();
#endif

	return 0;
}
0