結果
問題 | No.2011 Arbitrary Mod (Hidden) |
ユーザー |
![]() |
提出日時 | 2022-07-15 23:07:56 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 35 ms / 2,000 ms |
コード長 | 2,094 bytes |
コンパイル時間 | 2,152 ms |
コンパイル使用メモリ | 195,916 KB |
最終ジャッジ日時 | 2025-01-30 08:51:14 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 43 |
ソースコード
#define _USE_MATH_DEFINES#include <bits/stdc++.h>using namespace std;#define FOR(i,m,n) for(int i=(m);i<(n);++i)#define REP(i,n) FOR(i,0,n)#define ALL(v) (v).begin(),(v).end()using ll = long long;constexpr int INF = 0x3f3f3f3f;constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;constexpr double EPS = 1e-8;constexpr int MOD = 1000000007;// constexpr int MOD = 998244353;constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};template <typename T, typename U>inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }template <typename T, typename U>inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }struct IOSetup {IOSetup() {std::cin.tie(nullptr);std::ios_base::sync_with_stdio(false);std::cout << fixed << setprecision(20);}} iosetup;template <typename T>std::vector<std::pair<T, int>> prime_factorization(T n) {std::vector<std::pair<T, int>> res;for (T i = 2; i * i <= n; ++i) {if (n % i == 0) {int exponent = 0;for (; n % i == 0; n /= i) {++exponent;}res.emplace_back(i, exponent);}}if (n > 1) res.emplace_back(n, 1);return res;}long long euler_phi(long long n) {assert(n >= 1);long long res = n;for (long long i = 2; i * i <= n; ++i) {if (n % i == 0) {res -= res / i;while (n % i == 0) n /= i;}}return n > 1 ? res - res / n : res;}bool is_prime(const long long n) {if (n <= 1) return false;for (long long i = 2; i * i <= n; ++i) {if (n % i == 0) return false;}return true;}int main() {cout << 16711935 << '\n' << 1 << '\n';return 0;set<ll> bad_prime;for (int a = 20; a <= 5000; ++a) {for (const auto [p, _] : prime_factorization(a * a - 398)) bad_prime.emplace(p);}ll m = 1;for (ll p2 = 2; ; p2 *= 2) {if (!bad_prime.count(p2 + 1) && is_prime(p2 + 1)) {m *= p2 + 1;if (m >= 1000000) {cout << m << '\n';return 0;}}}}