結果
問題 | No.614 壊れたキャンパス |
ユーザー |
![]() |
提出日時 | 2019-04-19 16:29:52 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 534 ms / 2,000 ms |
コード長 | 13,713 bytes |
コンパイル時間 | 325 ms |
コンパイル使用メモリ | 35,200 KB |
実行使用メモリ | 43,468 KB |
最終ジャッジ日時 | 2024-09-22 11:32:20 |
合計ジャッジ時間 | 7,248 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <float.h>#include <limits.h>#include <math.h>#include <stdio.h>#include <stdlib.h>#include <string.h>#include <time.h>// 内部定数#define D_MOVE_MAX 200000 // 最大移動数#define D_TREE_MAX 400005 // 木の最大データ数#define D_TREE_WCNT 2 // 木の方向数#define D_TREE_LEFT 0 // 木の方向 - 左側#define D_TREE_RIGHT 1 // 木の方向 - 右側#define D_VTX_MAX D_TREE_MAX // 最大頂点数#define D_EDGE_MAX 500005 // 最大辺数#define D_HEAP_MAX D_VTX_MAX // 最大ヒープ数// 内部構造体 - 移動情報typedef struct Move {int miBNo; // 建物 1~int miFVal; // 移動元階int miTVal; // 移動先階} Move;// 内部構造体 - 木構造typedef struct Tree {int miBNo; // 建物 1~int miVal; // 階 1~int miVNo; // 頂点番号 0~int mi1Height[D_TREE_WCNT]; // 木の高さstruct Tree *mzp1Child[D_TREE_WCNT]; // 子} Tree;// 内部構造体 - 辺情報typedef struct Edge {int miVNo; // 接続先頂点int miVal; // 階差struct Edge *mzpNext; // 次の辺情報} Edge;// 内部構造体 - 頂点情報typedef struct Vtx {Edge *mzpEdge; // 辺long long mlSum; // 階合計} Vtx;// 内部構造体 - ヒープ情報typedef struct Heap {long long mlSum; // 階合計int miVNo; // 頂点番号 0~} Heap;// 内部変数static FILE *szpFpI; // 入力static int siBCnt; // 建物数static Move sz1Move[D_MOVE_MAX]; // 移動情報static int siMCnt; // 移動情報数static Tree sz1Tree[D_TREE_MAX]; // 木の実データstatic int siTCnt; // 木の実データ数static Tree *szpTop; // 先頭の木データstatic Tree *szpTPre; // 1つ前の木データstatic Vtx sz1Vtx[D_VTX_MAX]; // 頂点static int siVCnt; // 頂点数static Edge sz1Edge[D_EDGE_MAX * 2]; // 辺static int siECnt; // 辺数static Heap sz1Heap[D_HEAP_MAX]; // ヒープstatic int siHCnt; // ヒープ数// 内部変数 - テスト用#ifdef D_TESTstatic int siRes;static FILE *szpFpA;static int siTNo;#endif// 出力intfOut(char *pcpLine // <I> 1行){char lc1Buf[1024];#ifdef D_TESTfgets(lc1Buf, sizeof(lc1Buf), szpFpA);if (strcmp(lc1Buf, pcpLine)) {siRes = -1;}#elseprintf("%s", pcpLine);#endifreturn 0;}// 木データ - 作成Tree *fTreeMake(int piBNo // <I> 建物 1~, int piVal // <I> 階 1~, int piVNo // <I> 頂点番号 0~){// 対象の木データTree *lzpTree = &(sz1Tree[siTCnt]);siTCnt++;// データセットmemset(lzpTree, 0, sizeof(Tree)); // 初期化lzpTree->miBNo = piBNo; // 建物lzpTree->miVal = piVal; // 階lzpTree->miVNo = piVNo; // 頂点番号return lzpTree;}// 木データ - 比較 - 建物昇順 - 階昇順intfTreeCmp(Tree *pzpTree // <I> 木データ, int piBNo // <I> 建物 1~, int piVal // <I> 階 1~){// 建物昇順if (piBNo < pzpTree->miBNo) { // 左側return -1;}else if (piBNo > pzpTree->miBNo) { // 右側return 1;}// 階昇順if (piVal < pzpTree->miVal) { // 左側return -1;}else if (piVal > pzpTree->miVal) { // 右側return 1;}return 0;}// 木データ - 検索Tree *fTreeSrh(int piBNo // <I> 建物 1~, int piVal // <I> 階 1~){// 先頭の木データTree *lzpNow = szpTop;// 検索while (1) {// データ有無if (lzpNow == NULL) {return NULL;}// 比較int liRet = fTreeCmp(lzpNow, piBNo, piVal);if (liRet == 0) { // 一致return lzpNow;}// 子へ移動if (liRet < 0) { // 左側lzpNow = lzpNow->mzp1Child[D_TREE_LEFT];}else { // 右側lzpNow = lzpNow->mzp1Child[D_TREE_RIGHT];}}return NULL;}// 木データ - 高さ取得intfTreeGetHeight(Tree *pzpTree // <I> 対象の木情報){// データ有無if (pzpTree == NULL) {return 0;}if (pzpTree->mi1Height[D_TREE_LEFT] >= pzpTree->mi1Height[D_TREE_RIGHT]) {return pzpTree->mi1Height[D_TREE_LEFT] + 1;}else {return pzpTree->mi1Height[D_TREE_RIGHT] + 1;}}// 木データ - 右回転(親は子の右下へ・子は親の左上へ)intfTreeRttR(Tree **pzppTree // <I> 回転対象){// 現在の子Tree *lzpChild = (*pzppTree)->mzp1Child[D_TREE_LEFT];// 右回転(*pzppTree)->mzp1Child[D_TREE_LEFT] = lzpChild->mzp1Child[D_TREE_RIGHT]; // 親の左側 = 子の右側(*pzppTree)->mi1Height[D_TREE_LEFT] = lzpChild->mi1Height[D_TREE_RIGHT]; // 親の高さ(左) = 子の高さ(右)lzpChild->mzp1Child[D_TREE_RIGHT] = *pzppTree; // 子の右側 = 親lzpChild->mi1Height[D_TREE_RIGHT] = fTreeGetHeight(*pzppTree); // 子の高さ(右) - 親の高さ*pzppTree = lzpChild; // 親 = 子return 0;}// 木データ - 左回転(親は子の左下へ・子は親の右上へ)intfTreeRttL(Tree **pzppTree // <I> 回転対象){// 現在の子Tree *lzpChild = (*pzppTree)->mzp1Child[D_TREE_RIGHT];// 左回転(*pzppTree)->mzp1Child[D_TREE_RIGHT] = lzpChild->mzp1Child[D_TREE_LEFT]; // 親の右側 = 子の左側(*pzppTree)->mi1Height[D_TREE_RIGHT] = lzpChild->mi1Height[D_TREE_LEFT]; // 親の高さ(右) = 子の高さ(左)lzpChild->mzp1Child[D_TREE_LEFT] = *pzppTree; // 子の左側 = 親lzpChild->mi1Height[D_TREE_LEFT] = fTreeGetHeight(*pzppTree); // 子の高さ(左) - 親の高さ*pzppTree = lzpChild; // 親 = 子return 0;}// 木データ - 追加・削除の共通処理// 戻り値:[1]高さの変更あり [0]高さの変更なしintfTreeComAddDel(Tree **pzppNow // <I> 現在の木情報, int piWay // <I> 対象の方向){// 高さの変更があるかチェックint liNew = fTreeGetHeight((*pzppNow)->mzp1Child[piWay]);if ((*pzppNow)->mi1Height[piWay] == liNew) { // 変化なしreturn 0;}(*pzppNow)->mi1Height[piWay] = liNew; // 更新// 高さが離れている場合、回転if ((*pzppNow)->mi1Height[D_TREE_LEFT] - (*pzppNow)->mi1Height[D_TREE_RIGHT] > 1) { // 左側が高いfTreeRttR(pzppNow); // 右回転}else if ((*pzppNow)->mi1Height[D_TREE_RIGHT] - (*pzppNow)->mi1Height[D_TREE_LEFT] > 1) { // 右側が高いfTreeRttL(pzppNow); // 左回転}return 1;}// 木データ - 追加// 戻り値:[1]高さの変更あり [0]高さの変更なし [-1]追加なしintfTreeAdd(Tree **pzppNow // <I> 現在の木情報, int piBNo // <I> 建物 1~, int piVal // <I> 階 1~, int piVNo // <I> 頂点番号 0~){// 作成if (*pzppNow == NULL) {*pzppNow = fTreeMake(piBNo, piVal, piVNo);return 1;}// 比較int liRet = fTreeCmp(*pzppNow, piBNo, piVal);// 一致if (liRet == 0) {return -1;}// 方向の判別int liWay;if (liRet < 0) { // 左側liWay = D_TREE_LEFT;}else { // 右側liWay = D_TREE_RIGHT;}// 下位へliRet = fTreeAdd(&((*pzppNow)->mzp1Child[liWay]), piBNo, piVal, piVNo);if (liRet < 1) { // 高さの変更なし or 追加なしreturn liRet;}// 追加・削除の共通処理return fTreeComAddDel(pzppNow, liWay);}// 辺 - 追加intfAddEdge(int piVFNo // <I> 頂点 - 元 0~, int piVTNo // <I> 頂点 - 先 0~, int piVal // <I> 階差 0~){sz1Edge[siECnt].miVNo = piVTNo;sz1Edge[siECnt].miVal = piVal;sz1Edge[siECnt].mzpNext = sz1Vtx[piVFNo].mzpEdge;sz1Vtx[piVFNo].mzpEdge = &sz1Edge[siECnt];siECnt++;return 0;}// 同一建物内の頂点を辺で結合intfTreeVNoConv(Tree *pzpNow // <I> 現在の木情報){// データ有無if (pzpNow == NULL) {return 0;}// 左側fTreeVNoConv(pzpNow->mzp1Child[D_TREE_LEFT]);// 1つ前の木データif (szpTPre != NULL) {if (szpTPre->miBNo == pzpNow->miBNo) {int liVal = abs(szpTPre->miVal - pzpNow->miVal);fAddEdge(szpTPre->miVNo, pzpNow->miVNo, liVal);fAddEdge(pzpNow->miVNo, szpTPre->miVNo, liVal);}}szpTPre = pzpNow;// 右側fTreeVNoConv(pzpNow->mzp1Child[D_TREE_RIGHT]);return 0;}// ヒープ - 比較 - 階合計昇順intfHeapCmp(int piNo1 // <I> 配列番号1 0~, int piNo2 // <I> 配列番号2 0~){// 階合計昇順if (sz1Heap[piNo1].mlSum < sz1Heap[piNo2].mlSum) {return -1;}else if (sz1Heap[piNo1].mlSum > sz1Heap[piNo2].mlSum) {return 1;}return 0;}// ヒープ - 親子関係チェック// 戻り値:[>=0]:変更した子の配列番号 [-1]:変更なしintfHeapChk(int piPNo // <I> 親の配列番号 0~){int liRet;// 最小値int liMNo = piPNo;// 左の子と比較int liCNo = piPNo * 2 + 1;if (liCNo < siHCnt) {liRet = fHeapCmp(liMNo, liCNo);if (liRet == 1) {liMNo = liCNo;}}// 右の子と比較liCNo = piPNo * 2 + 2;if (liCNo < siHCnt) {liRet = fHeapCmp(liMNo, liCNo);if (liRet == 1) {liMNo = liCNo;}}// 自分が最小値であるかチェックif (piPNo == liMNo) {return -1;}// 値の交換Heap lzWork;memcpy(&lzWork, &sz1Heap[liMNo], sizeof(Heap));memcpy(&sz1Heap[liMNo], &sz1Heap[piPNo], sizeof(Heap));memcpy(&sz1Heap[piPNo], &lzWork, sizeof(Heap));return liMNo;}// ヒープ - キュー追加intfHeapEnqueue(long long plSum // <I> 階合計, int piVNo // <I> 頂点番号 0~){int liRet;// 階合計 - チェックif (sz1Vtx[piVNo].mlSum <= plSum) {return 0;}sz1Vtx[piVNo].mlSum = plSum;// 末尾に追加sz1Heap[siHCnt].mlSum = plSum;sz1Heap[siHCnt].miVNo = piVNo;siHCnt++;// 親子関係チェックint liNo = siHCnt - 1;while (1) {// 親の配列番号liNo = (liNo - 1) / 2;// 親子関係チェックliRet = fHeapChk(liNo);if (liRet < 0) {break;}}return 0;}// ヒープ - キュー取得intfHeapDequeue(Heap *pzpRet // <O> 取得先){// データ数if (siHCnt < 1) {return -1;}// 取得memcpy(pzpRet, &sz1Heap[0], sizeof(Heap));siHCnt--;// データ数if (siHCnt < 1) {return 0;}// 末尾を先頭へmemcpy(&sz1Heap[0], &sz1Heap[siHCnt], sizeof(Heap));// 親子関係チェックint liNo = 0;while (liNo >= 0) {liNo = fHeapChk(liNo);}return 0;}// 実行メインlong longfMain(){int i, liRet, liWork;char lc1Buf[1024];// 建物数・移動情報数・開始階・終了階 - 取得int liSVal, liEVal;fgets(lc1Buf, sizeof(lc1Buf), szpFpI);sscanf(lc1Buf, "%d%d%d%d%d", &siBCnt, &siMCnt, &liWork, &liSVal, &liEVal);// 開始階・終了階の頂点登録fTreeAdd(&szpTop, 1, liSVal, 0);fTreeAdd(&szpTop, siBCnt, liEVal, 1);siVCnt = 2;// 移動情報 - 取得for (i = 0; i < siMCnt; i++) {fgets(lc1Buf, sizeof(lc1Buf), szpFpI);sscanf(lc1Buf, "%d%d%d", &sz1Move[i].miBNo, &sz1Move[i].miFVal, &sz1Move[i].miTVal);// 頂点登録liRet = fTreeAdd(&szpTop, sz1Move[i].miBNo, sz1Move[i].miFVal, siVCnt);if (liRet >= 0) {siVCnt++;}liRet = fTreeAdd(&szpTop, sz1Move[i].miBNo + 1, sz1Move[i].miTVal, siVCnt);if (liRet >= 0) {siVCnt++;}}// 頂点 - 初期化for (i = 0; i < siVCnt; i++) {sz1Vtx[i].mlSum = LLONG_MAX;}// 同一建物内の頂点を辺で結合szpTPre = NULL;fTreeVNoConv(szpTop);// 移動情報で頂点を辺で結合for (i = 0; i < siMCnt; i++) {Tree *lzpFrom = fTreeSrh(sz1Move[i].miBNo, sz1Move[i].miFVal);Tree *lzpTo = fTreeSrh(sz1Move[i].miBNo + 1, sz1Move[i].miTVal);fAddEdge(lzpFrom->miVNo, lzpTo->miVNo, 0);}// ヒープ - 初期値fHeapEnqueue(0, 0);// 頂点 - 移動while (1) {// ヒープ - 取得Heap lzHeap;liRet = fHeapDequeue(&lzHeap);if (liRet != 0) {break;}// 辺でループEdge *lzpEdge = sz1Vtx[lzHeap.miVNo].mzpEdge;while (lzpEdge != NULL) {// ヒープ - 追加fHeapEnqueue(lzHeap.mlSum + (long long)lzpEdge->miVal, lzpEdge->miVNo);// 次の辺へlzpEdge = lzpEdge->mzpNext;}}// 到着チェックif (sz1Vtx[1].mlSum == LLONG_MAX) {return -1;}else {return sz1Vtx[1].mlSum;}}// 1回実行intfOne(){long long llRet;char lc1Buf[1024];// データ - 初期化siTCnt = 0; // 木の実データ数szpTop = NULL; // 先頭の木データmemset(sz1Vtx, 0, sizeof(sz1Vtx)); // 頂点siECnt = 0; // 辺数siHCnt = 0; // ヒープ数// 入力 - セット#ifdef D_TESTsprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo);szpFpI = fopen(lc1Buf, "r");sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo);szpFpA = fopen(lc1Buf, "r");siRes = 0;#elseszpFpI = stdin;#endif// 実行メインllRet = fMain();// 出力sprintf(lc1Buf, "%lld\n", llRet);fOut(lc1Buf);// 残データ有無#ifdef D_TESTlc1Buf[0] = '\0';fgets(lc1Buf, sizeof(lc1Buf), szpFpA);if (strcmp(lc1Buf, "")) {siRes = -1;}#endif// テストファイルクローズ#ifdef D_TESTfclose(szpFpI);fclose(szpFpA);#endif// テスト結果#ifdef D_TESTif (siRes == 0) {printf("OK %d\n", siTNo);}else {printf("NG %d\n", siTNo);}#endifreturn 0;}// プログラム開始intmain(){#ifdef D_TESTint i;for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {siTNo = i;fOne();}#elsefOne();#endifreturn 0;}