結果

問題 No.774 tatyamと素数大富豪
ユーザー ceni1055ceni1055
提出日時 2018-12-22 00:34:31
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,332 bytes
コンパイル時間 1,726 ms
コンパイル使用メモリ 171,416 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-25 09:48:17
合計ジャッジ時間 3,149 ms
ジャッジサーバーID
(参考情報)
judge4 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 2 ms
5,376 KB
testcase_02 WA -
testcase_03 AC 4 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 AC 1 ms
5,376 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#include <assert.h>

using namespace std;

long long powMod(long long x, long long k, long long m)
{
  if (k == 0)
    return 1;
  if (k % 2 == 0)
    return powMod(x * x % m, k / 2, m);
  else
    return x * powMod(x, k - 1, m) % m;
}

bool suspect(long long a, int s, long long d, long long n)
{
  long long x = powMod(a, d, n);
  if (x == 1)
    return true;
  for (int r = 0; r < s; ++r)
  {
    if (x == n - 1)
      return true;
    x = x * x % n;
  }
  return false;
}
// {2,7,61,-1}                 is for n < 4759123141 (= 2^32)
// {2,3,5,7,11,13,17,19,23,-1} is for n < 10^16 (at least)
bool isPrime(long long n)
{
  if (n <= 1 || (n > 2 && n % 2 == 0))
    return false;
  int test[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, -1};
  long long d = n - 1, s = 0;
  while (d % 2 == 0)
    ++s, d /= 2;
  for (int i = 0; test[i] < n && test[i] != -1; ++i)
    if (!suspect(test[i], s, d, n))
      return false;
  return true;
}

int main()
{
  int n;
  cin >> n;
  vector<string> a(n);
  for (int i = 0; i < n; i++)
    cin >> a[i];

  long long ans = 0;
  do
  {
    string tmp;
    for (int i = 0; i < n; i++)
      tmp += a[i];
    long long check = stoll(tmp);
    if (isPrime(check))
      ans = max(ans, check);
  } while (next_permutation(a.begin(), a.end()));
  cout << ans << endl;

  return 0;
}
0