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

提出ソース

結果

問題	No.308 素数は通れません
ユーザー	minami
提出日時	2019-03-30 18:43:12
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	WA
実行時間	-
コード長	3,407 bytes
コンパイル時間	7,089 ms
コンパイル使用メモリ	401,128 KB
実行使用メモリ	6,824 KB
最終ジャッジ日時	2024-11-15 22:25:07
合計ジャッジ時間	8,315 ms
ジャッジサーバーID （参考情報）	judge4 / judge2

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

テストケース

テストケース表示

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

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

ソースコード

raw source code

#include "bits/stdc++.h"
using namespace std;
#ifdef _DEBUG
#include "dump.hpp"
#else
#define dump(...)
#endif

#include <boost/multiprecision/cpp_int.hpp>
namespace mp = boost::multiprecision;
using Int = mp::cpp_int;

#define int long long
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
#define all(c) begin(c),end(c)
const int INF = sizeof(int) == sizeof(long long) ? 0x3f3f3f3f3f3f3f3fLL : 0x3f3f3f3f;
const int MOD = 1'000'000'007;
template<class T> bool chmax(T &a, const T &b) { if (a < b) { a = b; return true; } return false; }
template<class T> bool chmin(T &a, const T &b) { if (b < a) { a = b; return true; } return false; }


// エラトステネスの篩
vector<char> eratos(int n) {
	vector<char> is_prime(n + 1, true);
	is_prime[0] = is_prime[1] = false;
	for (int i = 2; i*i <= n; i++)
		if (is_prime[i]) {
			int j = i + i;
			while (j <= n) {
				is_prime[j] = false;
				j += i;
			}
		}
	return is_prime;
}

// (a*b) % mod 
template<typename Int>
Int modmul(Int a, Int b, Int mod) {
	Int x = 0, y = a % mod;
	while (b > 0) {
		if (b & 1)x = x + y % mod;
		y = y * 2 % mod;
		b >>= 1;
	}
	return x % mod;
}

// mod^2 が long long の最大値より大きければオーバーフローするので掛け算に modmul を使う
template<typename Int>
Int modpow(Int a, Int e, Int mod) {
	Int res = 1;
	while (e > 0) {
		if (e & 1)res = modmul(res, a, mod);
		a = modmul(a, a, mod);
		e >>= 1;
	}
	return res;
}

// 素数判定（Miller-Rabin primality test）
// 2^24程度から
// long long で2^32以上を判定する場合は modpow で modmul を使う
// millerRabinPrimalityTest(n, 5)
template<typename Int>
bool millerRabinPrimalityTest(Int x, int iteration) {
	if (x < 2)return false;
	if (x != 2 && x % 2 == 0)return false;
	Int s = x - 1;
	while (s % 2 == 0)s /= 2;
	for (int i = 0; i < iteration; i++) {
		Int a = rand() % (x - 1) + 1, tmp = s;
		Int mod = modpow(a, tmp, x);
		while (tmp != x - 1 && mod != 1 && mod != x - 1) {
			mod = modmul(mod, mod, x);
			tmp *= 2;
		}
		if (mod != x - 1 && tmp % 2 == 0)return false;
	}
	return true;
}

signed main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	int N; cin >> N;
	if (N < 50) {
		auto e = eratos(N);
		static const int di[] = { 1,0,-1,0 };
		static const int dj[] = { 0,1,0,-1 };
		using State = tuple<int, int, int>;
		rep(w, 1, N + 1) {
			dump(w);
			vector<bool> f(N + 1);
			int H = (N - 1) / w + 1, W = w;
			auto inrange = [&](int i, int j) { return i >= 0 && i < H && j >= 0 && j < W; };
			auto bfs = [&](int si, int sj, int gi, int gj) {
				queue<State> q;
				q.emplace(0, si, sj);
				while (q.size()) {
					int d, ci, cj; tie(d, ci, cj) = q.front(); q.pop();
					if (ci == gi && cj == gj)return d;
					rep(i, 0, 4) {
						int ni = ci + di[i], nj = cj + dj[i];
						int nx = ni * w + nj + 1;
						if (!inrange(ni, nj))continue;
						if (nx > N)continue;
						if (e[nx])continue;
						if (f[nx])continue;
						f[nx] = true;
						q.emplace(d + 1, ni, nj);
					}
				}
				return -1LL;
			};
			auto res = bfs(0, 0, (N - 1) / w, (N - 1) % w);
			dump(res);
			if (res != -1) {
				cout << w << endl;
				break;
			}
		}
	}
	else {
		for (int w = 8;; w += 2) {
			if (millerRabinPrimalityTest(Int(1 + w), 5))continue;
			if (((N - 1) % w) == 0 && millerRabinPrimalityTest(Int(N - w), 5))continue;
			cout << w << endl;
			break;
		}
	}
	return 0;
}

yukicoder

結果

テストケース

ソースコード