結果
問題 |
No.251 大きな桁の復習問題(1)
|
ユーザー |
|
提出日時 | 2015-11-08 21:14:24 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 7 ms / 5,000 ms |
コード長 | 1,002 bytes |
コンパイル時間 | 610 ms |
コンパイル使用メモリ | 60,852 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-17 17:30:28 |
合計ジャッジ時間 | 1,487 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#include <iostream> #include <cstdio> #include <vector> #include <algorithm> #include <utility> #include <string> using namespace std; const long long P = 129402307; // x^y % mod long long mod_pow(long long x, long long y, long long m) { long long ret = 1; while (y > 0) { if (y & 1) { ret = (ret * x) % m; } x = (x * x) % m; y >>= 1; } return ret; } int main() { string nstr, mstr; long long n, m; cin >> nstr >> mstr; n = 0; for (int i = 0; i < nstr.size(); i++) { n = (n * 10 + nstr[i] - '0') % P; } if (n == 0) { cout << (mstr == "0" ? 1 : 0) << endl; return 0; } // 今、n = n % pだから n < p // またpが素数 // なので n, pは互いに素 // フェルマーの小定理 // a^(p-1) ≡ 1 (mod p) m = 0; for (int i = 0; i < mstr.size(); i++) { m = (m * 10 + mstr[i] - '0') % (P - 1); } /* long long ans = 1; for (int i = 0; i < m; i++) { ans = (ans * n) % P; } */ long long ans = mod_pow(n, m, P); cout << ans << endl; return 0; }