#include #include #include #include #include #include #include // 内部定数 #define D_ON 1 // 汎用フラグ - ON #define D_OFF 0 // 汎用フラグ - OFF #define D_MOD (long long)1000000007 // 除数(10の9乗+7) #define D_ARRAY_MAX 1000 // 最大配列数 #define D_PRM_MAX 100000 // 最大素数数 // 内部変数 static FILE *szpFpI; // 入力 static int si1Array[D_ARRAY_MAX]; // 配列 static int siACnt; // 配列数 static int siMCnt; // 積数 static char sc1PWork[D_PRM_MAX + 5]; // 素数作成用ワーク static int si1Prm[D_PRM_MAX]; // 素数 static int siPCnt; // 素数数 static int si2PA[D_ARRAY_MAX][D_PRM_MAX]; // 素数化配列 // 内部変数 - テスト用 #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昇順 int fSortFnc( const void *pzpVal1 // 値1 , const void *pzpVal2 // 値2 ) { int *lipVal1 = (int *)pzpVal1; int *lipVal2 = (int *)pzpVal2; // int昇順 if (*lipVal1 > *lipVal2) { return 1; } else if (*lipVal1 < *lipVal2) { return -1; } return 0; } // 素数化配列 - 作成 int fPrm( ) { int i, j; // 開始値 int liNow = 2; // 作成 - 開始 while (1) { // 素数の検索 while (sc1PWork[liNow] != D_OFF) { liNow++; } // 上限チェック if (liNow > D_PRM_MAX) { break; } // 素数化 int liCnt = 0; for (i = 0; i < siACnt; i++) { while (si1Array[i] % liNow == 0) { si1Array[i] /= liNow; si2PA[i][siPCnt]++; liCnt++; } } if (liCnt > 0) { si1Prm[siPCnt] = liNow; siPCnt++; } // ワークに素数以外をセット int liVal = liNow; while (liVal <= D_PRM_MAX) { sc1PWork[liVal] = D_ON; liVal += liNow; } } // 残った素数 int li1Array[D_ARRAY_MAX]; memcpy(li1Array, si1Array, sizeof(li1Array)); qsort(li1Array, siACnt, sizeof(int), fSortFnc); for (i = 0; i < siACnt; i++) { if (li1Array[i] < 2) { continue; } // 素数化 int liCnt = 0; for (j = 0; j < siACnt; j++) { if (si1Array[j] % li1Array[i] == 0) { si2PA[j][siPCnt]++; liCnt++; } } if (liCnt > 0) { si1Prm[siPCnt] = li1Array[i]; siPCnt++; } } return 0; } // 実行メイン int fMain( ) { int i, j; char lc1Buf[1024]; // 配列数・積数 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d%d", &siACnt, &siMCnt); // 配列 - 取得 for (i = 0; i < siACnt; i++) { fscanf(szpFpI, "%d", &si1Array[i]); } fgets(lc1Buf, sizeof(lc1Buf), szpFpI); // 素数化配列 - 作成 fPrm(); // 最小公倍数 - 作成 long long llVal = 1; for (i = 0; i < siPCnt; i++) { // 最大個数 - 初期値 int liMax = 0; for (j = 0; j < siMCnt; j++) { liMax += si2PA[j][i]; } // 最大個数 - 更新 int liSum = liMax; for (j = 0; j < siACnt - 1; j++) { liSum -= si2PA[j][i]; liSum += si2PA[(siMCnt + j) % siACnt][i]; if (liMax < liSum) { liMax = liSum; } } // 値化 for (j = 0; j < liMax; j++) { llVal *= si1Prm[i]; llVal %= D_MOD; } } return (int)llVal; } // 1回実行 int fOne( ) { int liRet; char lc1Buf[1024]; // データ - 初期化 siPCnt = 0; // 素数数 memset(sc1PWork, D_OFF, sizeof(sc1PWork)); // 素数作成用ワーク memset(si2PA, 0, sizeof(si2PA)); // 素数化配列 // 入力 - セット #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(); // 結果 - セット sprintf(lc1Buf, "%d\n", liRet); // 結果 - 出力 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; }