#include #include #include #include #include #include #include // 内部定数 #define D_ON 1 // 汎用フラグ - ON #define D_OFF 0 // 汎用フラグ - OFF #define D_ARRAY_MAX 100000 // 最大配列数 #define D_PRM_MAX 1000 // 最大素数リスト数 // 内部構造体 - 値情報 typedef struct Val { int miVal; // 値 int miMCnt; // 倍数数 } Val; // 内部構造体 - 素数情報 typedef struct Prm { int miVal; // 素数 int miPCnt; // 素数数 } Prm; // 内部変数 static FILE *szpFpI; // 入力 static Val sz1Array[D_ARRAY_MAX]; // 数列 static int siACnt; // 数列数 static int si1Prm[D_PRM_MAX]; // 素数リスト static char sc1PWork[D_PRM_MAX + 5]; // 素数リスト作成用ワーク static int siPCnt; // 素数リスト数 static Prm sz1Div[D_PRM_MAX]; // 素因数分解 static int siDCnt; // 素因数分解数 static int siMax; // 最大倍数数 // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; #endif // ソート関数 - 値昇順 int fSortFnc( const void *pzpVal1 // 値1 , const void *pzpVal2 // 値2 ) { Val *lzpVal1 = (Val *)pzpVal1; Val *lzpVal2 = (Val *)pzpVal2; // 値昇順 if (lzpVal1->miVal > lzpVal2->miVal) { return 1; } else if (lzpVal1->miVal < lzpVal2->miVal) { return -1; } return 0; } // 数列 - 検索 // 戻り値:[>=0]配列番号 [-1]なし int fBSrh( int piVal // 値 ) { // 初期範囲 int liSNo = 0; int liENo = siACnt - 1; // 検索 while (1) { // 中間位置 int liMNo = (liSNo + liENo) / 2; // 一致チェック if (piVal == sz1Array[liMNo].miVal) { return liMNo; } // 範囲を絞る if (piVal < sz1Array[liMNo].miVal) { // 左側へ liENo = liMNo - 1; } else { // 右側へ liSNo = liMNo + 1; } // 範囲チェック if (liSNo > liENo) { return -1; } } return -1; } // 素数リストの作成 - エラトステネスの篩 int fMakePrm( int piMax // 最大値 ) { // 初期化 siPCnt = 0; memset(sc1PWork, D_OFF, sizeof(sc1PWork)); // 開始値 int liNow = 2; // 上限値 int liLimit = (int)sqrt((double)piMax); // 作成 - 開始 while (1) { // 素数の検索 while (sc1PWork[liNow] != D_OFF) { liNow++; } // 上限チェック if (liNow > liLimit) { break; } // 素数リスト - 追加 si1Prm[siPCnt] = liNow; siPCnt++; // ワークに素数以外をセット int liVal = liNow; while (liVal <= piMax) { sc1PWork[liVal] = D_ON; liVal += liNow; } } // 残った素数を追加 while (liNow <= piMax) { if (sc1PWork[liNow] == D_OFF) { si1Prm[siPCnt] = liNow; siPCnt++; } liNow++; } return 0; } // 素因数分解 int fDivPrm( int piVal // 値 ) { int i; // 初期化 memset(sz1Div, 0, sizeof(sz1Div)); siDCnt = 0; // 素因数分解 for (i = 0; i < siPCnt; i++) { while (piVal % si1Prm[i] == 0) { piVal /= si1Prm[i]; sz1Div[siDCnt].miPCnt++; } if (sz1Div[siDCnt].miPCnt > 0) { sz1Div[siDCnt].miVal = si1Prm[i]; siDCnt++; if (piVal == 1) { break; } } } // 平方以上の素数 if (piVal != 1) { sz1Div[siDCnt].miPCnt = 1; sz1Div[siDCnt].miVal = piVal; siDCnt++; } return 0; } // 素数 - 選択 int fSelPrm( int piANo // 数列の配列番号 , int piVal // 現在値 , int piPCnt // 素数数 ) { int i; // 選択終了 if (piPCnt == siDCnt) { // 値検索 int liNo = fBSrh(piVal); if (liNo < 0) { return 0; } if (liNo == piANo) { // 自分 return 0; } // 倍数数 int liMCnt = sz1Array[piANo].miMCnt + 1; // 倍数数 - 更新 if (sz1Array[liNo].miMCnt < liMCnt) { sz1Array[liNo].miMCnt = liMCnt; } // 最大倍数数 - 更新 if (siMax < liMCnt) { siMax = liMCnt; } return 0; } // 選択 for (i = 0; i <= sz1Div[piPCnt].miPCnt; i++) { fSelPrm(piANo, piVal, piPCnt + 1); piVal *= sz1Div[piPCnt].miVal; } return 0; } // 実行メイン int fMain( int piTNo // テスト番号 1~ ) { int i; char lc1Buf[1024], lc1Out[1024]; // データ - 初期化 memset(sz1Array, 0, sizeof(sz1Array)); // 数列 siMax = 0; // 最大倍数数 // 入力 - セット #ifdef D_TEST sprintf(lc1Buf, ".\\Test\\T%d.txt", piTNo); szpFpI = fopen(lc1Buf, "r"); sprintf(lc1Buf, ".\\Test\\A%d.txt", piTNo); szpFpA = fopen(lc1Buf, "r"); siRes = 0; #else szpFpI = stdin; #endif // 数列数 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d", &siACnt); // 数列 - 取得 for (i = 0; i < siACnt; i++) { fscanf(szpFpI, "%d", &sz1Array[i].miVal); } fgets(lc1Buf, sizeof(lc1Buf), szpFpI); // 数列 - ソート qsort(sz1Array, siACnt, sizeof(Val), fSortFnc); // 素数リストの作成 fMakePrm(D_PRM_MAX); // 倍数数 - セット for (i = siACnt - 1; i > 0; i--) { // 素因数分解 fDivPrm(sz1Array[i].miVal); // 素数 - 選択 fSelPrm(i, 1, 0); } // 結果 - セット sprintf(lc1Out, "%d\n", siMax + 1); // 結果 - 表示 #ifdef D_TEST fgets(lc1Buf, sizeof(lc1Buf), szpFpA); if (strcmp(lc1Buf, lc1Out)) { siRes = -1; } #else printf("%s", lc1Out); #endif // 残データ有無 #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", piTNo); } else { printf("NG %d\n", piTNo); } #endif return 0; } int main() { #ifdef D_TEST int i; for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) { fMain(i); } #else fMain(0); #endif return 0; }