#include #include #include #include #include #include #include // 内部定数 #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 // 1行 ) { 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 // 三角形頂点1X , double pdY1 // 三角形頂点1Y , double pdX2 // 三角形頂点2X , double pdY2 // 三角形頂点2Y , double pdX3 // 三角形頂点3X , double pdY3 // 三角形頂点3Y , double pdX // 対象点X , double pdY // 対象点Y ) { 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); } // 1・2点目 - 選択 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; } // 3点目 - 選択 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; } // 1直線 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; } // 1回実行 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; }