結果

問題 No.732 3PrimeCounting
コンテスト
ユーザー donkorin_
提出日時 2018-09-09 16:24:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 259 ms / 3,000 ms
コード長 1,429 bytes
コンパイル時間 1,437 ms
コンパイル使用メモリ 169,656 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-23 14:39:19
合計ジャッジ時間 6,090 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 89
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<(b);++i)
#define erep(i,a,b) for(int i=a;i<=(int)(b);++i)
#define per(i,a,b) for(int i=(a);i>(b);--i)
#define eper(i,a,b) for(int i=(a);i>=b;--i)
#define pb push_back
#define mp make_pair
#define INF (1<<30)-1
#define MOD 1000000007
#define all(x) (x).begin(),(x).end()
#define vii vector<int>
#define vll vector<long long>
using namespace std;
typedef long long ll;
typedef pair<int,int> Pii;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a = b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a = b; return 1; } return 0; }
int dy[]={0, 0, 1, -1};
int dx[]={1, -1, 0, 0};
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int lcm(int a,int b){return a/gcd(a, b)*b;}

vector<int> makePrimes(int n)
{
    vector<int> res, pr(n + 1, 1);
    pr[0] = pr[1] = 0;
    rep(p, 2, sqrt(n) + 2) if (pr[p]) for (int x = p * 2; x <= n; x += p) pr[x] = 0;
    rep(p, 2, n + 1) if (pr[p])
        res.push_back(p);
    return res;
}

int n, cnt[400001];
ll ans;
int main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
    cin >> n;
    auto primes = makePrimes(3 * n + 1);
    for (int i = 0; i < primes.size() && primes[i] <= n; ++i) {
        for (int j = primes.size() - 1; primes[j] > primes[i]; --j) ans += cnt[primes[j] - primes[i]];
        rep(k, 0, i) cnt[primes[i] + primes[k]]++;
    } 
    cout << ans << endl;
    return 0;
}
0