結果

問題 No.732 3PrimeCounting
コンテスト
ユーザー 梧桐
提出日時 2025-12-23 22:32:00
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,140 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 766 ms
コンパイル使用メモリ 71,732 KB
実行使用メモリ 7,852 KB
最終ジャッジ日時 2025-12-23 22:33:02
合計ジャッジ時間 52,926 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 89
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:33:12: warning: ignoring return value of ‘FILE* freopen(const char*, const char*, FILE*)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   33 |     freopen("psum.in", "r", stdin);
      |     ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~
main.cpp:34:12: warning: ignoring return value of ‘FILE* freopen(const char*, const char*, FILE*)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   34 |     freopen("psum.out", "w", stdout);
      |     ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~
main.cpp:37:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   37 |     scanf("%d", &n);
      |     ~~~~~^~~~~~~~~~

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

typedef long long LL;

const int N = 100010;

int n;
LL ans;
bool is_prime[N];
vector<int> primes;

void EulerSieve() {
    fill(is_prime, is_prime + N, true);
    is_prime[0] = is_prime[1] = false;
    for (int i = 2; i < N; ++i) {
        if (is_prime[i]) primes.push_back(i);
        for (auto p : primes) {
            if (1LL * i * p >= N) continue;
            is_prime[i * p] = false;
            if (i % p == 0) {
                break;
            }
        }
    }
}

// 1 - 100000有9527个质数
int main() {
    freopen("psum.in", "r", stdin);
    freopen("psum.out", "w", stdout);
    
    EulerSieve();
    scanf("%d", &n);

    for (int i = 0; i < primes.size(); ++i) {
        if (primes[i] > n) break;
        for (int j = i + 1; j < primes.size(); ++j) {
            if (primes[j] > n) break;
            for (int k = j + 1; k < primes.size(); ++k) {
                if (primes[k] > n) break;
                if (is_prime[primes[i] + primes[j] + primes[k]]) ++ans;
            }
        }
    }

    printf("%lld\n", ans);

    return 0;
}
0