結果
問題 | No.732 3PrimeCounting |
ユーザー |
![]() |
提出日時 | 2022-07-12 04:00:29 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 436 ms / 3,000 ms |
コード長 | 806 bytes |
コンパイル時間 | 4,625 ms |
コンパイル使用メモリ | 251,932 KB |
最終ジャッジ日時 | 2025-01-30 06:41:49 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 89 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; #define rep2(i, m, n) for (int i = (m); i < (n); ++i) #define rep(i, n) rep2(i, 0, n) #define drep2(i, m, n) for (int i = (m)-1; i >= (n); --i) #define drep(i, n) drep2(i, n, 0) typedef long long ll; int main(){ int n; cin >> n; vector<ll> f1(n+1); vector<ll> g2(2*n+1); vector<bool> dp(3*n+1,false); for (int i=2; i<=3*n; i++){ if (!dp[i]){ if (i<=n){ f1[i] += 1; g2[2*i] += 1; } for (int j=2*i; j<=3*n; j+=i){ dp[j] = true; } } } vector<ll> f2 = convolution_ll(f1, f1); vector<ll> f3 = convolution_ll(f2, f1); vector<ll> g3 = convolution_ll(g2, f1); ll ans = 0; for (int i=2; i<=3*n; i++){ if (!dp[i]){ ans += (f3[i] - 3*g3[i])/6; } } cout << ans << endl; }