#include #include #include #include #include #include #include // 内部定数 #define D_MOD 129402307 // 除数 #define D_CHAR_MAX 100005 // 最大文字数 #define D_VAL_UNIT 1000000000000000000 // 値 - 単位 #define D_VAL_CNT 10000 // 値 - 数 // 内部構造体 - 値情報 typedef struct Val { long long ml1Val[D_VAL_CNT]; // 値 [0]:最小桁 [0>]:上位桁 int miVCnt; // 値数 int miSign; // 符号 1 or -1 } Val; // 内部変数 static FILE *szpFpI; // 入力 static char sc1N[D_CHAR_MAX], sc1M[D_CHAR_MAX]; // N, M // 内部変数 - テスト用 #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 fValAdjK( Val *pzpVal // 値 ) { // 値数 if (pzpVal->ml1Val[pzpVal->miVCnt] != 0) { pzpVal->miVCnt++; } if (pzpVal->ml1Val[pzpVal->miVCnt - 1] == 0) { pzpVal->miVCnt--; } // 符号 if (pzpVal->ml1Val[pzpVal->miVCnt - 1] > 0) { pzpVal->miSign = 1; } else if (pzpVal->ml1Val[pzpVal->miVCnt - 1] < 0) { pzpVal->miSign = -1; } else { pzpVal->miSign = 0; } return 0; } int fValAdj( Val *pzpVal // 値 ) { int i; // 単位 for (i = 0; i < pzpVal->miVCnt; i++) { pzpVal->ml1Val[i + 1] += pzpVal->ml1Val[i] / D_VAL_UNIT; pzpVal->ml1Val[i] %= D_VAL_UNIT; fValAdjK(pzpVal); } // 符号統一 for (i = pzpVal->miVCnt - 2; i >= 0; i--) { int liSign; if (pzpVal->ml1Val[i] > 0) { liSign = 1; } else if (pzpVal->ml1Val[i] > 0) { liSign = -1; } else { continue; } if (liSign == pzpVal->miSign) { continue; } pzpVal->ml1Val[i] += D_VAL_UNIT * pzpVal->miSign; pzpVal->ml1Val[i + 1] -= pzpVal->miSign; fValAdjK(pzpVal); } return 0; } // 値 - 作成 int fValMakeS( char *pcpVal // 値の文字列 , Val *pzpVal // 作成値 ) { char lc1Buf[1024]; // 初期化 memset(pzpVal, 0, sizeof(Val)); // 取得文字数 int liCCnt = 0; long long llVal = D_VAL_UNIT; while (llVal >= 10) { llVal /= 10; liCCnt++; } // 文字列開始位置・符号 int liSNo, liSign; if (pcpVal[0] == '-') { liSNo = 1; liSign = -1; } else { liSNo = 0; liSign = 1; } // 文字列終了位置 int liENo = -1; for (liENo = liSNo; ; liENo++) { if (pcpVal[liENo] < '0' || '9' < pcpVal[liENo]) { liENo--; break; } } // 値 - 取得 while (liENo >= liSNo) { int liNo = liENo - liCCnt + 1; if (liNo < liSNo) { liCCnt += liNo - liSNo; liNo = liSNo; } liENo -= liCCnt; memcpy(lc1Buf, &pcpVal[liNo], liCCnt); lc1Buf[liCCnt] = '\0'; sscanf(lc1Buf, "%lld", &pzpVal->ml1Val[pzpVal->miVCnt]); pzpVal->ml1Val[pzpVal->miVCnt] *= liSign; pzpVal->miVCnt++; } // 値 - 調整 fValAdj(pzpVal); return 0; } // 値 - 剰余 long long fValMod( Val *pzpVal // 値 , long long plMod // 剰余値 ) { int i; long long llVal = 0; for (i = pzpVal->miVCnt - 1; i >= 0; i--) { llVal *= D_VAL_UNIT; llVal += pzpVal->ml1Val[i]; llVal %= plMod; } return llVal; } // べき乗 - 取得 int fGetModPower( int piBase // 基数 , int piIdx // 指数 ) { // 平方リスト - 作成 int li1Val[100]; li1Val[0] = piBase; int liCnt = 1; int liIdx = 1; while (piIdx > liIdx) { li1Val[liCnt] = (int)((long long)li1Val[liCnt - 1] * (long long)li1Val[liCnt - 1] % (long long)D_MOD); liCnt++; liIdx += liIdx; } // べき乗 - 取得 int liVal = 1; while (piIdx > 0) { if (piIdx >= liIdx) { piIdx -= liIdx; liVal = (int)((long long)liVal * (long long)li1Val[liCnt - 1] % (long long)D_MOD); } liCnt--; liIdx /= 2; } return liVal; } // 実行メイン int fMain( ) { // N - 取得 fgets(sc1N, sizeof(sc1N), szpFpI); if (sc1N[0] == '0') { return 0; } // M - 取得 fgets(sc1M, sizeof(sc1M), szpFpI); if (sc1M[0] == '0') { return 1; } // N - 数値化 Val lzVal; fValMakeS(sc1N, &lzVal); int liN = (int)fValMod(&lzVal, D_MOD); // M - 数値化 fValMakeS(sc1M, &lzVal); int liM = (int)fValMod(&lzVal, D_MOD - 1); // NのM乗 return fGetModPower(liN, liM); } // 1回実行 int fOne( ) { int liRet; char lc1Buf[1024]; // 入力 - セット #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; }