結果
問題 |
No.826 連絡網
|
ユーザー |
![]() |
提出日時 | 2019-05-03 22:16:49 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,150 bytes |
コンパイル時間 | 2,104 ms |
コンパイル使用メモリ | 176,796 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-12-31 18:06:03 |
合計ジャッジ時間 | 10,172 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | WA * 30 |
ソースコード
#include "bits/stdc++.h" #define in std::cin #define out std::cout #define rep(i,N) for(LL i=0;i<N;++i) typedef long long int LL; std::map<LL, LL> decomposit_prime(LL n) { LL d = 2; std::map<LL, LL>res; while (d * d <= n) { if (n % d == 0) { ++res[d]; n /= d; } else ++d; } ++res[n]; return res; } LL gcd(LL a, LL b) { if (a < b) std::swap(a, b); LL r = a % b; while (r != 0) { a = b; b = r; r = a % b; } return b; } LL lcm(LL a, LL b) { return a / gcd(a, b) * b; } int main() { LL N, P; in >> N >> P; std::vector<LL>primes, nums; for (LL i = 2; i <= N; ++i) { bool flag = true; for (LL j = 2; j * j <= i; ++j) { if (i % j == 0) { flag = false; break; } } if (flag) primes.push_back(i); } auto temp = decomposit_prime(P); for (auto t : temp) nums.push_back(t.first); std::vector<LL>ok; for (LL i : primes) { for (LL j : nums) { if (lcm(i, j) <= N) { ok.push_back(i); break; } } } LL alllcm = 1; for (LL i : ok) alllcm = lcm(alllcm, i); LL ans = 0; for (LL i = 2; i <= N; ++i) { if (gcd(alllcm, i) != 1) ++ans; } out << ans << std::endl; }