#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]; // 素数化配列A static int si2PB[D_ARRAY_MAX][D_PRM_MAX]; // 素数化配列B // 内部変数 - テスト用 #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; } // 素数化配列A - 作成 int fPrmA( ) { 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; } // 素数化配列B - 作成 int fPrmB( ) { int i, j, k; for (i = 0; i < siACnt; i++) { for (j = 0; j < siMCnt; j++) { int liNo = (i + j) % siACnt; for (k = 0; k < siPCnt; k++) { si2PB[i][k] += si2PA[liNo][k]; } } } 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); // 素数化配列A - 作成 fPrmA(); // 素数化配列B - 作成 fPrmB(); // 最小公倍数 - 作成 int li1Lcm[D_PRM_MAX]; memset(li1Lcm, 0, sizeof(li1Lcm)); for (i = 0; i < siACnt; i++) { for (j = 0; j < siPCnt; j++) { if (li1Lcm[j] < si2PB[i][j]) { li1Lcm[j] = si2PB[i][j]; } } } // 値化 long long llVal = 1; for (i = 0; i < siPCnt; i++) { for (j = 0; j < li1Lcm[i]; 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)); // 素数化配列A memset(si2PB, 0, sizeof(si2PB)); // 素数化配列B // 入力 - セット #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; }