結果

問題 No.308 素数は通れません
ユーザー moti
提出日時 2015-12-01 02:34:58
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,452 bytes
コンパイル時間 958 ms
コンパイル使用メモリ 104,036 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-14 06:18:20
合計ジャッジ時間 3,304 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19 WA * 88
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <complex>
#include <queue>
#include <deque>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <iomanip>
#include <assert.h>
#include <array>
#include <cstdio>
#include <cstring>
#include <random>
#include <functional>
#include <numeric>
#include <bitset>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)
#define all(c) (c).begin(), (c).end()
#define zero(a) memset(a, 0, sizeof a)
#define minus(a) memset(a, -1, sizeof a)
#define minimize(a, x) a = std::min(a, x)
#define maximize(a, x) a = std::max(a, x)

typedef long long ll;
int const inf = 1<<29;

typedef __int128_t INT;

INT pow(INT x, INT n, INT M){
  INT tmp = 1;

  if ( n > 0 ){
    tmp = pow(x, n/2, M);
    if ( n%2 == 0 ) tmp = (tmp*tmp)%M;
    else tmp = (((tmp*tmp)%M)*x)%M;
  }
  return tmp;
}

bool is_prime(INT const& x) {
  if(x == 2) { return 1; }
  if(x < 2 || (x%2 == 0)) { return 0; }
  return pow(2, x-1, x) == 1;
}

int main() {

  string s; cin >> s; INT N = 0;
  rep(i, s.size()) {
    N *= 10; N += s[i]-'0';
  }

  if(N < 7) {
    cout << (int)N-1 << endl; exit(0);
  }

  if((N % 7 == 1 || is_prime(N-1)) && is_prime(N-7)) {
    int cnt = 1;
    while((N % (7+cnt) == 1 || is_prime(N-1)) && is_prime(N-7-cnt)) cnt++;
    cout << 7 + cnt << endl;
  }
  else {
    cout << 7 << endl;
  }

  return 0;
}
0