#64256 (C++11) No.308 素数は通れません

提出ソース

結果

問題	No.308 素数は通れません
ユーザー	moti
提出日時	2015-12-09 20:04:05
言語	C++11 (gcc 11.4.0)
結果	WA
実行時間	-
コード長	5,041 bytes
コンパイル時間	1,058 ms
コンパイル使用メモリ	117,952 KB
実行使用メモリ	6,948 KB
最終ジャッジ日時	2024-09-15 05:58:17
合計ジャッジ時間	3,457 ms
ジャッジサーバーID （参考情報）	judge4 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,812 KB
testcase_01	AC	2 ms 6,940 KB
testcase_02	AC	2 ms 6,944 KB
testcase_03	AC	2 ms 6,940 KB
testcase_04	AC	1 ms 6,944 KB
testcase_05	AC	2 ms 6,944 KB
testcase_06	AC	2 ms 6,940 KB
testcase_07	AC	2 ms 6,944 KB
testcase_08	AC	1 ms 6,940 KB
testcase_09	AC	1 ms 6,940 KB
testcase_10	AC	2 ms 6,940 KB
testcase_11	AC	1 ms 6,940 KB
testcase_12	AC	2 ms 6,944 KB
testcase_13	AC	2 ms 6,944 KB
testcase_14	AC	2 ms 6,940 KB
testcase_15	AC	2 ms 6,940 KB
testcase_16	AC	2 ms 6,940 KB
testcase_17	AC	2 ms 6,944 KB
testcase_18	AC	2 ms 6,940 KB
testcase_19	AC	1 ms 6,944 KB
testcase_20	AC	2 ms 6,944 KB
testcase_21	AC	2 ms 6,940 KB
testcase_22	AC	2 ms 6,944 KB
testcase_23	AC	1 ms 6,940 KB
testcase_24	AC	2 ms 6,940 KB
testcase_25	AC	2 ms 6,944 KB
testcase_26	AC	1 ms 6,940 KB
testcase_27	AC	2 ms 6,940 KB
testcase_28	AC	1 ms 6,944 KB
testcase_29	AC	1 ms 6,940 KB
testcase_30	AC	1 ms 6,940 KB
testcase_31	AC	2 ms 6,940 KB
testcase_32	AC	1 ms 6,940 KB
testcase_33	AC	2 ms 6,940 KB
testcase_34	AC	2 ms 6,944 KB
testcase_35	AC	1 ms 6,940 KB
testcase_36	AC	1 ms 6,940 KB
testcase_37	AC	2 ms 6,940 KB
testcase_38	AC	2 ms 6,944 KB
testcase_39	AC	1 ms 6,940 KB
testcase_40	AC	2 ms 6,940 KB
testcase_41	AC	2 ms 6,940 KB
testcase_42	AC	1 ms 6,944 KB
testcase_43	AC	2 ms 6,940 KB
testcase_44	AC	1 ms 6,940 KB
testcase_45	AC	2 ms 6,944 KB
testcase_46	AC	2 ms 6,940 KB
testcase_47	AC	2 ms 6,944 KB
testcase_48	AC	2 ms 6,944 KB
testcase_49	AC	2 ms 6,944 KB
testcase_50	AC	2 ms 6,940 KB
testcase_51	AC	2 ms 6,940 KB
testcase_52	AC	2 ms 6,940 KB
testcase_53	AC	1 ms 6,944 KB
testcase_54	AC	1 ms 6,940 KB
testcase_55	AC	2 ms 6,940 KB
testcase_56	AC	2 ms 6,940 KB
testcase_57	AC	2 ms 6,944 KB
testcase_58	AC	2 ms 6,940 KB
testcase_59	AC	2 ms 6,940 KB
testcase_60	AC	2 ms 6,944 KB
testcase_61	AC	2 ms 6,944 KB
testcase_62	AC	2 ms 6,944 KB
testcase_63	AC	2 ms 6,940 KB
testcase_64	AC	2 ms 6,940 KB
testcase_65	AC	2 ms 6,940 KB
testcase_66	AC	2 ms 6,944 KB
testcase_67	AC	2 ms 6,940 KB
testcase_68	AC	1 ms 6,940 KB
testcase_69	AC	2 ms 6,944 KB
testcase_70	AC	2 ms 6,940 KB
testcase_71	AC	2 ms 6,944 KB
testcase_72	AC	2 ms 6,940 KB
testcase_73	AC	2 ms 6,940 KB
testcase_74	AC	1 ms 6,940 KB
testcase_75	AC	1 ms 6,944 KB
testcase_76	AC	2 ms 6,940 KB
testcase_77	AC	2 ms 6,940 KB
testcase_78	AC	2 ms 6,940 KB
testcase_79	AC	1 ms 6,944 KB
testcase_80	WA	-
testcase_81	WA	-
testcase_82	AC	1 ms 6,940 KB
testcase_83	AC	2 ms 6,944 KB
testcase_84	WA	-
testcase_85	AC	2 ms 6,940 KB
testcase_86	WA	-
testcase_87	WA	-
testcase_88	AC	2 ms 6,944 KB
testcase_89	AC	1 ms 6,944 KB
testcase_90	WA	-
testcase_91	AC	2 ms 6,944 KB
testcase_92	AC	2 ms 6,940 KB
testcase_93	AC	1 ms 6,940 KB
testcase_94	AC	2 ms 6,940 KB
testcase_95	AC	1 ms 6,944 KB
testcase_96	AC	2 ms 6,944 KB
testcase_97	AC	2 ms 6,944 KB
testcase_98	AC	2 ms 6,940 KB
testcase_99	AC	1 ms 6,944 KB
testcase_100	AC	1 ms 6,940 KB
testcase_101	WA	-
testcase_102	AC	2 ms 6,944 KB
testcase_103	AC	1 ms 6,940 KB
testcase_104	AC	2 ms 6,944 KB
testcase_105	AC	2 ms 6,944 KB
testcase_106	AC	2 ms 6,944 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#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 __int128_t ll;
int const inf = 1<<29;

namespace math {

ll mulmod(ll a, ll b, ll mod) {
  ll x = 0, y = a % mod;
  while(b > 0) {
    if(b & 1) {
      x = (x + y) % mod;
    }
    y = (y * 2) % mod;
    b /= 2;
  }
  return x % mod;
}

}

namespace math {

ll modular_exp(ll base, ll exp, ll mod) {
  ll x = 1;
  ll y = base;
  while(exp > 0) {
    if(exp & 1) {
      x *= y; x %= mod;
    }
    y *= y; y %= mod;
    exp >>= 1;
  }
  return x % mod;
}

}

struct xor128 {
  unsigned x,y,z,w;
  xor128(): x(89165727), y(157892372), z(7777777), w(757328) {}
  unsigned next() {
    unsigned t=x^(x<<11);
    x=y;y=z;z=w;
    return w=w^(w>>19)^t^(t>>8);
  }
  unsigned next(unsigned k) {
    return next()%k;
  }
} rndgen;

namespace math {
namespace prime {

bool miller_rabin(ll p, int iteration) {
  if(p < 2) { return false; }
  if(p != 2 && p & 1) { return false; }
  ll s = p - 1;
  while(s & 1) { s >>= 1; }
  rep(_, iteration) {
    ll a = rndgen.next(p - 1) + 1;
    ll temp = s;
    ll mod = modular_exp(a, temp, p);
    while(temp != p - 1 && mod != 1 && mod != p - 1) {
      mod = mulmod(mod, mod, p);
      temp <<= 1;
    }
    if(mod != p - 1 && temp & 1) {
      return false;
    }
  }
  return true;
}

}
}

set<ll> carmichael = {
  561, 1105, 1729, 2465, 2821, 6601, 8911, 10585, 15841, 29341, 41041, 46657, 52633, 62745, 63973, 75361, 101101, 115921, 126217, 162401, 172081, 188461, 252601, 278545, 294409, 314821, 334153, 340561, 399001, 410041, 449065, 488881, 512461, 530881, 552721, 656601, 658801, 670033, 748657, 825265, 838201, 852841, 997633, 1024651, 1033669, 1050985, 1082809, 1152271, 1193221, 1461241, 1569457, 1615681, 1773289, 1857241, 1909001, 2100901, 2113921, 2433601, 2455921, 2508013, 2531845, 2628073, 2704801, 3057601, 3146221, 3224065, 3581761, 3664585, 3828001, 4335241, 4463641, 4767841, 4903921, 4909177, 5031181, 5049001, 5148001, 5310721, 5444489, 5481451, 5632705, 5968873, 6049681, 6054985, 6189121, 6313681, 6733693, 6840001, 6868261, 7207201, 7519441, 7995169, 8134561, 8341201, 8355841, 8719309, 8719921, 8830801, 8927101, 9439201, 9494101, 9582145, 9585541, 9613297, 9890881,
  10877581, 12945745, 13992265, 16778881, 18162001, 27336673, 28787185, 31146661, 36121345, 37167361, 40280065, 41298985, 41341321, 41471521, 47006785, 67371265, 67994641, 69331969, 74165065,  
    321197185, 413631505, 417241045, 496050841, 509033161, 611397865, 612347905, 638959321, 672389641, 832060801, 834720601, 868234081, 945959365, 986088961, 1074363265, 1177800481, 1210178305, 1256855041, 1410833281, 1481619601,
5394826801, 6295936465, 12452890681, 13577445505, 15182481601, 20064165121, 22541365441, 24673060945, 26242929505, 26602340401, 27405110161, 28553256865, 33203881585, 38059298641, 39696166081, 40460634865 ,
  232250619601, 306177962545, 432207073585, 576480525985, 658567396081, 689702851201, 747941832001, 1013666981041, 1110495895201, 1111586883121, 1286317859905, 1292652236161, 1341323384401, 1471186523521, 1567214060545,
    9746347772161, 11537919313921, 11985185775745, 14292786468961, 23239986511105, 24465723528961, 26491881502801, 27607174936705, 30614445878401, 30912473358481, 34830684315505, 51620128928641,    
};

inline ll pow(ll x, ll n, ll mod){
  ll y = 1;
  if(n){
    y = pow(x, n >> 1, mod);
    if (!(n & 1)) { y = (y * y) % mod; }
    else { y = (((y * y) % mod) * x) % mod; }
  }
  return y;
}

bool fermat_test(ll const& x, ll a) {
  if(carmichael.count(x)) { return false; }
  if(x == a) { return 1; }
  if(x < a || (x % a == 0)) { return 0; }
  return pow(a, x-1, x) == 1;
}

bool fermat_test(ll const& x) {
  return fermat_test(x, 2) && fermat_test(x, 3);
}

map<ll, int> mp = {
  {4,3}, 
  {6,5}, 
  {8,7}, 
  {9,7}, 
  {10,7}, 
  {12,11}, 
  {14,13}, 
  {15,7}, 
  {16,7}, 
  {18,8}, 
  {20,19}, 
  {21,19}, 
  {22,7}, 
  {24,23}, 
  {25,23},
};

int main() {

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

  if(mp.find(N) != mp.end()) {
    cout << mp[N] << endl;
    exit(0);
  }

  if(N < 114514893) {
    if((N % 8 == 1 || fermat_test(N-1)) && fermat_test(N-8)) {
      cout << 14 << endl;
    }
    else {
      cout << 8 << endl;
    }
  }
  else {
    if((N % 8 == 1) && math::prime::miller_rabin(N-8, 100)) {
      cout << 14 << endl;
    }
    else {
      cout << 8 << endl;
    }
  }

  return 0;
}

yukicoder

結果

テストケース

ソースコード