結果
問題 | No.2011 Arbitrary Mod (Hidden) |
ユーザー | 👑 emthrm |
提出日時 | 2022-07-15 23:07:56 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,094 bytes |
コンパイル時間 | 2,115 ms |
コンパイル使用メモリ | 205,816 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-27 20:23:39 |
合計ジャッジ時間 | 5,240 ms |
ジャッジサーバーID (参考情報) |
judge5 / 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; } } } }