PKU ACM poj 1175 源代码

最新推荐文章于 2020-02-28 14:30:00 发布

原创最新推荐文章于 2020-02-28 14:30:00 发布 · 886 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#struct #allocation #input #c #matrix #存储

本文介绍了一个星空图中识别并标记不同星群的算法。通过深度优先搜索（DFS）找到每个星群，然后对其进行标准化处理及旋转，以便与其他星群进行比较。相似星群将被标记为相同的字母，实现了对星图的有效分类。

思路：

1.（使用DFS算法）搜索得到每一个cluster；

2.和已经存储的cluster进行比较；

3. 如果是一个新的cluster，将其进行旋转，并存储起来供后续比较使用；

4.为了进行比较，需要对得到的cluster进行一些处理：排序、坐标处理。

源代码：

/*Starry Night Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 1102 Accepted: 423 Description High up in the night sky, the shining stars appear in clusters of various shapes. A cluster is a non-empty group of neighbouring stars, adjacent in horizontal, vertical or diagonal direction. A cluster cannot be a part of a larger cluster. Clusters may be similar. Two clusters are similar if they have the same shape and number of stars, irrespective of their orientation. In general, the number of possible orientations for a cluster is eight, as Figure 1 exemplifies. The night sky is represented by a sky map, which is a two-dimensional matrix of 0's and 1's. A cell contains the digit 1 if it has a star, and the digit 0 otherwise. Given a sky map, mark all the clusters with lower case letters. Similar clusters must be marked with the same letter; non-similar clusters must be marked with different letters. You mark a cluster with a lower case letter by replacing every 1 in the cluster by that lower case letter. The order of allocation of letters depends on the first appearance of different clusters, from top to bottom and from left to right. Input Your program is to read from standard input. The first two lines contain, respectively, the width W and the height H of a sky map. The sky map is given in the following H lines, of W characters each. 0 <= W (width of the sky map) <= 100 0 <= H (height of the sky map) <= 100 0 <= Number of clusters <= 500 0 <= Number of non-similar clusters <= 26 (a..z) 1 <= Number of stars per cluster <= 160 Output Your program is to write to standard output. The output contains the same map as the input, except that the clusters are marked as described above. Sample Input 23 15 10001000000000010000000 01111100011111000101101 01000000010001000111111 00000000010101000101111 00000111010001000000000 00001001011111000000000 10000001000000000000000 00101000000111110010000 00001000000100010011111 00000001110101010100010 00000100110100010000000 00010001110111110000000 00100001110000000100000 00001000100001000100101 00000001110001000111000 Sample Output a000a0000000000b0000000 0aaaaa000ccccc000d0dd0d 0a0000000c000c000dddddd 000000000c0b0c000d0dddd 00000eee0c000c000000000 0000e00e0ccccc000000000 b000000e000000000000000 00b0f000000ccccc00a0000 0000f000000c000c00aaaaa 0000000ddd0c0b0c0a000a0 00000b00dd0c000c0000000 000g000ddd0ccccc0000000 00g0000ddd0000000e00000 0000b000d0000f000e00e0b 0000000ddd000f000eee000 Hint */ #include "stdlib.h" #include "stdio.h" /*0 <= W (width of the sky map) <= 100 0 <= H (height of the sky map) <= 100 0 <= Number of clusters <= 500 0 <= Number of non-similar clusters <= 26 (a..z) 1 <= Number of stars per cluster <= 160*/ #define CLUSTER_SHAPE_NUMBER 26 #define NUMBER_OF_STARS_PER_CLUSTER 160 #define SKY_MAP_WIDTH 100 #define SKY_MAP_HEIGHT 100 #define ORIENTATION_NUMBER 8 #define STAR_COLOR_WHITE 0 #define STAR_COLOR_GRAY 1 #define STAR_COLOR_BLACK 2 typedef struct _STAR_ST_ { int iPositionHight; int iPositionWidth; }STAR_ST; typedef struct _CLUSTER_STAR_ST_ { int iStarNum; STAR_ST astStar[NUMBER_OF_STARS_PER_CLUSTER]; }CLUSTER_STAR_ST; typedef struct _CLUSTER_SHAPE_ST_ { int iStarNum; char cLowerCaseLetter; struct { STAR_ST astStar[NUMBER_OF_STARS_PER_CLUSTER]; }astOrientation[ORIENTATION_NUMBER]; }CLUSTER_SHAPE_ST; typedef struct _CLUSTER_ST_ { int iShapeNum; CLUSTER_SHAPE_ST astClusterShape[CLUSTER_SHAPE_NUMBER]; }CLUSTER_ST; typedef struct _SKY_MAP_ST_ { int iWidth; int ihight; char acMapNode[SKY_MAP_HEIGHT][SKY_MAP_WIDTH]; }SKY_MAP_ST; CLUSTER_ST gstCluster; SKY_MAP_ST gstSkyMap; CLUSTER_STAR_ST gstClusterStar; unsigned char gaucStarColor[SKY_MAP_HEIGHT][SKY_MAP_WIDTH]; /************************************************************************/ /* */ /* */ /* */ /* */ /************************************************************************/ void SkyMapInput(void) { int iLoop_Width; int iLoop_Hight; unsigned char aucBuff[SKY_MAP_WIDTH+1]; scanf("%d",&gstSkyMap.iWidth); scanf("%d",&gstSkyMap.ihight); for(iLoop_Hight = 0; iLoop_Hight < gstSkyMap.ihight; iLoop_Hight++) { scanf("%s",aucBuff); for(iLoop_Width = 0; iLoop_Width < gstSkyMap.iWidth; iLoop_Width++) { if('1' == aucBuff[iLoop_Width]) { gstSkyMap.acMapNode[iLoop_Hight][iLoop_Width] = 1; } else { gstSkyMap.acMapNode[iLoop_Hight][iLoop_Width] = 0; } } } return; } void StarMapDFS(int viHightIndex,int viWidthIndex) { int iHightDelta; int iWidthDelta; int iNextHightIndex; int iNextWidthIndex; gaucStarColor[viHightIndex][viWidthIndex] = STAR_COLOR_GRAY; gstClusterStar.astStar[gstClusterStar.iStarNum].iPositionHight = viHightIndex; gstClusterStar.astStar[gstClusterStar.iStarNum].iPositionWidth = viWidthIndex; gstClusterStar.iStarNum++; for (iHightDelta = -1; iHightDelta < 2; iHightDelta++) { for (iWidthDelta = -1; iWidthDelta < 2; iWidthDelta++) { if (0 == iHightDelta&& 0 == iWidthDelta) { continue; } iNextHightIndex = viHightIndex+iHightDelta; iNextWidthIndex = viWidthIndex+iWidthDelta; if (0 > iNextHightIndex || gstSkyMap.ihight <= iNextHightIndex) { continue; } if (0 > iNextWidthIndex || gstSkyMap.iWidth <= iNextWidthIndex) { continue; } if (1 == gstSkyMap.acMapNode[iNextHightIndex][iNextWidthIndex] &&STAR_COLOR_WHITE == gaucStarColor[iNextHightIndex][iNextWidthIndex]) { StarMapDFS(iNextHightIndex,iNextWidthIndex); } } } gaucStarColor[viHightIndex][viWidthIndex] = STAR_COLOR_BLACK; return; } int StarCmp(const void *a,const void *b) { if (((STAR_ST *)a)->iPositionHight == ((STAR_ST *)b)->iPositionHight) { return ((STAR_ST *)a)->iPositionWidth -((STAR_ST *)b)->iPositionWidth; } else { return ((STAR_ST *)a)->iPositionHight -((STAR_ST *)b)->iPositionHight; } } void AdjustClusterStarIndex(STAR_ST *vpastStar,int viStarNum) { int iLoopStar; int iMinHight = SKY_MAP_HEIGHT; int iMinWidth = SKY_MAP_WIDTH; for(iLoopStar = 0; iLoopStar < viStarNum ; iLoopStar++) { if(iMinHight > vpastStar[iLoopStar].iPositionHight) { iMinHight = vpastStar[iLoopStar].iPositionHight; } if(iMinWidth > vpastStar[iLoopStar].iPositionWidth) { iMinWidth = vpastStar[iLoopStar].iPositionWidth; } } for(iLoopStar = 0; iLoopStar < viStarNum ; iLoopStar++) { vpastStar[iLoopStar].iPositionHight = vpastStar[iLoopStar].iPositionHight - iMinHight ; vpastStar[iLoopStar].iPositionWidth = vpastStar[iLoopStar].iPositionWidth - iMinWidth ; } return; } void CreatOtherOrientation(CLUSTER_SHAPE_ST *vpstClusterShape) { int iLoop; int iLoopStar; for(iLoop = 1; iLoop < 4; iLoop++) { for(iLoopStar = 0; iLoopStar < vpstClusterShape->iStarNum ; iLoopStar++) { vpstClusterShape->astOrientation[iLoop].astStar[iLoopStar].iPositionHight = SKY_MAP_WIDTH - vpstClusterShape->astOrientation[iLoop-1].astStar[iLoopStar].iPositionWidth; vpstClusterShape->astOrientation[iLoop].astStar[iLoopStar].iPositionWidth = vpstClusterShape->astOrientation[iLoop-1].astStar[iLoopStar].iPositionHight; } qsort(vpstClusterShape->astOrientation[iLoop].astStar,vpstClusterShape->iStarNum,sizeof(STAR_ST ),StarCmp); AdjustClusterStarIndex(vpstClusterShape->astOrientation[iLoop].astStar,vpstClusterShape->iStarNum); } for(iLoopStar = 0; iLoopStar < vpstClusterShape->iStarNum ; iLoopStar++) { vpstClusterShape->astOrientation[4].astStar[iLoopStar].iPositionHight = SKY_MAP_HEIGHT - vpstClusterShape->astOrientation[0].astStar[iLoopStar].iPositionHight; vpstClusterShape->astOrientation[4].astStar[iLoopStar].iPositionWidth = vpstClusterShape->astOrientation[0].astStar[iLoopStar].iPositionWidth; } qsort(vpstClusterShape->astOrientation[4].astStar,vpstClusterShape->iStarNum,sizeof(STAR_ST ),StarCmp); AdjustClusterStarIndex(vpstClusterShape->astOrientation[4].astStar,vpstClusterShape->iStarNum); for(iLoop = 5; iLoop < 8; iLoop++) { for(iLoopStar = 0; iLoopStar < vpstClusterShape->iStarNum ; iLoopStar++) { vpstClusterShape->astOrientation[iLoop].astStar[iLoopStar].iPositionHight = SKY_MAP_WIDTH - vpstClusterShape->astOrientation[iLoop-1].astStar[iLoopStar].iPositionWidth; vpstClusterShape->astOrientation[iLoop].astStar[iLoopStar].iPositionWidth = vpstClusterShape->astOrientation[iLoop-1].astStar[iLoopStar].iPositionHight; } qsort(vpstClusterShape->astOrientation[iLoop].astStar,vpstClusterShape->iStarNum,sizeof(STAR_ST ),StarCmp); AdjustClusterStarIndex(vpstClusterShape->astOrientation[iLoop].astStar,vpstClusterShape->iStarNum); } return; } void printfAllOrientation(CLUSTER_SHAPE_ST *vpstClusterShape) { SKY_MAP_ST stSkyMap; int iLoop_Width; int iLoop_Hight; int iLoop_Star; int iLoop_Shape; if(2 >= vpstClusterShape->iStarNum ) { return; } printf("==============================================================/n"); for (iLoop_Shape = 0; iLoop_Shape < 8; iLoop_Shape++) { memset(&stSkyMap,0,sizeof(SKY_MAP_ST)); for (iLoop_Star = 0; iLoop_Star < vpstClusterShape->iStarNum; iLoop_Star++) { stSkyMap.acMapNode[vpstClusterShape->astOrientation[iLoop_Shape].astStar[iLoop_Star].iPositionHight][vpstClusterShape->astOrientation[iLoop_Shape].astStar[iLoop_Star].iPositionWidth] = '*'; } for(iLoop_Hight = 0; iLoop_Hight < gstSkyMap.ihight; iLoop_Hight++) { for(iLoop_Width = 0; iLoop_Width < gstSkyMap.iWidth; iLoop_Width++) { if(0 != stSkyMap.acMapNode[iLoop_Hight][iLoop_Width]) { printf("%c",stSkyMap.acMapNode[iLoop_Hight][iLoop_Width]); } else { printf("0"); } } printf("/n"); } printf("/n"); } printf("==============================================================/n"); return; } int StarryNightMain(void)// { int iLoop_Width; int iLoop_Hight; int iLoop_Star; int iLoop_Shape; int iLoop_Orientation; char cLowerCaseLetterResource = 'a'; char cCurrentLowerCaseLetter; CLUSTER_STAR_ST stClusterStarTemp; SkyMapInput(); for(iLoop_Hight = 0; iLoop_Hight < gstSkyMap.ihight; iLoop_Hight++) { for(iLoop_Width = 0; iLoop_Width < gstSkyMap.iWidth; iLoop_Width++) { if (1 == gstSkyMap.acMapNode[iLoop_Hight][iLoop_Width] &&STAR_COLOR_WHITE == gaucStarColor[iLoop_Hight][iLoop_Width]) { gstClusterStar.iStarNum = 0; StarMapDFS(iLoop_Hight,iLoop_Width); qsort(gstClusterStar.astStar,gstClusterStar.iStarNum,sizeof(STAR_ST ),StarCmp); memcpy(&stClusterStarTemp,&gstClusterStar,sizeof(CLUSTER_STAR_ST)); AdjustClusterStarIndex(gstClusterStar.astStar,gstClusterStar.iStarNum); for (iLoop_Shape = 0; iLoop_Shape < gstCluster.iShapeNum; iLoop_Shape++) { if (gstClusterStar.iStarNum == gstCluster.astClusterShape[iLoop_Shape].iStarNum) { if(2 > gstClusterStar.iStarNum) { cCurrentLowerCaseLetter = gstCluster.astClusterShape[iLoop_Shape].cLowerCaseLetter; goto FOUND; } for (iLoop_Orientation = 0; iLoop_Orientation < ORIENTATION_NUMBER; iLoop_Orientation++) { if(0 == memcmp(gstCluster.astClusterShape[iLoop_Shape].astOrientation[iLoop_Orientation].astStar,gstClusterStar.astStar,gstClusterStar.iStarNum*sizeof(STAR_ST) )) { cCurrentLowerCaseLetter = gstCluster.astClusterShape[iLoop_Shape].cLowerCaseLetter; goto FOUND; } } } } if (iLoop_Shape >= gstCluster.iShapeNum) { cCurrentLowerCaseLetter= cLowerCaseLetterResource; cLowerCaseLetterResource++; gstCluster.astClusterShape[gstCluster.iShapeNum].iStarNum = gstClusterStar.iStarNum; gstCluster.astClusterShape[gstCluster.iShapeNum].cLowerCaseLetter = cCurrentLowerCaseLetter; memcpy(gstCluster.astClusterShape[gstCluster.iShapeNum].astOrientation[0].astStar,gstClusterStar.astStar,gstClusterStar.iStarNum*sizeof(STAR_ST) ); CreatOtherOrientation(&gstCluster.astClusterShape[gstCluster.iShapeNum]); //printfAllOrientation(&gstCluster.astClusterShape[gstCluster.iShapeNum]); gstCluster.iShapeNum++; } FOUND: for (iLoop_Star = 0; iLoop_Star < stClusterStarTemp.iStarNum; iLoop_Star++) { gstSkyMap.acMapNode[stClusterStarTemp.astStar[iLoop_Star].iPositionHight][stClusterStarTemp.astStar[iLoop_Star].iPositionWidth] = cCurrentLowerCaseLetter; } } } } //printf("/n"); for(iLoop_Hight = 0; iLoop_Hight < gstSkyMap.ihight; iLoop_Hight++) { for(iLoop_Width = 0; iLoop_Width < gstSkyMap.iWidth; iLoop_Width++) { if(0 != gstSkyMap.acMapNode[iLoop_Hight][iLoop_Width]) { printf("%c",gstSkyMap.acMapNode[iLoop_Hight][iLoop_Width]); } else { printf("0"); } } printf("/n"); } return 0; } /**************************END*******************************************/