結果
問題 |
No.2552 Not Coprime, Not Divisor
|
ユーザー |
|
提出日時 | 2023-11-25 14:58:00 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,038 bytes |
コンパイル時間 | 1,924 ms |
コンパイル使用メモリ | 197,260 KB |
最終ジャッジ日時 | 2025-02-18 00:51:55 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 TLE * 3 |
ソースコード
#include <bits/stdc++.h> #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") using namespace std; int main() { cin.tie(0); cout.tie(0); ios::sync_with_stdio(false); int N; cin >> N; long long ans = 0; for(int i = 1; i < N; i++) { int j = i; vector<int> p; for(int k = 2; k * k <= i; k++) { if(j % k == 0) { p.emplace_back(k); while(j % k == 0) { j /= k; } } if(j == 1) { break; } } if(j > 1) { p.emplace_back(j); } int sz = p.size(); for(int bit = 1; bit < 1 << sz; bit++) { int S = 1; for(int j = 0; j < sz; j++) { if(bit >> j & 1) { S *= p[j]; } } int add = N / S - i / S; if(__builtin_popcount(bit) % 2 == 0) { add *= -1; } //cout << i << ' ' << bit << ' ' << add << endl; ans += add; } if((int) p.size() > 0) { ans -= N / i - 1; } //cout<<i<<' '<<ans<<endl; } cout << ans << '\n'; return 0; }