結果

問題 No.1140 EXPotentiaLLL!
ユーザー kurotemkokurotemko
提出日時 2020-08-03 20:08:50
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,591 bytes
コンパイル時間 2,135 ms
コンパイル使用メモリ 205,040 KB
実行使用メモリ 43,392 KB
最終ジャッジ日時 2024-09-13 14:56:04
合計ジャッジ時間 8,963 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 284 ms
43,392 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 AC 132 ms
43,264 KB
testcase_09 AC 132 ms
43,264 KB
testcase_10 AC 131 ms
43,392 KB
testcase_11 WA -
testcase_12 AC 133 ms
43,264 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

long long mod = 1000000007;
vector<bool> vp;
vector<long long> vk(5000001);

// ------------------------------
// 素数列挙
// 参考:https://ei1333.github.io/luzhiled/snippets/math/prime-table.html
void prime_table(long long n) {
  vector<bool> prime(n+1, true);
  if(n >= 0) prime[0] = false;
  if(n >= 1) prime[1] = false;
  for(long long i = 2; i*i <= n; ++i) {
    if(!prime[i]) continue;
    for(long long j = i+i; j <= n; j += i) {
      prime[j] = false;
    }
  }
  vp = prime;
}
// ------------------------------

void kaijou() {
    vk[0] = 1;
    for(long long i = 1; i <= 5000000; ++i) vk[i] = (vk[i-1]*i)%mod;
}

// ------------------------------
// 繰り返し二乗法
// 参考:https://math.nakaken88.com/textbook/cp-binary-exponentiation-and-recursive-function/
long long modPow(long long a, long long n, long long p) {
  if (n == 0) return 1;
  if (n == 1) return a % p;
  if (n % 2 == 1) return (a * modPow(a, n - 1, p)) % p;
  long long t = modPow(a, n / 2, p);
  return (t * t) % p;
}
// ------------------------------

int main() {

    // あらかじめ素数テーブルと階乗テーブルを作成
    // ちゃんとmodを取っておく
    prime_table(5000000);
    kaijou();

    long long t;
    cin >> t;

    for(long long i = 0; i < t; ++i) {
        long long a;
        long long p;
        scanf("%lld %lld", &a, &p);
        if(vp[p]) {
            long long ans = modPow((a%p), vk[p], p);
            printf("%lld\n", ans);
        }
        else {
            printf("-1\n");
        }
    }
}
0