結果
問題 | No.1514 Squared Matching |
ユーザー |
![]() |
提出日時 | 2021-04-29 20:18:23 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2,307 ms / 4,000 ms |
コード長 | 1,587 bytes |
コンパイル時間 | 2,113 ms |
コンパイル使用メモリ | 199,600 KB |
最終ジャッジ日時 | 2025-01-21 02:12:57 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
コンパイルメッセージ
main.cpp: In function ‘int main()’: main.cpp:56:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 56 | scanf("%d", &n); | ~~~~~^~~~~~~~~~
ソースコード
#include <bits/stdc++.h> using namespace std; #pragma region Math Osa_k struct osa_k { private: std::vector<int> spf, pr; public: osa_k() = default; osa_k(int MAX) : spf(MAX + 1) { for(int i = 2; i <= MAX; i++) { if(spf[i] == 0) { spf[i] = i; pr.push_back(i); } for(int j = 0; j < pr.size() and pr[j] <= spf[i] and i * pr[j] <= MAX; j++) spf[i * pr[j]] = pr[j]; } } std::vector<std::pair<int, int>> prime_factorization(int n) { std::vector<std::pair<int, int>> divisor; if(n == 1) return divisor; int before = spf[n], cnt = 0; while(n > 1) { if(spf[n] == before) { cnt++; n /= spf[n]; } else { divisor.emplace_back(before, cnt); before = spf[n]; cnt = 1; n /= spf[n]; } } divisor.emplace_back(before, cnt); return divisor; } int smallestprimefactor(const int n) const { return spf[n]; } bool is_prime(const int n) const { return n == spf[n]; } }; #pragma endregion int embed[6] = {0, 99371160, 207162592, 318140852, 431159896, 545761220}; int n, res; int main() { scanf("%d", &n); osa_k pf(n); /* 愚直 for(int a = 1; a <= n; a++) { auto v = pf.prime_factorization(a); int b = 1; for(auto [val, cnt] : v) if(cnt & 1) b *= val; res += sqrt(n / b); } printf("%d\n", res); */ int base = n / 10000000 * 10000000; res = embed[n / 10000000]; for(int a = base + 1; a <= n; a++) { auto v = pf.prime_factorization(a); int b = 1; for(auto [val, cnt] : v) if(cnt & 1) b *= val; int d = sqrt(n / b); res += d + min(d, (int)sqrt(base / b)); } printf("%d\n", res); }